[POI2000]Repetitions

题目大意：
　　求多个字符串的LCS。

思路：
　　同SPOJ-LCS2，不过因为SPOJ上数据比较水，当时用错误的写法过掉了，这次用正确的写法重新过了一遍。
　　拓扑排序按照每个状态的len值，用计数排序实现。
　　每个状态往上更新时，应该对std::min(s[p].maxlen,s[q].len)取max（每个状态能匹配的最长子串长度不超过len）。

#include<cstdio>
#include<cstring>
#include<algorithm>
const int inf=0x7fffffff;
const int LEN=2001;
int n;
char s[LEN];
class SuffixAutomaton {
    private:
        static const int SIGMA_SIZE=26;
        struct State {
            int len,link,go[SIGMA_SIZE],maxlen,min;
        };
        int sz,root,last,len;
        int cnt[LEN],top[LEN<<1];
        State s[LEN<<1];
        int newState(const int l) {
            sz++;
            s[sz].len=l;
            s[sz].min=inf;
            return sz;
        }
        int idx(const char ch) {
            return ch-'a';
        }
        void extend(const char ch) {
            int w=idx(ch);
            int p=last;
            int new_p=newState(s[p].len+1);
            while(p&&!s[p].go[w]) {
                s[p].go[w]=new_p;
                p=s[p].link;
            }
            if(!p) {
                s[new_p].link=root;
                s[root].go[w]=new_p;
            } else {
                int q=s[p].go[w];
                if(s[q].len==s[p].len+1) {
                    s[new_p].link=q;
                } else {
                    int new_q=newState(s[p].len+1);
                    memcpy(s[new_q].go,s[q].go,sizeof s[q].go);
                    s[new_q].link=s[q].link;
                    s[q].link=s[new_p].link=new_q;
                    while(p&&s[p].go[w]==q) {
                        s[p].go[w]=new_q;
                        p=s[p].link;
                    }
                }
            }
            last=new_p;
        }
    public:
        void build(char str[]) {
            root=last=newState(0);
            len=strlen(str);
            for(int i=0;str[i];i++) extend(str[i]);
        }
        void top_sort() {
            for(int i=1;i<=sz;i++) cnt[s[i].len]++;
            for(int i=len;i;i--) cnt[i-1]+=cnt[i];
            for(int i=1;i<=sz;i++) top[cnt[s[i].len]--]=i;
        }
        void match(char str[]) {
            int p=root,tmp=0;
            for(int i=0;str[i];i++) {
                int w=idx(str[i]);
                if(s[p].go[w]) {
                    tmp++;
                    p=s[p].go[w];
                } else {
                    while(p&&!s[p].go[w]) {
                        p=s[p].link;
                    }
                    if(!p) {
                        tmp=0;
                        p=root;
                    } else {
                        tmp=s[p].len+1;
                        p=s[p].go[w];
                    }
                }
                s[p].maxlen=std::max(s[p].maxlen,tmp);
            }
            for(int i=1;i<=sz;i++) {
                int p=top[i],q=s[p].link;
                s[q].maxlen=std::max(s[q].maxlen,std::min(s[p].maxlen,s[q].len));
                s[p].min=std::min(s[p].min,s[p].maxlen);
                s[p].maxlen=0;
            }
        }
        int lcs() {
            int ret=0;
            for(int i=1;i<=sz;i++) ret=std::max(ret,s[i].min);
            return ret;
        }
};
SuffixAutomaton sam;
int main() {
    scanf("%d",&n);
    scanf("%s",s);
    sam.build(s);
    sam.top_sort();
    for(int i=2;i<=n;i++) {
        scanf("%s",s);
        sam.match(s);
    }
    printf("%d
",sam.lcs());
    return 0;
}