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)的更多相关文章

  1. hdu 5773 The All-purpose Zero 最长上升子序列+树状数组

    题目链接:hdu 5773 The All-purpose Zero 官方题解:0可以转化成任意整数,包括负数,显然求LIS时尽量把0都放进去必定是正确的. 因此我们可以把0拿出来,对剩下的做O(nl ...

  2. HDU 5773:The All-purpose Zero(贪心+LIS)

    http://acm.hdu.edu.cn/showproblem.php?pid=5773 The All-purpose Zero Problem Description   ?? gets an ...

  3. hdu 2865 Polya计数+(矩阵 or 找规律 求C)

    Birthday Toy Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Tota ...

  4. HDU 1394 Minimum Inversion Number(线段树求最小逆序数对)

    HDU 1394 Minimum Inversion Number(线段树求最小逆序数对) ACM 题目地址:HDU 1394 Minimum Inversion Number 题意:  给一个序列由 ...

  5. HDU 5773 The All-purpose Zero (变形LIS)

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5773 0可以改变成任何数,问你严格递增的子序列最长是多少. 猜测0一定在最长上升子序列中用到,比如2 ...

  6. 【动态规划】【二分】【最长上升子序列】HDU 5773 The All-purpose Zero

    题目链接: http://acm.hdu.edu.cn/showproblem.php?pid=5773 题目大意: T组数据,n个数(n<=100000),求最长上升子序列长度(0可以替代任何 ...

  7. HDU 5773 The All-purpose Zero(树状数组)

    [题目链接] http://acm.hdu.edu.cn/showproblem.php?pid=5773 [题目大意] 给出一个非负整数序列,其中的0可以替换成任意整数,问替换后的最长严格上升序列长 ...

  8. hdu 5773 The All-purpose Zero 线段树 dp

    The All-purpose Zero 题目连接: http://acm.hdu.edu.cn/showproblem.php?pid=5773 Description ?? gets an seq ...

  9. HDU 5773 最长上升子序列

    题意 给出一个序列 问它的最长严格上升子序列多长 这个序列中的0可以被替代为任何数 n的范围给出了1e5 所以平常的O(n*n)lis不能用了 在kuangbin的模板里有O(nlogn)的模板 套上 ...

随机推荐

  1. node.js + mongodb 做项目的详解(二)

    这次内容是结合bootstrap把登陆注册做好,还有就是express的中间件等问题. 看这篇博客之前建议先看我上篇写的那篇博客http://www.cnblogs.com/hubwiz/p/4118 ...

  2. Web项目管理工具精选(下)

    原文:Web项目管理工具精选(下) 我们在上篇中已推介『代码管理.任务管理.支付工具.数据记录.Dashboard Analytics.客户支持』六个方面的工具.本文将介绍剩下七类工具. A/B测试 ...

  3. PAE 分页模式详解

    2016-11-18 记得之前看windows内核原理与实现的时候,在内存管理部分,看到涉及到PAE模式的部分,提到此模式下可以让系统在虚拟地址还是32位宽的情况下,支持64GB的物理内存或者更多.当 ...

  4. Websocket - Websocket原理(握手、解密、加密)、基于Python实现简单示例

    一.Websocket原理(握手.解密.加密) WebSocket协议是基于TCP的一种新的协议.WebSocket最初在HTML5规范中被引用为TCP连接,作为基于TCP的套接字API的占位符.它实 ...

  5. shell_01

    定义变量: name='qwer' 不解析任何字符 name="qwer" 会解析$和\特殊字符 name1=123;name2=456 定义多个变量 now_date=`date ...

  6. mysql 正则表达式 regexp rlike not rlike

    regexp   rlike    not rlike 今天在上班的时候突然遇到了一个问题,C+组的同事要删除mysql里面的一些特殊数据,想要正则进行匹配.于是问我,没想到还真的把我难住了,问题如下 ...

  7. 图文讲解zabbix安装全过程(5)

    在了解<zabbix硬件.软件需求>之后,在你心里应该有备选的机器.今天开始安装zabbix.zabbix需要LNMP或者LAMP环境.环境的搭建不在本章范围内. LNMP环境配置Linu ...

  8. python全栈开发从入门到放弃之socket并发编程之协程

    一.为什么会有协程 本节的主题是基于单线程来实现并发,即只用一个主线程(很明显可利用的cpu只有一个)情况下实现并发,为此我们需要先回顾下并发的本质:切换+保存状态 cpu正在运行一个任务,会在两种情 ...

  9. C++ error C2064:

    error C2064: term does not evaluate to a function 错误原因,非函数,当成函数使用了

  10. 在Idea中连接数据库并生成实体类(mybatis逆向生成实体类)

    1.连接数据库 (1)按下图 ,  点击view-----选择tool windows----------选择database并点击 (2)弹出Database窗口 点击加号------------选 ...