[bzoj4722]由乃

身为10班人，就凭标题，这道题是一定要做的。

但是做了才发现有毒。。。。所以是信念和题解和大腿支撑了我！

先“假设”自己实力过硬，推出了结论：当区间过大时，必定存在一种方案可以输出”Yuno“。剩下就是这区区几个数的乱搞了。

数据结构学傻先不动声色地敲了一个线段树维护所有的立方操作，并且我发现自从noip2015挂过lazy tag以后就再也不会写错了。。。。这就是传说中教训的力量

然后我还是请教了一些更优秀的维护方法，不过万变不离其宗，这个常数说不要就不要了吧。。。。

接下来，我开始想当然：差不多了。。。再敲个快速幂。。。随便搞下。。。差不多了。。。

后来幸好看了下V，不是质数，那又要求phi啊？！

而且不一定互质，不仅要求phi而且要提取质因数啊？！

不想写啊，只能向倍增势力低头。

再后来，思考“区区几个数的乱搞”，好像并不容易。看了题解说得很含蓄，发现是个乱搞的思想。

对于区间【L，R】，以mid为界，左边枚举两个不相交子集，右边枚举两个不相交子集，判断它们的合法性。

于是就T了，原来我的复杂度和常数达到了十分可怕的境地。这就很不爽了。因为我不会卡常。

再回过头来看那个“假设”，有实力的人证一证，没实力的人写个暴力找下规律，发现可以把题解中的下界13改为10，这样就卡过去了。

（当然，你写9也能A）。。。。数据太随机了。

事情还没结束，time排名垫底了，原来是我乱搞写错了。四重循环，不挂才怪。正确做法还得写一个桶，把左边所有值扔进桶里，这样询问就飞起来了。

最后shy说：你这样也能过，说明常数实在太小了。

2333

改了一下。。。仍旧快垫底？

#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 200020
#define ll long long
],tree[N*];ll fa[N][];
],V[];
ll a[N],b[N];
void push_down(int p,int st,int ed){
    ){
        f[p+p]+=f[p];f[p+p+]+=f[p];;
        tree[p+p]+=f[p]*(mid-st+);tree[p+p+]+=f[p]*(ed-mid);
        f[p]=;
    }
}
void update(int L,int R,int st,int ed,int p){
    ;;
    )push_down(p,st,ed);
    if(R<=mid)update(L,R,st,mid,p+p);else
     ,ed,p+p+);else
      update(L,mid,st,mid,p+p),update(mid+,R,mid+,ed,p+p+);
}
int query(int x,int st,int ed,int p){
    if(st==x&&ed==x)return tree[p];
    )push_down(p,st,ed);
    ;
    if(x<=mid)return query(x,st,mid,p+p);else
               ,ed,p+p+);
}
ll quickmi(ll x,ll k){
    ll tmp=,a=x;
    ){
        ==){tmp=tmp*a%mo;}
        a=a*a%mo;k=k/;
    }return tmp;
}
void work(int st,int ed,int cc){
    ;
    <<(mid-st+))-,l2=(<<(ed-mid))-;
    ;u<=l1;u++)
     ;v<=l1-u;v++))||v!=){
            ,a2=,a3=,a4=;
            for(int i=st;i<=mid;i++){
                <<(i-st);
                )a1+=b[i]+;)a2+=b[i]+;
            }
            if(a1>a2)std::swap(a1,a2);
            ;return;}
            V[a2-a1]=;tot++;bb[tot]=a2-a1;
     }
       ;x<=l2;x++)
        ;y<=l2-x;y++))||y!=){
            ,a2=,a3=,a4=;
            ;i<=ed;i++){
                <<(i-mid-);
                )a3+=b[i]+;)a4+=b[i]+;
            }
            if(a3>a4)std::swap(a3,a4);
            ;return;}
            ){fool=;return;}
        } 

}
int main()
{
//freopen("yono.in","r",stdin);//freopen("my.out","w",stdout);
    scanf("%d%d%d",&n,&m,&mo);
    ;i<=n;i++)scanf("%lld",&a[i]);
    ;i<mo;i++) fa[i][]=i*i%mo*i%mo;
    ;j<=;j++)
      ;i<mo;i++)
        fa[i][j]=fa[fa[i][j-]][j-];
    ;cc<=m;cc++){
        int op,L,R;scanf("%d%d%d",&op,&L,&R);
        ){update(L,R,,n,);continue;}
        >)printf("Yuno\n");else{
            for(int i=L;i<=R;i++){
                ll tt=query(i,,n,),x=tt,tmp=,now=a[i];
                b[i]=a[i];
                ;j>=;j--)<<j)){
                   x-=<<j;b[i]=fa[b[i]][j]%mo;
                }
            }
            fool=;tot=;
            work(L,R,cc);
            ;i<=tot;i++)V[bb[i]]=;
            )printf("Yuno\n");else printf("Yuki\n");
        }
    }
} 

bzoj4722