BZOJ2946: [Poi2000]公共串

为啥没人写hash或者sa？不明觉厉一脸懵逼…………
于是我也跟风写一发SAM吧…>_<…
把第一个串建SAM，后面的每个放到SAM里跑，匹配到对应的节点，就与当前长度取max（len[p]=max(len[p],tmp)），表示这个串在这个节点能匹配的最大长度，然后每个串取min，意味着取出来的是公共字串（ans[i]=min(ans[i],len[i]))，最后所有ans取max，表示所有公共字串里最长的那个。注意每个串搞定之后要toposort一遍…>_<…
这个取max取min取max我也是无语了…………

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
//by:MirrorGray
using namespace std;
const int N=4111;
char s[N];
int n,cnt,root,last;
int gs[N],q[N],len[N],ans[N];
struct SAM{
    int nxt[26],fa,mx;
}tr[N];

int nwn(int l){
    tr[++cnt].mx=l;
    return cnt;
}

int extend(int x){
    int np=nwn(tr[last].mx+1);
    int p=last;last=np;
    for(;p && !tr[p].nxt[x];p=tr[p].fa)tr[p].nxt[x]=np;
    if(!p)return tr[np].fa=root,233;
    int q=tr[p].nxt[x];
    if(tr[q].mx==tr[p].mx+1)return tr[np].fa=q,233;
    int nq=nwn(tr[p].mx+1);
    tr[nq].fa=tr[q].fa;tr[q].fa=tr[np].fa=nq;
    memcpy(tr[nq].nxt,tr[q].nxt,sizeof(tr[q].nxt));
    for(;p && tr[p].nxt[x]==q;p=tr[p].fa)tr[p].nxt[x]=nq;
}

void toposort(){
    for(int i=1;i<=n;i++)gs[i]=0;
    for(int i=1;i<=cnt;i++)gs[ tr[i].mx ]++;
    for(int i=1;i<=n;i++)gs[i]+=gs[i-1];
    for(int i=cnt;i;i--)q[ gs[tr[i].mx]-- ]=i;
}

int main(){
    int m;scanf("%d",&m);m--;
    scanf("%s",s+1);n=strlen(s+1);
    root=last=nwn(0);
    for(int i=1;s[i];i++)extend(s[i]-'a');
    toposort();
    for(int i=1;i<=cnt;i++)ans[i]=tr[i].mx;
    while(m--){
        memset(len,0,sizeof(len));
        scanf("%s",s+1);
        int p=root,tmp=0;
        for(int i=1;s[i];i++){
            int c=s[i]-'a';
            while(p!=root && !tr[p].nxt[c])p=tr[p].fa,tmp=tr[p].mx;
            if(tr[p].nxt[c])tmp++,p=tr[p].nxt[c];
            len[p]=max(len[p],tmp);
        }
        for(int i=cnt;i;i--)len[ tr[q[i]].fa ]=max(len[ tr[q[i]].fa ],len[ q[i] ]);
        for(int i=1;i<=cnt;i++)ans[i]=min(ans[i],len[i]);
    }
    int pr=0;
    for(int i=1;i<=cnt;i++)pr=max(pr,ans[i]);
    printf("%d\n",pr);
    return 0;
}