BZOJ 3207 花神的嘲讽计划Ⅰ

title: ‘BZOJ 3207 花神的嘲讽计划Ⅰ’
categories: BZOJ
date: 2016-1-6 13:40:00
tags: [主席树,Hash]

---

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.txt
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

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

Solution

对每段长度为k的序列哈希一下，存到0～Hash_Max的权值线段树里面，加个可持久化。。。
之后就是随便搞搞了。。哈希我开的ULL自然溢出，这样(l+r)>>1那里可能会爆掉，需要换成(l>>1)+(r>>1)+(l&r&1)。然而，然而某神犇wyy离散化了哈希值，d了一波我的空间复杂度

Code

#include<cstdio>
#include<climits>
#include<cstdlib>
#include<cstring>
#include<algorithm>

#define maxn 2000000+5
#define maxm 8000000+5

using namespace std;

typedef unsigned long long ll;

struct Ha_Tree{
    int l,r;
    int sum;
}tr[maxm];

ll haha[maxn],ha[maxn];
int a[maxn],root[maxn];
int ind,n,m,k;
ll utk=1;

inline int in(){
    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=x*10+ch-'0',ch=getchar();
    return x*f;
}

void Update(ll l,ll r,int last,int &k,ll pos,int val){
    k=++ind;
    tr[k]=tr[last]; tr[k].sum+=val;
    if(l==r) return;
    ll mid=(l>>1)+(r>>1)+(l&r&1);
    if(pos<=mid) Update(l,mid,tr[last].l,tr[k].l,pos,val);
    else Update(mid+1,r,tr[last].r,tr[k].r,pos,val);
}

int Query(ll l,ll r,int pl,int pr,ll pos){
    if(l==r) return tr[pr].sum-tr[pl].sum;
    ll mid=(l>>1)+(r>>1)+(l&r&1);
    if(pos<=mid) return Query(l,mid,tr[pl].l,tr[pr].l,pos);
    else return Query(mid+1,r,tr[pl].r,tr[pr].r,pos);
}

int main(){
    freopen("3207.in","r",stdin);
    scanf("%d%d%d",&n,&m,&k);
    for(int i=1;i<=n;i++)
    a[i]=in();
    for(int i=1;i<=n;i++){
    haha[i]=a[i]+haha[i-1]*23333;
    if(i<=k) utk=utk*23333;
    }
    for(int i=k;i<=n;i++)
    ha[i]=haha[i]-haha[i-k]*utk;
    for(int i=k;i<=n;i++)
    Update(0,ULLONG_MAX,root[i-1],root[i],ha[i],1);
    for(int i=1;i<=m;i++){
    int l,r; ll moha=0;
    l=in(); r=in();
    for(int j=1;j<=k;j++) moha=moha*23333+in();
    if(Query(0,ULLONG_MAX,root[l+k-2],root[r],moha)) puts("No");
    else puts("Yes");
    }
    return 0;
}