51NOD最长公共子序列Lcs(动态规划+打印路径)

给出两个字符串A B，求A与B的最长公共子序列（子序列不要求是连续的）。
比如两个串为：

abcicba
abdkscab

ab是两个串的子序列，abc也是，abca也是，其中abca是这两个字符串最长的子序列。
Input
第1行：字符串A
第2行：字符串B
(A,B的长度 <= 1000)
Output
输出最长的子序列，如果有多个，随意输出1个。
Input示例
abcicba
abdkscab
Output示例
abca

求最长公共子序列直接用dp.
打印路径的话，
我们直接从后往前开始算,
如果s1[i]==s2[j],把字符压入vector,
接着就按照dp的思路你过来即可。
dp[i+1][j+1]=max(dp[i+1][j],dp[i][j+1])
那个比较大,那么就选择哪个的路径。
逆过来就是,dp[i+1][j]>dp[i][j+1],那么我们dp[i+1][j+1]就是选择dp[i+1][j]的路径，
逆过来就是--j,以此类推--i,若相等，--i或者--j都行。
知道i>=0&&j>=0不成立跳出，我们讲vector反着输出即可。

#include<iostream>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cctype>
#include<cmath>
#include<ctime>
#include<string>
#include<stack>
#include<deque>
#include<queue>
#include<list>
#include<set>
#include<map>
#include<cstdio>
#include<limits.h>
#define MOD 1000000007
#define fir first
#define sec second
#define fin freopen("/home/ostreambaba/文档/input.txt", "r", stdin)
#define fout freopen("/home/ostreambaba/文档/output.txt", "w", stdout)
#define mes(x, m) memset(x, m, sizeof(x))
#define Pii pair<int, int>
#define Pll pair<ll, ll>
#define INF 1e9+7
#define inf 0x3f3f3f3f
#define Pi 4.0*atan(1.0)

#define lowbit(x) (x&(-x))
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define max(a,b) a>b?a:b

typedef long long ll;
typedef unsigned long long ull;
const double eps = 1e-9;
const int maxn = 1000+5;
using namespace std;

inline int read(){
    int x(0),f(1);
    char ch=getchar();
    while (ch<'0'||ch>'9') {if (ch=='-') f=-1;ch=getchar();}
    while (ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
vector<char> vec;
int dp[maxn][maxn];
void solve(string s1,string s2){
    vec.clear();
    mes(dp,0);
    int n=s1.size();
    int m=s2.size();
    for(int i=0;i<n;++i){
        for(int j=0;j<m;++j){
            if(s1[i]==s2[j]){
                dp[i+1][j+1]=dp[i][j]+1;
            }else{
                dp[i+1][j+1]=max(dp[i][j+1],dp[i+1][j]);
            }
        }
    }
   //cout<<dp[n][m]<<endl;
    for(int i=n-1,j=m-1;j>=0&&i>=0;){
        if(s1[i]==s2[j]){
            vec.push_back(s1[i]);
            --i,--j;
        }else{
            if(dp[i][j+1]>dp[i+1][j]){
                --i;
            }else{
                --j;
            }
        }
    }
    //reverse(vec.begin(),vec.end());
    for(int i=vec.size()-1;i>=0;--i){
        printf("%c",vec[i]);
    }
    printf("\n");
}

int main(){
    string s1,s2;
    cin>>s1>>s2;
    solve(s1,s2);
}