LCS

LCS（Longest Common Subsequence），最长公共子序列，典型的dp问题。

问题：设字符串A、B，求A和B的LCS。

二维dp：用dp[i][j]数组表示以A[i]结尾、B[j]结尾的字符串的LCS长度，那么

(1)dp[i][j] = dp[i - 1][j - 1] + 1 if A[i] == B[j]

(2)dp[i][j] = max(dp[i - 1][j], dp[i][j - 1]  if A[i] != B[j]

这种方法的时间复杂度是O(nm)，空间复杂度也是O(nm)，记录路径时，可以记录dp[i][j]是由哪一步转化而来。

#define N 1005
#define DIAGONAL 0
#define LEFT 1
#define UP 2

char str1[N],str2[N];
int dp[N][N],path[N][N];
int dir[3][2] = {{-1,-1},{-1,0},{0,-1}};

void _print(int i,int j)
{
	if(i == 1)
	{
		printf("%c",str1[i]);
		return;
	}
	if(j == 1)
	{
		printf("%c",str2[j]);
		return;
	}
	int ni,nj;
	ni = i + dir[path[i][j]][0];
	nj = j + dir[path[i][j]][1];
	_print(ni,nj);
	if(path[i][j] == DIAGONAL)
		printf(" %c",str1[i]);
}

void print()
{
	_print(strlen(str1) - 1,strlen(str2) - 1);
	puts("");
}

int lcs()
{
	memset(dp,0,sizeof(dp));
	memset(path,0,sizeof(path));
	for(int i = 1;str1[i];i++)
		for(int j = 1;str2[j];j++)
			if(str1[i] == str2[j])
				dp[i][j] = dp[i - 1][j - 1] + 1;
			else
			{
				dp[i][j] = max(dp[i - 1][j],dp[i][j - 1]);
				if(dp[i][j] == dp[i - 1][j])
					path[i][j] = LEFT;
				else
					path[i][j] = UP;
			}
	return dp[strlen(str1) - 1][strlen(str2) - 1];
}