bzoj1597 土地购买

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

DP斜率优化的一道题

方程:f[i] = min(f[j - 1] +　f[j] * x[i]);

优化后(f[q[l] + 1] - f[q[l]]) / (y[q[l] + 1] - y[q[l] + 2])；

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
#define ll long long 
using namespace std;
struct Node
{
	ll x,y;
}a[50010];
int q[50010],n, tot;
ll x[50010],y[50010], f[50010];
int cmp(Node a, Node b)
{
	if(a.x == b.x) return a.y < b.y;
	else return a.x < b.x;
}

int main()
{
	scanf("%d",&n);
	for(int i = 1; i <= n; i ++)
		scanf("%d%d",&a[i].x, &a[i].y);
	
	sort(a + 1, a + n + 1, cmp);
	for(int i = 1; i <= n; i ++)
    {
        while(tot && a[i].y >= y[tot])tot --;
        x[++ tot] = a[i].x;
		y[tot] = a[i].y;
    }
	int l = 0, r = 0; 
	for(int i = 1; i <= n; i ++)
	{
		while(l < r && ((f[q[l] + 1] - f[q[l]]) / (y[q[l] + 1] - y[q[l] + 2])) < x[i]) l ++;
		f[i] = f[ q[l] ] + y[q[l] + 1] * x[i];
		while(l < r && ((f[i] - f[q[r]]) / (y[q[r] + 1] - y[i + 1]))< ((f[ q[r] ] - f[q[r - 1]]) / (y[q[r - 1] + 1]-y[q[r] + 1])) )r--;
		q[++ r] = i;
	}
	printf("%d",f[tot]);
	return 0;
}