BZOJ 2806 Cheat (后缀自动机+dp)

题意：

给出m个串组成的文章，在给出n个串，求n个串每个串的一个最优划分，划分：将串分成多分子串，L0表示其中能在文章中找到的子串的长度和，maxl表示其中子串的最大长度，在保证L0>=90%总长度前提下使得maxl最大，求出这个最大的maxl。

题解：

我们分析下这样的题，很明显要满足两种条件不易，那么我们可以枚举maxl，然后判断是否符合条件，对于如何判断我们可以用dp处理，dp[i]表示到i位置时出现在文章中的子串长度和，方程:dp[i]=max{dp[j]+i-j},那么我们要做的就是判断i-j是否出现在文章中！我们通过后缀自动机预处理每个位置i能匹配的最大长度，如果v[i]>i-j，那么明显这段是会出现在文章中的。由于复杂度)(n^2)所以考虑用优化（这个没想到），一直就对队列优化很盲，这题返程可以变形dp[i]-i=dp[j]-j，很明显存在单调递增性，然后就可以用队列优化，队列优化都是很简单，我就是学不会0_0.

#include<iostream>
#include<math.h>
#include<stdio.h>
#include<algorithm>
#include<string.h>
#include<vector>
#include<queue>
#include<map>
#include<set>
using namespace std;
#define B(x) (1<<(x))
void cmax(int& a,int b){ if(b>a)a=b; }
void cmin(int& a,int b){ if(b<a)a=b; }
typedef long long ll;
const int oo=0x3f3f3f3f;
const ll OO=1LL<<61;
const int MOD=2012;
const int maxn=2500005;
const int maxm=200005;
const int SIZE=maxn<<1;
const int type=3;
int Next[SIZE][type],fa[SIZE],len[SIZE];
int tol,last;
int dp[maxn],v[maxn];
int q[maxn],h,t;
char str[maxn];

int newNode(int x){

    len[tol]=x;
    fa[tol]=-1;
    for(int i=0;i<type;i++){
        Next[tol][i]=-1;
    }
    return tol++;
}

void Init(){
    tol=0;
    last=newNode(0);
}

void add(int k){

    int now=last;
    int end=newNode(len[now]+1);
    while(now!=-1&&Next[now][k]==-1){
        Next[now][k]=end;
        now=fa[now];
    }
    if(now==-1) fa[end]=0;
    else{
        int nxt=Next[now][k];
        if(len[nxt]==len[now]+1) fa[end]=nxt;
        else{
            int cnxt=newNode(len[now]+1);
            for(int i=0;i<type;i++) Next[cnxt][i]=Next[nxt][i];
            fa[cnxt]=fa[nxt];
            fa[nxt]=fa[end]=cnxt;
            while(now!=-1&&Next[now][k]==nxt){
                Next[now][k]=cnxt;
                now=fa[now];
            }
        }
    }
    last=end;
}

void Insert(char T[]){
    for(int i=0;T[i];i++)
        add(T[i]-'0');
}

void Search(char T[]){

    int now=0,l=0;
    for(int i=0;T[i];i++){

        int k=T[i]-'0';
        if(Next[now][k]!=-1){
            now=Next[now][k];
            l++;
        }else{
            while(now!=-1&&Next[now][k]==-1)now=fa[now];
            if(now==-1) l=now=0;
            else{
                l=len[now]+1;
                now=Next[now][k];
            }
        }
        v[i+1]=l;
    }
}

bool DP(int ml,int n){

    h=t=0;
    dp[0]=0;
    for(int i=1;i<=n;i++){

        dp[i]=dp[i-1];
        int L=i-v[i],R=i-ml;
        if(R>=0){
            while(h<t&&dp[R]-R>dp[q[h]]-q[h])t--;
            q[t++]=R;
        }
        while(h<t&&q[h]<L)h++;
        if(h<t)cmax(dp[i],dp[q[h]]+(i-q[h]));
    }
    return 10*dp[n]>=9*n;
}

int solve(int n){

    int l=0,r=n,mid;
    int ans=0;
    while(l<=r){

        mid=(l+r)>>1;
        if(DP(mid,n)){
            ans=mid;
            l=mid+1;
        }else r=mid-1;
    }
    return ans;
}

int main(){

    int N,M;
    while(scanf("%d %d",&N,&M)!=EOF){

        Init();
        for(int i=1;i<=M;i++){
            scanf("%s",str);
            Insert(str);
            add(2);
        }
        for(int i=1;i<=N;i++){
            scanf("%s",str);
            Search(str);
            printf("%d\n",solve(strlen(str)));
        }
    }
    return 0;
}
/**
1 2
10110
000001110
1011001100
*/