[SPOJ-LCS2]Longest Common Substring II

vjudge

题意

求多个串的最长公共子串

sol

首先拿第一个串建SAM。
然后对于每一个串，在SAM上跑匹配。
那么最终就是要把多个串跑出匹配的信息合并起来。
具体来说，对每个串记录一下在每个状态上的最长长度，然后同一个状态上对每个串的最长长度取个\(min\)，最后再把每个状态取个\(max\)。

code

#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
const int N = 2e5+5;
int n,last=1,tot=1,tr[N][26],fa[N],len[N],t[N],a[N],sum[N],ans[N],Ans;
char s[N];
void extend(int c)
{
	int v=last,u=++tot;last=u;
	len[u]=len[v]+1;
	while (v&&!tr[v][c]) tr[v][c]=u,v=fa[v];
	if (!v) fa[u]=1;
	else
	{
		int x=tr[v][c];
		if (len[x]==len[v]+1) fa[u]=x;
		else
		{
			int y=++tot;
			memcpy(tr[y],tr[x],sizeof(tr[y]));
			fa[y]=fa[x];fa[x]=fa[u]=y;len[y]=len[v]+1;
			while (v&&tr[v][c]==x) tr[v][c]=y,v=fa[v];
		}
	}
}
int main()
{
	scanf("%s",s+1);n=strlen(s+1);
	for (int i=1;i<=n;++i) extend(s[i]-'a');
	for (int i=1;i<=tot;++i) ++t[len[i]];
	for (int i=1;i<=tot;++i) t[i]+=t[i-1];
	for (int i=1;i<=tot;++i) a[t[len[i]]--]=i;
	for (int i=1;i<=tot;++i) ans[i]=len[i];
	while (scanf("%s",s+1)!=EOF)
	{
		n=strlen(s+1);
		memset(sum,0,sizeof(sum));
		for (int i=1,now=1,cnt=0;i<=n;++i)
		{
			int c=s[i]-'a';
			if (tr[now][c]) ++cnt,now=tr[now][c];
			else
			{
				while (now&&!tr[now][c]) now=fa[now];
				if (!now) cnt=0,now=1;
				else cnt=len[now]+1,now=tr[now][c];
			}
			sum[now]=max(sum[now],cnt);
		}
		for (int i=tot;i;--i) sum[fa[a[i]]]=max(sum[fa[a[i]]],sum[a[i]]);
		for (int i=1;i<=tot;++i) ans[i]=min(ans[i],sum[i]);
	}
	for (int i=1;i<=tot;++i) Ans=max(Ans,ans[i]);
	printf("%d\n",Ans);return 0;
}