[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
;
int next
[3],fail
,stp
,v
,Q
,f
;

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;
}