动态规划；最长公共子序列；LCS问题；

#include <iostream> using namespace std; class LCS { private: char *mString;//字符串m char *nString;//字符串n int **lcs;//降低空间复杂度，利用2*n维数组记录lcs[1:m][1:n]当前与前一行的状况 int **record;//为了能够从顶至下的构造出LCS进行的记录 int mLength;//m字符串的长度 int nLength;//n字符串的长度 public: //构造函数 LCS(int m,int n) { mLength=m; nLength=n; mString=new char[m+1]; nString=new char[n+1]; lcs=new int* [2]; for(int j=0;j<2;j++) { lcs[j]=new int[n+1]; } record=new int* [m+1]; for(int i=0;i<m+1;i++) { record[i]=new int[n+1]; } } //输入字符串 void input() { cout<<"输入字符串m"<<endl; cin>>mString; cout<<"输入字符串n"<<endl; cin>>nString; } //LCS-Length动态规划算法 void LCSLength() { //初始化lcs[0][j],lcs[i][0]均为0,即没有公共字串. for(int i=0;i<2;i++) { for(int l=0;l<nLength+1;l++) { lcs[i][l]=0; } } int zero=1;//用于处理仅用两行数组lcs的记录 for(int j=1;j<=mLength;j++) { for(int k=1;k<=nLength;k++) { if(mString[j]==nString[k]) { lcs[zero][k]=lcs[zeroReturn(zero)][k-1]+1; record[j][k]=0; } else if(lcs[zero][k-1]>=lcs[zeroReturn(zero)][k]) { lcs[zero][k]=lcs[zero][k-1]; record[j][k]=1; } else { lcs[zero][k]=lcs[zeroReturn(zero)][k]; record[j][k]=-1; } } zero=zeroReturn(zero); } cout<<lcs[1][nLength]<<endl; } //0,1取反操作 int zeroReturn(int i) { if(i==0) { return 1; } else { return 0; } } //调用printLCS void print() { printLCS(mLength,nLength); } //递归的根据record[][]打印LCS void printLCS(int m,int n) { if(m==0||n==0) { return; } else//根据record标记确定是否打印当前字符与判断向哪个方向前进 { if(record[m][n]==0) { printLCS(m-1,n-1); cout<<mString[m]<<" "; } else if(record[m][n]==1) { printLCS(m,n-1); } else { printLCS(m-1,n); } } } }; void main() { LCS test(5,3);//m字符串5个，n字符串3个 test.input();//输入字符串 test.LCSLength();//求LCS test.print();//打印结果 }

 

lcs[i][j]表示x字符串的前i个与y字符串的前j个的最长公共子序列

lcs[i][j]=0       i=0或j=0

lcs[i][j]=lcs[i-1][j-1]+1         x_i=y_i

lcs[i][j]=max(lcs[i-1][j],lcs[i][j-1])     x_i!=y_j

 

经过分析，i,j前的最长公共子序列由lcs[i-1][j-1],lcs[i-1][j],lcs[i][j-1],这些子问题决定，然后子问题也是最长的子序列，即这个递归式具有最优子结构性质，每一层递归都是最优的，并且由最优的子问题决定，并且子问题的最优值不受父问题影响，即无后效性，这就是能够自底向上的进行递推的重要依据。

另外一个使用动态规划的原因是由于递归过程中有许多重叠子问题，例如求lcs[i-1][j]时可能会计算到lcs[i-1][j-1],而求lcs[i][j-1]时也可能需要计算lcs[i-1][j-1].  所以，为了避免重复计算，我们应该制作一张表用来记录已经得到的结果，然后我们能够这样做的原因，就是因为递归的表达式具有最优子结构性质，即用到的都是最优的，计算出的也都是最优的，再次使用一个运算过的问题时，直接调用表中记录的这个最优值就可以了。