您的位置:首页 > 其它

最长公共子序列LCS(动态规划法)

2017-07-28 10:46 190 查看
//动态规划法
//最长公共子序列
//c数组指出子序列长度,b数组为具体路径
#include <iostream>
#define N 1000
#define UPLEFT 1
#define UP 2
#define LEFT 3

int c

,b

;

using namespace std;

void PRINT_LCS(int b[]
,char X[],int i,int j)  //打印路径
{
if(i==0 || j==0) return ;
if(b[i][j] == UPLEFT)
{
PRINT_LCS(b,X,i-1,j-1);
cout << X[i] << ends;
}
else if(b[i][j] == UP) PRINT_LCS(b,X,i-1,j);
else PRINT_LCS(b,X,i,j-1);
}

LCS(char X[],char Y[],int m,int n)
{

for(int i=0;i<=m;i++) c[i][0]=0;
for(int j=0;j<=n;j++) c[0][j]=0;
for(int i=0;i<=m;i++)
{
for(int j=0;j<=n;j++) b[i][j]=0;
}

for(int i=1;i<=m;i++)
{
for(int j=1;j<=n;j++)
{
if(X[i] == Y[j])
{
c[i][j]=c[i-1][j-1]+1;
b[i][j]=UPLEFT;
}
else if(c[i-1][j] >= c[i][j-1])
{
c[i][j]=c[i-1][j];
b[i][j]=UP;
}
else
{
c[i][j]=c[i][j-1];
b[i][j]=LEFT;
}
}
}
PRINT_LCS(b,X,m,n);
}

int main(int argc, char *argv[])
{
int m,n; cin >> m >> n;

char X
,Y
;

for(int i=1;i<=m;i++) cin >> X[i];
for(int i=1;i<=n;i++) cin >> Y[i];

LCS(X,Y,m,n);
cout << endl;
cout << c[m]
<<endl;
cout << "Hello World!" << endl;
return 0;
}
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: