NYoj 36 最长公共子序列[典型动态规划]

/*  状态转移方程:dp[i][j]={dp[i-1][j-1]+1,if(a[i]==b[j])                        max(dp[i-1][j],dp[i][j-1]])，if(a[i]!=b[j])}   表示a的长度为i，b的长度为j时，其最大公共子序列.   动态规划，要理解他的精髓还是要多做题练习啊.   细解:dp[i][j]={dp[i-1][j-1]+2(s1[i]==s2[j])                  max(dp[i-1][j],dp[i][j-1])(s1[i]!=s2[j])}     dp[i][j]表示s1前i和和s2的前j个的最大公共子序列，     他一来于其前一个状态。     如果s1[i]==s2[j],那么他的状态就是上一个状态+1；     否则他当前的最佳状态取决于其前一个状态的最佳.*/#include<stdio.h>#include<string.h>#include<algorithm>#define max(a,b) a>b?a:busing namespace std;int dp[1005][1005];int LCSlenth(char a[],char b[]){    int lena=strlen(a),lenb=strlen(b);    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]);        }    }    return dp[lena][lenb];}int main(){    int T;    scanf("%d",&T);    while(T--)    {        char A[1005],B[1005];        scanf("%s%s",A,B);        memset(dp,0,sizeof(dp));        printf("%d\n",LCSlenth(A,B));    }}