1046: [HAOI2007]上升序列 - BZOJ
Description
对于一个给定的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.
Input
第一行一个N,表示序列一共有N个元素第二行N个数,为a1,a2,…,an 第三行一个M,表示询问次数。下面接M行每行一个数L,表示要询问长度为L的上升序列。
Output
对于每个询问,如果对应的序列存在,则输出,否则打印Impossible.
Sample Input
6
3 4 1 2 3 6
3
6
4
5
Sample Output
Impossible
1 2 3 6
Impossible
数据范围
N<=10000
M<=1000
先nlogn倒着做完最长上升,然后用一个数组b存以i为开头的最长上升子序列的长度
每个询问,从前往后扫一遍,能用就用,就是字典序最小的了(MD,一开始以为是数字最小,还排了一个序,原来是下标)
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
1 const 2 maxn=10010; 3 var 4 f,a,b:array[0..maxn]of longint; 5 n,m,max:longint; 6 7 procedure find(x:longint); 8 var 9 l,r,mid:longint; 10 begin 11 l:=1; 12 r:=max+1; 13 while l<>r do 14 begin 15 mid:=(l+r)>>1; 16 if f[mid]>a[x] then l:=mid+1 17 else r:=mid; 18 end; 19 if l>max then max:=l; 20 if f[l]<a[x] then f[l]:=a[x]; 21 b[x]:=l; 22 end; 23 24 procedure init; 25 var 26 i:longint; 27 begin 28 read(n); 29 for i:=1 to n do 30 begin 31 f[i]:=-100000000; 32 read(a[i]); 33 end; 34 for i:=n downto 1 do 35 find(i); 36 end; 37 38 procedure work; 39 var 40 i,j,x,lasta:longint; 41 begin 42 read(m); 43 for i:=1 to m do 44 begin 45 read(x); 46 lasta:=-1000; 47 if x>max then write('Impossible') 48 else 49 for j:=1 to n do 50 if (b[j]>=x) and (a[j]>lasta) then 51 begin 52 dec(x); 53 lasta:=a[j]; 54 if x>0 then write(a[j],' ') 55 else write(a[j]); 56 if x=0 then break; 57 end; 58 writeln; 59 end; 60 end; 61 62 begin 63 init; 64 work; 65 end.