Longest Increasing Subsequence

描述

Given an unsorted array of integers, find the length of longest increasing subsequence.

For example,

Given [10, 9, 2, 5, 3, 7, 101, 18],

The longest increasing subsequence is [2, 3, 7, 101], therefore the length is 4. Note that there may be more than one LIS combination, it is only necessary for you to return the length.

Your algorithm should run in O(n^2) complexity.

Follow up: Could you improve it to O(n log n) time complexity?

解法1 动规

这是一个多阶段决策问题,求最长,是一个最优化问题,用 BFS, 贪心或DP。

如果用BFS,首先用数组中的所有元素作为根节点,形成n颗树,每棵树开始往下扩展,出现逆序则终止,最后计算每棵树的高度,取最大,就是最终结果。算法复杂度为O(n!)

本题中,一个节点往下扩展的时候,没有一个确定的准则,让你走哪些边,本题不具有贪心选择性质,因此不能用贪心法。

我们来尝试用DP来解决这题。最重要的是要定义出状态。首先从状态扩展这方面看,对于数组中的一个元素,它往后走,凡是比它大的元素,都可以作为下一步,因此这里找不到突破口。

我们换一个角度,从结果来入手,我们要求的最长递增子序列,一个递增子序列,肯定是有首尾两个端点的,假设我们定义 f[i] 为以第i个元素为起点的最长递增子序列,那么 f[i]f[j]之间没有必然联系,这个状态不好用。

假设定义f[i]为以第i个元素为终点的最长递增子序列,那么如果i<jnums[i]<nums[j],则f[i]一定是f[j]的前缀。这个状态能表示子问题之间的关系,可以接着深入下去。

现在正式开始定义,我们定义状态f[i]为第i个元素为终点的最长递增子序列的长度,那么状态转移方程是

f[j] = max{f[i], 0 <= i < j && f[i] < f[j]} + 1

有了状态和状态转移方程,代码就不难写了。

  1. // Longest Increasing Subsequence
  2. // 时间复杂度O(nlogn),空间复杂度O(n)
  3. class Solution {
  4. public:
  5. int lengthOfLIS(vector<int>& nums) {
  6. if (nums.empty()) return 0;
  7. // f[i]表示以i结尾的最长递增子序列的长度
  8. vector<int> f(nums.size(), 1);
  9. int global = 1;
  10. for (int j = 1; j < nums.size(); ++j) {
  11. for (int i = 0; i < j; ++i) {
  12. if (nums[i] < nums[j]) {
  13. f[j] = max(f[j], f[i] + 1);
  14. }
  15. }
  16. global = max(global, f[j]);
  17. }
  18. return global;
  19. }
  20. };

解法2 Insert Position

本题最后有一个进阶问题,能不能O(n log n) 解决?有。

维护一个单调递增序列,遍历数组,二分查找每一个数在单调序列中的位置,然后替换之。

  1. // Longest Increasing Subsequence
  2. // 时间复杂度O(nlogn),空间复杂度O(n)
  3. class Solution {
  4. public:
  5. int lengthOfLIS(vector<int>& nums) {
  6. vector<int> lis;
  7. for (auto x : nums) {
  8. int insertPos = lower_bound(lis, 0, lis.size(), x);
  9. if (insertPos >= lis.size()) {
  10. lis.push_back(x);
  11. } else {
  12. lis[insertPos] = x;
  13. }
  14. }
  15. return lis.size();
  16. }
  17. int lower_bound (const vector<int>& A, int first, int last, int target) {
  18. while (first != last) {
  19. int mid = first + (last - first) / 2;
  20. if (target > A[mid]) first = ++mid;
  21. else last = mid;
  22. }
  23. return first;
  24. }
  25. };

原文: https://soulmachine.gitbooks.io/algorithm-essentials/content/cpp/dp/longest-increasing-subsequence.html