对于一般的LCS问题,都属于NP问题。当数列的量为一定的时,都可以采用动态规划去解决。
动态规划的一个计算最长公共子序列的方法如下,以两个序列 X、Y 为例子:
设有二维数组 f[i][j] 表示 X 的 i 位和 Y 的 j 位之前的最长公共子序列的长度,则有:
f[1][1] = same(1,1)
f[i][j] = max{f[i-1][j-1] + same(i,j),f[i-1][j],f[i][j-1]}
其中,same(a,b)当 X 的第 a 位与 Y 的第 b 位完全相同时为“1”,否则为“0”。
此时,f[i][j]中最大的数便是 X 和 Y 的最长公共子序列的长度,依据该数组回溯,便可找出最长公共子序列。
该算法的空间、时间复杂度均为O(n^2)。
1 #include <iostream> 2 #include <fstream> 3 #include <cstring> 4 #include <string> 5 #include <cmath> 6 using namespace std; 7 const int MAX=1020; 8 int dp[MAX][MAX];//保存信息 9 string str_1,str_2; 10 int main() 11 { 12 #ifndef ONLINE_JUDGE 13 freopen("F://code//txt//25.txt","r",stdin); 14 #endif 15 while(getline(cin,str_1)&&str_1.size()) 16 { 17 getline(cin,str_2); 18 memset(dp,0,sizeof dp); 19 for(int i=1;i<=str_1.size();i++) 20 { 21 for(int j=1;j<=str_2.size();j++) 22 { 23 if(str_1[i-1]==str_2[j-1]) 24 dp[i][j]=dp[i-1][j-1]+1; 25 else 26 dp[i][j]=max(dp[i-1][j],dp[i][j-1]); 27 } 28 } 29 cout<<dp[str_1.size()][str_2.size()]<<endl; 30 } 31 32 return 0; 33 }