二分求LIS并打印结果

1275: God‘s ladder [DP]

时间限制: 1 Sec 内存限制: 128 MB  Special Judge

题目描述

天明来到神之宫殿,在他眼前出现了若干个石柱,每个石柱上有1枚金币,天明可以任意选择一个石柱开始,然后向前方的石柱瞬移,而且他所瞬移到的石柱的高度必须要大于现在所在石柱的高度。

求天明所能获得的最大金币数以及任意一种可以获得这么多金币的路线(每个石柱的高度)。

输入

第一行一个数n,表示石柱的个数。

然后2~n+1行,每行一个石柱的高度h[i],分别是1,2,,n石柱的高度。

2<n<200100;1<= h[i] < 15000000;

输出

第一行为一个数m,表示最大金币数。

2~m+1行,每行一个数,分别是每次所瞬移到的石柱的高度。

样例输入

7

3

1

2

5

9

6

7

样例输出

5

1

2

5

6

7

代码:

#include<iostream>

#include<cstdio>

#include<algorithm>

#include<cstring>

using namespace std;

#define  inf  0x3f3f3f3f

int dp[200220],list[200220],a[200220];//dp[i]保存lis为i时最小的元素,list保存每个元素的lis;

int main()

{

int n,i,j,k,p;

while (cin>>n){p=0;

memset(dp,inf,sizeof(dp));

for(i=0;i<n;i++){

scanf("%d",&a[i]);

*lower_bound(dp,dp+n,a[i])=a[i];          //记录更新长度对应的最大潜力元素时,顺便记录下该元素对应的LIS长度

list[p++]=lower_bound(dp,dp+n,a[i])-dp+1;

}

k=lower_bound(dp,dp+n,inf)-dp;     //总的LIS

cout<<k<<endl;

int m=k;

for(i=n-1;i>=0;i--){                          //从右至左扫描一遍直至找全LIS长度

if(k==list[i])

dp[k--]=a[i];

if(!k) break;

}

for(i=1;i<=m;i++)printf("%d\n",dp[i]);

}

return 0;

}

更快的解法,减小二分搜索的次数:

#include<iostream>

#include<cstdio>

#include<algorithm>

#include<cstring>

using namespace std;

#define  inf  0x3f3f3f3f

int dp[200220],list[200220],a[200220];

int main()

{

int n,i,j,k,p;

while (cin>>n){p=0;k=0;

memset(dp,inf,sizeof(dp));

for(i=0;i<n;i++){

scanf("%d",&a[i]);

if(a[i]>dp[k])  dp[++k]=a[i];      // 设置k为当前最长lis值

else*lower_bound(dp,dp+k,a[i])=a[i];

list[p++]=lower_bound(dp,dp+k,a[i])-dp+1;

}

//    k=lower_bound(dp,dp+n,inf)-dp;

k=k+1;

printf("%d\n",k);

int m=k;

for(i=n-1;i>=0;i--){

if(k==list[i])

dp[k--]=a[i];

if(!k) break;

}

for(i=1;i<=m;i++)printf("%d\n",dp[i]);

}

return 0;

}

从右至左扫描的必要性:

由于是求单调上升子序列,所以最大值定在右侧。

又具有相同lis的不同元素中,潜力最大的也在右侧。

即假设一个元素的lis为n(n>2),则在这个元素之前定有lis为(n-1)且小于此元素的原宿存在。

时间: 2024-10-05 23:08:54

二分求LIS并打印结果的相关文章

HDU ACM 1025 Constructing Roads In JGShining&amp;#39;s Kingdom-&amp;gt;二分求解LIS+O(NlogN)

#include<iostream> using namespace std; //BFS+优先队列(打印路径) #define N 500005 int c[N]; int dp[N]; //dp[i]保存的是长度为i的最长不降子序列的最小尾元素 int BS(int n,int x) //二分查找下标,当x比全部元素小时下标为1,比全部元素大时下标为n+1. { int low,high,mid; low=1,high=n; while(low<=high) { mid=(low+h

