P2900 [USACO08MAR]土地征用Land Acquisition

Description

农夫John准备扩大他的农场,他正在考虑N (1 <= N <= 50,000) 块长方形的土地. 每块土地的长宽满足(1 <= 宽 <= 1,000,000; 1 <= 长 <= 1,000,000). 每块土地的价格是它的面积,但FJ可以同时购买多快土地. 这些土地的价格是它们最大的长乘以它们最大的宽, 但是土地的长宽不能交换. 如果FJ买一块3×5的地和一块5×3的地,则他需要付5×5=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组买这些土地: 第一组:100×1, 第二组1×100, 第三组20×5 和 15×15 plot. 每组的价格分别为100,100,300, 总共500.

代码：

#include<iostream>#include<cstdio>#include<cstring>using namespace std;const int maxn=1000000+5;int read(){int x=0,f=1;char ch;ch=getchar();while(ch<'0'||ch>'9'){if(ch=='-') f=-1; ch=getchar();}while(ch>='0'&&ch<='9'){x=x*10+ch-'0'; ch=getchar();}return x*f;}int n,a,b,c;int w[maxn],q[maxn];long long s[maxn],f[maxn];long long g(int i){return a*s[i]*s[i]+b*s[i];}long long h(int j){return f[j]+a*s[j]*s[j]-b*s[j];}double slope(int j1,int j2){return (h(j1)-h(j2))/(2.0*a*(s[j1]-s[j2]));}long long js(int x){return a*x*x+b*x+c;}int main(){int i;n=read();a=read();b=read();c=read();for(i=1;i<=n;i++)w[i]=read();for(i=1;i<=n;i++)s[i]=s[i-1]+w[i];int h=1,p=0;q[++p]=0;for(i=1;i<=n;i++){while(p>h&&slope(q[h],q[h+1])<=s[i]) h++;int t=q[h];f[i]=f[t]+js(s[i]-s[t]);while(p>h&&slope(q[p],i)<slope(q[p-1],q[p])) p--;q[++p]=i;}printf("%lld",f[n]);return 0;}