动态规划———最长公共子序列（LCS）

最长公共子序列;
言简意赅就是在两个序列a，b里找同时是a,b两个序列的子序列并且相同的长度最长的那个。
我们用一个二维的dp[i][j]去表示以a[i]结尾和以b[j]结尾的最长的公共子序列的长度。
这样的话就可以很容易得到更新的dp[i][j]的方程。

我们来理解这个方程为什么是这个样子的，
dp[i][j]表示的是以a[i]结尾的子序列和以b[j]结尾的子序列的最长长度。
那么，我们来想一下dp[i][j]是从那些地方更新过来的。
一个二维数组要更新而且还是和他的前一个状态有关系则肯定只有三个地方更新得到dp[i-1][j-1],dp[i][j-1],dp[i-1][j]三个方向更新过来。
那三个状态分别表示什么呢？

dp[i-1][j-1]: 以a[i-1]结尾的子序列与b[j-1]结尾的子序列的最长公共子序列。
dp[i][j-1]: 以a[i]结尾的子序列与b[j-1]结尾的子序列的最长公共子序列。
dp[i-1][j]: 以a[i-1]结尾的子序列与b[j]结尾的子序列的最长公共子序列。

所以：
如果当前比较的两个元素是一样的则：当前的长度可以在不包含这两个元素的最长公共子序列的基础上长度+1。
但是如果当前比较元素不一样，则为了是的现在的长度任是最长的就要去继承之前的两个比较长的那一种情况。

模板题目 HDU-1159
ac code:

#include <iostream>
#include <cstdio>
#include <cstring>
#include <string>
#include <cmath>
#include <algorithm>
using namespace std;
const int maxn = 1e4+10;
int dp[maxn][maxn];
int main()
{
    string a,b;
    while(cin>>a>>b){
        dp[0][0]=0;
        int alen = a.size();
        int blen = b.size();
        for(int i=0;i<alen;i++){
            for(int j=0;j<blen;j++){
                if(a[i]==b[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[alen][blen]<<endl;
    }

}

打印出来最长的序列是什么 POJ - 2250
我们可以在确定最长长度之后往回退就好了 在相等的地方记录一下。
ac code:

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <string>
using namespace std;
const int maxn =110;
int dp[maxn][maxn];
char a[maxn][40];
char b[maxn][40];
string ans[maxn];
int main()
{
    char s[maxn];
    while(~scanf("%s",s)){
        memset(dp,0,sizeof(dp));
        dp[0][0]=0;
        int cnt_a=0,cnt_b=0,cnt=0;
        strcpy(a[cnt_a++],s);
        while(1){
            scanf("%s",s);
            if(strcmp(s,"#")==0)break;
            strcpy(a[cnt_a++],s);
        }
        while(1){
            scanf("%s",s);
            if(strcmp(s,"#")==0)break;
            strcpy(b[cnt_b++],s);
        }
        for(int i=0;i<cnt_a;i++){
            for(int j=0;j<cnt_b;j++){
                if(strcmp(a[i],b[j])==0){
                    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[cnt_a][cnt_b]<<endl;
        // for(int i=0;i<cnt_a;i++){
        //     cout<<a[i]<<endl;
        // }
        // for(int i=0;i<cnt_b;i++){
        //     cout<<b[i]<<endl;
        // }
        int i=cnt_a,j=cnt_b;
        while(dp[i][j]){
            if(strcmp(a[i-1],b[j-1])==0){
                ans[cnt++]=a[i-1];
                i--;j--;
            }
            else if(dp[i][j-1]>dp[i-1][j]){
                j--;
            }
            else{
                i--;
            }
        }
        for(int i=cnt-1;i>=0;i--){
            if(i==0)cout<<ans[i]<<endl;
            else cout<<ans[i]<<" ";
        }
    }


    return 0;
}

HDU - 1503
**题意:**给你两个字符串你要找到最短的包含这两个字符串，并且字符顺序不变。
思路：这就是最长公共子序列的变形，你想一下，要使的求的的字符串最短当然要最长公共子序列因为他们一个顶两啊，
ac code:

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <string>
#include <cmath>
#include <cstring>
using namespace std;
const int maxn =110;
int dp[maxn][maxn];
int main()
{
    string a,b,ans;
    while(cin>>a>>b){
        int alen = a.size();
        int blen = b.size();
        ans="";
        dp[0][0]=0;
        for(int i=0;i<alen;i++){
            for(int j=0;j<blen;j++){
                if(a[i]==b[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]);
            }
        }
        int i=alen,j=blen;
        while(dp[i][j]){
            if(a[i-1]==b[j-1]){
                ans+=a[i-1];
                i--,j--;
            }
            else if(dp[i][j-1]>dp[i-1][j]){
                ans+=b[j-1];
                j--;
            }
            else {
                ans+=a[i-1];
                i--;
            }
        }
        while(i!=0){
            ans+=a[i-1];
            i--;
        }
        while(j!=0){
            ans+=b[j-1];
            j--;
        }
        for(int i=ans.size()-1;i>=0;i--)
        cout<<ans[i];
        cout<<endl;
    }
    return 0;
}

如有错误，望评论或私信指正，感激不尽。