SPOJ 1811-LCS POJ 2774

SAM表示真的很难艹

SPOJ LCS

题目链接

推荐学习的博客博客地址

  题意：求两个字符串的最长公共子串

  解法：时间限制及其的严格，SAM可以在线性时间内解决

性质:每个子串都可以由S走到

~ SAM真的很强大，但是我太弱

另外，一份代码两次AC。POJ 2774

 
 
/************************************************************************* 
    > File Name: lcs.cpp 
    > Author: cy 
    > Mail: 1002@qq.com  
    > Created Time: 14/11/25 11:58:50 
 ************************************************************************/ 
 
#include<iostream> 
#include<cstring> 
#include <algorithm> 
#include<cstdlib> 
#include<vector> 
#include<cmath> 
#include<stdlib.h> 
#include<iomanip> 
#include<list> 
#include<deque> 
#include<map> 
#include <stdio.h> 
#include <queue> 
 
const int maxn=1000+5; 
 
#define inf 0x3f3f3f3f 
  #define INF 0x3FFFFFFFFFFFFFFFLL 
#define rep(i,n) for(i=0;i<n;i++) 
 #define reP(i,n) for(i=1;i<=n;i++) 
 
#define ull unsigned long long 
 #define ll long long 
 
#define cle(a) memset(a,0,sizeof(a)) 
#define MAXN 500010+5 
#define MAXC 26 
using namespace std; 
char ch[MAXN];int ws[MAXN],rk[MAXN],sz[MAXN],ans[MAXN]; 
 
struct Suffix_Automaton{ 
    int child[MAXN][MAXC]; 
 
    int f[MAXN];//指向最近的接受态 
    int l[MAXN];//此节点代表的最长串的长度 
    int head,tail,tot;//头节点,最后的接受态,总节点 
    void clear(){ 
        memset(child,0,sizeof(child)); 
        memset(f,0,sizeof(f)); 
        memset(l,0,sizeof(l)); 
        head=tail=tot=1; 
    } 
    void add(int c){ 
        int p=tail,np=++tot;l[np]=l[p]+1;tail=np; 
        for(;p&&!child[p][c];p=f[p]) child[p][c]=np; 
        if(!p) f[np]=head; 
        else if(l[child[p][c]]==l[p]+1) f[np]=child[p][c]; 
        else{ 
            int q=child[p][c],r=++tot; 
            memcpy(child[r],child[q],MAXC*sizeof(int)); 
            f[r]=f[q];l[r]=l[p]+1;f[q]=f[np]=r; 
            for(;p&&child[p][c]==q;p=f[p]) child[p][c]=r; 
        } 
    } 
    void insert(char*ch){ 
        int l=strlen(ch); 
        for(int i=0;i<l;i++) add(ch[i]-'a'); 
    } 
    int solve(char *ch) 
    { 
        int len=strlen(ch); 
        int i; 
        int ip; 
        int temp=0; 
        int ans=0; 
        int p=head; 
        for(i=0;i<len;i++){ 
            ip=ch[i]-'a'; 
            if(child[p][ip]){ 
                temp++,p=child[p][ip]; 
            } 
            else{ 
                while(true){ 
                    if(!p)break; 
                    if(child[p][ip])break; 
                    p=f[p]; 
                } 
                if(!p){ 
                    temp=0; 
                    p=head; 
                } 
                else{ 
                    temp=l[p]+1; 
                    p=child[p][ip]; 
                } 
            } 
            ans=max(ans,temp); 
        } 
        return ans; 
    } 
}SAM; 
int main() 
{ 
#ifndef ONLINE_JUDGE 
     freopen("in.txt","r",stdin); 
     //freopen("out.txt","w",stdout); 
#endif 
    SAM.clear(); 
    scanf("%s",&ch); 
    SAM.insert(ch); 
    scanf("%s",&ch); 
    printf("%d\\n",SAM.solve(ch)); 
    return 0; 
}