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