求最长公共子串

poj2774，codevs3160

题目描述 Description

给出两个由小写字母组成的字符串，求它们的最长公共子串的长度。

输入描述 Input Description

读入两个字符串

输出描述 Output Description

输出最长公共子串的长度

样例输入 Sample Input

yeshowmuchiloveyoumydearmotherreallyicannotbelieveit
yeaphowmuchiloveyoumydearmother

样例输出 Sample Output

27

数据范围及提示 Data Size & Hint

单个字符串的长度不超过100000

 

错误代码

        else{
            while(now&&!c[now][s2[i]-'a'])now=f[now],qwq;
                if(c[now][s2[i]-'a']){
                now=c[now][s2[i]-'a'];len=l[now]+1;if(len>ans)ans=len;
            }
        
        }    

一定要判断now是否跳到了0，否则妥妥地死循环

#include<cstring>
#include<cstdio>
const int N=200005;
#define qwq puts("*")
char s1[N>>1],s2[N>>1];
int pre=1,cnt=1,f[N],l[N],l1,l2,c[N][27],ans;

inline void ins(int x){
    int p=pre,np=++cnt;pre=np;l[np]=l[p]+1;
    for(;p&&!c[p][x];p=f[p])c[p][x]=np;
    if(!p)f[np]=1;
    else{
        int q=c[p][x];
        if(l[q]==l[p]+1)f[np]=q;
        else{
            int nq=++cnt;
            memcpy(c[nq],c[q],sizeof(c[q]));
            f[nq]=f[q];l[nq]=l[p]+1;
            f[q]=f[np]=nq;
            for(;c[p][x]==q;p=f[p])c[p][x]=nq;
        }
    }    
}
int main(){
    scanf("%s%s",s1,s2);
    l1=strlen(s1);l2=strlen(s2);//printf("%d %d\n",l1,l2);
    for(int i=0;i<l1;++i)ins(s1[i]-'a');
    int now=1,len=0;
    for(int i=0;i<l2;++i){
        if(c[now][s2[i]-'a'])++len,now=c[now][s2[i]-'a'];
        else{
            while(now&&!c[now][s2[i]-'a'])now=f[now];
            if(!now)len=0,now=1;
            else{
                len=l[now]+1;now=c[now][s2[i]-'a'];
            }
        }
        if(len>ans)ans=len;
    }
    printf("%d",ans);
    return 0;
}