[bzoj3207][可持久化线段树][哈希表]花神的嘲讽计划Ⅰ

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

题解

神题啊我还想着啥ac自动机kmp的
谁知道蹦出来了个我不会的哈希
对于哈希的学习我推荐这个
然后其实只要对这个串搞一搞哈希，因为给出的匹配子串长度稳定是k嘛
对于每个点，存一下以他为结尾的串的哈希值
每次给出的串搞一搞哈希，扔进可持久化线段树里找一下x~y这里面有没有这个哈希值
结果我把可持久化线段树写成了主席树？？？

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
const long long ha=107;
unsigned long long hash[110000];
long long a[110000];
struct Q
{
    unsigned long long c;int p;
}u[110000];
bool cmp(Q n1,Q n2){return n1.c<n2.c;}
int n,m,k;
struct node
{
    int lc,rc,c;
}tr[5110000];int tot;
void add(int &now,int l,int r,int p)
{
    if(now==0)now=++tot;
    tr[now].c++;
    if(l==r)return ;
    int mid=(l+r)/2;
    if(p<=mid)add(tr[now].lc,l,mid,p);
    else add(tr[now].rc,mid+1,r,p);
}
int rt[210000],Rank[210000];
void merge(int &x,int y)
{
    if(x==0){x=y;return ;}
    if(y==0)return ;
    tr[x].c+=tr[y].c;
    merge(tr[x].lc,tr[y].lc);
    merge(tr[x].rc,tr[y].rc);
}
bool fd_kth(int x,int y,int l,int r,int kk)
{
    if(tr[y].c-tr[x].c==0)return false;
    if(l==r)
    {
        if(tr[y].c-tr[x].c>0)return true;
        return false;
    }
    int mid=(l+r)/2,c=mid-l+1;
    if(kk<=c)return fd_kth(tr[x].lc,tr[y].lc,l,mid,kk);
    else return fd_kth(tr[x].rc,tr[y].rc,mid+1,r,kk-c);
}
long long b[110000];
int fd(long long p)
{
    int l=1,r=n;
    while(l<=r)
    {
        int mid=(l+r)/2;
        if(u[mid].c<p)l=mid+1;
        else if(u[mid].c>p)r=mid-1;
        else return mid-k+1;
    }
    return -1;
}
int main()
{
    scanf("%d%d%d",&n,&m,&k);
    for(int i=1;i<=n;i++)scanf("%lld",&a[i]);
    hash[0]=0;
    for(int i=1;i<=n;i++)hash[i]=hash[i-1]*ha+a[i];
    unsigned long long Mod=1;
    for(int i=1;i<=k;i++)Mod*=ha;
    for(int i=k;i<=n;i++)
    {
        u[i].c=hash[i]-hash[i-k]*Mod;
        u[i].p=i;
    }
    sort(u+1,u+1+n,cmp);
    for(int i=k;i<=n;i++)Rank[u[i].p]=(i-k+1);
    for(int i=k;i<=n;i++)
    {
        add(rt[i],1,n,Rank[i]);
        merge(rt[i],rt[i-1]);
    }
    while(m--)
    {
        int x,y;
        scanf("%d%d",&x,&y);
        for(int i=1;i<=k;i++)scanf("%lld",&b[i]);
        unsigned long long tmp=0;
        for(int i=1;i<=k;i++)tmp=tmp*ha+b[i];
        int P=fd(tmp);
        if(P==-1){printf("Yes\n");continue;}
        if(fd_kth(rt[x+k-2],rt[y],1,n,P))printf("No\n");
        else printf("Yes\n");
    }
    return 0;
}