spoj 1812 lcs2

題目大意：給你n個串，求它們的最長公共子串   n<=10,每個串<=100000

這題是後綴自動機的論文題

具體的做法是，先以第一個串做一個後綴自動機，然後用剩下的每個串去匹配就行了

具體來說，就是用每個串都像lcs那樣去做就行了，然後把每個節點可以擴展的最長長度記錄下來，再用兒子更新父親，就OK了

#include<cmath>
#include<cstdio>
#include<cstring>
#include<cassert>
#include<iostream>
#include<algorithm>
using namespace std;
const int maxn=100011;
char s[maxn]; int n;
void init(){scanf("%s",s); n=strlen(s);}
struct Tsam{
	struct Tnode{
		Tnode *ch[26],*f; int ml,nl;
		void clear(){memset(ch,0,sizeof(ch)); f=NULL; ml=nl=0;}
	}e[maxn<<1],*root,*last;
	int tot;
	Tnode *newnode(){e[tot].clear(); return e+(tot++);}
	void clear(){tot=0; root=last=newnode();}
	void add(int c){
		Tnode *np=newnode(),*p=last;
		for (np->ml=p->ml+1,last=np;p && !p->ch[c];p=p->f) p->ch[c]=np;
		if (!p) return np->f=root,void();
		Tnode *q=p->ch[c];
		if (q->ml==p->ml+1) return np->f=q,void();
		Tnode *r=newnode(); *r=*q; r->ml=p->ml+1;
		for (q->f=np->f=r;p && p->ch[c]==q;p=p->f) p->ch[c]=r;
	}
}sam;
void work(){
	static int ws[maxn]; static Tsam::Tnode *tp[maxn<<1];
	sam.clear(); for (int i=0;i<n;++i) sam.add(s[i]-'a');
	for (int i=0;i<=n;++i) ws[i]=0;
	for (int i=0;i<sam.tot;++i) ++ws[sam.e[i].ml];
	for (int i=1;i<=n;++i) ws[i]+=ws[i-1];
	for (int i=sam.tot-1;i>=0;--i) tp[--ws[sam.e[i].ml]]=sam.e+i;
	int ans;
	while (scanf("%s",s)!=EOF){
		Tsam::Tnode *now=sam.root;
		for (int i=0,res=0;s[i];++i){
			int c=s[i]-'a';
			if (now->ch[c]) ++res,now=now->ch[c];
			else{
				for (;now && !now->ch[c];now=now->f);
				if (!now) now=sam.root,res=0;
				else res=now->ml+1,now=now->ch[c];
			}
			if (res>now->nl) now->nl=res;
		}
		for (int i=sam.tot-1;i>=0;--i){
			Tsam::Tnode *x=tp[i];
			if (x->ml>x->nl) x->ml=x->nl;
			if (x->f && x->f->nl<x->nl) x->f->nl=x->nl;
			x->nl=0;
		}
	}
	ans=0;
	for (int i=1;i<sam.tot;++i) if (sam.e[i].ml>ans) ans=sam.e[i].ml;
	cout<<ans<<endl;
}
int main(){
	init();
	work();
	return 0;
}

還有一點要注意，如果不用ml來存每個擴展最大的最小值，在用兒子更新父親的時候，就要判斷有沒有超過父親的ml，WA了好久。。。