我们发现(看题解)知道当区间长度大于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)即可
接下来是怎么查询的问题
证明:长度为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;
}