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;
- }