题目描述
一个数的序列bi,当b1 < b2 < … < bS的时候,我们称这个序列是上升的。对于给定的一个序列(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).
你的任务,就是对于给定的序列,求出最长上升子序列的长度。
输入
输入的第一行是序列的长度N (1 <= N <= 1000)。第二行给出序列中的N个整数,这些整数的取值范围都在0到10000。
输出
最长上升子序列的长度。
样例输入
7
1 7 3 5 9 4 8
样例输出
4
思路
解法一
首先介绍O(n^2)的方法。将问题变成一个个子问题。列如给定的数组 a[7] = {1,6,4,2,3,9,8}。
这里用到的是一个自底向上的寻找最优子结构的的思想。粗俗来说:如果你想要得到七个数里面的最长子序列,你可以先找前6个数里面的最长子序列,同理,你又必须得找前5个数里面的最长子序列,直到子序列为1
状态转移方程
maxLen(k)表示以ak做为“终点”的最长上升子序列的长度那么:
#include <bits/stdc++.h>
const int N = 1e4 + 1;
using namespace std;
int a[N];
int maxLen[N];
int L = 1;
int main() {
int n;
cin >> n;
for (int i = 1;i <= n;i++) {
cin >> a[i];
maxLen[i] = 1;
}
for (int i = 2;i <= n;i++) {
for (int j = 1;j < i;j++)
if (a[i] > a[j]) maxLen[i] = max(maxLen[i],maxLen[j] + 1);
L = max(L,maxLen[i]);
}
cout << L;
return 0;
}
解法二
二分思想,用这个方法可以把时间复杂度缩小为nlog2(n),我们现在利用一个数组来存放最长子序列 d[i](i 为当前长度),这个数组有这样的规定:不管怎么样,d[i]总的保存已遍历过最大的数,如果遇到比他小的,则在的d[i] 前面的数中,找到第一个比他大的值,然后替换(这里必须明白:替换并不会改变d数组的长度,所以在最后得到的结果是最长子序列长度的排好序的最小序列 ),如果比d[i]大直接接在后面即可,具体算法请看代码
#include <bits/stdc++.h>
const int N = 1e4 + 1;
using namespace std;
int a[N];
int d[N];
int main() {
int n;
cin >> n;
for (int i = 1;i <= n;i++)
cin >> a[i];
d[1] = a[1];//将第一个元素存入
int len = 1;
for (int i = 2;i <= n;i++) {
if (a[i] > d[len]) d[++len] = a[i];
else {
int j = lower_bound(d+1,d + len + 1,a[i]) - d;
d[j] = a[i];//替换序列中比他第一个大的元素
}
}
cout << len;
return 0;
}