bzoj 3207 花神的嘲讽计划Ⅰ

Description

背景

花神是神，一大癖好就是嘲讽大J，举例如下：

“哎你傻不傻的！【hqz：大笨J】”

“这道题又被J屎过了！！”

“J这程序怎么跑这么快！J要逆袭了！”

……

描述

这一天DJ在给吾等众蒟蒻讲题，花神在一边做题无聊，就跑到了一边跟吾等众蒟蒻一起听。以下是部分摘录：

“J你在讲什么！”

“我在讲XXX！”

“哎你傻不傻的！这么麻烦，直接XXX再XXX就好了！”

“……”

“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

Sample Output

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
");
55       else if (query(root[l+k-2],root[r],0,inf,s)) printf("No
");
56       else printf("Yes
");
57     }
58 }