poj 1226 Substrings

题意：就是求n个字符串的最长公共子串，子串是可以反转的

思路：strstr，kmp爆搜，还可以用后缀数组

实测kmp比strstr快，16ms那个是后缀数组的

8389284	Euler_M	1226	Accepted	368K	141MS	G++	1911B	2011-03-28 16:44:59
8389280	Euler_M	1226	Accepted	364K	63MS	G++	1907B	2011-03-28 16:44:19
7957084	Euler_M	1226	Accepted	508K	16MS	C++	2443B	2010-12-02 20:07:26

 

kmp(strstr只用改一下就可以了)：

/* * File: main.cpp * Author: Mi * * Created on 2011年3月28日, 下午4:09 */ #include <cstdlib> #include <stdio.h> #include <string.h> #include <algorithm> #define N 105 using namespace std; /* * */ char str[N][N]; int P[N],P1[N]; void next(char *B,int m,int *P) { int i,j=0; P[1]=0; for(i=2;i<=m;i++) { while(j>0&&B[j+1]!=B[i]) j=P[j]; if(B[j+1]==B[i]) j++; P[i]=j; } } int kmp(char *A,char *B,int n,int m,int *P) { int i,j=0; for(i=1;i<=n;i++) { while(j>0&&B[j+1]!=A[i]) j=P[j]; if(B[j+1]==A[i]) j++; if(j==m) return 1; } return 0; } int main(int argc, char** argv) { int t; scanf("%d",&t); while(t--) { int n,ans=0; scanf("%d",&n); for(int i=1;i<=n;i++) scanf("%s",str[i]+1); for(int j=1;j<=strlen(str[1]+1);j++) { for(int k=j;k<=strlen(str[1]+1);k++) { char temp[N],temp1[N]; strncpy(temp+1,str[1]+j,k-j+1); int l=k-j+1,cnt=0; temp[l+1]='/0'; strcpy(temp1+1,temp+1); reverse(temp1+1,temp1+1+l); next(temp,l,P); next(temp1,l,P1); for(int h=2;h<=n;h++) if(kmp(str[h],temp,strlen(str[h]+1),l,P)||kmp(str[h],temp1,strlen(str[h]+1),l,P1)) // if(strstr(str[h]+1,temp+1)!=NULL||strstr(str[h]+1,temp1+1)!=NULL) cnt++; else break; if(cnt==n-1&&l>ans) ans=l; } } printf("%d/n",ans); } return 0; }

后缀数组代码：

#include "stdio.h" #include "string.h" #define maxn 20201 int wa[maxn],wb[maxn],wv[maxn],ws[maxn]; int cmp(int *r,int a,int b,int l) {return r[a]==r[b]&&r[a+l]==r[b+l];} void da(int *r,int *sa,int n,int m) { int i,j,p,*x=wa,*y=wb,*t; for(i=0;i<m;i++) ws[i]=0; for(i=0;i<n;i++) ws[x[i]=r[i]]++; for(i=1;i<m;i++) ws[i]+=ws[i-1]; for(i=n-1;i>=0;i--) sa[--ws[x[i]]]=i; for(j=1,p=1;p<n;j*=2,m=p) { for(p=0,i=n-j;i<n;i++) y[p++]=i; for(i=0;i<n;i++) if(sa[i]>=j) y[p++]=sa[i]-j; for(i=0;i<n;i++) wv[i]=x[y[i]]; for(i=0;i<m;i++) ws[i]=0; for(i=0;i<n;i++) ws[wv[i]]++; for(i=1;i<m;i++) ws[i]+=ws[i-1]; for(i=n-1;i>=0;i--) sa[--ws[wv[i]]]=y[i]; for(t=x,x=y,y=t,p=1,x[sa[0]]=0,i=1;i<n;i++) x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++; } return; } int rank[maxn],height[maxn]; void calheight(int *r,int *sa,int n) { int i,j,k=0; for(i=1;i<=n;i++) rank[sa[i]]=i; for(i=0;i<n;height[rank[i++]]=k) for(k?k--:0,j=sa[rank[i]-1];r[i+k]==r[j+k];k++); return; } int r[maxn],sa[maxn]; int who[maxn],yes[101]={0},ii=0; int len,n; int check(int mid) { int i,j,k,t,s; for(i=2;i<=len;i=j+1) { for(;height[i]<mid && i<=len;i++); for(j=i;height[j]>=mid;j++); if(j-i+1<n) continue; ii++;s=0; for(k=i-1;k<j;k++) if((t=who[sa[k]])!=0) if(yes[t]!=ii) yes[t]=ii,s++; if(s>=n) return 1; } return 0; } char st[110]; int main() { int i,j,k,min,mid,max,nn; scanf("%d",&nn); while(nn-->0) { scanf("%d",&n); len=0; for(i=1;i<=n;i++) { scanf("%s",st); k=strlen(st); for(j=0;j<k;j++) { r[j+len]=st[j]+200; who[j+len]=i; } r[len+k]=2*i-1; who[len+k]=0; len+=k+1; for(j=0;j<k;j++) { r[j+len]=st[k-1-j]+200; who[j+len]=i; } r[len+k]=2*i; who[len+k]=0; len+=k+1; } len--; r[len]=0; da(r,sa,len+1,328); calheight(r,sa,len); height[len+1]=-1; min=1;max=100; while(min<=max) { mid=(min+max)>>1; if(check(mid)) min=mid+1; else max=mid-1; } if(n==1) max=len/2; printf("%d/n",max); } return 0; }