bzoj 2946: [Poi2000]公共串 后缀自动机

题目大意:

给出几个由小写字母构成的单词，求它们最长的公共子串的长度。

题解:

我们首先根据一个串构建后缀自动机
然后我们依次让每个串上去跑。
在每个节点记录一下当前已经匹配的最大的长度
然后在每次跑不同的串的过程中相应节点取min值

然后取所有的min的max输出就好了..

其实可以说成:ans是parent树中子节点的max中相应点不同串的min的max值.

...好晕啊...

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long ll;
inline void read(int &x){
    x=0;char ch;bool flag = false;
    while(ch=getchar(),ch<'!');if(ch == '-') ch=getchar(),flag = true;
    while(x=10*x+ch-'0',ch=getchar(),ch>'!');if(flag) x=-x;
}
const int maxn = 2048;
struct Node{
    int nx[26];
    int len,fa;
}T[maxn<<1];
int last,nodecnt;
inline void init(){
    last = nodecnt = 0;
    T[0].len = 0;T[0].fa = -1;
    for(int k=0;k<26;++k) T[0].nx[k] = 0;
}
inline void insert(int c){
    int cur = ++ nodecnt,p;
    T[cur].len = T[last].len + 1;
    for(p = last;p != -1 && !T[p].nx[c];p = T[p].fa) T[p].nx[c] = cur;
    if(p == -1) T[cur].fa = 0;
    else{
        int q = T[p].nx[c];
        if(T[q].len == T[p].len + 1) T[cur].fa = q;
        else{
            int co = ++ nodecnt;
            T[co].len = T[p].len + 1;T[co].fa = T[q].fa;
            for(int i=0;i<26;++i) T[co].nx[i] = T[q].nx[i];
            for(;p != -1 && T[p].nx[c] == q;p = T[p].fa) T[p].nx[c] = co;
            T[cur].fa = T[q].fa = co;
        }
    }last = cur;
}
inline void build(char *s){
    init();int len = strlen(s);
    for(int i=0;i<len;++i) insert(s[i]-'a');
}
int q[maxn<<1],sum[maxn<<1];
inline void Sort(){
    for(int i=1;i<=nodecnt;++i) sum[T[i].len] ++ ;
    for(int i=1;i<=4000;++i) sum[i] += sum[i-1];
    for(int i=1;i<=nodecnt;++i) q[sum[T[i].len]--] = i;
}
int f[maxn<<1],g[maxn<<1];
inline void find(char *s){
    memset(g,0,sizeof g);
    int p = 0,len = strlen(s);
    for(int i=0,l = 0;i<len;++i){
        int c = s[i] - 'a';
        if(T[p].nx[c]) ++l,p =T[p].nx[c];
        else{
            while(p && !T[p].nx[c]) p = T[p].fa;
            if(!p) l = 0;
            else l = T[p].len+1,p = T[p].nx[c];

        }g[p] = max(g[p],l); 
    }
    for (int i=nodecnt;i;i--) g[T[q[i]].fa] = max(g[T[q[i]].fa],g[q[i]]);
    for(int i=0;i<=nodecnt;++i) f[i] = min(f[i],g[i]);
}
char s[maxn<<1];
int main(){
    int n;read(n);scanf("%s",s);build(s);
    for(int i=1;i<=nodecnt;++i) f[i] = T[i].len;
    Sort();
    for(int i=2;i<=n;++i){
        scanf("%s",s);
        find(s);
    }
    int ans = 0;
    for(int i=1;i<=nodecnt;++i){
        ans = max(ans,f[i]);
    }printf("%d\n",ans);
    getchar();getchar();
    return 0;
}

转载于:https://www.cnblogs.com/Skyminer/p/6522900.html