最长公共子序列Common Subsequence POJ1458

#include <stdio.h>
#include <string.h>
int max(int a,int b)
{
    return a<b?b:a;
}
int main()
{
    int dp[1005][1005];
    char a[1005],b[1005];
    int i,j,k;
    int lena,lenb;
    int maxn;
    while(scanf("%s%s",a+1,b+1)!=-1)
    {
        maxn=0;
        //printf("%s %s\n",a+1,b+1);
        lena=strlen(a+1);
        lenb=strlen(b+1);
        memset(dp,0,sizeof(dp));
        for(j=1;j<=lena;j++)
        {
            for(k=1;k<=lenb;k++)
            {
                if(a[j]==b[k])
                    dp[j][k]=dp[j-1][k-1]+1;
                else
                    dp[j][k]=max(dp[j-1][k-1],max(dp[j-1][k],dp[j][k-1]));
                if(dp[j][k]>maxn)
                    maxn=dp[j][k];
            }
        //printf("1\n");
        }
        printf("%d\n",maxn);
    }
    return 0;
}