[bzoj4722] 由乃
题意：两种操作。1、将所给区间划分成两个集合，集合中的元素的贡献为a[i]+1，求是否能找出两个集合使两个集合的总贡献相等。2、区间立方（有模数）
题解：
倍增+线段树+搜索（meet in the middle）
区间立方：由于模数比较小，我们可以用倍增预处理某个数的2^j次立方，然后用线段树维护查询的次数，单次复杂度：O（logn）
划分集合：
首先有个结论，由于所有数的最大范围只有1000，所以超过13个数一定能通过加减等于0
对于范围小于等于13的区间，首先查出每个数操作后得到的数，然后用meet in the middle暴搜出解，单次复杂度：O（logn+3^7）
总复杂大概是：O（m*(logn+3^7)）
然后这题我犯了一大堆傻逼错误，调了一个下午，我也是醉了
1 #include<iostream> 2 #include<cstdio> 3 #include<cstdlib> 4 #include<cstring> 5 #include<algorithm> 6 #include<cmath> 7 #define ll long long 8 #define ls x<<1 9 #define rs x<<1|1 10 #define N 100010 11 #define M 15000 12 using namespace std; 13 14 int n,m,v,top; 15 int a[N],tr[N*4],lazy[N*4],f[1010][22],st[N],val[20]; 16 bool flg,vis[N]; 17 18 int gi() { 19 int x=0,o=1; char ch=getchar(); 20 while(ch!='-' && (ch<'0' || ch>'9')) ch=getchar(); 21 if(ch=='-') o=-1,ch=getchar(); 22 while(ch>='0' && ch<='9') x=x*10+ch-'0',ch=getchar(); 23 return o*x; 24 } 25 26 void pre() { 27 for(int i=0; i<v; i++) f[i][0]=i*i*i%v;//错点1：初始化对象错误 28 for(int j=1; j<=20; j++) 29 for(int i=0; i<v; i++) 30 f[i][j]=f[f[i][j-1]][j-1]; 31 } 32 33 void pushdown(int x) { 34 if(!lazy[x]) return; 35 tr[ls]+=lazy[x],tr[rs]+=lazy[x]; 36 lazy[ls]+=lazy[x],lazy[rs]+=lazy[x]; 37 lazy[x]=0; 38 } 39 40 void update(int x, int l, int r, int ql, int qr) { 41 if(ql<=l && r<=qr) { 42 tr[x]++,lazy[x]++;//其实只要开一个就够了，反正每次只会查询叶子= =、 43 return; 44 } 45 pushdown(x); 46 int mid=(l+r)>>1; 47 if(qr<=mid) update(ls,l,mid,ql,qr); 48 else if(ql>mid) update(rs,mid+1,r,ql,qr); 49 else update(ls,l,mid,ql,mid),update(rs,mid+1,r,mid+1,qr); 50 } 51 52 int query(int x, int l, int r, int qx) { 53 if(l==r) return tr[x]; 54 pushdown(x);//错点3：查询时没有pushdown 55 int mid=(l+r)>>1; 56 if(qx<=mid) return query(ls,l,mid,qx); 57 else return query(rs,mid+1,r,qx); 58 } 59 60 int get(int x) { 61 int y=query(1,1,n,x),ret=a[x]; 62 for(int i=20; i>=0; i--) { 63 if(y&(1<<i)) ret=f[ret][i]; 64 } 65 return ret; 66 } 67 68 void dfs1(int dep, int k, int sum, int siz) { 69 if(dep>k) { 70 vis[sum+M]=1,st[++top]=sum; 71 if(!sum && siz) flg=1; 72 return; 73 } 74 dfs1(dep+1,k,sum,siz); 75 if(flg) return; 76 dfs1(dep+1,k,sum+val[dep]+1,siz+1); 77 if(flg) return; 78 dfs1(dep+1,k,sum-val[dep]-1,siz+1); 79 } 80 81 void dfs2(int dep, int k, int sum, int siz) { 82 if(dep>k) { 83 if((vis[-sum+M] && sum!=0) || (!sum && siz)) flg=1; 84 return; 85 } 86 dfs2(dep+1,k,sum,siz); 87 if(flg) return;//错点5：没打分号 88 dfs2(dep+1,k,sum+val[dep]+1,siz+1); 89 if(flg) return; 90 dfs2(dep+1,k,sum-val[dep]-1,siz+1); 91 } 92 93 bool check(int x) { 94 for(int i=1; i<=top; i++) vis[st[i]+M]=0; 95 top=0,flg=0; 96 dfs1(1,x/2,0,0); 97 if(flg) return true; 98 dfs2(x/2+1,x,0,0); 99 if(flg) return true; 100 return false;//错点4：返回false时没有清空数组 101 } 102 103 int main() { 104 n=gi(),m=gi(),v=gi(); 105 for(int i=1; i<=n; i++) a[i]=gi(); 106 pre(); 107 for(int i=1; i<=m; i++) { 108 int op=gi(),l=gi(),r=gi(); 109 if(op==2) update(1,1,n,l,r); 110 else { 111 if(r-l+1>13) puts("Yuno"); 112 else { 113 for(int j=l; j<=r; j++) val[j-l+1]=get(j);//错点2：重复用i循环 114 if(check(r-l+1)) puts("Yuno"); 115 else puts("Yuki"); 116 } 117 } 118 } 119 return 0; 120 }
