【洛谷2215】【HNOI2007】上升序列

题目描述

对于一个给定的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.

输入输出格式

输入格式：

第一行一个N，表示序列一共有N个元素第二行N个数，为a1,a2,…,an第三行一个M，表示询问次数。下面接M行每行一个数L，表示要询问长度为L的上升序列。

输出格式：

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

输入输出样例

输入样例#1：

6
3 4 1 2 3 6
3
6
4
5

输出样例#1：

Impossible
1 2 3 6
Impossible

说明

数据范围  N<=10000  M<=1000

题解

由于数据范围是10000，所以n^2的求LIS不行了，可以用树状数组优化，但这题有m次询问，所以和以前的树状数组求LIS不同，这次f[i]表示以i为开头的LIS长度（由于是字典序最小（标号）），发现树状数组不能求了，因为这次求的是比它大的max,用了一个小技巧，把树状数组倒过来，i这个数的位置变成了n-i+1,这样大的数就到了前面，而小的数到了后面，这样就可以求前缀max了（非常巧妙）。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=10010;
int n,m,k,a[maxn],b[maxn],c[maxn],ans[maxn],num;
int t[maxn],f[maxn];
int lowbit(int x){
	return x&(-x);
}
void add(int x,int y){
	x=n-x+1;
	for(int i=x;i<=n;i+=lowbit(i))
	t[i]=max(t[i],y);
}
int find(int x){
	x=n-x+1;//数倒过来储存，前面大后面小，前缀和就是比他大的数中的max
	int ans=0;
	for(int i=x;i>=1;i-=lowbit(i))
	ans=max(t[i],ans);
	return ans;
}
int main(){
	scanf("%d",&n);
	for(int i=1;i<=n;i++){
		scanf("%d",&a[i]);
		c[i]=b[i]=a[i];
    }
	sort(b+1,b+n+1);
	num=unique(b+1,b+n+1)-b-1;
	for(int i=1;i<=n;i++){
		c[i]=lower_bound(b+1,b+num+1,c[i])-b;
	}

/* for(int i=1;i<=n;i++)
	cout<<c[i]<<" ";*/
	int len=0,last=0;
	for(int i=n;i>=1;i--){
		f[i]=find(c[i]+1)+1;
		add(c[i],f[i]);
		len=max(len,f[i]);
	}
	scanf("%d",&m);
	for(int i=1;i<=m;i++){
		scanf("%d",&k);
		if(k>len){
			printf("Impossible\n");
		}
		else{
			last=0;
			for(int i=1;i<=n;i++)
			if(f[i]>=k&&a[i]>last){
				if(k==0) break;
				else{
					printf("%d ",a[i]);
					k--;
					last=a[i];
				}
			}
			printf("\n");
		}
	}

	return 0;
}