zoukankan      html  css  js  c++  java
  • 最长上升子序列(LIS)nlogn模板

    参考https://www.cnblogs.com/yuelian/p/8745807.html

    注意最长上升子序列用lower_bound,最长不下降子序列用upper_bound
    比如123458, 加入了5
    假设求最长上升子序列
    这个时候只能替换5,不能替换8(严格上升)
    虽然没有用,但是这样不会错,写upper_bound就错了。
    假设求最长不下降子序列
    这样应该替换8,替换5并不是最优的
    所以用upper_bound

    最长上升子序列(LIS)nlogn模板

    #include<cstdio>
    #include<algorithm>
    #define REP(i, a, b) for(int i = (a); i < (b); i++)
    using namespace std;
    
    const int MAXN = 51234;
    int a[MAXN], f[MAXN], n; //a数组从0开始,f数组从1开始 
    
    int main()
    {
    	scanf("%d", &n);
    	REP(i, 0, n) scanf("%d", &a[i]);
    	int len = 1;
    	f[1] = a[0]; //初始化 
    	REP(i, 1, n)
    	{
    		if(a[i] > f[len]) f[++len] = a[i]; //这里是++len 若是不下降就改为>= 
    		else f[lower_bound(f + 1, f + len + 1, a[i]) - f] = a[i]; //注意f数组是从1开始 
    	}
    	printf("%d
    ", len);
    	return 0;
    }

    最长不下降子序列(LIS)nlogn模板

    #include<cstdio>
    #include<algorithm>
    #define REP(i, a, b) for(int i = (a); i < (b); i++)
    using namespace std;
    
    const int MAXN = 51234;
    int a[MAXN], f[MAXN], n; 
    
    int main()
    {
    	scanf("%d", &n);
    	REP(i, 0, n) scanf("%d", &a[i]);
    	int len = 1;
    	f[1] = a[0]; 
    	REP(i, 1, n)
    	{
    		if(a[i] >= f[len]) f[++len] = a[i]; //>改成>= 
    		else f[upper_bound(f + 1, f + len + 1, a[i]) - f] = a[i];  //lower_bound改成upper_bound
    	}
    	printf("%d
    ", len);
    	return 0;
    }

    如果要求最长下降子序列或者最长不上升子序列符号改变,同时二分加上cmp即可

    另外有个神奇的定理
    如果是求一个数组最少分成几组最长不上升子序列的话
    答案就是最长上升子序列(上升改成下降也成立)
    导弹拦截那题要用到

    输出路径的版本,见https://blog.csdn.net/lxcxingc/article/details/81238008

    #include<cstdio>
    #include<algorithm>
    #define REP(i, a, b) for(int i = (a); i < (b); i++)
    using namespace std;
    
    const int MAXN = 51234;
    int a[MAXN], f[MAXN];
    int ans[MAXN], pos[MAXN], n;
    
    int main()
    {
    	scanf("%d", &n);
    	REP(i, 0, n) scanf("%d", &a[i]);
    	int len = 1;
    	f[1] = pos[1] = a[0];
    	REP(i, 1, n)
    	{
    		if(a[i] > f[len]) f[++len] = a[i], pos[i] = len;
    		else f[pos[i] = lower_bound(f + 1, f + len + 1, a[i]) - f] = a[i];
    	}
    	printf("%d
    ", len);
    	
    	int maxx = 1e9, t = len;
    	for(int i = n - 1; i >= 0; i--)
    	{
    		if(t == 0) break;
    		if(pos[i] == t && maxx > a[i])
    		{
    			maxx = a[i];
    			ans[t--] = a[i];
    		}
    	}
    	REP(i, 1, len + 1) printf("%d ", ans[i]);
    	puts("");
    	
    	return 0;	
    } 
  • 相关阅读:
    [转]责任链模式
    spring中常见注解描述
    [转]外观模式
    [转]策略模式
    分布式事务实践
    SpringBoot入门
    服务器性能调优
    kvm qemu内幕介绍
    xen 保存快照的实现之 —— device model 状态保存
    xen hypercall 的应用层实现
  • 原文地址:https://www.cnblogs.com/sugewud/p/9819456.html
Copyright © 2011-2022 走看看