LCS 最长公共子串

在研究sphinx的时候，看到sphinx的排序算法是LCS
于是小研究了一下LCS，原理：

LCS问题就是求两个字符串最长公共子串的问题。解法就是用一个矩阵来记录两个字符串中所有位置的两个字符之间的匹配情况，若是匹配则为1，否则为0。然后求出对角线最长的1序列，其对应的位置就是最长匹配子串的位置.
下面是字符串21232523311324和字符串312123223445的匹配矩阵，前者为X方向的，后者为Y方向的。不难找到，红色部分是最长的匹配子串。通过查找位置我们得到最长的匹配子串为：21232
0 0 0 1 0 0 0 1 1 0 0 1 0 0 0
　　 0 1 0 0 0 0 0 0 0 1 1 0 0 0 0
　 　1 0 1 0 1 0 1 0 0 0 0 0 1 0 0
　 　0 1 0 0 0 0 0 0 0 1 1 0 0 0 0
　 　1 0 1 0 1 0 1 0 0 0 0 0 1 0 0
　 　0 0 0 1 0 0 0 1 1 0 0 1 0 0 0
　　 1 0 1 0 1 0 1 0 0 0 0 0 1 0 0
　 　1 0 1 0 1 0 1 0 0 0 0 0 1 0 0
　 　0 0 0 1 0 0 0 1 1 0 0 1 0 0 0
　 　0 0 0 0 0 0 0 0 0 0 0 0 0 1 0
　 　0 0 0 0 0 0 0 0 0 0 0 0 0 1 0
　 　0 0 0 0 0 1 0 0 0 0 0 0 0 0 0
　 　0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
但是在0和1的矩阵中找最长的1对角线序列又要花去一定的时间。通过改进矩阵的生成方式和设置标记变量，可以省去这部分时间。下面是新的矩阵生成方式：
0 0 0 1 0 0 0 1 1 0 0 1 0 0 0
　　 0 1 0 0 0 0 0 0 0 2 1 0 0 0 0
　 　1 0 2 0 1 0 1 0 0 0 0 0 1 0 0
　 　0 2 0 0 0 0 0 0 0 1 1 0 0 0 0
　 　1 0 3 0 1 0 1 0 0 0 0 0 1 0 0
　 　0 0 0 4 0 0 0 2 1 0 0 1 0 0 0
　　 1 0 1 0 5 0 1 0 0 0 0 0 2 0 0
　 　1 0 1 0 1 0 1 0 0 0 0 0 1 0 0
　 　0 0 0 2 0 0 0 2 1 0 0 1 0 0 0
　 　0 0 0 0 0 0 0 0 0 0 0 0 0 1 0
　　 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0
　 　0 0 0 0 0 1 0 0 0 0 0 0 0 0 0
　 　0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
当字符匹配的时候，我们并不是简单的给相应元素赋上1，而是赋上其左上角元素的值加一。我们用两个标记变量来标记矩阵中值最大的元素的位置，在矩阵生成的过程中来判断当前生成的元素的值是不是最大的，据此来改变标记变量的值，那么到矩阵完成的时候，最长匹配子串的位置和长度就已经出来了。


代码如下：

def LCS(x, y):
	print x,y
	if (len(x) == 0 or len(y) == 0):
		return ""
	else:
		a = x[0]
		b = y[0]
		if (a == b):
			return LCS(x[1:], y[1:]) + a
		else:
			return cxMax(LCS(x[1:], y), LCS(x, y[1:]))


def cxMax(a, b):
	if (a >= b):
		return a
	else:
		return b

if  __name__ == '__main__':
	cxStra = "123456789"
	cxStrb = "123z"

	lista = []
	listb = []

	for i in range(0, len(cxStra)):
		lista.append(cxStra[i])
	for  i in range(0, len(cxStrb)):
		listb.append(cxStrb[i])

	d = LCS(lista, listb)

	print 'LCS is:%s' % d[::-1]