BZOJ3207 花神的嘲讽计划I

Time Limit: 10 Sec Memory Limit: 128 MB

Summary

给你一个模式串P,q个询问,对每个询问回答从Pl到Pr是否存在与给定串相同的子串,同时有所有的给定串长度均为k(k≤20) n,m≤200000

Description

背景

花神是神,一大癖好就是嘲讽大J,举例如下:

“哎你傻不傻的!【hqz:大笨J】”

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

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

……

描述

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

“J你在讲什么!”

“我在讲XXX!”

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

“……”

“J你XXX讲过了没?”

“……”

“那个都不讲你就讲这个了?哎你傻不傻的!”

“……”

DJ对这种情景表示非常无语,每每出现这种情况,DJ都是非常尴尬的。

经过众蒟蒻研究,DJ在讲课之前会有一个长度为N方案,我们可以把它看作一个数列;

同样,花神在听课之前也会有一个嘲讽方案,有M个,每次会在x到y的这段时间开始嘲讽,为了减少题目难度,每次嘲讽方案的长度是一定的,为K。

花神嘲讽DJ让DJ尴尬需要的条件:

在xy的时间内DJ没有讲到花神的嘲讽方案,即J的讲课方案中的xy没有花神的嘲讽方案【这样花神会嘲讽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不会尴尬

//2018.7.6新加数据一组.未重测.

Solution

对原串做一遍字符串哈希,取出每一个K长度的串的哈希值,建主席树。查询的时候,先判断区间有没有K大,然后判断是给定串的哈希值是否出现过,然后再在对应区间主席树里面查找是否存在即可。主席树入门题吧

#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#define lowbit(x) ((x)&(-(x)))
using namespace std;
#define REP(i,a,n) for(register int i=(a);i<=(n);++i)
#define PER(i,a,n) for(register int i=(a);i>=(n);--i)
#define FEC(i,x) for(register int i=head[x];i;i=g[i].ne)
template<typename A>inline void read(A&a){a=0;char c=0;int f=1;while(c<'0'||c>'9')(c=getchar())=='-'?f=-1:0;while(c>='0'&&c<='9')a=(a<<3)+(a<<1)+c-'0',c=getchar();f==-1?a=-a:0;}
char buf[30];template<typename A>inline void write(A a){if(a<0)putchar('-'),a=-a;int top=0;if(!a)buf[top=1]='0';while(a)buf[++top]=a%10+'0',a/=10;while(top)putchar(buf[top--]);}
typedef long long ll;typedef unsigned long long ull;
template<typename A,typename B>inline bool SMAX(A&x,const B&y){return y>x?x=y,1:0;}
template<typename A,typename B>inline bool SMIN(A&x,const B&y){return y<x?x=y,1:0;}

const int N=2e5+7,base=1997;
int n,m,K,a[N],dis,x,y,z;
ull h[N],bin[N],f[N],b[N],s;

inline void Make_Hash(){
	bin[0]=1;for(register int i=1;i<=n;++i)bin[i]=bin[i-1]*base;
	for(register int i=1;i<=n;++i)h[i]=h[i-1]*base+a[i];
}
inline ull GetHash(int l,int r){return h[r]-h[l-1]*bin[r-l+1];}

struct Node{int lc,rc,val;}t[N*20];int nod,RT[N];
inline void Insert(int &o,int p,int L,int R,int x,int k){t[o=++nod]=t[p];t[o].val+=k;if(L==R)return;int M=(L+R)>>1;x<=M?Insert(t[o].lc,t[p].lc,L,M,x,k):Insert(t[o].rc,t[p].rc,M+1,R,x,k);}
inline int Sum(int o,int p,int L,int R,int x){if(L==R)return t[o].val-t[p].val;int M=(L+R)>>1;return x<=M?Sum(t[o].lc,t[p].lc,L,M,x):Sum(t[o].rc,t[p].rc,M+1,R,x);}

inline int Find(ull x){int l=1,r=dis;while(l<r){int mid=(l+r)>>1;if(x<=b[mid])r=mid;else l=mid+1;}return b[l]==x?l:0;}
inline void LSH(){
	sort(b+1,b+n+1);for(register int i=1;i<=n;++i)if(i==1||b[i]!=b[i-1])b[++dis]=b[i];
	for(register int i=1;i<=n;++i)f[i]=Find(f[i]);
}
inline void Preprocess(){
	Make_Hash();for(register int i=1;i<=n-K+1;++i)f[i]=b[i]=GetHash(i,i+K-1);
	LSH();for(register int i=1;i<=n-K+1;++i)Insert(RT[i],RT[i-1],1,dis,f[i],1);
}

int main(){
	read(n),read(m);read(K);
	for(register int i=1;i<=n;++i)read(a[i]);
	Make_Hash();Preprocess();
	for(register int i=1;i<=m;++i){
		read(x),read(y);s=0;for(register int i=1;i<=K;++i)read(z),s=s*base+z;
		s=Find(s);if(y-x+1<K||!s||!Sum(RT[y-K+1],RT[x-1],1,dis,s))puts("Yes");else puts("No");
	}
}
原文地址:https://www.cnblogs.com/hankeke/p/BZOJ3207.html