BZOJ1597: [Usaco2008 Mar]土地购买

Description

农夫John准备扩大他的农场,他正在考虑N (1 <= N <= 50,000) 块长方形的土地. 每块土地的长宽满足(1 <= 宽 <= 1,000,000; 1 <= 长 <= 1,000,000). 每块土地的价格是它的面积,但FJ可以同时购买多快土地. 这些土地的价格是它们最大的长乘以它们最大的宽, 但是土地的长宽不能交换. 如果FJ买一块3x5的地和一块5x3的地,则他需要付5x5=25. FJ希望买下所有的土地,但是他发现分组来买这些土地可以节省经费. 他需要你帮助他找到最小的经费.

Input

* 第1行: 一个数: N

* 第2..N+1行: 第i+1行包含两个数,分别为第i块土地的长和宽

Output

* 第一行: 最小的可行费用.

Sample Input

4
100 1
15 15
20 5
1 100

输入解释:

共有4块土地.

Sample Output

500

HINT

FJ分3组买这些土地: 第一组:100x1, 第二组1x100, 第三组20x5 和 15x15 plot. 每组的价格分别为100,100,300, 总共500.

记长为x, 宽为y，我们以X（或Y）排序，一个比较容易想到的dp方程为(记f[i]为前i个矩形的最小可行费用)，f[i] = min(f[i], f[j] + x[i] * max(y[j + 1], y[j + 2]... y[i]))

这个DP是N^2的，而且由于MAX的存在我们完全无法优化。

把长宽抽象成平面上的点，则一个点左下角的点可以被这个点覆盖，也就是说我们可以删除一个点左下角的所有点，这不影响答案。

然后我们可以去掉DP中的MAX，f[i] = min(f[i], f[j] + x[i] * y[j + 1]) (此时y[j + 1]必然是(y[j + 1], y[j + 2]... y[i])中的最大值)。

这个DP像什么？

我们考虑对于i的2个决策状态j和k，假设f[j] + x[i] * y[j + 1] <= f[k] + x[i] * y[k + 1],我们对这个式子进行化简f[j] - f[k] <= x[i](y[k + 1] - y[j + 1])    
(f[j] - f[k])/(y[k + 1] - y[j + 1]) <= x[i]这东西是不是很像斜率？

那么我们就可以用一个队列维护斜率，保证队头2个点的直线斜率小于等于x[i]，否则弹出队头，然后f[i]的最优的那个j实际就是q[队头]，然后我们考虑加入一个点，因为这题的斜率满足单调递增，若是不满足则删除队尾，最后把i加入队尾。

附上代码：

#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int MAXX = 50005;
int sum[MAXX * 40];
int n, i, j, k, l, q[MAXX], x[MAXX], y[MAXX], len;
long long f[MAXX];
bool b[MAXX];
struct sb{
	int x, y;
};
sb a[MAXX];
inline int get()
{
	char c;
	while ((c = getchar()) < 48 || c > 57);
	int res = c - 48;
	while ((c = getchar()) >= 48 && c <= 57)
	res = res * 10 + c - 48;
	return res;
}
inline bool rule(const sb &a, const sb &b)
{
	return (a.x < b.x || (a.x == b.x && a.y < b.y));
}
inline void insert(int x)
{
	for(int i = x; i <= 1000001; i += (i & -i))
		sum[i] ++;
}
inline int find(int x)
{
	int tot = 0;
	for(int i = x; i >= 1; i -= (i & -i))
		tot += sum[i];
	return tot;
}
inline double xie(int i, int j)
{
	return (double)(f[i] - f[j]) / (double)(y[j + 1] - y[i + 1]);
}
int main()
{
	cin >> n;
	for(i = 1; i <= n; i ++)
		a[i].x = get(), a[i].y = get();
	sort(a + 1, a + 1 + n, rule);
	for(i = n; i >= 1; i --)
	{
		if (find(1000001) - find(a[i].y - 1) != 0) b[i] = 1;
		insert(a[i].y);
	}
	for(i = 1; i <= n; i ++)
		if (!b[i]) x[++len] = a[i].x, y[len] = a[i].y;
	int t = 1, w = 1;
	for(i = 1; i <= len; i ++)
	{
		while (t < w && xie(q[t + 1], q[t]) <= (double)x[i]) t ++;
		f[i] = f[q[t]] + x[i] * (long long)y[q[t] + 1];
		while (t < w && xie(q[w], q[w - 1]) >= xie(i, q[w])) w --;
		q[++w] = i;
	}
	cout << f[len] << endl;
}