bzoj 2806: [Ctsc2012]Cheat

传送门

　　好久没刷bzoj惹……

　　题意不说可以嘛。

　　首先二分答案。

　　SAM的事情搞完以后就是dp辣。

　　我们已经对于每个位置i，找到了最小的一个k，使得[k,i]这个子串在模版串中出现过。那么我们需要做的是把f[i]给min上f[k]到f[i-x]，直接搞是$n^2logn$的，套个数据结构也是两个log的。然而如果一个位置j不在合法的区间中，那么以后也不会进入，那么直接用一个单调队列维护就好了。

#include<cstdio>
#include<cstring>
#include<algorithm>
#define MN 1200000
using namespace std;
int read_p,read_ca,read_f;
inline int read(){
    read_p=;read_ca=getchar();read_f=;
    while(read_ca<''||read_ca>'') read_f=read_ca=='-'?-:read_f,read_ca=getchar();
    while(read_ca>=''&&read_ca<='') read_p=read_p*+read_ca-,read_ca=getchar();
    return read_p*read_f;
}
struct na{int l,f,t[];na(){t[]=t[]=;f=-;}}t[MN<<];
int n,m,la,num=,mi[MN],st[MN];
char s[MN];
inline void add(int x){
    int p=++num;t[p].l=t[la].l+;
    while (la!=-&&!t[la].t[x]) t[la].t[x]=p,la=t[la].f;
    if (la==-) t[p].f=;else{
        int o=t[la].t[x];
        if (t[o].l==t[la].l+) t[p].f=o;else{
            int np=++num;
            t[np]=t[o];
            t[np].l=t[la].l+;
            t[o].f=t[p].f=np;
            while (la!=-&&t[la].t[x]==o) t[la].t[x]=np,la=t[la].f;
        }
    }
    la=p;
}
inline bool ju(int x){
    mi[]=;
    int i,p,l,L=,R=;
    for (i=;s[i-];i++) mi[i]=1e9;
    for (i=,p=,l=;s[i];i++){
        while (p&&!t[p].t[s[i]-'']) p=t[p].f,l=t[p].l;
        if (t[p].t[s[i]-'']) p=t[p].t[s[i]-''],l++;
        if (i+>=x){
            while(L<=R&&mi[i+-x]<=mi[st[R]]) R--;
            st[++R]=i+-x;
        }
        while (L<=R&&st[L]<=i-l) L++;
        if (L<=R) if (mi[i+]>mi[st[L]]) mi[i+]=mi[st[L]];
        if (mi[i+]>mi[i]+) mi[i+]=mi[i]+;
    }
    return mi[i]*<=i;
}
int main(){
    n=read();m=read();
    for (int i=;i<=m;i++){
        la=;scanf("%s",s);
        for (int j=;s[j];j++) add(s[j]-'');
    }
    for (int i=;i<=n;i++){
        scanf("%s",s);
        int l=,r=strlen(s),mid;
        while(l<r) if (ju(mid=l+r+>>)) l=mid;else r=mid-;
        printf("%d\n",l);
    }
}