HDU 1159 Common Subsequence 公共子序列 简单DP

传送门：HDU 1159 Common Subsequence

分析见注释，入门DP。很经典。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
#define INF 0x3f3f3f3f
#define ms(x, y) memset(x, y, sizeof(x))
const double PI = acos(-1.0), eps = 1e-8;
// 不得不说是很经典的入门DP
int m[2001][2001],i,j,l1,l2;
int main(){
    char s1[2001],s2[2001];
    // 从 数组下标为1开始输入字符串 ---- 即&s[1]为数组的首地址
    // 为什么一定要从1开始输入呢？ 防止后面的状态转移方程越界
    // 统一成逻辑下标
    while(scanf("%s%s",s1+1,s2+1) != EOF) {
        l1=strlen(s1+1),l2=strlen(s2+1);
        ms(m, 0);
        for(i=1;i<=l1;i++) {
            for(j=1;j<=l2;j++) {
                // S1[i]!= s2[j]时，MaxLen(S1,S2)不会比MaxLen(S1,S2 j-1 )和MaxLen(S1 i-1 ,S2)两者之中任何一个小，也不会比两者都大。
                if(s1[i]==s2[j])
                    m[i][j]=m[i-1][j-1]+1; //如果是从0开始就越界了！
                else
                    // 对于邻近的状态 选择较大的 转移状态给现有状态
                    // 本身m[i][j]=0 需要从邻近的状态 转移状态 更新状态
                    m[i][j]=max(m[i-1][j],m[i][j-1]); // m[i][j] 状态可从 m[i-1][j]  m[i][j-1]这两种状态转移过来
            }
        }
        printf("%d\n",m[l1][l2]);
    }
    return 0;
}