问题:
https://www.acwing.com/activity/content/problem/content/1004/
思路:
朴素做法中的每一个状态的最大值都是由前面所有状态的最大值得到的。时间复杂度是状态的数量和状态的计算量的乘积,为 n2 。在本题中会超时,所以需要对状态数的计算进行优化。
使用一个数组储存 f[i] 状态时已经存在的最长子序列。比如:1 3 4 6 5 2 ,a[4] 的数值为 5,已经存在的最长子序列为 1 3 4 6 . 通过查找该子序列中第一个大于等于 5 的值可以得到为 6 将该值替换为 5 可以得到新的最长子序列为 1 3 4 5 ,该序列的长度以 f[4] 结尾的最长子序列。
查找方式选择二分查找,因为逼近的数据是第一个大于等于的数,所以使用 f[mid] < a[i] 条件来更新 l = mid 。
优化完成后状态数量为 n 每个状态对应的计算量为二分查找的数量为 log(n) ,则总共的时间复杂度为 nlog(n) 。
代码:
#include<iostream>
using namespace std ;
const int N = 1e5 + 10 ;
int n ;
// a[N]为储存的数据,q[N]为维护的最长子序列数组
int a[N], q[N] ;
int main()
{
cin >> n ;
for(int i = 1 ; i <= n ; i ++)
scanf("%d", &a[i]) ;
// 只需要知道a[n]状态对应的最长子序列即可,使用len储存即可。
int len = 0 ;
for(int i = 1 ; i <= n ; i ++)
{
// 二分法模板
int l = 0, r = len ;
while(l < r)
{
int mid = l + r + 1 >> 1 ;
if(q[mid] < a[i]) l = mid ;
else r = mid - 1 ;
}
// 此时的len就是a[i]数据对应的最长子序列。
len = max(len, r + 1) ;
q[r + 1] = a[i] ;
}
cout << len << endl ;
return 0 ;
}