hdu 5406 CRB and Apple（树状数组+dp）

题目链接：hdu 5406 CRB and Apple

将树按照高度从大到小排序，高度相同的美味值小的放前面。然后高度值就没有用了，等于说对得到的美味值序列求两个不相交的LIS，使得长度和尽量长。对序列元素离散化后，dp[i][j]表示说两个序列的终止位置分别为i,j的最大长度。从左向右考虑每个美味值，复杂度o(n)。对于一个值c，更新dp数组，枚举第一维状态i，复杂度o(n)。dp[i][c] = dp[c][i] = max(dp[i][j] | j < c)，o(n)，但是这里可以用树状数组或是其他数据结构，优化到log(n)。注意不能边枚举变修改，会对答案造成影响。

#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;
#define lowbit(x) ((x)&(-x))
const int maxn = 1005;

struct Item {
	int h, d;
	bool operator < (const Item& u) const { return h > u.h || (h == u.h && d < u.d); }
}P[maxn];

int N, M, idx[maxn], top, fenw[maxn][maxn];

void add (int* f, int x, int w) {
	while (x < maxn) {
		f[x] = max(f[x], w);
		x += lowbit(x);
	}
}

int find(int* f, int x) {
	int ret = 0;
	while (x) {
		ret = max(ret, f[x]);
		x -= lowbit(x);
	}
	return ret;
}

void init () {
	scanf("%d", &N);
	memset(fenw, 0, sizeof(fenw));
	for (int i = 0; i < N; i++) {
		scanf("%d%d", &P[i].h, &P[i].d);
		idx[i] = P[i].d;
	}
	sort(P, P + N);

	top = 1;
	sort(idx, idx + N);
	for (int i = 1; i < N; i++) {
		if (idx[i] == idx[top-1]) continue;
		idx[top++] = idx[i];
	}
	for (int i = 0; i < N; i++)
		P[i].d = lower_bound(idx, idx + top, P[i].d) - idx + 1;
}

int solve () {
	int ans = 0, tmp[maxn];
	for (int i = 0; i < N; i++) {
		int v = P[i].d;
		for (int j = 0; j <= top; j++)
			tmp[j] = find(fenw[j], v + 1) + 1;
		for (int j = 0; j <= top; j++) {
			ans = max(ans, tmp[j]);
			add(fenw[j], v + 1, tmp[j]);
			add(fenw[v], j + 1, tmp[j]);
		}
	}
	return ans;
}

int main () {
	int cas;
	scanf("%d", &cas);
	while (cas--) {
		init ();
		printf("%d\n", solve());
	}
	return 0;
}