动态规划 最长上升子序列 nlogn

来源:互联网 发布:linux执行文件 编辑:程序博客网 时间:2024/05/19 00:12

题目描述
LIS问题是最经典的动态规划基础问题之一。如果要求一个满足一定条件的最长上升子序列,你还能解决吗?
给出一个长度为N整数序列,请求出它的包含第K个元素的最长上升子序列。
例如:对于长度为6的序列<2,7,3,4,8,5>,它的最长上升子序列为<2,3,4,5>,但如果限制一定要包含第2个元素,那么满足此要求的最长上升子序列就只能是<2,7,8>了。
输入
第一行为两个整数N,K,如上所述。
接下来是N个整数,描述一个序列。
80%的数据,满足0 < n <=1000,0 < k<=n
100%的数据,满足0 < n <=200000,0 < k<=n
输出
请输出两个整数,即包含第K个元素的最长上升子序列长度。

样例输入
8 6
65 158 170 299 300 155 207 389

样例输出
4

首先把k前比k大,k后比k小的元素都去掉(他men没有什么用)
f[0]存答案,b[0]代表搞后数列变成了几个数
于是我们呢先把b[1] 插入f[1]中
若b[i]>f[f[0]] 直接插入不解释
如果f[1]>b[i] f[1]=b[i]
如果其他呢,我们在1到f[0] 这段区间里二分个最先大于b[i] 的位置,把这个位置替换成b[i]
然后就没了
哈哈哈
上代码

#include <cstdio>#include <iostream>#include <cstring>using namespace std;int n,k;int a[200005];int b[200005];int f[200005];inline int dng(int l,int r,int val){    while(l<=r){        int mid=(l+r)>>1;        if(f[mid]==val) return mid;        if(f[mid]>val) r=mid-1;        else l=mid+1;    }    return l;}int main(){    scanf("%d%d",&n,&k);    for(int i=1;i<=n;i++){        scanf("%d",&a[i]);        f[i]=1;    }    int maxn=0;    for(int i=1;i<k;i++){        if(a[i]<=a[k])         b[++b[0]]=a[i];    }    b[++b[0]]=a[k];    for(int i=k+1;i<=n;i++){        if(a[i]>=a[k])        b[++b[0]]=a[i];    }    f[0]=1;    f[1]=b[1];    for(int i=2;i<=b[0];i++){        if(b[i]>f[f[0]]) f[++f[0]]=b[i];        else if(f[1]>b[i]) f[1]=b[i];        else f[dng(1,f[0],b[i])]=b[i];    }    printf("%d",f[0]);    return 0;}