bzoj 2806: [Ctsc2012]Cheat （后缀自动机+dp+单调队列）

2806: [Ctsc2012]Cheat

Time Limit: 20 Sec   Memory Limit: 256 MB
Submit: 1056   Solved: 544
[ Submit][ Status][ Discuss]

Description

Input

第一行两个整数N，M表示待检查的作文数量，和小强的标准作文库
的行数
接下来M行的01串，表示标准作文库
接下来N行的01串，表示N篇作文

Output

N行，每行一个整数，表示这篇作文的Lo 值。

Sample Input

1 2
10110
000001110
1011001100

Sample Output

4

HINT

输入文件不超过1100000字节

注意：题目有改动,可识别的长度不小于90%即可，而不是大于90%

Source

[ Submit][ Status][ Discuss]

题解：后缀自动机+dp+单调队列

对标准作文库建立广义的后缀自动机，然后处理作文的时候先让其在后缀自动机上进行匹配。求出对于作文中的每个位置最长匹配的长度ans[i]（就是从这个点向前最多多长是标准作文库的子串）

然后我们可以二分L，然后判定。

f[i]表示的到i最长的匹配长度 f[i]=max(f[i-1],f[j]+i-j) i-ans[i]<=j<=i-L

这个dp是O(n^2)的 ，但是我们必须要降为O(n)

ans[i]+1>=ans[i+1]

i+1-ans[i]-1<=i+1-ans[i+1]

i-ans[i]<=i+1-ans[i+1]

所以i-ans[i]是单调不降的，我们可以用单调队列维护区间f[i]-i的最值。

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
#define N 1100003
#define eps 1e-7
using namespace std;
int l[N],ch[N][3],fa[N],ans[N],st[N],top,f[N];
int n,m,cnt,root,last,size,p,q,np,nq,pos[N],cn;
char s[N],s1[N];
struct data{
	int l,r,len;
}a[N];
int cmp(data a,data b)
{
	return a.l<b.l||a.l==b.l&&a.r<b.r;
}
void extend(int i)
{
    int x=s1[i]-'0';
	p=last; last=np=++size;
	l[np]=l[p]+1;
	for (;!ch[p][x]&&p;p=fa[p]) ch[p][x]=np;
	if (!p) fa[np]=root;
	else {
		q=ch[p][x];
		if (l[p]+1==l[q]) fa[np]=q;
		else {
			nq=++size; l[nq]=l[p]+1;
			memcpy(ch[nq],ch[q],sizeof(ch[nq]));
			fa[nq]=fa[q];
			fa[q]=fa[np]=nq;
			for (;ch[p][x]==q;p=fa[p]) ch[p][x]=nq;
		}
	}
}
void solve()
{
	p=root; int tmp=0;
	for (int i=1;i<=cnt;i++) {
		int x=s[i]-'0';
		if (ch[p][x]!=0) 
		  tmp++,p=ch[p][x];
		else {
			while (!ch[p][x]&&p) p=fa[p];
			if (!p) p=root,tmp=0;
			else tmp=l[p]+1,p=ch[p][x]; 
		}
		ans[i]=tmp;
	}
}
bool pd(int len)
{
	if (len==0) return true;
	int head=1,tail=0;
	f[0]=0; 
	for (int i=1;i<=cnt;i++) {
		f[i]=f[i-1];
		int pos=i-len;
		if (pos>=0){		 
		  while (f[st[tail]]-st[tail]<f[pos]-pos&&head<=tail)  tail--;
		  st[++tail]=pos;
		}
		while (head<=tail&&st[head]<i-ans[i]) head++;
		if (head<=tail) f[i]=max(f[i],f[st[head]]-st[head]+i);
	}
	if ((double)f[cnt]*1.0>=(double)cnt*0.9-eps) return true;
	else return false;
}
int main()
{
   freopen("a.in","r",stdin);
   freopen("my.out","w",stdout);
   scanf("%d%d",&n,&m);
   root=size=1;
   for (int i=1;i<=m;i++) {
   	 scanf("%s",s1+1);
   	 last=root;
   	 int len=strlen(s1+1);
   	 for (int j=1;j<=len;j++) 
		 extend(j);
   }
   for (int i=1;i<=n;i++) {
   	 scanf("%s",s+1);
   	 cnt=strlen(s+1);
     solve(); 
     cn=0; int l=0; int r=cnt;
   	 int ans1=0;
   	 while (l<=r) {
   	   int mid=(l+r)/2;
	   if (pd(mid)) ans1=max(ans1,mid),l=mid+1;
	   else r=mid-1;
	 }
	 printf("%d\n",ans1);
   }
}