上升序列

【问题描述】  
  
　　对于一个给定的S={a[1],a[2],a[3],…,a[N]}，若有P={a[x1],a[x2],a[x3],…,a[xm]}，满足（x1 < x2 < … < xm）且（a[x1] < a[x2] < … < a[xm]）。那么就称 P 为 S 的一个上升子序列。如果有多个P满足条件，那么我们想求字典序最小的那个。


　　给出S序列，给出若干询问。对于第i个询问，求出长度为Li的上升序列，如有多个，求出字典序最小的那个（即首先x1最小，如果不唯一，再看x2最小……），如果不存在长度为Li的上升序列，则打印Impossible. 
 
    
 【输入格式】  
  
　　第一行一个N，表示序列一共有N个元素
　　第二行N个数，为a1,a2,…,an
　　第三行一个M，表示询问次数。下面接M行每行一个数Li，表示要询问长度为Li的上升序列。


 
    
 【输出格式】  
   
　　对于每个询问，如果对应的序列存在，则输出，否则打印Impossible.
 
    
 【输入样例】   
   
6
3 4 1 2 3 6
3
6
4
5


 
    
 【输出样例】  
   
Impossible
1 2 3 6
Impossible 


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

原题是n<=100000的,我们题库借鉴(chao xi)过来的时候把数据改成了n<=10000的,实际上O(n^2+n*m)加手写输入也能过。

要求规定长度的上升子序列不难,这个题难在要求所输出的子序列字典序最小。(也就是下标小的要排在前)

观察求d[i]的过程,在求d[i]的过程中,我们不只是单纯的求到了每一个a[i]结尾(或者说是开头)的最长上升子序列的长度,我们也得到了以a[i]结尾(或开头的)每一个最长上升子序列,

所以对于整个问题,我们只需要先求出以a[i]开头的上升子序列的长度d[i],然后从左到右依次扫描一次,如果存在d[i]==Li,那么所求序列的第Li个元素就是a[i],并且Li-1,重复以上操作直到Li=0(子问题)。然后依次输出即可。

技巧:因为a[i]开头的情况不好计算,所以将原序列颠倒过后求最长下降子序列。

注意用nlog2n的算法生成d[i]的时候g数组应该是-inf,此时二分要查找的是g数组中第一个小于等于a[i]的元素的下标。

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cctype>
#include<iostream>
#define maxn 10010
#define oo 1000000010
using namespace std;
/*
	g(x)表示将a翻转后长度为x的下降子序列的最后一个元素的最大值 
*/
int n,m,cnt;
int a[maxn],g[maxn],d[maxn];

void initial()//g[i]初始化为-inf
{
	for(int i=1;i<=n;i++)
	g[i]=-oo;
}

void out(int x)//输出次数多可以手写减少时间
{
	if(x==0)  putchar('0');
	int s[15];
	int m=1;
	while(x!=0)
	{
		s[m]=x%10;
		x=x/10;
		m++;
	}
	for(int i=m-1;i>=1;i--)
	putchar(s[i]+'0');
}

int find(int v)<span style="font-family: Arial, Helvetica, sans-serif;">//查找g[1]...g[cnt]中第一个小于等于v的元素的下标,如果没找到就返回最后一个元素的下标</span>
{
	int A=1,B=cnt,ans=0;
	while(A<=B)
	{
		int mid=(A+B)>>1;
		if(g[mid]>v)ans=mid,A=mid+1;
		else B=mid-1;
	}
	return ans;
}


int main()
{
	//freopen("my.in","r",stdin);
	//freopen("my.out","w",stdout);
	
	scanf("%d",&n);
	
	for(int i=1;i<=n;i++)
	scanf("%d",&a[i]);
	
	initial();
	cnt=0;
	for(int i=n;i;i--)//O(nlog2n算法求d[i])
	{
		int t=find(a[i]);
		d[i]=t+1;
		if(a[i]>g[t+1])g[t+1]=a[i];
		cnt=max(t+1,cnt);
	}
	
	/*int cnt=0;
	memset(d,0,sizeof(d));
	d[n]=1;
	for(int i=n-1;i>=1;i--)
	{
		int t=0;
		for(int j=n;j>i;j--)
		if(a[i]<a[j])  t=max(t,d[j]);
		d[i]=t+1;
		cnt=max(cnt,d[i]);
	}*/ //O(n^2)算法 会超时  
	
	scanf("%d",&m);
	for(int k=1;k<=m;k++)
	{
		int x;
		scanf("%d",&x);
		
		int s=0;
		if(x>cnt)  printf("Impossible\n");//如果输入的Li比整个序列的LIS还要长显然不行
		
		else
		{
			for(int i=1;i<=n;i++)if(d[i]>=x && a[i]>s)
			{
				//printf("%d",a[i]);
				out(a[i]);
				if(x!=1)putchar(' ');//最后一个数后不要打空格
				x--;
				s=a[i];
				if(!x)break;//Li==0
			}
			printf("\n");
		}
	}	
	
	return 0;
}