最长公共子序LCS代码【Longest Common Subsequenct】

//作者：bzljn, liu。
//根据《算法导论》编写。
//头文件

#include "stdafx.h"
#include <string.h>
#include <stdio.h>
#include <iostream>
using namespace std;

//2012-01-21，liu.
//输入数组长度为s1[]=s1[m]，s2[]=s2[n]，
//但是c[][]=c[m+1][n+1]，b[][]=b[m+1][n+1]，参考《算法导论》，liu。
//注意：需要按顺序“左上角，上边，左边”，才能正确的返回c[m][n]。
//自底向上，构造LCS的c[i][j]和b[i][j]，两次for语句。
int  LCS_Length_liu(int *s1, int *s2, int m, int n, int **c, int **b);
//由上往下，输出LCS的c[i][j]和b[i][j]，递归。
void LCS_Print_liu(int *s1, int i, int j, int **b);
void LCS_test_liu();
// directions of LCS generation  
enum decreaseDir {kInit = 0, kLeft, kUp, kLeftUp};

//实现代码
int LCS_Length_liu(int *s1, int *s2, int m, int n, int **c, int **b)
{
	//assert(length(c)==length(b)==[m+1][n+1]);
	if(m<1 || n<1) { return 0; }
	//初始化。
	int i=0, j=0, p=0, k=0;
	for (i=0; i<m+1; i++)
	{
		for (j=0; j<n+1; j++) { c[i][j]=0; b[i][j]=0; }
	}
	for (i=1; i<m+1; i++)
	{
		for (j=1; j<n+1; j++) 
		{ 
			p=i-1, k=j-1;
			if(s1[p]==s2[k])//左上角
			{ c[i][j]=c[i-1][j-1]+1; b[i][j]=kLeftUp; } //记得+1。
			else if(c[i-1][j]>c[i][j-1])//上边
			{ c[i][j]=c[i-1][j]; b[i][j]=kUp; }
			else //if(c[i-1][j]<=c[i][j-1])//左边
			{ c[i][j]=c[i][j-1]; b[i][j]=kLeft; }
		}//End of: for(i)
	}//End of: for(j)
	//test.
	printf("Print out the c[i][j] and b[i][j].\n"); 
	for (i=0; i<m+1; i++)
	{
		for (j=0; j<n+1; j++) { printf("%d ", c[i][j]); }  
		//for (j=0; j<n+1; j++) { printf("%d ", b[i][j]); }
		printf("\n"); 
	}
	//printf("\n"); 
	//注意：需要按顺序“左上角，上边，左边”，才能正确的返回c[m][n]。
	printf("LCS_Length=c[m][n]=%d.\n", c[m][n]); 
	return c[m][n];
}
void LCS_Print_liu(int *s1, int i, int j, int **b)
{
	//assert(length(c)==length(b)==[m+1][n+1]);
	if(i<1 || j<1) { return; }
	if(b[i][j]==kLeftUp) //左上角
	{ 
		LCS_Print_liu(s1, i-1, j-1, b); //记得递归。
		printf("%d,", s1[i-1]);
	}
	else if(b[i][j]==kUp) //上边
	{ LCS_Print_liu(s1, i-1, j, b); }
	else //if(b[i][j]==kLeft) //左边
	{ LCS_Print_liu(s1, i, j-1, b); }
}
void LCS_test_liu()
{
	int s1[]={5,3,4,2,9};
	int s2[]={5,6,4,2,1};
	int m=sizeof(s1)/sizeof(s1[0]);
	int n=sizeof(s2)/sizeof(s2[0]);
	int **c=new int*[m+1]; //c[][]=c[m+1][n+1].
	int **b=new int*[m+1]; //b[][]=b[m+1][n+1].
	for (int i=0; i<m+1; i++)
	{ c[i]=new int[n+1]; b[i]=new int[n+1]; }
	//assert(length(c)==length(b)==[m+1][n+1]);
	//c[2][2]=5; printf("c[2][2]=%d.\n", c[2][2]);
	printf("The input data is: \n");
	for(int i=0; i<m; i++) { printf("%d, ", s1[i]); }
	printf("\n");
	for(int i=0; i<n; i++) { printf("%d, ", s2[i]); }
	printf("\n");
	printf("Result is: \n");
	LCS_Length_liu(s1, s2, m, n, c, b);
	LCS_Print_liu(s1, m, n, b);
	printf("\n");
}

参考文献：

/*
程序员编程艺术第十一章：最长公共子序
http://blog.csdn.net/v_july_v/article/details/6695482
三、动态规划算法解最长公共子序列LCS问题（2011.12.13重写）
http://blog.csdn.net/v_JULY_v/article/details/6110269
《算法导论》
*/
bzljn
完成日期：2012-01-21
链接：http://blog.csdn.net/bzljn/article/details/7213564