動態規劃法求解最長公共子序列問題
阿新 • • 發佈:2020-12-20
技術標籤:演算法
狀態轉移方程如下:
#include<stdio.h>
#include<vector>
#include<string.h>
#include<iostream>
using namespace std;
#define max(x,y) (x)>(y)?(x):(y)
#define MAX 51
int m,n;
string a,b;
int dp[MAX][MAX];
vector<char>subs;
void LCSlength(){
for(int i=0;i<=m;i++) {
dp[i][0]=0;
}
for(int i=0;i<=n;i++){
dp[0][i]=0;
}
for(int i=1;i<=m;i++){
for(int j=1;j<=n;j++){
if(a[i-1]==b[j-1]) {
dp[i][j]=dp[i-1][j-1]+1;
}
else {
dp[i][j]=max(dp[i][j-1],dp[i-1][j]);
}
}
}
}
void Buildsubs(){
int k=dp[m][n];
int i=m;
int j=n;
while(k> 0){
if(dp[i][j]==dp[i-1][j]) i--;
else if(dp[i][j]==dp[i][j-1]) j--;
else {
subs.push_back(a[i-1]);
i--;j--;k--;
}
}
}
int main(){
a="cccccffaagdsfsfcccc";
b="wwwdfsifdfosi";
m=a.length();
n=b.length();
LCSlength();
Buildsubs();
}