[BZOJ 1046] [HAOI2007] 上升序列 【DP】

题目链接：BZOJ - 1046

题目分析

先倒着做最长下降子序列，求出 f[i]，即以 i 为起点向后的最长上升子序列长度。

注意题目要求的是 xi 的字典序最小，不是数值！

如果输入的 l 大于最长上升子序列长度，输出 Impossible。

否则，从 1 向 n 枚举，贪心，如果 f[i] >= l，就选取 a[i]，同时 --l，然后继续向后找比 a[i] 大的第一个数判断是否 f[i] >= l （这时l已经减小了1）。

代码

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <algorithm>

using namespace std;

const int MaxN = 10000 + 5;

int n, m, l, MaxL;
int T[MaxN];

struct ES 
{
	int Num, Pos, w, v;
} E[MaxN];

inline bool CmpNum(ES e1, ES e2) {
	if (e1.Num == e2.Num) return e1.Pos < e2.Pos;
	return e1.Num < e2.Num;
}
inline bool CmpPos(ES e1, ES e2) {
	return e1.Pos < e2.Pos;
}

inline void Add(int x, int Num) {
	for (int i = x; i <= n; i += i & -i)
		T[i] = max(T[i], Num);
}
inline int Get(int x) {
	if (x == 0) return 0;
	int ret = 0;
	for (int i = x; i; i -= i & -i) 
		ret = max(ret, T[i]);
	return ret;
}

int main() 
{
	scanf("%d", &n);
	for (int i = 1; i <= n; ++i) {
		scanf("%d", &E[i].Num);
		E[i].Pos = i;
	}
	sort(E + 1, E + n + 1, CmpNum);
	int v_Index = 0;
	for (int i = 1; i <= n; ++i) {
		if (i == 1 || E[i].Num > E[i - 1].Num) ++v_Index;
		E[i].v = v_Index;
	}
	sort(E + 1, E + n + 1, CmpPos);
	MaxL = 0;
	for (int i = n; i >= 1; --i) {
		E[i].w = Get(n - (E[i].v + 1) + 1) + 1;
		MaxL = max(MaxL, E[i].w);
		Add(n - E[i].v + 1, E[i].w);
	}
	scanf("%d", &m);
	for (int i = 1; i <= m; ++i) {
		scanf("%d", &l);
		if (l > MaxL) {
			printf("Impossible
");
			continue;
		}
		int x = 0;
		for (int j = 1; j <= n; ++j) {
			if (E[j].v <= x) continue;
			if (E[j].w >= l) {
				x = E[j].v;
				--l;
				printf("%d", E[j].Num);
				if (l > 0) printf(" ");
				else break;
			} 
		}
		printf("
");
	}
	return 0;
}