poj 1226(转)

 

某大牛写的，方便个人今后参考使用(⊙o⊙)…

 

 

题意：就是求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只用改一下就可以了)：

[cpp] view plain copy

2. #include <cstdlib>
3. #include <stdio.h>
4. #include <string.h>
5. #include <algorithm>
6. #define N 105
7. using namespace std;
9. char str[N][N];
10. int P[N],P1[N];
11. void next(char *B,int m,int *P)
12. {
13. int i,j=0;
14. P[1]=0;
15. for(i=2;i<=m;i++)
16. {
17. while(j>0&&B[j+1]!=B[i])
18. j=P[j];
19. if(B[j+1]==B[i])
20. j++;
21. P[i]=j;
22. }
23. }
24. int kmp(char *A,char *B,int n,int m,int *P)
25. {
26. int i,j=0;
27. for(i=1;i<=n;i++)
28. {
29. while(j>0&&B[j+1]!=A[i])
30. j=P[j];
31. if(B[j+1]==A[i])
32. j++;
33. if(j==m)
34. return 1;
35. }
36. return 0;
37. }
38. int main(int argc, char** argv)
39. {
40. int t;
41. scanf("%d",&t);
42. while(t--)
43. {
44. int n,ans=0;
45. scanf("%d",&n);
46. for(int i=1;i<=n;i++)
47. scanf("%s",str[i]+1);
48. for(int j=1;j<=strlen(str[1]+1);j++)
49. {
50. for(int k=j;k<=strlen(str[1]+1);k++)
51. {
52. char temp[N],temp1[N];
53. strncpy(temp+1,str[1]+j,k-j+1);
54. int l=k-j+1,cnt=0;
55. temp[l+1]='/0';
56. strcpy(temp1+1,temp+1);
57. reverse(temp1+1,temp1+1+l);
58. next(temp,l,P);
59. next(temp1,l,P1);
60. for(int h=2;h<=n;h++)
61. if(kmp(str[h],temp,strlen(str[h]+1),l,P)||kmp(str[h],temp1,strlen(str[h]+1),l,P1))
62. // if(strstr(str[h]+1,temp+1)!=NULL||strstr(str[h]+1,temp1+1)!=NULL)
63. cnt++;
64. else
65. break;
66. if(cnt==n-1&&l>ans)
67. ans=l;
68. }
69. }
70. printf("%d/n",ans);
71. }
72. return 0;
73. }

后缀数组代码：

[cpp] view plain copy

1. #include "stdio.h"
2. #include "string.h"
3. #define maxn 20201
4. int wa[maxn],wb[maxn],wv[maxn],ws[maxn];
5. int cmp(int *r,int a,int b,int l)
6. {return r[a]==r[b]&&r[a+l]==r[b+l];}
7. void da(int *r,int *sa,int n,int m)
8. {
9. int i,j,p,*x=wa,*y=wb,*t;
10. for(i=0;i<m;i++) ws[i]=0;
11. for(i=0;i<n;i++) ws[x[i]=r[i]]++;
12. for(i=1;i<m;i++) ws[i]+=ws[i-1];
13. for(i=n-1;i>=0;i--) sa[--ws[x[i]]]=i;
14. for(j=1,p=1;p<n;j*=2,m=p)
15. {
16. for(p=0,i=n-j;i<n;i++) y[p++]=i;
17. for(i=0;i<n;i++) if(sa[i]>=j) y[p++]=sa[i]-j;
18. for(i=0;i<n;i++) wv[i]=x[y[i]];
19. for(i=0;i<m;i++) ws[i]=0;
20. for(i=0;i<n;i++) ws[wv[i]]++;
21. for(i=1;i<m;i++) ws[i]+=ws[i-1];
22. for(i=n-1;i>=0;i--) sa[--ws[wv[i]]]=y[i];
23. for(t=x,x=y,y=t,p=1,x[sa[0]]=0,i=1;i<n;i++)
24. x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++;
25. }
26. return;
27. }
28. int rank[maxn],height[maxn];
29. void calheight(int *r,int *sa,int n)
30. {
31. int i,j,k=0;
32. for(i=1;i<=n;i++) rank[sa[i]]=i;
33. for(i=0;i<n;height[rank[i++]]=k)
34. for(k?k--:0,j=sa[rank[i]-1];r[i+k]==r[j+k];k++);
35. return;
36. }
37. int r[maxn],sa[maxn];
38. int who[maxn],yes[101]={0},ii=0;
39. int len,n;
40. int check(int mid)
41. {
42. int i,j,k,t,s;
43. for(i=2;i<=len;i=j+1)
44. {
45. for(;height[i]<mid && i<=len;i++);
46. for(j=i;height[j]>=mid;j++);
47. if(j-i+1<n) continue;
48. ii++;s=0;
49. for(k=i-1;k<j;k++)
50. if((t=who[sa[k]])!=0)
51. if(yes[t]!=ii) yes[t]=ii,s++;
52. if(s>=n) return 1;
53. }
54. return 0;
55. }
56. char st[110];
57. int main()
58. {
59. int i,j,k,min,mid,max,nn;
60. scanf("%d",&nn);
61. while(nn-->0)
62. {
63. scanf("%d",&n);
64. len=0;
65. for(i=1;i<=n;i++)
66. {
67. scanf("%s",st);
68. k=strlen(st);
69. for(j=0;j<k;j++)
70. {
71. r[j+len]=st[j]+200;
72. who[j+len]=i;
73. }
74. r[len+k]=2*i-1;
75. who[len+k]=0;
76. len+=k+1;
77. for(j=0;j<k;j++)
78. {
79. r[j+len]=st[k-1-j]+200;
80. who[j+len]=i;
81. }
82. r[len+k]=2*i;
83. who[len+k]=0;
84. len+=k+1;
85. }
86. len--;
87. r[len]=0;
88. da(r,sa,len+1,328);
89. calheight(r,sa,len);
90. height[len+1]=-1;
91. min=1;max=100;
92. while(min<=max)
93. {
94. mid=(min+max)>>1;
95. if(check(mid)) min=mid+1;
96. else max=mid-1;
97. }
98. if(n==1) max=len/2;
99. printf("%d/n",max);
100. }
101. return 0;
102. }