bzoj 3207 花神的嘲讽计划Ⅰ
Description
背景
花神是神，一大癖好就是嘲讽大J，举例如下：
“哎你傻不傻的！【hqz：大笨J】”
“这道题又被J屎过了！！”
“J这程序怎么跑这么快！J要逆袭了！”
……
描述
这一天DJ在给吾等众蒟蒻讲题，花神在一边做题无聊，就跑到了一边跟吾等众蒟蒻一起听。以下是部分摘录：
1.
“J你在讲什么！”
“我在讲XXX！”
“哎你傻不傻的！这么麻烦，直接XXX再XXX就好了！”
“……”
2.
“J你XXX讲过了没？”
“……”
“那个都不讲你就讲这个了？哎你傻不傻的！”
“……”
DJ对这种情景表示非常无语，每每出现这种情况，DJ都是非常尴尬的。
经过众蒟蒻研究，DJ在讲课之前会有一个长度为N方案，我们可以把它看作一个数列；
同样，花神在听课之前也会有一个嘲讽方案，有M个，每次会在x到y的这段时间开始嘲讽，为了减少题目难度，每次嘲讽方案的长度是一定的，为K。
花神嘲讽DJ让DJ尴尬需要的条件：
在x~y的时间内DJ没有讲到花神的嘲讽方案，即J的讲课方案中的x~y没有花神的嘲讽方案【这样花神会嘲讽J不会所以不讲】。
经过众蒟蒻努力，在一次讲课之前得到了花神嘲讽的各次方案，DJ得知了这个消息以后欣喜不已，DJ想知道花神的每次嘲讽是否会让DJ尴尬【说不出话来】。
Input
第1行3个数N，M，K；
第2行N个数，意义如上；
第3行到第3+M-1行，每行K+2个数，前两个数为x,y,然后K个数，意义如上；
Output
对于每一个嘲讽做出一个回答会尴尬输出‘Yes’，否则输出‘No’
Sample Input
8 5 3
1 2 3 4 5 6 7 8
2 5 2 3 4
1 8 3 2 1
5 7 4 5 6
2 5 1 2 3
1 7 3 4 5
1 2 3 4 5 6 7 8
2 5 2 3 4
1 8 3 2 1
5 7 4 5 6
2 5 1 2 3
1 7 3 4 5
Sample Output
No
Yes
Yes
Yes
No
Yes
Yes
Yes
No
HINT
题中所有数据不超过2*10^9;保证方案序列的每个数字<=N
2~5中有2 3 4的方案，输出No，表示DJ不会尴尬
1~8中没有3 2 1的方案，输出Yes，表示DJ会尴尬
5~7中没有4 5 6的方案，输出Yes，表示DJ会尴尬
2~5中没有1 2 3的方案，输出Yes，表示DJ会尴尬
1~7中有3 4 5的方案，输出No，表示DJ不会尴尬
先字符串哈希，将所有长度为k的段哈希值按顺序存入主席树
询问时在对应主席树区间查找是否有这个哈希值
但此题比较玄学，试了很多模数都不行，要用unsigned类型自然溢出(相当于模一个很大的数)
这样就不用取模了，虽然会变慢
常数很大，要用一个技巧：
(l+r)/2写成(l>>1)+(r>>1)+(l&r&1)
1 #include<iostream> 2 #include<cstdio> 3 #include<cstring> 4 #include<algorithm> 5 #include<cmath> 6 using namespace std; 7 typedef unsigned long long ulol; 8 ulol inf=18446744073709551615ull; 9 ulol p=29; 10 ulol M=1; 11 int n,m,k,pos,ch[10000001][2],sum[10000001],a[100001],root[100001]; 12 ulol hash[100001]; 13 void update(int x,int &y,ulol l,ulol r,ulol k) 14 { 15 y=++pos; 16 ch[y][0]=ch[x][0];ch[y][1]=ch[x][1]; 17 sum[y]=sum[x]+1; 18 if (l==r) return; 19 ulol mid=(l>>1)+(r>>1)+(l&r&1); 20 if (k<=mid) update(ch[x][0],ch[y][0],l,mid,k); 21 else update(ch[x][1],ch[y][1],mid+1,r,k); 22 } 23 int query(int x,int y,ulol l,ulol r,ulol k) 24 { 25 if (l==r) return sum[y]-sum[x]; 26 ulol mid=(l>>1)+(r>>1)+(l&r&1); 27 if (k<=mid) return query(ch[x][0],ch[y][0],l,mid,k); 28 else return query(ch[x][1],ch[y][1],mid+1,r,k); 29 } 30 int main() 31 {int i,l,r,j; 32 ulol x; 33 ulol s; 34 cin>>n>>m>>k; 35 for (i=1;i<=n;i++) 36 scanf("%d",&a[i]); 37 for (i=1;i<=k;i++) 38 M=M*p; 39 hash[0]=1; 40 for (i=1;i<=n;i++) 41 hash[i]=hash[i-1]*p+a[i]; 42 for (i=k;i<=n;i++) 43 update(root[i-1],root[i],0,inf,hash[i]-hash[i-k]*M); 44 for (i=1;i<=m;i++) 45 { 46 scanf("%d%d",&l,&r); 47 s=1; 48 for (j=1;j<=k;j++) 49 { 50 scanf("%lld",&x); 51 s=s*p+x; 52 } 53 s=s-M; 54 if (r-l+1<k) printf("Yes\n"); 55 else if (query(root[l+k-2],root[r],0,inf,s)) printf("No\n"); 56 else printf("Yes\n"); 57 } 58 }
