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

2/22/2017来源:ASP.NET技巧人气:4072

题意略……

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