poj 3356 AGTC(dp,最小编辑距离)

poj 3356 AGTC

#include<stdio.h>
#include<string.h>
#define MAXN 1005
char s1[MAXN],s2[MAXN];
int dp[MAXN][MAXN];
int min(int a,int b) {return a<b?a:b;}
int main()
{
    int x,y;
    for(int i=0;i<MAXN;i++) dp[i][0]=dp[0][i]=i;
    while(scanf("%d%s",&x,s1+1)!=EOF)
    {
        scanf("%d%s",&y,s2+1);
        for(int i=1;i<=x;i++)
            for(int j=1;j<=y;j++)
            {
                dp[i][j]=min(dp[i-1][j]+1,dp[i][j-1]+1);
                dp[i][j]=min(dp[i][j],dp[i-1][j-1]+(s1[i]==s2[j]?0:1));
            }
        printf("%d\n",dp[x][y]);
    }
    return 0;
}