复习线性DP,少不了经典的LIS问题。进一步思考如何用二分来优化LIS、二分的实现。
时间复杂度分析:
暴力的写法两重循环,纯纯O(n^2),数据超1e5就被卡了。可以考虑用二分将一个
替换成
,时间复杂度降到O(
)。
思路:
暴力两重循环中,会产生很多的冗余,我们来想一想,如何去除这些冗余。
#include <iostream> #include <vector> #include <algorithm> using namespace std; int main() { int n;cin>>n; vector<int> a(n),dp(n); for(int i=0;i<n;i++)cin>>a[i]; dp[0]=1; for(int i=1;i<n;i++){ int max1=1; for(int j=0;j<i;j++){ if(a[i]>a[j])max1=max(max1,dp[j]+1); } dp[i]=max1; } cout<<*max_element(dp.begin(),dp.end()); return 0; }
1.假设,当前的下标为i,我们将1~i-1的上升子序列按照长度进行分类。
2.用一个数组q,该数组下标表示上升子序列长度,数组存的值是,上升子序列长度为x的选择中,结尾数字最小的结尾数字。例如x=2,有两个选择 1 3 和1 4,那我就存3。所以q[2]=3.
3.在目前得到的最大长度中,二分找到一个比当前数字小的结尾数字,那么,我这个数一定可以拼接在这个结尾数字的后面。假如这个结尾数字对应的长度为x,那么拼接后的长度为x+1,所以我们需要判断q[x+1]是否要更新成当前这个数字。
4.假如二分没找到比当前数字小的结尾数字,那么说明1~i-1的每一个数都大于等于当前的数。直接把这个数和q[1]比较就行。
这样的做法能减少很多冗余。比如说在1~i之间:
长度为1的:1,2,4,5,那我只保留1,因为1能和更多的数拼接。长度为2的:1 3,1 4,那我只保留1 3。以3结尾能从4开始选,以4结尾只能从5开始选。
可选的数少了,长度就短了,我们尽可能让他可选的数多一些(一个贪心的想法)。
每个长度的结尾数字是单调递增的(不再证明),所以可以使用二分。
代码实现:
#include <iostream> #include <vector> #include <algorithm> using namespace std; const int N=1e5+10; int q[N]; int main() { int n;cin>>n; vector<int> a(n); for(int i=0;i<n;i++)cin>>a[i];//数据输入 int len=0;//最终结果 for(int i=0;i<n;i++){//扫描整个数组 int l=0,r=len; while(l<r){//二分结尾数字 int mid=l+r+1>>1; if(a[i]>q[mid]){ l=mid; } else r=mid-1; } len=max(r+1,len); q[r+1]=a[i]; } cout<<len; return 0; }