【BZOJ3207】花神的嘲讽计划I 可持久化线段树/莫队

时间:2023-12-26 13:58:37

看到题目就可以想到hash

然后很自然的联想到可持久化权值线段树

WA:base取了偶数

这道题还可以用莫队做,比线段树快一些

可持久化线段树:

 #include<bits/stdc++.h>
#define ll long long
#define uint unsigned int
#define ull unsigned long long
#define inf 4294967295
#define N 100005
#define M 100005
#define base 100003
using namespace std;
inline int read(){
int x=,f=;char ch=getchar();
while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
while(ch>=''&&ch<=''){x=*x+ch-'';ch=getchar();}
return x*f;
}
int rt[N],n,m,k,tot,ls[M*],rs[M*],sz[M*];
uint hash[N],bas;
void add(int pre,int &x,uint l,uint r,uint pos){
x=++tot;ls[x]=ls[pre],rs[x]=rs[pre];sz[x]=sz[pre]+;
if(l==r)return;
uint mid=l/+r/;if((l&)&&(r&))mid++;
if(pos<=mid)add(ls[pre],ls[x],l,mid,pos);
else add(rs[pre],rs[x],mid+,r,pos);
}
int query(int L,int R,uint l,uint r,uint pos){
if(l==r)return sz[R]-sz[L];
uint mid=l/+r/;if((l&)&&(r&))mid++;
if(pos<=mid)query(ls[L],ls[R],l,mid,pos);
else query(rs[L],rs[R],mid+,r,pos);
}
int main(){
n=read();m=read();k=read();
hash[]=;
for(int i=;i<=n;i++)hash[i]=hash[i-]*base+(uint)read();
bas=;
for(int i=;i<=k;i++)bas*=base;
for(int i=;i<=n-k+;i++)add(rt[i-],rt[i],,inf,hash[i+k-]-hash[i-]*bas);
while(m--){
int x=read(),y=read();
uint ha=;
for(int i=;i<=k;i++)ha=ha*base+(uint)read();
query(rt[x-],rt[y-k+],,inf,ha)?puts("No"):puts("Yes");
}
return ;
}

莫队算法:

 #include<bits/stdc++.h>
#define ll long long
#define N 100005
#define base 100003
using namespace std;
inline int read(){
int x=,f=;char ch=getchar();
while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
while(ch>=''&&ch<=''){x=*x+ch-'';ch=getchar();}
return x*f;
}
int n,m,k,a[N],tot,pos[N],cnt[N];
bool ans[N];
ll hash[N],h[N],bas;
struct Node{
int l,r,v,id;
}q[N];
bool cmp(Node a,Node b){
if(pos[a.l]==pos[b.l])return a.r<b.r;
else return a.l<b.l;
}
bool find(int x){
if(x==-)return ;
return cnt[x]?:;
}
int main(){
n=read();m=read();k=read();
for(int i=;i<=n;i++)hash[i]=hash[i-]*base+(ll)read();
bas=;
for(int i=;i<=k;i++)bas*=base;
for(int i=;i<=n-k+;i++)h[++tot]=hash[i+k-]-hash[i-]*bas;
sort(h+,h++tot);
for(int i=;i<=n-k+;i++)a[i]=lower_bound(h+,h++tot,hash[i+k-]-hash[i-]*bas)-h;
int block=sqrt(n);
for(int i=;i<=n;i++)pos[i]=(i-)/block+;
for(int i=;i<=m;i++){
int x=read(),y=read();y=y-k+;
ll tmp=;
for(int j=;j<=k;j++)tmp=tmp*base+(ll)read();
int t=lower_bound(h+,h++tot,tmp)-h;
if(t>n||h[t]!=tmp)t=-;
q[i].l=x,q[i].r=y;q[i].id=i;q[i].v=t;
}
sort(q+,q++m,cmp);
int ql=,qr=;
for(int i=;i<=m;i++){
while(qr<q[i].r)cnt[a[++qr]]++;
while(ql>q[i].l)cnt[a[--ql]]++;
while(qr>q[i].r)cnt[a[qr--]]--;
while(ql<q[i].l)cnt[a[ql++]]--;
ans[q[i].id]=find(q[i].v);
}
for(int i=;i<=m;i++)puts(ans[i]?"No":"Yes");
return ;
}

3207: 花神的嘲讽计划Ⅰ

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 1627  Solved: 587
[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

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