BZOJ2806: [Ctsc2012]Cheat

这道题思路很好出来，但是不是很好写（自己太弱的原因）。

将那M个串连起来，建一个后缀自动机。

二分L，变成判定问题。

然后对每个询问串，在后缀自动机上跑出以i位置结尾的最长长度，设为c[i]。

在之后我们要做的是在这种条件下，求出最大覆盖长度。

之后我思路出了点问题，我以为就是用[i-c[i]+1，i]这样的线段做线段覆盖，wa了。

翻了翻题解，发现线段是可以在[L,c[i]]之间变化的。

列状态转移式：f[i]=max(f[i-1],max(f[i-j]+j))  j∈[L,c[i]]

可以看出i-c[i]和i-L都是单调不减的，搞一个单调队列即可。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<string>
#include<cstring>
#include<cmath>
#include<ctime>
#include<algorithm>
#include<iomanip>
#include<set>
#include<map>
#include<queue>
using namespace std;
#define mem1(i,j) memset(i,j,sizeof(i))
#define mem2(i,j) memcpy(i,j,sizeof(i))
#define LL long long
#define up(i,j,n) for(int i=(j);i<=(n);i++)
#define FILE "dealing"
#define poi vec
#define eps 1e-10
#define db double 
const int maxn=1010000,inf=1000000000,mod=1000000007;
int read(){
	int x=0,f=1,ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch<='9'&&ch>='0'){x=(x<<1)+(x<<3)+ch-'0',ch=getchar();}
	return f*x;
}
bool cmax(int& a,int b){return a<b?a=b,true:false;}
bool cmin(int& a,int b){return a>b?a=b,true:false;}
struct SAM{
	int pre[maxn],c[maxn][3],step[maxn],sa[maxn],cou[maxn],val[maxn],cnt,now,Len;
	SAM(){mem1(pre,0);mem1(c,0);mem1(step,0);cnt=now=1;}
	int extend(int x){
		int np,nq,q,p;
		p=now;now=np=++cnt;step[np]=step[p]+1;val[np]++;
		while(p&&!c[p][x])c[p][x]=np,p=pre[p];
		if(!p)pre[np]=1;
		else {
			q=c[p][x];
			if(step[q]==step[p]+1)pre[np]=q;
			else {
				step[nq=++cnt]=step[p]+1;
				mem2(c[nq],c[q]);
				pre[nq]=pre[q];
				pre[q]=pre[np]=nq;
				while(p&&c[p][x]==q)c[p][x]=nq,p=pre[p];
			}
		}
	}
	int getsort(){
		up(i,1,cnt)cou[step[i]]++;
		up(i,1,cnt)cou[i]+=cou[i-1];
		for(int i=cnt;i>=1;i--)sa[cou[step[i]]--]=i;
	}
	int walkprepare(){now=1,Len=0;}
	int walk(int x){
		while(pre[now]&&!c[now][x])now=pre[now],Len=step[now];
		if(!c[now][x])return 0;
		Len++;now=c[now][x];return Len;
	}
	int prepare(){cnt=now=1;}
	int build(char* s){
		int n=strlen(s+1);
		up(i,1,n)extend(s[i]);
	}
}a;
char s[maxn];
int c[maxn],f[maxn],q[maxn],head=0,tail=0;
int dcmp(db a){if(fabs(a)<=eps)return 0;return a<0?-1:1;}
bool check(int L){
	head=1;tail=0;
	int n=strlen(s+1),ans=0;
	up(i,1,n)f[i]=0;mem1(q,0);
	up(i,1,n){
		while(i>=L&&head<=tail&&f[i-L]+n-(i-L)>f[q[tail]]+n-q[tail])tail--;
		if(i>=L)q[++tail]=i-L;
		f[i]=f[i-1];
		while(head<=tail&&i-c[i]>q[head])head++;
		if(head<=tail&&i-q[head]>=L&&i-q[head]<=c[i])cmax(f[i],f[q[head]]+i-q[head]);
	}
	up(i,1,n)cmax(ans,f[i]);
	return dcmp(ans*1.0/n-0.90000)>=0;
}
int main(){
	freopen(FILE".in","r",stdin);
	freopen(FILE".out","w",stdout);
	int Q=read(),m=read();
	a.prepare();
	up(i,1,m){
		scanf("%s",s+1);int mm=strlen(s+1);
		up(j,1,mm)a.extend(s[j]-'0');
		a.extend(2);
	}
	while(Q--){
		scanf("%s",s+1);
		int left=0,right=strlen(s+1),mid,ans=0;
		up(i,1,right)c[i]=a.walk(s[i]-'0');
		while(left<=right){
			mid=(left+right)>>1;
			if(check(mid))left=mid+1,ans=mid;
			else right=mid-1;
		}
		printf("%d
",ans);
	}
	return 0;
}