mat eight def span tro con 常數 lin 條件

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

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

bzoj 3207 花神的嘲諷計劃Ⅰ