最长单调上升子序列问题

问题表述：

给出一个由n个数组成的序列A[1..n]，找出它的最长单调上升子序列。即求最大的m和a1,a2……,am，使得a1 < a2 < …… < am且A[a1] <= A[a2] <= …… <= A[am]。

 

基本方法：动态规划

设f(i)为以A[i]结尾的最长单调上升子序列的长度，那么所求m即为f = max{f(i), 1 <= i <= n}。

f(i) = 1 + max{0, f(j) | A[j] <= A[i], 1 <= j < i}                                                                            ---  (1)

 

基本性质：如果A[i] < A[j]，但是f(i) > f(j)，那么在计算f(k), k > i, k > j的过程中再也不必用到A[j]，即在（1）的解不可能是

f(k) = f(j) + 1，A[j]可以删除了。

证明：反证，如果（1）的解为f(k) = f(j) + 1，那么A[j] <= A[k]，又A[i]  < A[j]，那么A[i] <= A[k]，

且f(k') = 1 + f(i) > 1 + f(j) = f(k)，这与f(k) = f(j) + 1为（1）的解矛盾。

 

基本性质可以表叙为：在保留下来的元素中，如果A[i] < A[j]，那么f(i) <= f(j)。显然的如果i < j，那么f(i) < f(j)，

因为f(j) >= 1 + f(i) > f(i)。

 

得到算法

 

static void max_nondeac(const int *A, const int el_num) { int idx_kept[el_num];/*index set kept*/ int path[el_num]; /*predecessor of each element*/ int idx_kept_num = 0; /*the number of index kept*/ int i, larger; idx_kept[idx_kept_num++] = 0; path[0] = -1; for (i = 1; i < el_num; i++){ /*larger: idx of idx_kept*/ if ((larger = search_first_larger(idx_kept, A, A[i], idx_kept_num)) < idx_kept_num){ /*del and insert*/