动态规划-最长公共子序列
2012-11-11 14:40
211 查看
问题描述
某个序列的子序列是从最初序列通过去除某些元素但不破坏余下元素的相对位置而形成的新序列。
给定两个序列X和Y,如果Z既是X的一个子序列又是Y的一个子序列,则称Z是X,Y的公共子序列。
给定两个序列X={x1,x2...,xm}和Y={y1,y2...yn},找出Z={z0,z1,…,zk}为它们的最长公共子序列.
最优子结构
设X={x1,x2...,xm}和Y={y1,y2...yn}为两个序列,Z={z0,z1,…,zk}为任意一个LCS。
1)如果xm=ym,那么zk=xm=ym,且Zk-1为Xm-1和Yn-1的一个LCS。
2)如果xm!=ym&&zk!=xm,那么Z是Xm-1和Y的一个LCS。
3)如果xm!=ym&&zk!=yn,那么Z是X和Yn-1的一个LCS。
状态转移方程
c[i,j]为序列Xi和Yi的一个LCS的长度。
代码实现
//LCS
#include<iostream>
#include<string>
using namespace std;
void LCS_Length(string x,string y,int (*c)[100],int (*b)[100])
{
int m,n;
m=x.length();
n=x.length();
int i,j;
//如果i或j等于0则c[i][j]=0;
for(i=1;i<=m;i++)
{
c[i][0]=0;
}
for(i=1;i<=n;i++)
{
c[0][i]=0;
}
//遍历两个字符串,依次标记c[i][j],c[i][j]标记了从x的开始到第i个元素与从
//y开始到第j个元素中LCS的长度
//数组b用来标记最长公共子串所要走的路线,该路线为两个字符串组成的矩阵中的对应的字母
for(i=1;i<=m;i++)
{
for(j=1;j<=n;j++)
{
if(x[i-1]==y[j-1])
{
c[i][j]=c[i-1][j-1]+1;
//代替箭头指向左上
b[i][j]=1;
}
else
{
if(c[i][j-1]>=c[i-1][j])
{
c[i][j]=c[i][j-1];
//代替箭头指左
b[i][j]=2;
}
else
{
c[i][j]=c[i-1][j];
//代替箭头指向上
b[i][j]=3;
}
}
}
}
}
void PrintAnswer(string x,int(*b)[100],int i,int j)
{
if(i==0||j==0)
{
return ;
}
else
{
if(b[i][j]==1)
{
PrintAnswer(x,b,i-1,j-1);
cout<<x[i-1]<<" ";
}
else if(b[i][j]==2)
{
PrintAnswer(x,b,i,j-1);
}
else
{
PrintAnswer(x,b,i-1,j);
}
}
}
int main()
{
string x="abcbda";
string y="bdcaba";
int c[100][100]={0};
int b[100][100]={0};
LCS_Length(x,y,c,b);
cout<<"the LCS is: "<<c[x.length()][y.length()]<<endl;
PrintAnswer(x,b,x.length(),y.length());
return 0;
}
参考:最长公共子串LCS代码实现-http://blog.csdn.net/liuzhanchen1987/article/details/7851413
某个序列的子序列是从最初序列通过去除某些元素但不破坏余下元素的相对位置而形成的新序列。
给定两个序列X和Y,如果Z既是X的一个子序列又是Y的一个子序列,则称Z是X,Y的公共子序列。
给定两个序列X={x1,x2...,xm}和Y={y1,y2...yn},找出Z={z0,z1,…,zk}为它们的最长公共子序列.
最优子结构
设X={x1,x2...,xm}和Y={y1,y2...yn}为两个序列,Z={z0,z1,…,zk}为任意一个LCS。
1)如果xm=ym,那么zk=xm=ym,且Zk-1为Xm-1和Yn-1的一个LCS。
2)如果xm!=ym&&zk!=xm,那么Z是Xm-1和Y的一个LCS。
3)如果xm!=ym&&zk!=yn,那么Z是X和Yn-1的一个LCS。
状态转移方程
c[i,j]为序列Xi和Yi的一个LCS的长度。
代码实现
//LCS
#include<iostream>
#include<string>
using namespace std;
void LCS_Length(string x,string y,int (*c)[100],int (*b)[100])
{
int m,n;
m=x.length();
n=x.length();
int i,j;
//如果i或j等于0则c[i][j]=0;
for(i=1;i<=m;i++)
{
c[i][0]=0;
}
for(i=1;i<=n;i++)
{
c[0][i]=0;
}
//遍历两个字符串,依次标记c[i][j],c[i][j]标记了从x的开始到第i个元素与从
//y开始到第j个元素中LCS的长度
//数组b用来标记最长公共子串所要走的路线,该路线为两个字符串组成的矩阵中的对应的字母
for(i=1;i<=m;i++)
{
for(j=1;j<=n;j++)
{
if(x[i-1]==y[j-1])
{
c[i][j]=c[i-1][j-1]+1;
//代替箭头指向左上
b[i][j]=1;
}
else
{
if(c[i][j-1]>=c[i-1][j])
{
c[i][j]=c[i][j-1];
//代替箭头指左
b[i][j]=2;
}
else
{
c[i][j]=c[i-1][j];
//代替箭头指向上
b[i][j]=3;
}
}
}
}
}
void PrintAnswer(string x,int(*b)[100],int i,int j)
{
if(i==0||j==0)
{
return ;
}
else
{
if(b[i][j]==1)
{
PrintAnswer(x,b,i-1,j-1);
cout<<x[i-1]<<" ";
}
else if(b[i][j]==2)
{
PrintAnswer(x,b,i,j-1);
}
else
{
PrintAnswer(x,b,i-1,j);
}
}
}
int main()
{
string x="abcbda";
string y="bdcaba";
int c[100][100]={0};
int b[100][100]={0};
LCS_Length(x,y,c,b);
cout<<"the LCS is: "<<c[x.length()][y.length()]<<endl;
PrintAnswer(x,b,x.length(),y.length());
return 0;
}
参考:最长公共子串LCS代码实现-http://blog.csdn.net/liuzhanchen1987/article/details/7851413
相关文章推荐
- 算法知识之最长公共子序列问题(动态规划)
- 使用动态规划解决最长公共子序列问题
- 动态规划之----最长公共子序列
- 【算法】 动态规划 最长公共子序列
- 动态规划——《3》最长公共子序列
- 动态规划——最长公共子序列问题(LCS)
- 动态规划-最长公共子序列计算
- 10635 - Prince and Princess(最长公共子序列优化,动态规划)
- 动态规划求解最长公共子序列问题
- 【动态规划】 之最长公共子序列LCS
- 动态规划——最长公共子序列长度
- 【动态规划】【最长公共子序列】Vijos P1111 小胖的水果
- 【动态规划】【最长公共子序列】Vijos P1111 小胖的水果
- LCS问题(最长公共子序列)-动态规划实现
- POJ 1458 Common Subsequence(最长公共子序列)__动态规划
- 动态规划——最长公共子序列(算法设计课题)
- 最长公共子序列问题(Java)-动态规划
- 最长公共子序列(动态规划求解)
- 动态规划之最长公共子序列
- 《算法导论》学习笔记——最长公共子序列(动态规划)