字符串的相似度

许多程序会大量使用字符串。对于不同的字符串，我们希望能够有办法判断其相似程序。我们定义一套操作方法来把两个不相同的字符串变得相同，具体的操作方法为：

　　1.修改一个字符（如把“a”替换为“b”）;

　　2.增加一个字符（如把“abdd”变为“aebdd”）;

　　3.删除一个字符（如把“travelling”变为“traveling”）;

　　比如，对于“abcdefg”和“abcdef”两个字符串来说，我们认为可以通过增加/减少一个“g”的方式来达到目的。上面的两种方案，都仅需要一次 。把这个操作所需要的次数定义为两个字符串的距离，而相似度等于“距离+1”的倒数。也就是说，“abcdefg”和“abcdef”的距离为1，相似度为1/2=0.5。

给定任意两个字符串，你是否能写出一个算法来计算它们的相似度呢？

思路：经典的动态规划方法

一：从前往后

记Ai为字符串A(a1a2a3 … am)的前i个字符（即为a1,a2,a3 … ai）

记Bj为字符串B(b1b2b3 … bn)的前j个字符（即为b1,b2,b3 … bj）

记 f(i , j)为使两个字符串和Ai和Bj相等的最小操作次数，则

当ai==bj时，f(i, j)=f(i-1, j-1)

当ai!=bj时， 

    若将他们修改为相等，则至少还需要f(i-1, j-1)次操作；

    若删除ai或在Bj后添加ai，则至少还要操作f(i-1, j)次；

    若删除bj或在Ai后添加bj，则至少还要操作f(i, j-1)次；

    也即，f(i, j) = min(f(i-1, j-1), f(i-1, j), f(i, j-1))+1

#include <iostream>
#include <string>
using namespace std;

const int N = 100;

int minValue(int va, int vb, int vc)
{
	int temp = va;

	if(vb < temp)
		temp = vb;
	if(vc < temp)
		temp = vc;

	return temp;
}

int main()
{
	string A;
	string B;

	cin>>A;
	cin>>B;

	int lenA = A.length();
	int lenB = B.length();

	int dp[N][N];

	for (int i=0; i<=lenA; i++)
		dp[i][0] = i;
	for (int j=0; j<=lenB; j++)
		dp[0][j] = j;

	for (int i=1; i<=lenA; i++)
	{
		for (int j=1; j<=lenB; j++)
		{
			if(A[i-1]==B[j-1])
				dp[i][j] = dp[i-1][j-1];
			else
				dp[i][j] = minValue(dp[i-1][j], dp[i][j-1], dp[i-1][j-1])+1;
		}
	}

	for (int i=0; i<=lenA; i++)
	{
		for (int j=0; j<=lenB; j++)
		{
			cout<<dp[i][j]<<" ";
		}
		cout<<endl;
	}

	return 0;
}

优化：由于f(i,j)的计算只需要前一列和当前列的数据，所以可以进行进一步的优化

#include <iostream>
#include <string>
#include <vector>
using namespace std;

const int N = 100;

int minValue(int va, int vb, int vc)
{
	int temp = va;

	if(vb < temp)
		temp = vb;
	if(vc < temp)
		temp = vc;

	return temp;
}

int main()
{
	string A;
	string B;

	cin>>A;
	cin>>B;

	int lenA = A.length();
	int lenB = B.length();

	vector<int> dp(lenA+2);
	for (int i=0; i<=lenA; i++)
		dp[i+1] = i;

	for (int i=1; i<=lenB; i++)
	{
		dp[0] = i;

		for (int j=1; j<=lenA; j++)
		{
			if (A[j-1]!=B[i-1])
				dp[j] = minValue(dp[j-1], dp[j], dp[j+1])+1;
		}

		for (int j=lenA; j>=0; j--)
			dp[j+1] = dp[j];
	}
	
	cout<<dp[lenA+1];

	return 0;
}

二：从后往前

若ai==bj， 则f(i, j)=f(i+1, j+1)

若ai!=bj，则f(i, j)=min{f(i+1, j+1), f(i+1, j), f(i, j+1)}+1

#include <iostream>
#include <string>
using namespace std;

const int N = 100;

int minValue(int va, int vb, int vc)
{
	int temp = va;

	if(vb < temp)
		temp = vb;
	if(vc < temp)
		temp = vc;

	return temp;
}

int main()
{
	string A;
	string B;

	cin>>A;
	cin>>B;

	int lenA = A.length();
	int lenB = B.length();

	int dp[N][N];

	for(int i=lenB; i>=0; i--)
		dp[lenA][i] = lenB - i;
	for(int i=lenA; i>=0; i--)
		dp[i][lenB] = lenA - i;

	for (int i=lenA-1; i>=0; i--)
	{
		for (int j=lenB-1; j>=0; j--)
		{
			if(A[i]==B[j])
			dp[i][j] = dp[i+1][j+1];
			else
			dp[i][j] = minValue(dp[i+1][j], dp[i][j+1], dp[i+1][j+1])+1;

		}
	}

	

	for (int i=0; i<=lenA; i++)
	{
		for (int j=0; j<=lenB; j++)
		{
			cout<<dp[i][j]<<" ";
		}
		cout<<endl;
	}

	return 0;
}