[BZOJ2806][Ctsc2012][后缀自动机][队列优化][DP]Cheat

题意略……

因为要让L的长度最大，可以二分L，然后check
要求在L的限制下最大的匹配长度，可以DP。
令f[i]表示前i个字符的最大匹配长度，g[i]表示1~i可以匹配的最长后缀的长度，则f[i]=max{f[j]+i-j | i-g[i]<=j<=1-L}
i-g[i]是单调不降的，所以这个DP可以用单调队列优化

#include <cstdio>
#include <iostream>
#include <string>
#include <cstring>
#include <algorithm>
#define N 400010

using namespace std;

int n,m,head,tail,p=1,cnt=1,len;
char A[N];
int next[N][3],fail[N],stp[N],v[N],Q[N],f[N];

inline void extend(int x){
  x-='0';
  int np=++cnt; stp[np]=stp[p]+1;
  while(p&&!next[p][x]) next[p][x]=np,p=fail[p];
  if(!p) fail[np]=1;
  else{
    int q=next[p][x];
    if(stp[q]==stp[p]+1) fail[np]=q;
    else{
      int nq=++cnt; stp[nq]=stp[p]+1;
      memcpy(next[nq],next[q],sizeof(next[q]));
      fail[nq]=fail[q];
      fail[q]=fail[np]=nq;
      while(p&&next[p][x]==q) next[p][x]=nq,p=fail[p];
    }
  }
  p=np;
}

inline void Getval(int len){
  int x=0,k=1;
  for(int i=1;i<=len;i++){
    if(next[k][A[i]-'0']) x++,k=next[k][A[i]-'0'];
    else{
      while(k&&!next[k][A[i]-'0']) k=fail[k];
      if(!k) k=1,x=0; else x=stp[k]+1,k=next[k][A[i]-'0'];
    }
    v[i]=x;
  }
}

inline bool check(int L){
  head=1; tail=0; f[0]=0;
  for(int i=1;i<=len;i++){
    f[i]=f[i-1];
    if(i-L>=0){
      while(head<=tail&&f[i-L]-i+L>f[Q[tail]]-Q[tail]) tail--;
      Q[++tail]=i-L;
    }
    while(head<=tail&&i-v[i]>Q[head]) head++;
    if(head<=tail)f[i]=max(f[i],f[Q[head]]-Q[head]+i);
  }
  return 9*len<=10*f[len];
}

int main(){
  freopen("1.in","r",stdin);
  freopen("1.out","w",stdout);
  scanf("%d%d",&n,&m);
  for(int i=1;i<=m;i++){
    scanf("%s",A+1);
    int len=strlen(A+1); A[++len]='2';
    for(int i=1;i<=len;i++) extend(A[i]);
  }
  for(int i=1;i<=n;i++){
    scanf("%s",A+1);
    len=strlen(A+1);
    int l=1,r=len,mid,Ans;
    Getval(len);
    while(l<=r) check(mid=l+r>>1)?(Ans=mid,l=mid+1):r=mid-1;
    printf("%d\n",Ans);
  }
  return 0;
}