bzoj2946 [Poi2000]公共串

SAM模板题。

主要思路是,先对一个串建SAM,然后拿其他串在上面匹配,对单个串每个状态取最大的步数,对所有的串每个状态取最小步数。

循环的同时需要对fail树按topsort序更新最大值。

为什么要从fail树更新最大值?我的理解是,单个节点上记录的Max,Min等是结束位置属于此节点right集的最大公共长度,fail树上父亲的right集是儿子的并,因此每次需要dp一下。

你可能会认为最大值不可能被破坏掉,但实际上是可能的。

#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 "pow"
#define poi vec
#define eps 1e-10
#define db double 
const int maxn=4010,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;}
int pre[maxn],c[maxn][26],len[maxn],Max[maxn],Min[maxn],now,cnt,np,p,q,nq,Len;
void clear(){mem1(c,0);mem1(len,0);mem1(pre,0);cnt=1,now=1;}
void expend(int x){
	np=++cnt;len[np]=len[now]+1;p=now;now=np;
	while(p&&!c[p][x])c[p][x]=np,p=pre[p];
	if(!p)pre[np]=1;
	else {
		q=c[p][x];
		if(len[q]==len[p]+1)pre[np]=q;
		else {
			len[nq=++cnt]=len[p]+1;
			pre[nq]=pre[q];
			pre[q]=pre[np]=nq;
			mem2(c[nq],c[q]);
			while(p&&c[p][x]==q)c[p][x]=nq,p=pre[p];
		}
	}
}
void walkprepare(){now=1,Len=0,mem1(Min,0);;}
int walk(int x){
	while(!c[now][x]&&pre[now])now=pre[now],Len=len[now];
	if(!c[now][x])return 0;
	now=c[now][x];Len++;return Len;
}
void build(char* s){
	clear();
	int n=strlen(s+1);
	up(i,1,n)expend(s[i]-'a');
}
int n,que[maxn],head,tail,ru[maxn];
char s[maxn];
int temp[maxn],cchu[maxn];
int main(){
	freopen(FILE".in","r",stdin);
	freopen(FILE".out","w",stdout);
	n=read();scanf("%s",s+1);build(s);
	mem1(Max,10);
	up(i,1,cnt)cchu[len[i]]++;
	up(i,1,cnt)cchu[i]+=cchu[i-1];//这样写的主要依据是max(pre s)=min(s)-1,以及从后往前进行更新。就是一个快速求top序的过程。
	for(int i=cnt;i>=1;i--)temp[cchu[len[i]]--]=i;//这样做稍快点,写topsort也没什么问题。
	up(i,1,n-1){
		walkprepare();
		scanf("%s",s+1);int m=strlen(s+1),v;
		up(j,1,m){
			v=walk(s[j]-'a');
			cmax(Min[now],v);
		}
		up(j,1,cnt)cmin(Max[j],Min[j]);
		for(int i=cnt;i>=1;i--)Max[pre[temp[i]]]=max(Max[temp[i]],Max[pre[temp[i]]]);
	}
	printf("%d\n",Max[1]);
	return 0;
}

  

posted @ 2017-02-20 15:56  CHADLZX  阅读(347)  评论(0编辑  收藏  举报