BZOJ 2946: [Poi2000]公共串

2946: [Poi2000]公共串

Time Limit: 3 Sec  Memory Limit: 128 MB
Submit: 787  Solved: 342
[Submit][Status][Discuss]

Description

 

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

任务：

l        读入单词

l        计算最长公共子串的长度

l        输出结果

 

Input

 

文件的第一行是整数 n，1<=n<=5，表示单词的数量。接下来n行每行一个单词，只由小写字母组成，单词的长度至少为1，最大为2000。

 

Output

仅一行，一个整数，最长公共子串的长度。

 

Sample Input

3
abcb
bca
acbc

Sample Output

２

HINT

Source

[Submit][Status][Discuss]

求若干字符串的最长公共子串。

和SPOJ的那道题简直一模一样，后缀自动机。

 #include <bits/stdc++.h>

 const int maxn = 5e5 + ;

 /* AUTOMATON */

 int last;

 int tail;

 int mini[maxn];

 int maxi[maxn];

 int step[maxn];

 int fail[maxn];

 int next[maxn][];

 inline void buildAutomaton(char *s)

 {

     last = ;

     tail = ;

     while (*s)

     {

         int c = *s++ - 'a';

         int p = last;

         int t = tail++;

         step[t] = step[p] + ;

         mini[t] = maxi[t] = step[t];

         while (p && !next[p][c])

             next[p][c] = t, p = fail[p];

         if (p)

         {

             int q = next[p][c];

             if (step[q] == step[p] + )

                 fail[t] = q;

             else

             {

                 int k = tail++;

                 fail[k] = fail[q];

                 fail[q] = fail[t] = k;

                 step[k] = step[p] + ;

                 mini[k] = maxi[k] = step[k];

                 for (int i = ; i < ; ++i)

                     next[k][i] = next[q][i];

                 while (p && next[p][c] == q)

                     next[p][c] = k, p = fail[p];

             }

         }

         else

             fail[t] = ;

         last = t;

     }

 }

 inline void searchAutomaton(char *s)

 {

     memset(maxi, , sizeof(maxi));

     for (int t = , k = ; *s; )

     {

         int c = *s++ - 'a';

         if (next[t][c])

             ++k, t = next[t][c];

         else

         {

             while (t && !next[t][c])

                 t = fail[t];

             if (t)

                 k = step[t] + , t = next[t][c];

             else

                 k = , t = ;

         }

         if (maxi[t] < k)

             maxi[t] = k;

     }

     for (int i = tail - ; i; --i)

         if (maxi[fail[i]] < maxi[i])

             maxi[fail[i]] = maxi[i];

     for (int i = ; i < tail; ++i)

         if (mini[i] > maxi[i])

             mini[i] = maxi[i];

 }

 inline int calculateAutomaton(void)

 {

     register int ret = ;

     for (int i = ; i < tail; ++i)

         if (ret < mini[i])

             ret = mini[i];

     return ret;

 }

 /* MAIN FUNC */

 char str[maxn];

 signed main(void)

 {

     int n;

     scanf("%d", &n);

     scanf("%s", str);

     buildAutomaton(str);

     for (int i = ; i < n; ++i)

         scanf("%s", str), searchAutomaton(str);

     printf("%d\n", calculateAutomaton());

 }

@Author: YouSiki