【BZOJ 1046】 1046: [HAOI2007]上升序列
1046: [HAOI2007]上升序列
Description
对于一个给定的S={a1,a2,a3,…,an},若有P={ax1,ax2,ax3,…,axm},满足(x1 < x2 < … < xm)且( ax1 < ax
2 < … < axm)。那么就称P为S的一个上升序列。如果有多个P满足条件,那么我们想求字典序最小的那个。任务给
出S序列,给出若干询问。对于第i个询问,求出长度为Li的上升序列,如有多个,求出字典序最小的那个(即首先
x1最小,如果不唯一,再看x2最小……),如果不存在长度为Li的上升序列,则打印Impossible.Input
第一行一个N,表示序列一共有N个元素第二行N个数,为a1,a2,…,an 第三行一个M,表示询问次数。下面接M
行每行一个数L,表示要询问长度为L的上升序列。N<=10000,M<=1000Output
对于每个询问,如果对应的序列存在,则输出,否则打印Impossible.
Sample Input
6
3 4 1 2 3 6
3
6
4
5Sample Output
Impossible
1 2 3 6
ImpossibleHINT
Source
【分析】
好坑啊,字典序最小指的是下标字典序。。。
那么,我们首先求出以每个数为开头上升序列长度,就是倒着做最长下降子序列。
然后,遇到第一个f[i]>=x的,就输出,然后更新条件,继续在后面找,时间是每次询问是线性的。
==个人觉得数值字典序最小更难吧。。。我打了一个数值字典序最小的,不知道对不对。。
LIS做法:
d[i]表示长度为i的IS中,最后一个数字最小的是什么(因为最后一个数字最小最有利于后面的)
显然d数组是单调递增的。显然我们每次插一个数进去只会改变一个d值或者让LIS增加一位。
若a[i]>a[d[mx]] ,那么可以吧a[i]放在目前的LIS后面,那么mx+1。
否则,二分查找一个x,使得x满足a[i]>d[x]&&a[i]<=d[x+1],那么可以把a[i]放在长度为x的IS后面,替换长度为x+1的IS。
下降序列大于小于号反一下就好了。
本题code:
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
#define Maxn 10010
#define INF 0xfffffff int a[Maxn],b[Maxn],f[Maxn];
int lt[Maxn]; int mymax(int x,int y) {return x>y?x:y;} int n,mx; int tdiv(int x,int l,int r)
{
while(l<r)
{
int mid=(l+r+)>>;
if(x<a[lt[mid]]) l=mid;
else r=mid-;
}
return l;
} void ffind()
{
mx=;
memset(lt,,sizeof(lt));
a[]=INF;
for(int i=n;i>=;i--)
{
if(mx==||a[i]<a[lt[mx]])
{
mx++;
lt[mx]=i;
f[i]=mx;
}
else
{
int x=tdiv(a[i],,mx);
f[i]=x+;
lt[x+]=i;
}
}
} void output(int x)
{
int now=-INF;
for(int i=;i<=n;i++)
{
if(f[i]>=x&&a[i]>now)
{
if(now==-INF) printf("%d",a[i]);
else printf(" %d",a[i]);
now=a[i];
x--;
}
if(x==) break;
}
} int main()
{
freopen("a.in","r",stdin);
freopen("a.out","w",stdout);
scanf("%d",&n);
for(int i=;i<=n;i++) scanf("%d",&a[i]);
ffind();
int q;
scanf("%d",&q);
while(q--)
{
int x;
scanf("%d",&x);
if(x>mx) printf("Impossible\n");
else
{
output(x);
printf("\n");
}
}
return ;
}
数值最小(不知道对不对的)
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
#define Maxn 10010
#define INF 0xfffffff int a[Maxn],b[Maxn],f[Maxn];
int lt[Maxn],nt[Maxn]; int mymax(int x,int y) {return x>y?x:y;} int n,mx; int tdiv(int x,int l,int r)
{
while(l<r)
{
int mid=(l+r+)>>;
if(x>a[lt[mid]]) l=mid;
else r=mid-;
}
return l;
} void ffind()
{
mx=;
memset(nt,,sizeof(nt));
memset(lt,,sizeof(lt));
a[]=-INF;
for(int i=;i<=n;i++)
{
if(mx==||a[i]>a[lt[mx]])
{
mx++;
lt[mx]=i;
nt[i]=lt[mx-];
}
else
{
int x=tdiv(a[i],,mx);
lt[x+]=i;
nt[i]=lt[x];
}
}
} void output(int x)
{
if(x==) return;
output(nt[x]);
printf("%d ",a[x]);
} int main()
{
scanf("%d",&n);
for(int i=;i<=n;i++) scanf("%d",&a[i]);
ffind();
int q;
scanf("%d",&q);
while(q--)
{
int x;
scanf("%d",&x);
if(x>mx) printf("Impossible\n");
else
{
output(lt[x]);
printf("\n");
}
}
return ;
}
2016-12-15 18:29:39
【BZOJ 1046】 1046: [HAOI2007]上升序列的更多相关文章
- BZOJ 1046: [HAOI2007]上升序列 LIS -dp
1046: [HAOI2007]上升序列 Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 3438 Solved: 1171[Submit][Stat ...
- 【BZOJ】1046 : [HAOI2007]上升序列
1046: [HAOI2007]上升序列 题意:给定S={a1,a2,a3,…,an}问是否存在P={ax1,ax2,ax3,…,axm},满足(x1 < x2 < … < xm)且 ...
- Bzoj 1046: [HAOI2007]上升序列 二分,递推
1046: [HAOI2007]上升序列 Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 3671 Solved: 1255[Submit][Stat ...
- BZOJ 1046: [HAOI2007]上升序列(LIS)
题目挺坑的..但是不难.先反向做一次最长下降子序列.然后得到了d(i),以i为起点的最长上升子序列,接下来贪心,得到字典序最小. ----------------------------------- ...
- bzoj 1046 : [HAOI2007]上升序列 dp
题目链接 1046: [HAOI2007]上升序列 Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 3620 Solved: 1236[Submit] ...
- BZOJ 1046: [HAOI2007]上升序列【贪心+二分状态+dp+递归】
1046: [HAOI2007]上升序列 Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 4987 Solved: 1732[Submit][Stat ...
- 1046: [HAOI2007]上升序列(dp)
1046: [HAOI2007]上升序列 Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 4999 Solved: 1738[Submit][Stat ...
- BZOJ1046 [HAOI2007]上升序列 【LIS + 字典序最小】
1046: [HAOI2007]上升序列 Time Limit: 10 Sec Memory Limit: 162 MB Submit: 5410 Solved: 1877 [Submit][St ...
- 【BZOJ1046】[HAOI2007]上升序列
[BZOJ1046][HAOI2007]上升序列 题面 bzoj 洛谷 题解 \(dp\)完之后随便搞一下即可,注意不要看错题 代码 #include <iostream> #includ ...
随机推荐
- (转)MongoDb的十个使用要点
从 mongodb 阶段性技术总结 中抽取并整理了对大家有帮助的十个要点: 1.mongodb 表名和字段名统一用小写字母 mongodb 是默认区分大小写的,为了避免以前在 mysql 下遇到的 ...
- What are Upgrade, Product and Package Codes used for? By pusu
Following content is reprinted from here, please go to the original website for more information. Au ...
- Navicat for mysql linux 破解方法
安装方法 进入下载页面:http://www.navicat.com.cn/download/navicat-for-mysql ,选择Linux版本进行下载,下载后解压安装包,运行 start_ ...
- Unity3D设置字体颜色大小,用于游戏分数显示设置等,
最近在学unity3d,慢慢的学会了许多unity的东西,今天记录下unity3d的Label字体大小及颜色的代码,下面是显示游戏中分数的代码,, public static int Score = ...
- linux 常用软件安装-目录
nginx apache php mysql oracle tomcat memcached mongodb sqlserver
- NSS_01 开始
工作中一直使用asp.net webform, 最近有一个新的小项目, 决定用asp.net mvc3, 边学习边工作吧,记录一下开发过程中的问题,因为工作嘛,只记录问题,可能不会很详细. 准备使用a ...
- NPOI Excel导入 导出
添加引用 using NPOI.HSSF.UserModel; using NPOI.SS.UserModel; using System; using System.Collections.Gene ...
- 在PHP中如何使用消息列队
/** * 消息列队服务 * @author zhou.tingze * @example * -----------------------------------Create----------- ...
- Iso8601 日期格式
unit Iso8601Unit; interface type TIso8601 = class(TObject) public class function DateTimeFromIso8601 ...
- Django1.6添加comments应用的简单过程
今天尝试为自己开发的博客加上评论功能,查看Django的文档,发现1.6之后Django不再自带comments这个app了,下面是官方文档上的说明: Django’s comment framewo ...