DP求两个字符串的编辑距离

给出两个字符串代表源串和目标串,可以对源串执行插入一个字符到任意位置,删除任意位置上的一个字符,替换任意位置上一个字符等3中操作,求最少需要执行多少次这样的步骤.

如:

AGCTAAGGCCTT

AGCTAAGGCCT

只需要在最后插入一个T就可以了.

再如:

snowy

sunny

最少需要3步:插入一个u,把o替换成n,删除w.

可以用dp[i][j]表示源串1 - i的子串转换到目标串子串1 - j所需要的最小步骤.

边界值:

dp[i][0] = i;//目标串是空串,要删光源串,最少需要i步

dp[0][j] = j;//源串是空串,要增加到目标串长度,最少需要j步

则方程可以表示为:

if (str1[i] == str2[j]){
          dp[i][j] = dp[i - 1][j - 1];
}else{
          dp[i][j] = min(dp[i][j - 1] + 1, min(dp[i - 1][j] + 1, dp[i - 1][j - 1] + 1));
}

str1[i]等于str2[j],则不需要任何操作就是dp[i - 1][j - 1];

否则dp[i][j - 1] + 1表示插入一个字符串需要的最少步骤,//首先把1 - i的字符串变成1 - j - 1的字符串,然后插入一个字符就成了1 - j长度的字符串

dp[i - 1][j]表示删除一个字符串需要的最少步骤,              //同上相反

dp[i - 1][j - 1]表示替换一个字符需要的最少步骤.            //这个好理解,就是把源串某个字符变一下

#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;

const int maxn = 10001;
int dp[maxn][maxn];
char str1[maxn], str2[maxn];
int main(){
	int n;
	scanf("%d", &n);
	while (n--){
		scanf("%s %s", str1 + 1, str2 + 1);
		int len1 = strlen(str1 + 1);
		int len2 = strlen(str2 + 1);
		for (int i = 1; i <= len1; ++i){
			dp[i][0] = i;
		}
		for (int i = 1; i <= len2; ++i){
			dp[0][i] = i;
		}
		for (int i = 1; i <= len1; ++i){
			for (int j = 1; j <= len2; ++j){
				if (str1[i] == str2[j]){
					dp[i][j] = dp[i - 1][j - 1];
				}else{
					dp[i][j] = min(dp[i][j - 1] + 1, min(dp[i - 1][j] + 1, dp[i - 1][j - 1] + 1));
				}
			}
		}
		printf("%d\n", dp[len1][len2]);
	}
	return 0;
}