51nod1006最长公共子序列LCS

51nod1006最长公共子序列LCS_第1张图片

输出路径的过程
结合dp数组和代码一看就懂了
  a b d k s c a b
a 1 1 1 1 1 1 1 1 
b 1 2 2 2 2 2 2 2 
c 1 2 2 2 2 3 3 3 
i 1 2 2 2 2 3 3 3 
c 1 2 2 2 2 3 3 3 
b 1 2 2 2 2 3 3 4 
a 1 2 2 2 2 3 4 4 

经典的dp题,关键题目让输出最后的公共序列,可以通过判断dp数组来从后往前找到状态变化的字母,最后反序输出

状态转移方程:
d p [ i ] [ j ] = m a x ( d p [ i − 1 ] [ j ] , d p [ i ] [ j − 1 ] , m a t [ i − 1 ] [ j − 1 ] + ( a [ i ] = = b [ j ] ? 1 : 0 ) ) ; dp[i][j] = max(dp[i-1][j], dp[i][j-1], mat[i-1][j-1] + (a[i]==b[j] ? 1 : 0)); dp[i][j]=max(dp[i1][j],dp[i][j1],mat[i1][j1]+(a[i]==b[j]?1:0));

#include 
#define INF 0x3f3f3f3f
#define d(x) cout << (x) << endl
using namespace std;
const int mod = 1e9 + 7;
typedef long long ll;
const int N = 1e3 + 10;
const int M = 1e4 + 10;

int dp[N][N];
int main()
{
    string a, b, c = "";
    cin >> a >> b;
    for (int i = 1; i <= a.size(); i++) {
        for (int j = 1; j <= b.size(); 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]);
            }
        }
    }
    int i = a.size() - 1;
    int j = b.size() - 1;
    while(i >= 0 && j >= 0){
        if(a[i] == b[j]){
            c += a[i];
            i--;
            j--;
        }else if(dp[i][j+1] > dp[i+1][j]){
            i--;
        }else{
            j--;
        }
    }
    reverse(c.begin(), c.end());
    cout << c << endl;
    return 0;
}

你可能感兴趣的:(解题报告,51NOD,-----动态规划-----,基础dp)