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;
}