动态规划基础篇之最长公共子序列问题

状态转移方程：dp[i][j] = (a[i-1]==b[j-1] ? dp[i-1][j-1] + 1 : max(dp[i-1][j] , dp[i][j-1])/记录最长公共子序列的长度；

另用一个数组c[]根据dp[][]记录的最大子序列的路径从a[]或者b[]中存下属于最长公共子序列中的每个字符，最后再反向输出。

题目链接：http://www.51nod.com/Tutorial/Chapter.html#chapterId=4

#include<iostream>
#include<stdio.h>
#include<algorithm>
#include<string.h>
using namespace std;
char a[1005],b[1005];
int dp[1005][1005];
long long lena,lenb;

void lcs()
{
    for(int i=1;i<=lena;i++)
    {
        for(int j=1;j<=lenb;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]);
        }
    }
}
void print()
{
    char c[1005];
    memset(c,0,sizeof(c));
    long long i=lena,j=lenb,z=0;
    while(i!=0&&j!=0)
    {
        if(a[i-1]==b[j-1])
        {
            c[z++]=a[--i];
            j--;
        }
        else if(dp[i-1][j]<dp[i][j-1]) j--;
        else if(dp[i][j-1]<=dp[i-1][j]) i--;
    }
    for(long long k=z-1;k>=0;k--)
    {
        printf("%c",c[k]);
    }
    printf("\n");
}
int main()
{
    memset(dp,0,sizeof(dp));
    scanf("%s",a);
    scanf("%s",b);
    lena = strlen(a);
    lenb = strlen(b);
    lcs();
    print();
    //cout<<dp[lena][lenb]<<endl;
    return 0;
}