SPOJ 1812 LCS2(SAM)

转载请注明出处，谢谢http://blog.csdn.net/ACM_cxlove?viewmode=contents    by---cxlove  

题目：求多个串的LCS

http://www.spoj.pl/problems/LCS2/ 

先将一个串建SAM，然后用后面的串去匹配，对于每一个串，保存最大值，对于不同的串，更新最小值。

SAM结点多两个值，ml表示多个串的最小值，nl表示当前串匹配的最大值。

对于SAM初学，要深刻理解出现次数向父亲传递，接收串数从儿子获取这句话。

#include<iostream>#include<cstdio>#include<map>#include<cstring>#include<cmath>#include<vector>#include<algorithm>#include<set>#include<string>#include<queue>#define inf 100000005#define M 40#define N 200015#define maxn 300005#define eps 1e-10#define zero(a) fabs(a)<eps#define Min(a,b) ((a)<(b)?(a):(b))#define Max(a,b) ((a)>(b)?(a):(b))#define pb(a) push_back(a)#define mp(a,b) make_pair(a,b)#define mem(a,b) memset(a,b,sizeof(a))#define LL unsigned long long#define MOD 1000000007#define lson step<<1#define rson step<<1|1#define sqr(a) ((a)*(a))#define Key_value ch[ch[root][1]][0]#define test puts("OK");#pragma comment(linker, "/STACK:1024000000,1024000000")using namespace std;struct SAM{    SAM *pre,*son[26];    int len,nl,ml;}*root,*tail,que[N],*b[N];int tot=0;char str[N/2];int cnt[N/2];void add(int c,int l){    SAM *p=tail,*np=&que[tot++];    np->len=np->ml=l;tail=np;    while(p&&p->son[c]==NULL) p->son[c]=np,p=p->pre;    if(p==NULL) np->pre=root;    else    {        SAM *q=p->son[c];        if(p->len+1==q->len) np->pre=q;        else        {            SAM *nq=&que[tot++];            *nq=*q;            nq->len=nq->ml=p->len+1;            np->pre=q->pre=nq;            while(p&&p->son[c]==q) p->son[c]=nq,p=p->pre;        }    }}int main(){   // freopen("1.in","r",stdin);    tail=root=&que[tot++];    scanf("%s",str);int l=strlen(str);    for(int i=0;str[i];i++) add(str[i]-'a',i+1);    for(int i=0;i<tot;i++) cnt[que[i].len]++;    for(int i=1;i<=l;i++) cnt[i]+=cnt[i-1];    for(int i=0;i<tot;i++) b[--cnt[que[i].len]]=&que[i];    while(scanf("%s",str)!=EOF)    {        SAM *p=root;        for(int i=0,l=0;str[i];i++)        {            int c=str[i]-'a';            if(p->son[c]) {l++;p=p->son[c];}            else            {                while(p&&p->son[c]==NULL) p=p->pre;                if(p==NULL) {l=0;p=root;}                else {l=p->len+1;p=p->son[c];}            }            if(l>p->nl) p->nl=l;        }        for(int i=tot-1;i>=0;i--)        {            p=b[i];            if(p->nl<p->ml) p->ml=p->nl;            if(p->pre&&p->pre->nl<p->nl) p->pre->nl=p->nl;            p->nl=0;        }    }    int ans=0;    for(int i=0;i<tot;i++) if(que[i].ml>ans) ans=que[i].ml;    printf("%d\n",ans);    return 0;}