最长公共子序列长度以及打印（动态规划入门）

POJ1458原题链接

北大郭炜老师在mooc上讲的 算法基础的例题

代码也是按照讲的写的

方法呢就是

串 A 以 i 为结尾和串 B 以 j 为结尾的最长公共子序列长度  maxlen[i][j]=max（maxlen[i][j-1]，maxlen[i-1][j]）

当A的第i-1字符个和B的第j-1个字符一样时   maxlen[i][j] = maxlen[i-1][j-1] + 1;

 

核心代码就是双重for循环里面的状态转移

if(str1[i-1]==str2[j-1])
    maxlen[i][j]=maxlen[i-1][j-1]+1;
else
    maxlen[i][j]=max(maxlen[i-1][j],maxlen[i][j-1]);

ac代码

#include<stdio.h>
#include<iostream>
#include<algorithm>
#include<string.h>
#define maxn 1001
using namespace std;
char str1[maxn],str2[maxn];
int maxlen[maxn][maxn];
int main(){
    while(cin>>str1>>str2){
        int l1=strlen(str1);
        int l2=strlen(str2);
        for(int i=0;i<=l2;i++)maxlen[0][i]=0;
        for(int i=0;i<=l1;i++)maxlen[i][0]=0;
        for(int i=1;i<=l1;i++){
            for(int j=1;j<=l2;j++){
                if(str1[i-1]==str2[j-1])
                    maxlen[i][j]=maxlen[i-1][j-1]+1;
                else
                    maxlen[i][j]=max(maxlen[i-1][j],maxlen[i][j-1]);
            }
        }
        cout<<maxlen[l1][l2]<<endl;
    }
    return 0;
}

另外 如果让你打印最长公共子序列

需要要一个回溯去寻找并且打印

比如题目 51nod1006   链接51nod1006

ac代码

#include<stdio.h>
#include<iostream>
#include<algorithm>
#include<string.h>
using namespace std;
const int maxn=1e3+7;
char str1[maxn],str2[maxn];
int maxlen[maxn][maxn];
int main(){
    scanf("%s%s",str1,str2);
    int l1=strlen(str1);
    int l2=strlen(str2);
    for(int i=1;i<=l1;i++){
        for(int j=1;j<=l2;j++){
            if(str1[i-1]==str2[j-1])maxlen[i][j]=maxlen[i-1][j-1]+1;
            else maxlen[i][j]=max(maxlen[i-1][j],maxlen[i][j-1]);
        }
    }
    int i=l1,j=l2,w=0;;
    char q[maxn];
    while(maxlen[i][j]){//回溯寻找
        if(maxlen[i-1][j]==maxlen[i][j])i--;
        else if(maxlen[i][j-1]==maxlen[i][j])j--;
        else {q[w++]=str1[i-1];j--;i--;};
    }
    for(int i=w-1;i>=0;i--)
    printf("%c",q[i]);
    return 0;
}