最长上升子序列问题

蒜白

已于 2023-08-21 10:26:02 修改

阅读量205

点赞数 3

文章标签：算法 c++ 动态规划

于 2023-08-21 10:11:26 首次发布

本文链接：https://blog.csdn.net/2302_77698668/article/details/132389868

版权

目录

一、问题详情

输入格式

输出格式

数据范围

输入样例：

输出样例：

二、问题分析

1.状态标识

2. 状态转移

三、代码

1.未优化版本

2.优化版本

思路如下:

2.1状态表示

2.2状态转移

小结：

一、问题详情

给定一个长度为 N的数列，求数值严格单调递增的子序列的长度最长是多少。

输入格式

第一行包含整数 N。

第二行包含 N 个整数，表示完整序列。

输出格式

输出一个整数，表示最大长度。

数据范围

1≤N≤1000，
−1e9≤数列中的数≤1e9

输入样例：

7
3 1 2 1 8 5 6

输出样例：

二、问题分析

1.状态标识：

用 f [ i ] 表示从第一个数字开始遍历到第 i 个数字结尾的子字符串最大长度

需要初始化 f [ i ] 为 1 ，因为每个子字符串最小长度为1.，用数组 a [ N ] 来存储每个数字

2. 状态转移：

需要对 i 之前的 f [ 0 ~ i - 1 ] 进行条件筛选，为什么不是 f [ i - 1 ]呢，因为

0~ i - 1 当中每个状态 f [ j ] 只要满足 a [ i ] > a [ j ] ,都可以由 f [ j ] 转移到 f [ i ] ;

三、代码

1.未优化版本

代码如下（示例）：

#include<iostream>
using namespace std;
const int N =1010;
int n,a[N],f[N];

int main(){
    cin>>n;
    for(int i=1;i<=n;i++){
        cin>>a[i];
    }
    
    int ans=1;
    for(int i=1;i<=n;i++){
        f[i]=1;
        for(int j=1;j<i;j++){
            if(a[j]<a[i]) f[i]=max(f[i],f[j]+1);
        } 
        ans=max(ans,f[i]);
    }
    cout<<ans;
}

2.优化版本

未优化版本的时间复杂度是O(n^2) ，但是我们还有更好的方法，可以将时间复杂度简化到

O(n * log n);

思路如下:

2.1状态表示

用 f [ i ] 表示长度为 i 的上升子字符串最后结尾的数字【也就是该字符串最大值】，那么我们最后输出的就是该子字符串的长度计量 i

2.2状态转移

假设有 a [ i + 1 ] > = a [ i ] ,那么 a [ i + 1 ] 一定无法纳入 f [ i + 1 ] ,但是我们找到 a [ i + 1 ] 在 f [ i ] 当中本该存在的位置，令该位置的编号为 j ,该位置是原 f [ i ] 字符串当中第一个大于等于 a [ i + 1 ]的数，如此更新下来，f [ i ] 数组当中每一位储存都是长度为 i 字符串对应的末端数字当中的最小值。（这里涵盖了贪心的思想）。那么我们最后取 i 的最大值就是所求的答案

这里我们用 a [ i + 1 ]的值替换 a [ j ]的值,不会导致 j 长度本身发生改变，相反，根据贪心的思路，我们每个长度子字符串最大值越小越有利于形成下一个更大长度子字符串。终上所述，该猜想成立。

代码如下：

#include<iostream>
#include<vector>
using namespace std;
const int N =1e5+7;
int a[N];
int main(){
    int n;cin>>n;
    for(int i=1;i<=n;i++ ) cin>>a[i];
    vector<int> sta;//我们在这里用 堆来代替f [ N ]
    sta.push_back(a[1]);
    for(int i=2;i<=n;i++){
        if(a[i]>sta.back()) sta.push_back(a[i]);
        else *lower_bound(sta.begin(),sta.end(),a[i])=a[i];
        //这一段是找到 a[i] 在 sta 中比a[i]大的第一个数并且替换成a[i]
    }
    cout<<sta.size()<<endl;
    return 0;
}