[bzoj3207]: 花神的嘲讽计划Ⅰ

来自FallDream的博客，未经允许，请勿转载，谢谢，

---

题意太...我来简化一下。

给一个n个数的序列，m次询问，每次询问[l,r]这个区间内是否有一个给定的长度为k(固定)的子串si

n,m<=10^5 k<=20 

k是固定的，考虑直接哈希+主席树，每个点求出以它为开头的子串的哈希值插入主席树。

k不大 所以哈希值可以暴力算233

数据有点鬼畜 还是上两个膜数比较稳

复杂度k(n+m)+mlogn

#include<cstdio>
#include<iostream>
#include<map>
#define MN 100000
#define mod 999911659
#define mod2 9875321
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 = x * 10 + ch - '0';ch = getchar();}
    return x * f;
}

struct Tree{int l,r,x;}T[MN*50];
int n,m,k,a[MN+5],cc=0,cc2=0,cnt=0,H[MN+5],H2[MN+5],rt[MN+5],rt2[MN+5];
map<int,int> mp,mp2;

int query(int x,int nx,int v,int r)
{
    int l=1,mid;
    while(l<r)
    {
        mid=l+r>>1;
        if(v<=mid) r=mid,x=T[x].l,nx=T[nx].l;
        else l=mid+1,x=T[x].r,nx=T[nx].r;
    }
    return T[nx].x-T[x].x;
}

void ins(int x,int nx,int v,int r)
{
    int l=1,mid;T[nx].x=T[x].x+1;
    while(l<r)
    {
        mid=l+r>>1;
        if(v<=mid) T[nx].r=T[x].r,T[nx].l=++cnt,nx=T[nx].l,x=T[x].l,r=mid;
        else T[nx].l=T[x].l,T[nx].r=++cnt,nx=T[nx].r,x=T[x].r,l=mid+1;
        T[nx].x=T[x].x+1;
    }
}

int main()
{
    n=read();m=read();k=read();
    for(int i=1;i<=n;++i) a[i]=read();
    for(int i=1;i+k-1<=n;++i)
    {
        int ha=0,ha2=0;
        for(int j=0;j<k;++j)
            ha=(1LL*ha*193+a[i+j])%mod,ha2=(1LL*ha2*73+a[i+j])%mod2;
        !mp[ha]?mp[ha]=++cc:0;H[i]=mp[ha];
        !mp2[ha2]?mp2[ha2]=++cc2:0;H2[i]=mp2[ha2];
    }
    for(int i=1;i+k-1<=n;++i) ins(rt[i-1],rt[i]=++cnt,H[i],cc),ins(rt2[i-1],rt2[i]=++cnt,H2[i],cc2);
    for(int i=1;i<=m;++i)
    {
        int l=read(),r=read(),ha=0,ha2=0;
        for(int i=1,j;i<=k;++i)
            ha=(1LL*ha*193+(j=read()))%mod,ha2=(1LL*ha2*73+j)%mod2;
        ha=mp[ha];ha2=mp2[ha2];
        if(!ha||!ha2||(!query(rt[l-1],rt[r-k+1],ha,cc)&&!query(rt2[l-1],rt2[r-k+1],ha2,cc2))) puts("Yes");
        else puts("No");
    }
    return 0;
}