题目描述
这是一个简单的动规板子题。
给出一个由 n(n≤5000) 个不超过 1e6 的正整数组成的序列。请输出这个序列的最长上升子序列的长度。
最长上升子序列是指,从原序列中按顺序取出一些数字排在一起,这些数字是逐渐增大的。
输入格式
第一行,一个整数 n,表示序列长度。
第二行有 n 个整数,表示这个序列。
输出格式
一个整数表示答案。
思路:
5000个数n^2肯定不现实
一眼动规
将长度为i(1<=i<=n)的子序列最后一位数的最小值存为f[i]
由于是上升子序列
序列中后面的数必然大于前面的
故f[i]必然大于f[i-1]
即f数组大体如下图
那么每输入一个数a
只需找到f[i]<a<=f[i+1]
随后将f[i+1]更新为a即可
(f[i]意为目前长度为i的子序列内尾数的最小值,若a>f[i],由于a是后输入的(每输入一个数,查找更新一遍),则可将a放置原子序列末端形成长度为i+1的子序列,此时若a<=f[i+1],即可直接更新f[i+1])
如果这里选择遍历一遍f数组,则时间复杂度为O(n^2),不可取
注意到f数组是个单调递增的有序数组
可对其进行二分找到合适位置
最后输出f数组存到的最大长度即可(意味着每次更新后都要判断并记录f数组的最大长度是否被更新,如被更新,加一即可)
理论可行,实践开始
输入:
预处理:
二分查找:
这里要注意mid要写为(l+r+1)/2而不可是(l+r)/2
否则当r=l+1时mid==l
下文又将l赋值为mid
便会陷入死递归
其次,这里的r记录的是f[i]<a<=f[i+1]时的f[i]而非f[i+1]
所以二分出来更新的是f[r+1]而非f[r]
输出:
最后贴上完整代码
#include<bits/stdc++.h>
using namespace std;
int n;
int a[5005];//存输入的n个数
int f[5005];//长度为i的子序列最后一位数最小为f[i]
int main()
{
cin>>n;
for(int i=1;i<=n;i++)
{
cin>>a[i];
}
f[1]=a[1];//只输入一个数时
int len=1;//f数组的最大长度
for(int i=2;i<=n;i++)
{
int l=0,r=len,mid;
while(l<r)
{
mid=(l+r+1)/2;
if(f[mid]<a[i]) l=mid;
else r=mid-1;
}
len=max(len,r+1);//更新f数组最大长度
f[r+1]=a[i];
}
cout<<len;
return 0;
}