LCS模板,求长度,并记录子串
1 //LCS模板,求长度,并记录子串 2 //亦可使用注释掉的那些代码,但所用空间会变大 3 #include<iostream> 4 #include<cstring> 5 #include<cmath> 6 #include<cstdlib> 7 #include<cstdio> 8 using namespace std; 9 #define N 5005 10 11 int len[N][N]; 12 char str1[N],str2[N],str3[N]; 13 int k; 14 15 int lcsLen(char *s1,int n1,char *s2,int n2)//求长度 16 { 17 for(int i=0;i<=n1;i++) 18 len[i][0]=0; 19 for(int j=0;j<=n2;j++) 20 len[0][j]=0; 21 for(int i=1;i<=n1;i++) 22 { 23 for(int j=1;j<=n2;j++) 24 { 25 if(s1[i-1]==s2[j-1]) 26 len[i][j]=len[i-1][j-1]+1; 27 else 28 len[i][j]=max(len[i-1][j],len[i][j-1]); 29 } 30 } 31 return len[n1][n2]; 32 } 33 34 void LCS(char *s1,char *s2,int i,int j)//递归求字串,存在tr3中 35 { 36 if(i==0||j==0) 37 return ; 38 if(s1[i-1]==s2[j-1]) 39 { 40 LCS(s1,s2,i-1,j-1); 41 str3[k++]=s1[i-1]; 42 } 43 else if(len[i-1][j]>=len[i][j-1]) 44 LCS(s1,s2,i-1,j); 45 else 46 LCS(s1,s2,i,j-1); 47 } 48 49 void lcs(char *s1,int n1,char *s2,int n2)//求字串 50 { 51 memset(str3,0,sizeof(str3)); 52 lcsLen(s1,n1,s2,n2); 53 k=0; 54 LCS(s1,s2,n1,n2); 55 } 56 int main() 57 { 58 int n1,n2; 59 while(cin >> str1 >> str2) 60 { 61 n1=strlen(str1); 62 n2=strlen(str2); 63 lcs(str1,n1,str2,n2); 64 cout << str3 << endl; 65 cout << lcsLen(str1,n1,str2,n2) << endl; 66 } 67 return 0; 68 }