DP—求最长公共子序列LCS（代码）

#include <iostream>
#include <cstring>
int ans[1003][1003];//ans[i][j]表示s1长度为i，s2长度为j时的最长子序列长度
int b[1003][1003];//方向数组记录ans[i][j]的来源
using namespace std;
void print_LCS(int m,int n,string s){
    if(m == 0 || n == 0)
        return ;
    if(b[m-1][n-1] == 0) {
        print_LCS(m - 1, n - 1,s);
        cout << s[m-1];
    }
    if(b[m-1][n-1] == 1){
        print_LCS(m - 1, n, s);
    }
    if(b[m-1][n-1] == -1)
        print_LCS(m, n - 1, s);
}
int main()
{
    int m,n;//s1,s2长度
    while(cin >> m >> n){
        string s1,s2;//两个字符串s1，s2长度分别m，n
        cin >> s1 >> s2;
        for(int i = 0;i <= m;i++)//当s1，s2其中有一个长度为零时对应的最长子序列也为0
            ans[i][0]=0;
        for(int j = 0;j <= n;j++)
            ans[0][j]=0;
        for(int i = 1;i <= m;i++)
            for(int j = 1;j <= n;j++){
                if(s1[i-1] == s2[j-1]){//如果s1，s2最后的一个字符相同那么最长子序列长度加1
                    ans[i][j] = ans[i-1][j-1] + 1;
                    b[i-1][j-1] = 0;
                }
                else {
                    if(ans[i-1][j] >= ans[i][j-1]){//如果s1，s2最后的一个字符不同，ans[i][j]等于ans[i-1][j]，ans[i][j-1]其中较大的
                        ans[i][j] = ans[i-1][j];
                        b[i-1][j-1] = 1;
                    }
                    else {
                        ans[i][j] = ans[i][j-1];
                        b[i-1][j-1] = -1;
                    }
                }
            }
        cout << ans[m][n] << endl;
        print_LCS(m, n, s1);
    }
}