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. vue系列

    1.vue基础 2.vue请求和跨域 3.自定义全局组件 4.vue-router 5.Vuex

  2. excel判断单元格包含指定内容的函数用=IF(COUNTIF(A1,"*内容*"),"0","1")

    前面我们聊过怎样将Excel包含某字符的单元格填充颜色,这边我们用另外一种方法来实现:excel判断单元格包含指定内容的函数 选中需要显示结果的单元格,假设我们要判断第一行第一列的单元格A1是否含有“ ...

  3. nodejs实战《一起学 Node.js》 使用 Express + MongoDB 搭建多人博客

    GitHub: https://github.com/nswbmw/N-blog N-blog 使用 Express + MongoDB 搭建多人博客 开发环境 Node.js: 6.9.1 Mong ...

  4. jQuery异步加载数据并添加事件示例

    当时项目是通过树形栏进行权限控制的,管理员可以对从数据库去的数据动态生成树形栏进行增删改查操作,可是用$(".XX").click();方法是不行的. 1.之前用的是jq1.4.3 ...

  5. nginx配置ThinkPHP5二级目录访问

    可以通过 http://www.mracale.com/项目名/模块名/方法名 进行访问 第一步 首先,你要确保在不配置二级目录的情况下,可以通过浏览器访问到.例如:http://www.mracal ...

  6. PAT 1151 LCA in a Binary Tree[难][二叉树]

    1151 LCA in a Binary Tree (30 分) The lowest common ancestor (LCA) of two nodes U and V in a tree is ...

  7. 20165324 实验二《Java面向对象程序设计》实验报告

    20165324 实验二<Java面向对象程序设计>实验报告 一.实验报告封面 课程:Java程序设计 班级:1653班 姓名:何春江 学号:20165324 指导教师:娄嘉鹏 实验日期: ...

  8. 推荐系统第5周--- 基于内容的推荐,隐语义模型LFM

    基于内容的推荐

  9. vmware下安装centos7

    下载vmware http://down-www.newasp.net/pcdown/big/wm_pro_14_win.rar 下载centos7 https://www.centos.org/do ...

  10. CodeForces - 343C Read Time (二分+贪心)

    题意:有N个指针头,M个标记,用这N个针头扫描所有的标记,针头之间互不影响,求扫描完M个标记的最短时间 分析:二分搜答案,mid为时间限制,则只要所有的点在mid秒内被扫描到即可. 对于每个指针,若其 ...