[Ynoi2012]NOIP2016人生巅峰

链接：https://www.luogu.com.cn/problem/P5527

题目描述：题目已经说的很清楚了。

题解：操作\(2\)容易让人往那个奇怪的线段树操作那一块想，但假设这么做的话还是没法做操作\(1\)。

考虑第一种操作的性质，我们发现：

\(1\).假如\([l,r]\)中有相等的元素，那么他一定能被表示出来。

\(2\).假如\([l,r]\)中没有相等的元素，那么他可以表示\(2\)的指数级别的数。

我们发现\(2\)在区间长度大于\(10\)时，就能表示出\(1024\)个数左右，由于要对小于\(1000\)的数取模，那么一定存在相同的数。那么我们可以将区间长度大于\(10\)时的直接判断他一定能被表示出来。

然后这道题就变成了一道搜索题。

但是还有\(2\)操作，由于只需要单点查询，不需要用到那个奇怪的线段树操作，我们可以开一个树状数组记以下每个数被开的幂，然后用扩展欧拉定理求这个数就行了。

#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
long long c[1000001],n,q,ans,s[1000001],phi_v,v,t[1000001];
bool cmp;
inline int lowbit(register int x)
{
   return x&(-x);
}
inline void add(register int x,register int y)
{
    for (;x<=n;x+=lowbit(x))
    	c[x]+=y;
  return;
}
inline int sum(register int x)
{
	ans=0;
	for (;x>=1;x-=lowbit(x))
		ans+=c[x];
	return ans;
}
inline long long fast_pow(register long long a,register long long b,register long long mod)
{
	if (b==0)
    	return 1;
	if (b&1)
		return fast_pow(a*a%mod,b/2,mod)*a%mod;
	else
		return fast_pow(a*a%mod,b/2,mod);
}
inline void dfs(register int x,register int y,register int sum1,register int sum2,register bool used)
{
	if (x==y+1)
	{
    	if (sum1==sum2&&used)
			cmp=1;
    	return;
  	}
	dfs(x+1,y,sum1+s[x]+1,sum2,1);
  	if (cmp)
		return;
  	dfs(x+1,y,sum1,sum2+s[x]+1,1);
  	if (cmp)
   		return;
  	dfs(x+1,y,sum1,sum2,used);
 	if (cmp)
    	return;
  	return;
}
int main()
{
	register int l,r,que;
    scanf("%lld%lld%lld",&n,&q,&v);
    for (int i=1;i<=n;++i)
    	scanf("%lld",&t[i]);
    phi_v=v;
    int tmp=v;
    for (register int i=2;i*i<=v;++i)
    	if (tmp%i==0)
    	{
    		while (tmp%i==0)
    			tmp/=i;
    		phi_v=phi_v/i*(i-1);
		}
	if (tmp>1)
		phi_v=phi_v/tmp*(tmp-1);
  	for (register int i=1;i<=q;++i)
    {
    	scanf("%d%d%d",&que,&l,&r);
      	if (que==2)
        {
	 		add(l,1);
	  		add(r+1,-1);
        }
    	if (que==1)
    	{
	  		if (r-l>11)
	      		puts("Yuno");
	  		else
	   		{
	    		cmp=0;
	      		for (register int j=l;j<=r;++j)
	      		{
	      			if (sum(j)>log(phi_v)/log(3))
						s[j]=fast_pow(t[j],fast_pow(3,sum(j),phi_v)+phi_v,v);
					else
						s[j]=fast_pow(t[j],fast_pow(3,sum(j),phi_v),v);
				}
	     		dfs(l,r,0,0,0);
	      		if (cmp==0)
					puts("Yuki");
	      		else
					puts("Yuno");
	    	}
		}
    }
	return 0;
}