bzoj4722: 由乃

我们发现（看题解）知道当区间长度大于13时一定可以找到两集合总贡献值相等
证明：长度为L的区间共有2^L种集合，而每个元素贡献小于等于1000（v<=1000），所以最多有1000*L种集合值，当L>13时2^L>1000*L
于是乎解决了第一个问题，接下来发现区间问题，那就线段树搞一搞233
……可修改操作是将原数立方，那单重不断标记叠加就炸了……
由于v很小，然后发现倍增可以很好解决这一问题，
设f[a][b]数组，表示a的2^b次立方叠乘的结果，那每个点最多立方100000次，所以只要预处理到2^16次（即b<=16）即可
接下来是怎么查询的问题

我们发现可以二分搜索，将整个区间分为两部分，在每个部分中枚举两个不相交的集合并算出两集合的差，若两部分算出来的差有相等则证明可以找到

搜素时要注意判断两集合皆空的情况

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<string>
#include<queue>
using namespace std;
int a[100005];
struct data{int l,r,lzy;}d[400005];
int n,m,v;
int f[1005][20];
int sum,op;
int e[100005];
void build(int k,int lc,int rc)
{
     int mid=(lc+rc)>>1;
     d[k].l=lc;
     d[k].r=rc;
     if(lc==rc)return ;
     build(k<<1,lc,mid);
     build((k<<1)+1,mid+1,rc);
}
void change(int lc,int rc,int k)
{
     if(lc==d[k].l&&rc==d[k].r){d[k].lzy++;return ;}
     if(rc<=d[k<<1].r){change(lc,rc,k<<1);return ;}
     if(lc>=d[(k<<1)+1].l){change(lc,rc,(k<<1)+1);return ;}
     change(lc,d[k<<1].r,k<<1);change(d[(k<<1)+1].l,rc,(k<<1)+1);
}
void dream(int x,int k)
{
     int s=0,t=a[x];
     while(k!=0){
         if(k&1)t=f[t][s];
         k/=2;
         s++;
         }
     a[x]=t;
}
void find(int lc,int rc,int k)
{
     if(d[k].l==d[k].r){dream(lc,d[k].lzy);d[k].lzy=0;return ;}
     if(d[k].lzy!=0){d[k<<1].lzy+=d[k].lzy;d[(k<<1)+1].lzy+=d[k].lzy;d[k].lzy=0;}
     if(lc==d[k].l&&rc==d[k].r){find(d[k<<1].l,d[k<<1].r,k<<1);find(d[(k<<1)+1].l,d[(k<<1)+1].r,(k<<1)+1);return ;}
     if(rc<=d[k<<1].r){find(lc,rc,k<<1);return ;}
     if(lc>=d[(k<<1)+1].l){find(lc,rc,(k<<1)+1);return ;}
     find(lc,d[k<<1].r,k<<1);find(d[(k<<1)+1].l,rc,(k<<1)+1);
}
void lcr(int lc,int rc,int ans,int ti)
{
     int k;
     if(lc>rc){
        if(op!=0){
        k=ans;
        if(k<0)k=-k;
        e[k]=ti;
        }
        return ;
        }
     ans+=(a[lc]+1);
     op++;
     lcr(lc+1,rc,ans,ti);
     op--;
     ans-=(a[lc]+1);
     lcr(lc+1,rc,ans,ti);
     op++;
     ans-=(a[lc]+1);
     lcr(lc+1,rc,ans,ti);
     op--;
}
int drop(int lc,int rc,int ans,int ti)
{
     int k,tx;
     if(lc>rc){
        if(op!=0){
        if(ans==0)return 1;
        k=ans;
        if(k<0)k=-k;
        if(e[k]==ti)return 1;
        else return 0;
        }
        if(e[0]==ti)return 1;
        else return 0;
        }
     ans+=(a[lc]+1);
     op++;
     tx=drop(lc+1,rc,ans,ti);
     if(tx==1)return 1;
     op--;
     ans-=(a[lc]+1);
     tx=drop(lc+1,rc,ans,ti);
     if(tx==1)return 1;
     op++;
     ans-=(a[lc]+1);
     tx=drop(lc+1,rc,ans,ti);
     if(tx==1)return 1;
     op--;
}
int main()
{
    int i,j,k,opt,lc,rc,lef=0;
    scanf("%d%d%d",&n,&m,&v);
    for(i=1;i<=n;i++)scanf("%d",&a[i]);
    build(1,1,n);
    for(i=0;i<v;i++)f[i][0]=(i*i*i)%v;
    for(i=1;i<=16;i++)
        for(j=0;j<v;j++)f[j][i]=f[f[j][i-1]][i-1];
    for(i=1;i<=m;i++){
        scanf("%d%d%d",&opt,&lc,&rc);
        if(opt==1){
           if((rc-lc)>13)printf("Yuno\n");
           else{ 
                 find(lc,rc,1);
                 op=0;
                 lcr(((lc+rc)>>1)+1,rc,0,i);
                 op=0;
                 if(drop(lc,(lc+rc)>>1,0,i)==1)printf("Yuno\n");
                 else printf("Yuki\n");
                 }
           }
        else change(lc,rc,1);
        }
    return 0;
}