Dynamic Programming 学习笔记(二) LIS vs LCS

 

LIS (Longest Increasing Subsequence) 最长上升子序列。

题目：给出一个序列a1,a2,a3,a4,a5,a6,a7….an,求它的一个子序列（设为s1,s2,…sn），使得这个子序列满足这样的性质，s1<s2<s3<…<sn并且这个子序列的长度最长。输出这个最长的长度。

LCS (Longest Common Subsequence) 最长公共子序列。

题目：一个数列 S，如果其分别是两个或多个已知数列的子序列，且是最长的序列，则S 称为已知序列的最长公共子序列。

LIS与LCS都是典型的DP问题。

 

对于LIS，有两种基本的算法。

1.  依次遍历整个序列，每一次求出从第一个数到当前这个数的最长上升子序列，直至遍历到最后一个数字为止，然后再取dp数组里最大的那个即为整个序列的最长上升子序列。我们用dp[i]来存放序列1-i的最长上升子序列的长度，那么dp[i]=max(dp[j])+1,(j∈[1, i-1]); 显然dp[1]=1,我们从i=2开始遍历后面的元素即可。O(N^2)复杂度。

int LIS(int arr[1000], int n)
{
	for(int i=1; i<=n; ++i)
		dp[i] = 0;
	int max;
	dp[1] = 1;
	for(int i = 2; i <= n; ++i)
	{
		for(int j=1; j < i; ++j)
		{
			if(arr[i] > arr[j] && dp[j] + 1 > dp[i])
				dp[i] = dp[j] + 1;
		}
		if (dp[i] > max)
			max = dp[i];
	}

	return max;
} 

 

2.  用f[k]表示长度为k的上升子序列最后一个数最小是多少。易知数组f是递增的。读到一个新的数x后，找到某个i使得x>f[i]且x<=f[i+1]，于是用x去更新f[i+1]；特别地，如果所有的f[i]都小于x，则增加f的长度。f的长度即为所求。由于f是递增的，因此可以采用二分查找，时间复杂度为O(nlgn)。

 

int bSearch(int num, int k)  
{  
    int low=1, high=k;  
    while(low<=high)  
    {  
        int mid=(low+high)/2;  
        if(num>=b[mid])  
            low=mid+1;  
        else   
            high=mid-1;  
    }  
    return low;  
}  
 
int LIS_bak()
{
	int low = 1, high = n;
	int k = 1;
	b[1] = p[1];
	for(int i=2; i<=n; ++i)
	{
		if(p[i]>=b[k])
			b[++k] = p[i];
		else
		{
			int pos = bSearch(p[i], k);
			b[pos] = p[i];
		}
	}
	return k;
} 

而对于LCS。

子问题的递归结构为

分析：设序列X=<x₁, x₂, …, x_m>和Y=<y₁, y₂, …, y_n>的一个最长公共子序列Z=<z₁, z₂, …, z_k>，则：

1.  若x_m=y_n，则z_k=x_m=y_n且Z_k-1是X_m-1和Y_n-1的最长公共子序列；

2.  若x_m≠y_n且z_k≠x_{m ，}则Z是X_m-1和Y的最长公共子序列；

3.  若x_m≠y_n且z_k≠y_n ，则Z是X和Y_n-1的最长公共子序列。

其中X_m-1=<x₁, x₂, …, x_m-1>，Y_n-1=<y₁, y₂, …, y_n-1>，Z_k-1=<z₁, z₂, …, z_k-1>。

用c[i,j]记录序列X_i和Y_j的最长公共子序列的长度。其中X_i=<x₁, x₂, …, x_i>，Y_j=<y₁, y₂, …, y_j>。当i=0或j=0时，空序列是X_i和Y_j的最长公共子序列，故c[i,j]=0。

下图能够说明这个问题。

 

注意：最长公共子串（Longest CommonSubstring）和最长公共子序列（LongestCommon Subsequence, LCS）的区别：子串（Substring）是串的一个连续的部分，子序列（Subsequence）则是从不改变序列的顺序，而从序列中去掉任意的元素而获得的新序列

  

参考文献：

http://www.wutianqi.com/?p=1850

http://www.matrix67.com/blog/archives/112

http://blog.csdn.net/v_july_v/article/details/6695482