【dp】关于LIS最长上升子序列

51nod1134 最长递增子序列 

题目链接

 

最近在刷dp的题， 感觉很奇妙。玄学递推22333
本题是求最长上升子序列， 不难得到转移方程

可以写出这样的代码：

#include <bits/stdc++.h>
using namespace std;

int n, a[50010], dp[50010];//dp[i]代表前i项的最长上升子序列

int LIS()
{
    int tmp_max = -1010;
    for(int i = 1; i<=n; i++)
    {
        dp[i] = 1;
        for(int j = 1; j<i; j++)
        {
            if(a[j] < a[i])
                dp[i] = max(dp[i], dp[j] + 1);
        }
        tmp_max = max(tmp_max, dp[i]);
    }
    return tmp_max;
}

int main()
{
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
    cin >> n;
    for(int i = 1; i<=n; i++)
    {
        cin >> a[i];
    }
    cout<<LIS();
    return 0;
}

但是 这样写的LIS的时间复杂度是n^2的  数据范围是50000  肯定会TLE的  所以我们得优化这种算法，达到nlogn的级别  

有了这样的想法后，就可以写出nlogn的做法了！
 

#include <bits/stdc++.h>
using namespace std;

int n, a[50010], dp[50010];//dp[i]代表前i项的最长上升子序列

int LIS()
{
    int len = 1;
    dp[1] = a[1];
    for(int i = 2; i<=n; i++)
    {
        if(a[i] > dp[len])//当前值大于最大的值时  表示序列还可以上升
        {
            dp[++len] = a[i];
        }
        else
        {
            int pos = lower_bound(dp+1, dp+len+1, a[i]) - dp;
            dp[pos] = a[i];
        }
    }
    return len;
}

int main()
{
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
    cin >> n;
    for(int i = 1; i<=n; i++)
    {
        cin >> a[i];
    }
    cout<<LIS();
    return 0;
}