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[..] = ,可以看出来它的LIS长度为5。n
下面一步一步试着找出它。
我们定义一个序列B,然后令 i = to 逐个考察这个序列。
此外,我们用一个变量Len来记录现在最长算到多少了
首先,把d[]有序地放到B里,令B[] = ,就是说当只有1一个数字2的时候,长度为1的LIS的最小末尾是2。这时Len=
然后,把d[]有序地放到B里,令B[] = ,就是说长度为1的LIS的最小末尾是1,d[]=2已经没用了,很容易理解吧。这时Len=
接着,d[] = ,d[]>B[],所以令B[+]=B[]=d[]=,就是说长度为2的LIS的最小末尾是5,很容易理解吧。这时候B[..] = , ,Len=
再来,d[] = ,它正好加在1,5之间,放在1的位置显然不合适,因为1小于3,长度为1的LIS最小末尾应该是1,这样很容易推知,长度为2的LIS最小末尾是3,于是可以把5淘汰掉,这时候B[..] = , ,Len =
继续,d[] = ,它在3后面,因为B[] = , 而6在3后面,于是很容易可以推知B[] = , 这时B[..] = , , ,还是很容易理解吧? Len = 了噢。
第6个, d[] = ,你看它在3和6之间,于是我们就可以把6替换掉,得到B[] = 。B[..] = , , , Len继续等于3
第7个, d[] = ,它很大,比4大,嗯。于是B[] = 。Len变成4了
第8个, d[] = ,得到B[] = ,嗯。Len继续增大,到5了。
最后一个, d[] = ,它在B[] = 4和B[] = 8之间,所以我们知道,最新的B[] =,B[..] = , , , , ,Len = 。
于是我们知道了LIS的长度为5。
!!!!! 注意。这个1,,,,9不是LIS,它只是存储的对应长度LIS的最小末尾。有了这个末尾,我们就可以一个一个地插入数据。虽然最后一个d[] = 7更新进去对于这组数据没有什么意义,但是如果后面再出现两个数字 和 ,那么就可以把8更新到d[], 9更新到d[],得出LIS的长度为6。
然后应该发现一件事情了:在B中插入数据是有序的,而且是进行替换而不需要挪动——也就是说,我们可以使用二分查找,将每一个数字的插入时间优化到O(logN)~~~~~于是算法的时间复杂度就降低到了O(NlogN)~! 转自:http://blog.csdn.net/shuangde800/article/details/7474903
原理
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<vector>
#include<stack>
#include<queue>
#include<cmath>
#include<map>
using namespace std; const int maxn=+; int a[maxn];
int d[maxn];
int n; int Search(int x,int left,int right)
{
while(left<=right)
{
int mid=left+(right-left+)/;
if(x<d[mid]) right=mid-;
else left=mid+;
}
return left;
} int dp(int n)
{
int len=;
d[]=a[];
for(int i=;i<n;i++)
{
if(a[i]>d[len])
{
len=len+;
d[len]=a[i];
}
else
{
int pos=Search(a[i],,len);
d[pos]=a[i];
}
}
return len;
} int main()
{
//freopen("D:\\input.txt","r",stdin);
int T;
scanf("%d",&T);
for(int kase=;kase<=T;kase++)
{
scanf("%d",&n);
int num=;
int cnt=;
for(int i=;i<n;i++)
{
int x;
scanf("%d",&x);
if(x==) num++;
else a[cnt++]=x-num;
}
if(num==n) {printf("Case #%d: %d\n",kase,num);continue;}
int ans=dp(cnt);
printf("Case #%d: %d\n",kase,ans+num);
}
return ;
}
HDU 5773 The All-purpose Zero(O(nlgn)求LIS)的更多相关文章
- hdu 5773 The All-purpose Zero 最长上升子序列+树状数组
题目链接:hdu 5773 The All-purpose Zero 官方题解:0可以转化成任意整数,包括负数,显然求LIS时尽量把0都放进去必定是正确的. 因此我们可以把0拿出来,对剩下的做O(nl ...
- HDU 5773:The All-purpose Zero(贪心+LIS)
http://acm.hdu.edu.cn/showproblem.php?pid=5773 The All-purpose Zero Problem Description ?? gets an ...
- hdu 2865 Polya计数+(矩阵 or 找规律 求C)
Birthday Toy Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Tota ...
- HDU 1394 Minimum Inversion Number(线段树求最小逆序数对)
HDU 1394 Minimum Inversion Number(线段树求最小逆序数对) ACM 题目地址:HDU 1394 Minimum Inversion Number 题意: 给一个序列由 ...
- HDU 5773 The All-purpose Zero (变形LIS)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5773 0可以改变成任何数,问你严格递增的子序列最长是多少. 猜测0一定在最长上升子序列中用到,比如2 ...
- 【动态规划】【二分】【最长上升子序列】HDU 5773 The All-purpose Zero
题目链接: http://acm.hdu.edu.cn/showproblem.php?pid=5773 题目大意: T组数据,n个数(n<=100000),求最长上升子序列长度(0可以替代任何 ...
- HDU 5773 The All-purpose Zero(树状数组)
[题目链接] http://acm.hdu.edu.cn/showproblem.php?pid=5773 [题目大意] 给出一个非负整数序列,其中的0可以替换成任意整数,问替换后的最长严格上升序列长 ...
- hdu 5773 The All-purpose Zero 线段树 dp
The All-purpose Zero 题目连接: http://acm.hdu.edu.cn/showproblem.php?pid=5773 Description ?? gets an seq ...
- HDU 5773 最长上升子序列
题意 给出一个序列 问它的最长严格上升子序列多长 这个序列中的0可以被替代为任何数 n的范围给出了1e5 所以平常的O(n*n)lis不能用了 在kuangbin的模板里有O(nlogn)的模板 套上 ...
随机推荐
- Xception
Xception(Deep Learning with Depth-wise Separable convolutions)——google Inception-V3 Xception 并不是真正意义 ...
- Go学习笔记一:解析toml配置文件
本文系作者原创,转载请注明出处https://www.cnblogs.com/sonofelice/p/9085291.html . 一些mysql或者日志路径的信息需要放在配置文件中.那么本博文主要 ...
- Linux命令:nl
全称:number lines of files 用途:显示的时候添加行号. 格式:nl [OPTION]... [FILE]... 类型:nl is /usr/bin/nl 说明: 该命令主要就是针 ...
- windows系统的对象管理
windows中的对象和高级编程语言中所说的对象还欧区别,准确来讲,windows中的对象其实指的是一种数据结构并且是一种带着“对象头(object head)” 的数据结构! 所以windows中 ...
- ssm框架集成Quartz定时器
第一步:添加依赖 <dependency> <groupId>org.quartz-scheduler</groupId> <artifactId>qu ...
- linux下dubbo调试 ---telnet命令
linux下启动dubbo服务端, 怎么调试? 方法有二: 1. 自己写简单消费者功能,进行各种情况测试.(这确实是有必要的) 2. 使用telnet直接连接上dubbo,使用命令调用,然后调试.(这 ...
- sap 动态字段的使用
1.sap中的动态字段其实类似c#里面的指针,都是指向的地址. 只是写法不同,sap 中的field-symbols就是类似指针. , b TYPE i. FIELD-SYMBOLS <fs&g ...
- 如何保护自己的GitHub代码不被别人覆盖
我们在自己的github上创建了免费的公开代码,为了防止别人通过git push upstream master 覆盖了自己原有的代码,需要作一下设置:Settings->Branchs,然后在 ...
- PAT 1033 To Fill or Not to Fill[dp]
1033 To Fill or Not to Fill(25 分) With highways available, driving a car from Hangzhou to any other ...
- 使用 Task 简化异步编程
.Net 传统异步编程概述 .NET Framework 提供以下两种执行 I/O 绑定和计算绑定异步操作的标准模式: 异步编程模型 (APM),在该模型中异步操作由一对 Begin/End 方法(如 ...