【bzoj1046】上升序列

题意

对于一个给定的S={a1,a2,a3,…,an},若有P={ax1,ax2,ax3,…,axm},满足(x1 < x2 < … < xm)且( ax1 < ax2 < … < axm)。那么就称P为S的一个上升序列。如果有多个P满足条件,那么我们想求字典序最小的那个。任务给出S序列,给出若干询问。对于第i个询问,求出长度为Li的上升序列,如有多个,求出字典序最小的那个(即首先x1最小,如果不唯一,再看x2最小……),如果不存在长度为Li的上升序列,则打印Impossible.

\(1\leq n\leq 10000,1\leq m\leq 1000\)

分析

观察数据范围和时限:10s,支持\(O(n^2)\)

假如不要求最小字典序,那么就是dp傻题了。

最小字典序的处理方法,就是从后往前,倒序处理\(f\)数组。

\(f[i]\)表示以\(i\)为终点的从\(i\)之后出发的最长上升序列的最大长度,\(pre[i]\)表示以\(i\)为\(f[i]\)从哪个位置推过来的。

处理的时候尽可能选后面的。

然后对于每个询问,貌似...

找到第一个\(x\)满足\(f[x]>=Li\),然后根据\(pre\)输出\(Li\)个就行了?

尝试能不能找到错误。

发现还可能真的有错误。

因为我们的\(f\)是在子序列尽可能长的情况下满足最小字典序的,子序列尽可能长,这可能会存在更短的子序列,而字典序更小。

所以只能从头到尾重新扫一遍就好啦。

遇到一个\(f[x]\geq Li\),就说明取当前这个有解,且当前这个\(x\)最小,\(Li--\)。

不断输出。

这时候pre根本不需要了。

Extended:

如果要求\(a_i\)字典序怎么办?

方案:按照两个关键字\((a_i,i)\)从小到大排序,和上面的思路相同,要求在后面,逐个尝试。

小结

(1)关于字典序

字典序最小/最大的处理方法:逆序处理。

然后输出就按照自己的记录,或者逐个贪心输出。

代码

#include <cstdio>
#include <cctype>
#include <climits>
#include <algorithm>
using namespace std;

#define rep(i,a,b) for (int i=(a);i<=(b);i++)
#define per(i,a,b) for (int i=(a);i>=(b);i--)

const int N=16384;
const int MIN=INT_MIN;

int n;
int a[N];

int f[N];
int cnt;

int m;

int rd(void) {
    int x=0,f=1; char c=getchar();
    for (;!isdigit(c);c=getchar()) if (c=='-') f=-1;
    for (;isdigit(c);c=getchar()) x=x*10+c-'0';
    return x*f;
}

void Solve(int c) {
    int lastF=MIN;
    rep(i,1,n)
        if (lastF<a[i]&&f[i]>=c) {
            printf("%d",a[i]);
            if (c==1) printf("\n"); else printf(" ");
            lastF=a[i],c--;
            if (!c) break;
        }
}

int main(void) {
    #ifndef ONLINE_JUDGE
    freopen("bzoj1046.in","r",stdin);
    freopen("bzoj1046.out","w",stdout);
    #endif

    n=rd();
    rep(i,1,n) a[i]=rd();

    per(i,n,1) {
        f[i]=1;
        rep(j,i+1,n)
            if (a[i]<a[j])
                f[i]=max(f[i],f[j]+1);
    }
    cnt=*max_element(f+1,f+n+1);

    m=rd();
    rep(i,1,m) {
        int x=rd();
        if (x>cnt)
            printf("Impossible\n");
        else Solve(x);
    }

    return 0;
}

