算法设计与分析 SCAU8596 最长上升子序列（优先做）

8596 最长上升子序列（优先做）

时间限制:300MS 代码长度限制:10KB
提交次数:255 通过次数:118

题型: 编程题 语言: G++;GCC;VC;JAVA

Description
当元素 ai1 < ai2 < … < aiK. 就说这个序列是有序上升的。

给定序列(a1, a2, …, aN)，存在许多这样的子序列(ai1, ai2, …, aiK)，
其中1 <= i1 < i2 < … < iK <= N.
也就是说，子序列是原序列允许挑选若干不连续的元素形成的序列。

举个例子，序列 (1, 7, 3, 5, 9, 4, 8) 就有许多个上上子序列，比如(1, 7), (3, 4, 8) 等。
所有这些上升子序列中最长的长度为4，比如 (1, 3, 5, 8).

你来编程实现，当给定一个初始序列，寻找这个序列的最长上升子序列。

输入格式

此例包含多个测试cases（少于100组test cases）。
每一组test case包含2行。
第一行是这组case的序列长度 N。（N的范围0~10000）
第二行包含 N个整数的一个序列，用空格间隔这N个整数， 1 <= N <= 10000。

当N为0时，表示测试结束。

输出格式

输出必须对每个test case，都有一个整数结果，表示这组case的最长上升子序列的长度。

输入样例
7
1 7 3 5 9 4 8
6
1 8 3 6 5 9
5
1 2 3 4 5
0

输出样例
4
4
5

解题思路

动态规划

• 状态定义：dp[i] 的值代表 nums 以 nums[i] 结尾的最长子序列长度。

• 转移方程： 设 j∈[0,i)，考虑每轮计算新 dp[i] 时，遍历 [0,i) 列表区间，做以下判断：

  1. 当 nums[i] > nums[j] 时： nums[i] 可以接在 nums[j] 之后（此题要求严格递增），此情况下最长上升子序列长度为 dp[j] + 1 ；
  2. 当 nums[i] <= nums[j] 时： nums[i] 无法接在 nums[j] 之后，此情况上升子序列不成立，跳过。
  • 上述所有 1. 情况 下计算出的 dp[j] + 1 的最大值，为直到 i 的最长上升子序列长度（即 dp[i] ）。实现方式为遍历 j 时，每轮执行 dp[i] = max(dp[i], dp[j] + 1)。
  • 转移方程： dp[i] = max(dp[i], dp[j] + 1) for j in [0, i)。

• 初始状态：dp[i] 所有元素置 1，含义是每个元素都至少可以单独成为子序列，此时长度都为 1。

• 返回值：返回 dp 列表最大值，即可得到全局最长上升子序列长度。

复杂度分析：

• 时间复杂度 O(N^2) ： 遍历计算 dp 列表需 O(N)，计算每个 dp[i] 需 O(N)。

• 空间复杂度 O(N) ： dp 列表占用线性大小额外空间。

代码如下

#include <iostream>
#include <stdlib.h>
#include <stdio.h>
#include <math.h>

using namespace std;

int main()
{
    int n, i, j;
    while(true) {
        scanf("%d", &n);
        if(n == 0) break;

        int nums[200005];
        for (int i = 0; i < n; i++) {
            scanf("%d", &nums[i]);
        }

        int dp[200005]; // 动态规划数组

        for(i = 0; i < n; i++)
            dp[i] = 1;

        int res = 1;
        for (i = 0; i < n; i++) {
            for (j = 0; j < i; j++) {
                // 说明是递增，可以接上去
                if(nums[j] < nums[i]) {
                    dp[i] = max(dp[i], dp[j] + 1);
                    if(dp[i] > res)
                        res = dp[i];
                }
            }
        }

        printf("%d\n", res);
    }

    return 0;
}