当前位置: 代码迷 >> 综合 >> [bzoj1046][DP][乱搞]上升序列
  详细解决方案

[bzoj1046][DP][乱搞]上升序列

热度:94   发布时间:2023-12-19 05:09:50.0

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

题解

啊一开始只会 n 2 l o g n^2log n2log的…
预处理一个向前向后的LIS
从前往后填数
设往前的LIS是f1[i],往后的LIS是f2[i],当前填了的长度是len
如果len+f2[i]-1>=L,那么这个数显然可填
填上就好了…

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
#define zero 10005
using namespace std;
inline int read()
{int f=1,x=0;char ch=getchar();while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}return x*f;
}
inline void write(int x)
{if(x<0)putchar('-'),x=-x;if(x>9)write(x/10);putchar(x%10+'0');
}
inline void print(int x){write(x);printf(" ");}
struct LSnode{int y,p;}w[11000];int tt;
bool LScmp(LSnode n1,LSnode n2){return n1.y<n2.y;}
int nxt[11000],f1[11000],f2[11000],cal[11000],n,m,L;
int s[41000],fac[41000];
int lowbit(int x){return x&-x;}
void change(int x,int c){while(x<=40000)s[x]=max(s[x],c),x+=lowbit(x);}
int findmax(int x){int ret=0;while(x>=1)ret=max(ret,s[x]),x-=lowbit(x);return ret;}
int a[41000],alen;
int main()
{n=read();for(int i=1;i<=n;i++)w[i].y=read(),w[i].p=i;sort(w+1,w+1+n,LScmp);tt=1;cal[w[1].p]=tt;fac[tt]=w[1].y;for(int i=2;i<=n;i++){if(w[i].y!=w[i-1].y)tt++;cal[w[i].p]=tt;fac[tt]=w[i].y;}for(int i=1;i<=n;i++)f1[i]=findmax(cal[i]-1)+1,change(cal[i],f1[i]);memset(s,0,sizeof(s));for(int i=1;i<=n;i++)cal[i]=-cal[i];for(int i=n;i>=1;i--)f2[i]=findmax(cal[i]+zero-1)+1,change(cal[i]+zero,f2[i]);for(int i=1;i<=n;i++)cal[i]=-cal[i];m=read();while(m--){L=read();bool bk=false;for(int i=1;i<=n;i++)if(f1[i]>=L){bk=true;break;}if(!bk)puts("Impossible");else{alen=1;while(alen<=L){int la=cal[a[alen-1]];for(int i=a[alen-1]+1;i<=n;i++)if(cal[i]>la&&alen+f2[i]-1>=L){a[alen]=i;break;}alen++;}for(int i=1;i<alen;i++)print(fac[cal[a[i]]]);puts("");}}return 0;
}