【BZOJ】【3207】花神的嘲讽计划 I
字符串Hash+可持久化线段树
好神奇的转化……
蒟蒻一开始还去想AC自动机去了……然而由于a[i]的范围是小于等于n，怎么也想不出一个时间复杂度合理的方法
膜拜了题解0。0原来是字符串Hash！
首先每个询问的长度都是固定的！K！那么我们就可以看作一个长度为K的数！然后就可以利用字符串Hash将一个长度为n的数列转化成一个长度为n-k+1的数列！查询一个特定的序列在一个区间中是否出现过，就变成了在一段区间中，某个数是否出现过！然后就可以用可持久化线段树来维护了……
长度为K这个条件好神啊……居然可以这样用！
（裸的可持久化线段树好像写起来顺手了一些……
1 /************************************************************** 2 Problem: 3207 3 User: Tunix 4 Language: C++ 5 Result: Accepted 6 Time:3164 ms 7 Memory:79404 kb 8 ****************************************************************/ 9 10 //BZOJ 3207 11 #include<vector> 12 #include<cstdio> 13 #include<cstring> 14 #include<cstdlib> 15 #include<iostream> 16 #include<algorithm> 17 #define rep(i,n) for(int i=0;i<n;++i) 18 #define F(i,j,n) for(int i=j;i<=n;++i) 19 #define D(i,j,n) for(int i=j;i>=n;--i) 20 #define pb push_back 21 using namespace std; 22 inline int getint(){ 23 int v=0,sign=1; char ch=getchar(); 24 while(ch<'0'||ch>'9'){ if (ch=='-') sign=-1; ch=getchar();} 25 while(ch>='0'&&ch<='9'){ v=v*10+ch-'0'; ch=getchar();} 26 return v*sign; 27 } 28 const int N=2e5+10,INF=~0u>>2; 29 typedef long long LL; 30 typedef unsigned long long u64; 31 /******************tamplate*********************/ 32 33 int n,m,k,a[N],x[N],y[N]; 34 u64 s[N],ask[N],c[N]; 35 int rt[N],cnt,tot; 36 struct node{ 37 int cnt,l,r; 38 }t[N*30]; 39 #define mid (l+r>>1) 40 void update(int &o,int l,int r,int pos){ 41 t[++tot]=t[o], o=tot, t[o].cnt++; 42 if (l==r) return; 43 if (pos<=mid) update(t[o].l,l,mid,pos); 44 else update(t[o].r,mid+1,r,pos); 45 } 46 int query(int i,int j,int val){ 47 i=rt[i],j=rt[j]; 48 int l=1,r=cnt; 49 while(l<r){ 50 if (t[j].cnt-t[i].cnt==0) return 0; 51 if (val<=mid) r=mid,i=t[i].l,j=t[j].l; 52 else l=mid+1,i=t[i].r,j=t[j].r; 53 } 54 return t[j].cnt-t[i].cnt; 55 } 56 int main(){ 57 #ifndef ONLINE_JUDGE 58 freopen("3207.in","r",stdin); 59 freopen("3207.out","w",stdout); 60 #endif 61 n=getint(); m=getint(); k=getint(); 62 F(i,1,n) a[i]=getint(); 63 n=n-k+1; 64 F(i,1,n){ 65 F(j,i,i+k-1) s[i]=s[i]*131+a[j]; 66 c[++cnt]=s[i]; 67 } 68 F(i,1,m){ 69 x[i]=getint(); y[i]=getint()-k+1; 70 F(j,1,k) ask[i]=ask[i]*131+getint(); 71 c[++cnt]=ask[i]; 72 } 73 sort(c+1,c+cnt+1); 74 cnt=unique(c+1,c+cnt+1)-c-1; 75 F(i,1,n) s[i]=lower_bound(c+1,c+cnt+1,s[i])-c; 76 F(i,1,m) ask[i]=lower_bound(c+1,c+cnt+1,ask[i])-c; 77 F(i,1,n){ 78 rt[i]=rt[i-1]; 79 update(rt[i],1,cnt,s[i]); 80 } 81 F(i,1,m){ 82 if (query(x[i]-1,y[i],ask[i])) puts("No"); 83 else puts("Yes"); 84 } 85 return 0; 86 }
3207: 花神的嘲讽计划Ⅰ
Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 816 Solved: 311
[Submit][Status][Discuss]
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不会尴尬
Source
