BZOJ 2806: [Ctsc2012]Cheat(单调队列优化dp+后缀自动机)

传送门

解题思路

　　肯定先要建出来广义后缀自动机。刚开始以为是个二分＋贪心，写了一下结果(20)分。说一下正解，首先显然(L_0)具有单调性，是可以二分的。考虑二分后怎样判合法，对于分割序列很容易想到(dp)，设(f_i)表示前(i)个字符匹配成功数量，那么有转移方程(f_i=max(f_j+i-j)(i-j>=L) 且 (j)到(i)可以匹配 ())，(L)是二分出来的限制，判断是否能匹配可以预处理，预处理出(mth_i)表示(i)最多能与往前(mth_i)位匹配成功，那么第二个条件就变成了(j>=i-mth_i)。如果这样做是(O(n^2logn))的，实测可以拿到(75)分2333。考虑优化，发现(i-mth_i)具有单调性，因为每移动一格(i)会(+1)，而(mth_i)最多(+1)。那么可以用一个单调递减队列来优化，每次将(i-lim)入队，取出队头更新答案，时间复杂度为(O(nlogn))。

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>

using namespace std;
const int N=1100005<<1;

int n,m,res,ans,f[N],mth[N],q[N],hd,tl;
char s[N];

struct SAM{
	int ch[N][2],fa[N],l[N],lst,cnt;
	void Insert(int c){
		int p=lst,np=++cnt; l[np]=l[p]+1; lst=cnt;
		for(;p && !ch[p][c];p=fa[p]) ch[p][c]=np;
		if(!p) fa[np]=1;
		else {
			int q=ch[p][c]; 
			if(l[p]+1==l[q]) fa[np]=q;
			else {
				int nq=++cnt; 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][c]==q;p=fa[p]) ch[p][c]=nq;
			}
		}
	}
	void prework(int len){
		int now=1,num=0;
		for(int i=1;i<=len;i++){
			if(ch[now][s[i]-'0']) now=ch[now][s[i]-'0'],num++;
			else {
				for(;now && !ch[now][s[i]-'0'];now=fa[now]);
				if(!now) now=1,num=0;
				else num=l[now]+1,now=ch[now][s[i]-'0'];
			}
			mth[i]=num;
		}
	}
	bool check(int lim,int len){	
		int tmp; hd=1; tl=0; 
		for(int i=1;i<=len;i++){
			f[i]=f[i-1]; if(i<lim) continue; tmp=i-mth[i];
			while(hd<=tl && f[q[tl]]-q[tl]<=f[i-lim]-i+lim) tl--;
			q[++tl]=i-lim;
			while(hd<=tl && q[hd]<tmp) hd++;
			if(hd<=tl) f[i]=max(f[i],f[q[hd]]-q[hd]+i);
		}
//		for(int i=1;i<=len;i++){
//			f[i]=f[i-1];
//			for(int j=max(0,i-mth[i]);j+lim<=i;j++)
//				f[i]=max(f[i],f[j]+i-j);
//		}
		return len-f[len]>res?0:1;
	}
	void solve(int len){	
		prework(len);
		int L=1,R=len,mid;
		while(L<=R){
			mid=(L+R)>>1;
			if(check(mid,len)) L=mid+1,ans=mid;
			else R=mid-1;
		}
		printf("%d
",ans);
	}
}sam;

int main(){
	scanf("%d%d",&n,&m); sam.cnt=1; int len;
	for(int i=1;i<=m;i++){
		scanf("%s",s+1); sam.lst=1;
		len=strlen(s+1);
		for(int j=1;j<=len;j++) sam.Insert(s[j]-'0');
	}
	for(int i=1;i<=n;i++){
		scanf("%s",s+1); len=strlen(s+1);
		res=(len*9+9)/10; res=len-res;
		sam.solve(len); ans=0;
	}
	return 0;
}