bzoj 3207 花神的嘲讽计划Ⅰ 主席树+hash

花神的嘲讽计划Ⅰ

Time Limit: 10 Sec Memory Limit: 128 MB
Submit: 3112 Solved: 1086
[Submit][Status][Discuss]

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不会尴尬

Source

原创 Memphis

题解：

　　主席树吧，两棵线段树形式相同，满足前缀和的形式是主席树的基本操作

　　这道题k确定，这不就随便hash就可以瞎搞了。

 1 #include<cmath>
 2 #include<cstdio>
 3 #include<cstring>
 4 #include<iostream>
 5 #include<algorithm>
 6 
 7 #define ull unsigned ll
 8 #define ll long long
 9 #define inf 18446744073709551615UL
10 #define N 200007
11 #define MX 10000007
12 using namespace std;
13 int read()
14 {
15     int x=0,f=1;char ch=getchar();
16     while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
17     while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
18     return x*f;
19 }
20 
21 int n,m,K,cnt;
22 int val[N],b[N],rt[N];
23 int sum[MX],a[MX][2];
24 ull H[N];
25 
26 void ins(int x,int &y,unsigned ll val)
27 {
28     y=++cnt;int tmp=y;
29     sum[y]=sum[x]+1;
30     unsigned ll l=0,r=inf,mid;
31     while(l!=r)
32     {
33         mid=l/2+r/2;
34         if((l&1)&&(r&1))mid++;
35         int t=0;
36         if(val>mid)t=1,l=mid+1;
37         else r=mid-1;
38         a[y][t]=++cnt;
39         a[y][t^1]=a[x][t^1];
40         x=a[x][t];y=a[y][t];sum[y]=sum[x]+1;
41     }
42     y=tmp;
43 }
44 bool query(int x,int y,unsigned ll val)
45 {
46     unsigned ll l=0,r=inf,mid;
47     while(l!=r)
48     {
49         mid=l/2+r/2;
50         if((l&1)&&(r&1))mid++;
51         int t=0;
52         if(val>mid) t=1,l=mid+1;
53         else r=mid-1;
54         x=a[x][t],y=a[y][t];
55     }
56     return (sum[y]-sum[x])!=0;
57 }
58 int main()
59 {
60     n=read(),m=read(),K=read();
61     for(int i=1;i<=n;i++)
62         val[i]=read();
63     for(int i=1;i<=n;i++)
64         H[i]=H[i-1]*107+val[i];
65     unsigned ll M=1;
66     for(int i=1;i<=K;i++)M*=107;
67         for(int i=K;i<=n;i++)
68             ins(rt[i-1],rt[i],H[i]-H[i-K]*M);
69     while(m--)
70     {
71         int x=read(),y=read();
72         for(int i=1;i<=K;i++)
73             b[i]=read();
74         unsigned ll now=0;
75         for(int i=1;i<=K;i++)
76             now=now*107+b[i];
77         if(query(rt[x+K-2],rt[y],now))puts("No");
78         else puts("Yes");
79     }
80 }

o(￣ヘ￣o＃)，hzw的代码。