BZOJ1046 [HAOI2007]上升序列 【LIS + 字典序最小】

1046: [HAOI2007]上升序列

Time Limit: 10 Sec  Memory Limit: 162 MB
Submit: 5410  Solved: 1877
[Submit][Status][Discuss]

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<=1000

Output

　　对于每个询问，如果对应的序列存在，则输出，否则打印Impossible.

Sample Input

6
3 4 1 2 3 6
3
6
4
5

Sample Output

Impossible
1 2 3 6
Impossible

LIS的nlogn算法又用上了，但还是很不熟练

问题要我们算出字典序最小的方案

我们可以根据f[i]用O(n)的复杂度直接扫一遍，当前f[i]还在所求范围内而且A[i]满足条件就输出，保证了字典序最小

总的O(nlogn + nm)不会爆

#include<iostream>#include<cstdio>#include<cstring>#include<algorithm>#define LL long long int#define REP(i,n) for (int i = 1; i <= (n); i++)#define Redge(u) for (int k = head[u]; k != -1; k = edge[k].next)using namespace std;const int maxn = 10005,maxm = 100005,INF = 1000000000;inline int RD(){int out = 0,flag = 1; char c = getchar();while (c < 48 || c > 57) {if (c == '-') flag = -1; c = getchar();}while (c >= 48 && c <= 57) {out = (out << 1) + (out << 3) + c - '0'; c = getchar();}return out * flag;}int A[maxn],f[maxn],bac[maxn],pos[maxn],pre[maxn],ans[maxn],n,len = 0;int main(){n = RD();REP(i,n) A[i] = RD();for (int i = n; i > 0; i--){int l = 0,r = len,mid;while (l < r){mid = l + r + 1 >> 1;if (bac[mid] && A[bac[mid]] > A[i]) l = mid;else r = mid - 1;}f[i] = l + 1; pre[i] = bac[l];if (!bac[f[i]] || A[i] > A[bac[f[i]]]) bac[f[i]] = i;len = max(len,f[i]);}int m = RD(),v,last,first;while (m--){v = RD();if (v > len) printf("Impossible\n");else {last = 0; first = true;for (int i = 1; i <= n; i++)if (f[i] >= v && A[i] > last){if (first) first = false; else printf(" ");printf("%d",A[i]);last = A[i]; v--;if (!v) break;}printf("\n");}}return 0;}