★【动态规划】【单调性优化】【NOI2011】NOI嘉年华

【问题描述】
NOI2011 在吉林大学开始啦!为了迎接来自全国各地最优秀的信息学选手,吉林大学决定举办两场盛大的 NOI 嘉年华活动,分在两个不同的地点举办。每个嘉年华可能包含很多个活动,而每个活动只能在一个嘉年华中举办。
现在嘉年华活动的组织者小安一共收到了 n 个活动的举办申请,其中第 i 个活动的起始时间为 Si,活动的持续时间为 Ti。这些活动都可以安排到任意一个嘉年华的会场,也可以不安排。
小安通过广泛的调查发现,如果某个时刻,两个嘉年华会场同时有活动在进行(不包括活动的开始瞬间和结束瞬间),那么有的选手就会纠结于到底去哪个会场,从而变得不开心。所以,为了避免这样不开心的事情发生,小安要求不能有两个活动在两个会场同时进行(同一会场内的活动可以任意进行)。
另外,可以想象,如果某一个嘉年华会场的活动太少,那么这个嘉年华的吸引力就会不足,容易导致场面冷清。所以小安希望通过合理的安排,使得活动相对较少的嘉年华的活动数量最大。
此外,有一些活动非常有意义,小安希望能举办,他希望知道,如果第 i 个活动必须举办(可以安排在两场嘉年华中的任何一个),活动相对较少的嘉年华的活动数量的最大值。
【输入格式】
从文件 show.in 中读入数据。
输入的第一行包含一个整数 n,表示申请的活动个数。
接下来 n 行描述所有活动,其中第 i 行包含两个整数 Si、Ti,表示第 i 个活动从时刻 Si 开始,持续 Ti 的时间。
【输出格式】
输出到文件 show.out 中。
输出的第一行包含一个整数,表示在没有任何限制的情况下,活动较少的嘉年华的活动数的最大值。
接下来 n 行每行一个整数,其中第 i 行的整数表示在必须选择第 i 个活动的前提下,活动较少的嘉年华的活动数的最大值。
【评分标准】
对于一个测试点:
如果输出格式不正确(比如输出不足 n+1 行),得 0 分;
如果输出文件第一行不正确,而且后 n 行至少有一行不正确,得 0 分;
如果输出文件第一行正确,但后 n 行至少有一行不正确,得 4 分;
如果输出文件第一行不正确,但后 n 行均正确,得 6 分;
如果输出文件中的 n+1 行均正确,得 10 分。
【样例输入】
5
8 2
1 5
5 3
3 2
5 3

【样例输出】
2
2
1
2
2
2

【样例说明】
在没有任何限制的情况下,最优安排可以在一个嘉年华安排活动 1, 4,而在另一个嘉年华安排活动 3, 5,活动 2 不安排。
【数据规模与约定】
所有测试数据的范围和特点如下表所示

此题考察动态规划及其单调性优化。
区间离散化，设A = {嘉年华1的活动}, B = {嘉年华2的活动}，C = {未安排的活动}。
设num[i][j]为包含在区间[i, j]中的区间个数；
pre[i][j]为在区间[0, i]中放入j个区间到B后，最多能放入A的个数；
suf[i][j]为在区间[i, +∞)中放入j个区间到B后，最多能放入A的个数。

那么有以下两个转移方程：



当然程序中用了向后递推的方法。

第一问就很简单地算出来了，即：

或

对于第二问，先考虑朴素做法：
设g为满足条件的最大方案。
那么：

所以最终答案为：


现在考虑优化：
对于确定的i, j，记

那么则有：随着x的增大，使得g(i, j)最大的y是逐渐减小的。
于是对于g的计算可以优化到O (n ^ 3)
Accode：

#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <cstring>
#include <string>
using namespace std;
const int maxN = 410, INF = 0x3f3f3f3f;
int pre[maxN][maxN], suf[maxN][maxN], g[maxN][maxN];
int num[maxN][maxN], tab[maxN], L[maxN], R[maxN], n, m;

int main()
{
	freopen("show.in", "r", stdin);
	freopen("show.out", "w", stdout);
	scanf("%d", &n);
	for (int i = 0; i < n; ++i)
	{
		scanf("%d%d", L + i, R + i);
		tab[i << 1] = L[i];
		tab[(i << 1) + 1] = R[i] += L[i];
	}
	sort(tab, tab + (n << 1));
	int cnt = unique(tab, tab + (n << 1)) - tab;
	//unique函数，可以将有序的列表去重。
	//这里的cnt为最终离散化过后的区间右界。
	for (int i = 0; i < n; ++i)
		L[i] = lower_bound(tab, tab + cnt, L[i]) - tab,
		R[i] = lower_bound(tab, tab + cnt, R[i]) - tab;
	//离散化过程。
	for (int i = 0; i < cnt; ++i)
	{
		for (int j = 0; j < n; ++j)
			if (L[j] >= i) ++num[i][R[j]];
		for (int j = i + 1; j < cnt; ++j)
			num[i][j] += num[i][j - 1];
	}
	memset(pre, ~0x3f, sizeof pre); pre[0][0] = 0;
	memset(suf, ~0x3f, sizeof suf); suf[cnt - 1][0] = 0;
	for (int i = 0; i < cnt; ++i)
	{
		for (int j = 0; j < n + 1; ++j)
		if (pre[i][j] > ~INF)
			pre[i][pre[i][j]] = max(pre[i][pre[i][j]], j);
		//这里是将pre设一个初始值。
		for (int j = n - 1; j > -1; --j)
			pre[i][j] = max(pre[i][j], pre[i][j + 1]);
		for (int j = 0; j < n + 1; ++j)
		for (int k = i + 1; k < cnt; ++k)
		if (pre[i][j] > ~INF)
			pre[k][pre[i][j]] = max(pre[k][pre[i][j]], j + num[i][k]);
		//向后递推，相当于在B中放入了pre[i][j]个区间，
		//在A中放入了j + num[i][k]个区间。
	}
	for (int i = cnt - 1; i > -1; --i)
	{
		for (int j = 0; j < n + 1; ++j)
		if (suf[i][j] > ~INF)
			suf[i][suf[i][j]] = max(suf[i][suf[i][j]], j);
		for (int j = n - 1; j > -1; --j)
			suf[i][j] = max(suf[i][j], suf[i][j + 1]);
		for (int j = 0; j < n + 1; ++j)
		for (int k = i - 1; k > -1; --k)
		if (suf[i][j] > ~INF)
			suf[k][suf[i][j]] = max(suf[k][suf[i][j]], j + num[k][i]);
	}
	for (int i = 0; i < cnt; ++i)
	for (int j = i; j < cnt; ++j)
	{
		g[i][j] = ~INF;
		for (int x = 0, y = n; x < n + 1; ++x)
		{
			while (y > -1 && x + y > num[i][j] + pre[i][x] + suf[j][y]) --y;
			//让x + y（即放入B的个数）刚刚小于
			//num[i][j] + pre[i][x] + suf[j][y]（即放入A的个数），
			//这时的g值最大。
			if (y > -1) g[i][j] = max(g[i][j], x + y); //
		}
	}
	int ans = 0;
	for (int i = 0; i < n + 1; ++i)
		ans = max(ans, min(i, suf[0][i])); //
	printf("%d\n", ans);
	for (int i = 0; i < n; ++i)
	{
		int ans = 0;
		for (int j = 0; j < L[i] + 1; ++j)
		for (int k = R[i]; k < cnt; ++k)
			ans = max(ans, g[j][k]);
		printf("%d\n", ans);
	}
	return 0;
}