BZOJ1597 [USACO2008]土地购买

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.

正解：斜率优化DP

解题报告：这到题目暴力应该好想，就是f[i]=min(f[j]+max(x[k])*max(y[k]))，直接n^2转移

但是n<=50000，所以我们要想怎么优化它，我们看见有个x[k]*y[k]，是不是很像斜率优化的样子

但是因为里面有两个取MAX，导致它没有单调性，x[k]我们可以在外面sort让它单调递增

当x递增时，如果到了点k，前面有一个点j的y也比它小，那么j的费用可以直接并在k这一组里面来

所以把无用的点删完之后(蒟蒻的我用的是并查集，大佬们直接删点)y也是一个单调递减的，然后直接上斜率优化

#include <iostream>
#include <iomanip>
#include <cstdlib>
#include <cstdio>
#include <cmath>
#include <string>
#include <cstring>
#include <algorithm>
#define RG register
#define int long long
const int N = 100000;
const int inf = 2147483641;

using namespace std;

int gi(){
	char ch=getchar();int x=0;
	while(ch<‘0‘ || ch>‘9‘) ch=getchar();
	while(ch>=‘0‘ && ch<=‘9‘) x=x*10+ch-‘0‘,ch=getchar();
	return x;
}

struct date{
	int a,b;
	bool operator < (const date c) const{
		if (a==c.a) return b>c.b;
		return a<c.a;
	}
}f[N];

int st[N],fa[N],dp[N],l=1,r=1,sa[N],t;

int find(int a){
	return a==fa[a]?a:fa[a]=find(fa[a]);
}

long double cal(int k,int y){
	if (f[find(k+1)].b==f[find(y+1)].b) return inf;
	return (long double)(dp[y]-dp[k])/(long double)(f[find(k+1)].b-f[find(y+1)].b);
}

main(){
	int n=gi();
	for (RG int i=1; i<=n; ++i) f[i]=(date){gi(),gi()},fa[i]=i;
	st[1]=fa[0]=0,sort(f+1,f+n+1);
	for (RG int i=1; i<=n; ++i){
		while(t && f[sa[t]].b<=f[i].b) fa[sa[t]]=i,--t;
		sa[++t]=i;
		while(l<r && cal(st[l],st[l+1])<=f[i].a) ++l;
		dp[i]=dp[st[l]]+f[i].a*f[find(st[l]+1)].b;
		while(l<r && cal(st[r-1],st[r])>=cal(st[r],i)) --r;
		st[++r]=i;
	}
	printf("%lld",dp[n]);
	return 0;
}