【算法学习笔记】72.LCS 最大公公子序列 动态规划 SJTU OJ 1065 小M的生物实验1

#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;

//最长公共子序列长度 LCS dp

char dna1[1000+10];
char dna2[1000+10];
int len1 , len2;

int dp[1000+10][1000+10];
//表示dna1 从第1位到i 和 dna2 从第一位到j的 相似度


void init(){
    cin>>dna1;
    cin>>dna2;
    len1 = strlen(dna1);
    len2 = strlen(dna2);
}

int build(){
    memset(dp,0,sizeof(dp));
    for (int i = 1; i <= len1; ++i){
        for (int j= 1; j <= len2; ++j){
            bool ok = (dna1[i-1]==dna2[j-1]);
            dp[i][j]= max(
                dp[i-1][j-1] + ok,
                max(
                dp[i-1][j],
                dp[i][j-1]
                )
                );

        }
    }
    return dp[len1][len2];
}


int main(int argc, char const *argv[])
{
    init();
    cout<<build()<<endl;
    return 0;
}