poj 3356 AGTC 经典动态规划题

详细讲解可见

挑战编程 11.2编辑距离所讲解的，非常充分完整。

状态转移方程dp[i][j]定义为使第一个串前i长度与第二串前j相同的最小带价

如果X[i] == Y[i] 则dp[i][j] =  min3(dp[i-1][j-1],dp[i][j-1]+1,dp[i-1][j]+1)

不然dp[i][j] = min3(dp[i-1][j-1]+1,dp[i][j-1]+1,dp[i-1][j]+1)
dp[i-1][j-1]+1 ，1是通过替换是第i,j字符相等的带价

,dp[i][j-1]+1/dp[i-1][j]+1 是通过插入/删除使第i,j个字符相等的带价

Source Code
Problem: 3356	
User: csjiaxin
Memory: 4668K	
Time: 47MS
Language: G++	
Result: Accepted
Source Code
#include <stdio.h>
#include <iostream>

#include <cstring>

using namespace std;

#define MAX_N 1000
char X[MAX_N+1],Y[MAX_N+1];
int dp[MAX_N+1][MAX_N+1];

int min3(int a,int b,int c)
{
	int ret = a;
	if(ret > b)
		ret = b;
	if(ret > c)
		ret =c;

	return ret;
}
#define printf //
int main()
{
	int m,n;
	while(cin>>m)
	{
		memset(dp,0,sizeof(dp));
		for(int i=1;i<=m;i++)
		{
			cin>>X[i];
		}
		cin>>n;
		for(int i=1;i<=n;i++)
			cin>>Y[i];
	
		for(int i=0;i<=m;i++)
			dp[i][0] = i;
	
		for(int j=0;j<=n;j++)
			dp[0][j] =j;
	
		for(int i=1;i<=m;i++)
		{
			for(int j=1;j<=n;j++)
			{
				dp[i][j] = min3(dp[i-1][j-1]+1,dp[i-1][j]+1,dp[i][j-1]+1);
	
				if(X[i] == Y[j])
					dp[i][j] = min(dp[i][j],dp[i-1][j-1]);
	
				printf("dp[%d][%d] =%d\t",i,j,dp[i][j]);
			}
			printf("\n");
		}
	
		cout<<dp[m][n]<<endl;
	}
}