UVa 103 - Stacking Boxes (LIS,打印路径)

链接:UVa 103 题意:给n维图形,它们的边长是{d1,d2,d3...dn},  对于两个n维图形,求满足其中一个的所有边长 按照任意顺序都一一对应小于另一个的边长,这样的最长序列的个数,并且打印任意一个最长子串的路径, 例如:a(9,5,7,3),b(6,10,8,2),c(9,7,5,1),a和b不满足,但c和b满足 分析:首先对没组边长从小到大排序,再对各组图形按最小边排序,再求最大子串, 对于打印路径,可以逆序循环,也可递归求解 #include<cstdio> #include

HDU 5773 The All-purpose Zero(O(nlgn)求LIS)

http://acm.hdu.edu.cn/showproblem.php?pid=5773 题意: 求LIS,其中的0可以看做任何数. 思路: 因为0可以看做任何数,所以我们可以先不管0,先求一遍LIS,最后再加上0的个数就可以了.当然,每个数需要减去它前面0的个数. 还有这题如果用dp求LIS是要超时的,从别人那里学习了更快的求LIS的方法. 假设存在一个序列d[1..9] = 2 1 5 3 6 4 8 9 7,可以看出来它的LIS长度为5.n 下面一步一步试着找出它. 我们定义一个序列B

二分求幂,快速求解a的b次幂

一个引子 如何求得a的b次幂呢,那还不简单,一个for循环就可以实现! void main(void) { int a, b; int ans = 1; cin >> a >> b; for (int i = 1; i <= b; i++) { ans *= a; } cout << ans; } 那么如何快速的求得a的b次幂呢?上面的代码还可以优化吗? 当然是ok的!下面就介绍一种方法-二分求幂. 二分求幂 所谓二分求幂,即是将b次幂用二进制表示,当二进制位k位

HDU - 1588 Gauss Fibonacci (矩阵快速幂+二分求等比数列和)

Description Without expecting, Angel replied quickly.She says: "I'v heard that you'r a very clever boy. So if you wanna me be your GF, you should solve the problem called GF~. " How good an opportunity that Gardon can not give up! The "Prob

hdu 3641 数论 二分求符合条件的最小值数学杂题

http://acm.hdu.edu.cn/showproblem.php?pid=3641 学到: 1.二分求符合条件的最小值 /*==================================================== 二分查找符合条件的最小值 ======================================================*/ ll solve() { __int64 low = 0, high = INF, mid ; while(low <=

二分求幂(快速求幂,二进制求幂)

二分求幂, 非递归求法(二进制求法): 比如 2^5就是5个2相乘,按照5的二进制求 3^10就是8个3相乘,再2个3相乘. 处理幂的二进制,具体实现代码如下: long long quickmulti(long long a,long long b) { long long res=1; while(b) { if(b&1) //如果最后一位为1,则res*=a; res*=a; a*=a; //a*=a b>>=1; //b%=2 } return res; }

bzoj1692: [Usaco2007 Dec]队列变换(hash+二分求LCP)

以前一直用SA求LCP,今天学习了一波hash+二分求LCP的姿势,也是nlogn而且常数更小了. hash+二分可以logn比较两个后缀的字典序大小,求出LCP然后比较LCP后一个字符的字典序 #include<iostream> #include<cstdlib> #include<cstring> #include<cstdio> #define ull unsigned long long using namespace std; ull mul[1

线段树求LIS并统计最长子序列个数

以下面的题目为例(题目和代码在最后面),给定一个数列(长度最大为10000),求出最长的先增后减子序列长度及个数.做法是先求出以每一个位置结尾的最长单增子序列长度以及以该位置开头的最长单减子序列长度,然后遍历所有位置找出最大先增后减子序列长度. 以最长单增序列(LIS)为例,由于不仅需要整个序列LIS的长度,还要保存以每个位置为结尾位置的LIS长度.记以a[i]结尾的LIS长度为dp[i],则 dp[i] = max{dp[j] | a[j] < a[i]} + 1 这就是一个RMQ问题,涉及单