最长公共子串

POJ1458……

http://acm.pku.edu.cn/JudgeOnline/problem?id=1458

 

算法思想如下：

 

设，A={a1,a2,a3...ax}和B={b1,b2,b3...by}是两串，B'={bi,bj...bk};且B'是串B的子集，B'串是a串的子串。则，B'中每个元素必出现在A串中，即，对某bi有bi==aj（i<=x&&j<=y）。且，A中与B'中对应的元素出现的先后顺序与之在B'中出现的先后顺序相同。也就是说，B'中的每一个元素均出现在A中，且在A中出现的先后顺序与其在B'中出现的先后顺序相同。

这样，就可以先把A中的每个元素出现的位置保存起来。因为A中可能存在重复的元素，所以使用了二维数行代表每个字母，列中的数字是其在串A中出现的位置。对某一字母在A中出现的位置，用链表连接起来，表头是A【0】，其值是当前字母在A中出现的第一个位置。若A[0]==0，则表示该字母不出现于A中。

读入B。对B中的每个字母，查看其是否出现在A中。若出现，则取其第一次在A中出现的位置，与B中前一个出现在A中的字母的位置i相比，判断其位置是否合法（该位置在i之后）；若合法，则说明其是A,B的公共字串。否则，应该在关于A的二维数组中去掉当前字母的位置i（因为其非法i后可能还有其他位置）。

 

代码：

 

#include <stdio.h> #include <stdlib.h> #define MAXROWS 26 //行为字母 #define MAXCOLS 20 //列为该字母在所给字符串中出现的位置 /*START：POJ-1458(最长公共字串)*/ int main(void) { int buff[MAXROWS][MAXCOLS]; int m,n; int i,cur,pre,len; char tmp[MAXCOLS]; while(scanf("%s",&tmp[1])!=EOF) { for(m=0;m<MAXROWS;m++) buff[m][0]=0; for(i=1;tmp[i]!='/0';i++){ m=tmp[i]-'a'; for(n=0;buff[m][n]!=0;n=buff[m][n]) ; buff[m][n]=i; buff[m][i]=0; }//for scanf("%s",&tmp[1]); len=0; pre=0; for(i=1;tmp[i]!='/0';i++){ m=tmp[i]-'a'; if(buff[m][0]!=0){ cur=buff[m][0]; if(cur>pre){ len++; buff[m][0]=buff[m][cur]; pre=cur; } else buff[m][0]=0; } } printf("%d",len); }//while return EXIT_SUCCESS; } /*END*/

 

该算法是线性算法。时间复杂度为O(m+n)。程序要做的是读入两个串并对其线性处理。处理过程中不需要回溯等。

比较初级的算法是利用回溯的算法。http://blog.csdn.net/jqandjq/archive/2009/03/24/4021302.aspx时间复杂度应该是O(m*n)

 

有类似的题目，如：最大子序列和等。如下：

 

#include <stdio.h> #include <stdlib.h> #include <limits.h> /*START：求最大子序列和*/ int main(void) { FILE *fp1,*fp2; int num,thisum,maxsum; fp1=fopen("input.txt","r"); if(fp1==NULL) perror("input.txt"); fp2=fopen("output.txt","w"); if(fp2==NULL) perror("output.txt"); thisum=0; maxsum=INT_MIN; while(fscanf(fp1,"%d",&num)!=EOF) { thisum+=num; if(thisum>maxsum) maxsum=thisum; if(thisum<0) thisum=0; } fprintf(fp2,"%d",maxsum); fclose(fp1); fclose(fp2); return EXIT_SUCCESS; } /*END*/

 

这是清华的一道复试题。曾经在weiss的书上看到过，当时列出了4个算法，第三个是个递归。这个号称完美的线性联机算法。原来死活看不懂……其实想想也很简单，依照上面的方法可以很容易证明其正确性。

补一句：weiss书上约定，若和为负数，则最大和为0；我这里去掉了改限制。