动态规划之数组dp简单学习（以51nod 1006 最长公共子序列Lcs为例）

说明

参考博文：https://blog.csdn.net/godjing007/article/details/89048554

状态转移方程
dp[i][j]=max(dp[i-1][j],dp[i][j-1],dp[i-1][j-1]+(a[i]==b[j]?1:0)
关键是找到矩阵中第一次增加的地方就是想要的点

代码

#include<iostream>
#include<string>
#include<stack>
#include<algorithm>
using namespace std;
#define N 1000
int dp[N][N];//动态矩阵
int main(){
string a,b;
stack<char> s;
cin>>a>>b;
int i,j;
for(i=1;i<=a.length();i++){
    for(j=1;j<=b.length();j++){
        if(a[i-1]==b[j-1])
            dp[i][j]=dp[i-1][j-1]+1;
        else
            dp[i][j]=max(dp[i-1][j],dp[i][j-1]);

    }
}
i-=2,j-=2;
while(i>=0&&j>=0){
    if(a[i]==b[j]){
            s.push(a[i]);
    i--;
    j--;}
    else if(dp[i][j+1]>dp[i+1][j]) i--;//往变化最快的地方走
    else j--;
}
char q;
while(!s.empty()){
q=s.top();
cout<<q;
s.pop();
}
}