编程之美3.3计算字符串的相似度

推荐http://www.cnblogs.com/yujunyong/articles/2004724.html

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

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

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

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

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

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

　　原文的分析与解法　　

　　不难看出，两个字符串的距离肯定不超过它们的长度之和（我们可以通过删除操作把两个串都转化为空串）。虽然这个结论对结果没有帮助，但至少可以知道，任意两个字符串的距离都是有限的。

　　我们还是就住集中考虑如何才能把这个问题转化成规模较小的同样的子问题。如果有两个串A=xabcdae和B=xfdfa，它们的第一个字符是 相同的，只要计算A[2,...,7]=abcdae和B[2,...,5]=fdfa的距离就可以了。但是如果两个串的第一个字符不相同，那么可以进行 如下的操作（lenA和lenB分别是A串和B串的长度）。

　1.删除A串的第一个字符，然后计算A[2,...,lenA]和B[1,...,lenB]的距离。

　　2.删除B串的第一个字符，然后计算A[1,...,lenA]和B[2,...,lenB]的距离。

　　3.修改A串的第一个字符为B串的第一个字符，然后计算A[2,...,lenA]和B[2,...,lenB]的距离。

　　4.修改B串的第一个字符为A串的第一个字符，然后计算A[2,...,lenA]和B[2,...,lenB]的距离。

　　5.增加B串的第一个字符到A串的第一个字符之前，然后计算A[1,...,lenA]和B[2,...,lenB]的距离。

　　6.增加A串的第一个字符到B串的第一个字符之前，然后计算A[2,...,lenA]和B[1,...,lenB]的距离。

　　在这个题目中，我们并不在乎两个字符串变得相等之后的字符串是怎样的。所以，可以将上面的6个操作合并为：

　　1.一步操作之后，再将A[2,...,lenA]和B[1,...,lenB]变成相字符串。

　　2.一步操作之后，再将A[2,...,lenA]和B[2,...,lenB]变成相字符串。

　　3.一步操作之后，再将A[1,...,lenA]和B[2,...,lenB]变成相字符串。

　　这样，很快就可以完成一个递归程序。

#include <iostream>
#include <string>
using namespace std;
int min(int a,int b,int c){
	int t = a>b?b:a;
	t = t>c?c:t;
	return t;
}
int cal(string stra,int abegin,int aend,string strb,int bbegin,int bend){
	if(abegin>aend){
		if(bbegin>bend)
			return 0;
		else
			return bend-bbegin+1;
	}
	if(bbegin>bend){
		if(abegin>aend)
			return 0;
		else
			return aend-abegin+1;
	}
	if(stra[abegin]==strb[bbegin])
		return cal(stra,abegin+1,aend,strb,bbegin+1,bend);
	int t1 = cal(stra,abegin,aend,strb,bbegin+1,bend);
	int t2 = cal(stra,abegin+1,aend,strb,bbegin,bend);
	int t3 = cal(stra,abegin+1,aend,strb,bbegin+1,bend);
	return min(t1,t2,t3)+1;
}
int main(){
	string stra = "hello";
	string strb = "selww";
	int lena = stra.size();
	int lenb = strb.size();
	cout<<cal(stra,0,lena-1,strb,0,lenb-1)<<endl;
	return 0;
}

动态规划DP,分析可以到推荐链接中查看

#include <iostream>
#include <string>
using namespace std;
int min(int a,int b,int c){
	int t = a>b?b:a;
	t = t>c?c:t;
	return t;
}
int cal(){
	string stra = "abdd";
	string strb = "aebdd";
	int lena = stra.length();
	int lenb = strb.length();
	
	int c[10][10];
	for(int i=0;i<lena;i++)c[i][lenb]=lena-i;
	for(int j=0;j<lenb;j++)c[lena][j]=lenb-j;
	c[lena][lenb] = 0;
	for(i=lena-1;i>=0;i--)
		for(j=lenb-1;j>=0;j--){
			cout<<i<<" "<<j<<endl;
			if(stra[i]==strb[j])
				c[i][j]=c[i+1][j+1];
			else{
				c[i][j]=min(c[i][j+1],c[i+1][j],c[i+1][j+1])+1;
			}
		}
	return c[0][0];

}
int main(){
	cout<<cal()<<endl;
	return 0;
}