最长公共子序列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;
}