BZOJ_3207_花神的嘲讽计划Ⅰ

BZOJ_3207_花神的嘲讽计划Ⅰ_哈希+主席树

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

可以对每个询问的序列哈希，转化为询问区间里是否有x，主席树维护即可。

注意如果用ull自然溢出的哈希的话，可能l+r爆ull，所以要先除后加。

代码:

#include <stdio.h>
#include <string.h>
#include <algorithm>
using namespace std;
inline char nc() {
    static char buf[100000],*p1,*p2;
    return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline int rd() {
    register int x=0;
    register char s=nc();
    while(s<'0'||s>'9') s=nc();
    while(s>='0'&&s<='9') x=(x<<3)+(x<<1)+s-'0',s=nc();
    return x;
}
typedef unsigned long long ll;
#define N 200001
const ll maxn=-1;
ll mi[N],h[N],base=131;
int root[N],n,K,cnt,siz[N*50],ls[N*50],rs[N*50];
ll get_h(int l,int r) {
    return h[r]-h[l-1]*mi[r-l+1];
}
void insert(int x,int &y,ll l,ll r,ll v) {
    y=++cnt; 
    if(l==r) {
        siz[y]=siz[x]+1; return ;
    }
    ll mid=(l/2+r/2+(l&r&1));
    if(v<=mid) rs[y]=rs[x],insert(ls[x],ls[y],l,mid,v);
    else ls[y]=ls[x],insert(rs[x],rs[y],mid+1,r,v);
    siz[y]=siz[ls[y]]+siz[rs[y]];
}
bool find(int x,int y,ll l,ll r,ll v) {
    if(siz[y]-siz[x]<=0) return 0;
    if(l==r) return siz[y]-siz[x];
    ll mid=(l/2+r/2+(l&r&1));
    if(v<=mid) return find(ls[x],ls[y],l,mid,v);
    else return find(rs[x],rs[y],mid+1,r,v);
}
int main() {
    int m;
    n=rd(); m=rd(); K=rd();
    int i,x,y,z;
    for(mi[0]=1,i=1;i<=n;i++) {
        x=rd();
        mi[i]=mi[i-1]*base;
        h[i]=h[i-1]*base+x;
    }
    for(i=1;i<=n-K+1;i++) {
        insert(root[i-1],root[i],0,maxn,get_h(i,i+K-1));
    }
    while(m--) {
        ll has=0;
        x=rd(); y=rd()-K+1;
        for(i=1;i<=K;i++) {
            z=rd();
            has=has*base+z;
        }
        puts(find(root[x-1],root[y],0,maxn,has)?"No":"Yes");
    }
}