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

看到题目就可以想到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=0,f=1;char ch=getchar();
  while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
  while(ch>='0'&&ch<='9'){x=10*x+ch-'0';ch=getchar();}
  return x*f;
}
int rt[N],n,m,k,tot,ls[M*40],rs[M*40],sz[M*40];
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]+1;
  if(l==r)return;
  uint mid=l/2+r/2;if((l&1)&&(r&1))mid++; 
  if(pos<=mid)add(ls[pre],ls[x],l,mid,pos);
  else add(rs[pre],rs[x],mid+1,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/2+r/2;if((l&1)&&(r&1))mid++;
  if(pos<=mid)query(ls[L],ls[R],l,mid,pos);
  else query(rs[L],rs[R],mid+1,r,pos);
}
int main(){
  n=read();m=read();k=read();
  hash[0]=0;
  for(int i=1;i<=n;i++)hash[i]=hash[i-1]*base+(uint)read();
  bas=1;
  for(int i=1;i<=k;i++)bas*=base;
  for(int i=1;i<=n-k+1;i++)add(rt[i-1],rt[i],1,inf,hash[i+k-1]-hash[i-1]*bas);
  while(m--){
    int x=read(),y=read();
    uint ha=0;
    for(int i=1;i<=k;i++)ha=ha*base+(uint)read();
    query(rt[x-1],rt[y-k+1],1,inf,ha)?puts("No"):puts("Yes"); 
  }
  return 0;
}

莫队算法：

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

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