SPOJ - LCS 后缀自动机入门

LCS - Longest Common Substring

A string is finite sequence of characters over a non-empty finite set Σ.

In this problem, Σ is the set of lowercase letters.

Substring, also called factor, is a consecutive sequence of characters occurrences at least once in a string.

Now your task is simple, for two given strings, find the length of the longest common substring of them.

Here common substring means a substring of two or more strings.

Input

The input contains exactly two lines, each line consists of no more than 250000 lowercase letters, representing a string.

Output

The length of the longest common substring. If such string doesn't exist, print "0" instead.

Example

Input:
alsdfkjfjkdsal
fdjskalajfkdsla

Output:
3

Notice: new testcases added

题意：

　　求两个串的最长公共子串。

题解：

　　先用第一个串构造出后缀自动机，然后逐个的匹配第二个串

　　如果当前节点失配，利用Suffix Links 往后跳就可以了。

#include <bits/stdc++.h>
inline long long read(){long long x=,f=;char ch=getchar();while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}while(ch>=''&&ch<=''){x=x*+ch-'';ch=getchar();}return x*f;}
using namespace std;

const int N = 3e5+;

const long long mod = ;

int isPlus[N * ],endpos[N * ];int d[N * ];
int tot,slink[*N],trans[*N][],minlen[*N],maxlen[*N],pre;
int newstate(int _maxlen,int _minlen,int* _trans,int _slink){
    maxlen[++tot]=_maxlen;minlen[tot]=_minlen;
    slink[tot]=_slink;
    if(_trans)for(int i=;i<;i++)trans[tot][i]=_trans[i],d[_trans[i]]+=;
    return tot;
}
int add_char(char ch,int u){
    int c=ch-'a',v=u;
    int z=newstate(maxlen[u]+,-,NULL,);
    isPlus[z] = ;
    while(v&&!trans[v][c]){trans[v][c]=z;d[z]+=;v=slink[v];}
    if(!v){ minlen[z]=;slink[z]=;return z;}
    int x=trans[v][c];
    if(maxlen[v]+==maxlen[x]){slink[z]=x;minlen[z]=maxlen[x]+;return z;}
    int y=newstate(maxlen[v]+,-,trans[x],slink[x]);
    slink[z]=slink[x]=y;minlen[x]=minlen[z]=maxlen[y]+;
    while(v&&trans[v][c]==x){trans[v][c]=y;d[x]--,d[y]++;v=slink[v];}
    minlen[y]=maxlen[slink[y]]+;
    return z;
}
void init_sam() {
    for(int i = ; i <= tot; ++i)
        for(int j = ; j < ; ++j) trans[i][j] = ;
    pre = tot = ;
}
char a[N],b[N];
int main() {
    scanf("%s%s",a,b);
    init_sam();
    int n = strlen(a);
    int m = strlen(b);
    for(int i = ; i < n; ++i)
        pre = add_char(a[i],pre);
    int ans = ;
    int now = ;
    int sum = ;
    for(int i = ; i < m; ++i) {
        if(trans[now][b[i]-'a']) {
            now = trans[now][b[i]-'a'];
        } else {
            while(now) {
                now = slink[now];
                if(trans[now][b[i]-'a']) break;
            }
            if(!now) now = ,sum = ;
            else {
                if(trans[now][b[i]-'a']) {
                    sum = maxlen[now] + ;
                    now = trans[now][b[i]-'a'];
                }
            }
        }
        ans = max(ans,sum);
    }
    printf("%d\n",ans);
}