【扩展KMP】POJ_3450｜ HDU_2328 Corporate Identity

原题直通车：POJ_3450 Corporate Identity HDU_2328
Corporate Identity

题意概述：找出N个串中最长公共子串

分析：

一、可以直接枚举其中一个串的所有字串，跟所有串进行匹配找到结果。

二、用其中一个串的每一个后缀和其它所有串分别求一次扩展KMP，并找到后缀中能在所有串中匹配到的最长前缀，最后在每个后缀的最长前缀中找出最长且字典序最小的。

代码参考：

KMP：

#include<iostream>
#include<cstring>
#include<string>
#include<cstdio>
using namespace std;
const int maxn=4444;

char x[222],ans[222], f[maxn][222];
int next[222], n;

void strcpy(char *y,int i,int len){
for(int j=0;j<len;++i,++j) x[j]=y[i];
x[len]='\0';
}
void Get_Next(char *s){
int len=strlen(s), i = 0, j = -1;
next[0]=-1;
while(i<len){
if(j==-1||s[i]==s[j]){
if(s[++i]==s[++j]) next[i]=next[j];
else next[i]=j;
}
else j=next[j];
}
}
bool KMP(char *s,char *a){
int ls=strlen(s), la=strlen(a);
Get_Next(a);
int i=0, j=0;
while(i<ls&&j<la){
if(j==-1||s[i]==a[j])
++i,++j;
else j=next[j];
}
if(j==la) return true;
return false;
}
int main(){
while(~scanf("%d",&n),n){
for(int i=0; i<n; ++i) scanf("%s",f[i]);
int len = strlen(f[0]);
bool flag = false;
ans[0]='\0';
for(int i=len; i && !flag; --i)
for(int j=0,k; j<=len-i; ++j) {
strcpy(f[0], j, i);
for(k=1; k<n; ++k) if( !KMP(f[k],x) ) break;
if(k==n && (!flag || strcmp(ans,x)>0)) {
strcpy(ans, x);
flag=true;
}
}
if(flag) puts(ans);
else puts("IDENTITY LOST");
}
return 0;
}

扩展KMP：

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
char s[4005][205], f[4005][205], ans[205];
int next[205], extend[205];
void get_next(char *g) {
int lg = strlen(g), j = 0;
next[0] = lg;
while(j+1<lg && g[j] == g[j+1]) ++j;
next[1] = j;
for(int i=2, k=1; i<lg; ++i) {
int p = k+next[k]-1;
int l = next[i-k];
if(i+l < p+1) next[i] = l;
else {
j = (p-i+1>0 ? p-i+1 : 0);
while(i+j<lg && g[j] == g[i+j]) ++j;
next[i] = j;
k = i;
}
}
}
int ExKMP(char *s, char *g) {
int ls = strlen(s), lg = strlen(g), j = 0, max_len;
while(j<lg && j<ls && s[j] == g[j]) ++j;
max_len = extend[0] = j;
for(int i=1, k=0; i<ls; ++i) {
int p = k+extend[k]-1;
int l = next[i-k];
if(i+l < p+1) extend[i] = l;
else {
j = (p-i+1>0 ? p-i+1 : 0);
while(i+j<ls && j<lg && s[i+j] == g[j]) ++j;
extend[i] = j;
k = i;
}
if(extend[i] > max_len) max_len = extend[i];
}
return max_len;
}
int main() {
int n;
while(~scanf("%d", &n) && n) {
for(int i=0; i<n; ++i) scanf("%s", s[i]);
int len = strlen(s[0]), max_len = -1;
for(int i=0; i<len; ++i) {
get_next(s[0]+i);
int mx = -1;
for(int j=1; j<n; ++j) {
int ml = ExKMP(s[j], s[0]+i);
if(mx == -1 || ml < mx) mx = ml;
}
if(mx > max_len || (mx == max_len&&strncmp(ans, s[0]+i, mx)>0)) {
strncpy(ans, s[0]+i, mx);
max_len = mx;
ans[mx] = '\0';
}
}
if(max_len>0) puts(ans);
else puts("IDENTITY LOST");
}
return 0;
}