【bzoj1046】上升序列的更多相关文章

  1. [bzoj1046]上升序列

    以i为开头的最长上升子序列,那么就是反过来以i为结尾的最长下降子序列,预处理出来后,不断向后找到下一个数即可 1 #include<bits/stdc++.h> 2 using names ...

  2. Noip前的大抱佛脚----赛前任务

    赛前任务 tags:任务清单 前言 现在xzy太弱了,而且他最近越来越弱了,天天被爆踩,天天被爆踩 题单不会在作业部落发布,所以可(yi)能(ding)会不及时更新 省选前的练习莫名其妙地成为了Noi ...

  3. 【bzoj1046】 HAOI2007—上升序列

    http://www.lydsy.com/JudgeOnline/problem.php?id=1046 (题目链接) 题意 给出一个数列,求数列中长度为L的下标字典序最小的上升子序列. Soluti ...

  4. 【BZOJ1046】[HAOI2007]上升序列

    [BZOJ1046][HAOI2007]上升序列 题面 bzoj 洛谷 题解 \(dp\)完之后随便搞一下即可,注意不要看错题 代码 #include <iostream> #includ ...

  5. 【BZOJ1046】上升序列(动态规划,贪心)

    [BZOJ1046]上升序列(动态规划,贪心) 题面 BZOJ 洛谷 题解 我一开始看错题了,一度以为是字典序最小的序列. 最后发现它要求的字典序是位置的字典序最小. 那就很好办了. 设\(f[i]\ ...

  6. 【BZOJ-1046】上升序列 DP + 贪心

    1046: [HAOI2007]上升序列 Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 3723  Solved: 1271[Submit][Stat ...

  7. BZOJ1046 [HAOI2007]上升序列

    Description 对于一个给定的S={a1,a2,a3,…,an},若有P={ax1,ax2,ax3,…,axm},满足(x1 < x2 < … < xm)且( ax1 < ...

  8. 2014.8.15模拟赛【公主的工作】&&bzoj1046[HAOI2007]上升序列

    bzoj题目是这样的 Description 对于一个给定的S={a1,a2,a3,…,an},若有P={ax1,ax2,ax3,…,axm},满足(x1 < x2 < … < xm ...

  9. [BZOJ1046] [HAOI2007] 上升序列 (dp)

    Description 对于一个给定的S={a1,a2,a3,…,an},若有P={ax1,ax2,ax3,…,axm},满足(x1 < x2 < … < xm)且( ax1 < ...

随机推荐

  1. 网络编程(一)——InetAddress

    网络编程(一)--InetAddress InetAddress类在java中代表的是IP地址,它有两个子类分别是Inet4Address和Inet6Address,其中Inet4Address代表的 ...

  2. 真正意义上下一代 Windows Embedded:有关 Windows 10 "Athens" 的事

    早在去年 BUILD 大会上,Joe Belfiore 就在一台巨大的 "脚踩钢琴" 上简短展示过 "Windows on Device"(又名 Windows ...

  3. jquery之clone()方法详解

    clone()函数用于克隆当前匹配元素集合的一个副本,并以jQuery对象的形式返回.你也可以简单地理解为:克隆当前jQuery对象. 你还可以指定是否复制这些匹配元素(甚至它们的子元素)的附加数据( ...

  4. XUnit的使用

    安装: 通过NuGet获取XUnit的类库 实战使用 https://github.com/kerryjiang/SuperSocket.ProtoBase/tree/master/test

  5. 第五章 consul key/value

    1.key/value作用 动态修改配置文件 支持服务协同 建立leader选举 提供服务发现 集成健康检查 2.使用 2.1.查看全部key/value 说明: 使用?recurse参数来指定查看多 ...

  6. source命令

    source命令用法: source  文件名作用:在当前bash环境下读取并执行文件中的命令.注:该命令通常用命令“.”来替代.如:source .bash_rc 与 . .bash_rc 是等效的 ...

  7. win8下安装matlab7.0

    在win8下安装matlab7.0会出现一些兼容性的问题,需要设置系统环境变量,修改方式如下. 1.设置环境变量,方法:在你的安装目录的\MATLAB7\bin\win32有一个叫做atlas_Ath ...

  8. hdu 5927 Auxiliary Set 贪心

    Auxiliary Set Time Limit: 9000/4500 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others) Pr ...

  9. Android内存泄露测试

    Android性能测试过程中的一些常用命令: CPU: adb shell top -n | grep "+PackageName 内存: adb shell dumpsys meminfo ...

  10. [转]-Android Studio 快捷键整理分享-SadieYu

    文章编辑整理:Android Studio 中文组 - SadieYu Alt+回车 导入包,自动修正 Ctrl+N   查找类 Ctrl+Shift+N 查找文件 Ctrl+Alt+L  格式化代码 ...