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

来源：互联网发布：车机软件编辑：程序博客网时间：2024/05/01 02:07

题意略……

因为要让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 400010using 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;}

0 0