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希望买下所有的土地,但是他发现分组来买这些土地可以节省经费. 他需要你帮助他找到最小的经费

Solution

---

容易发现如果一块土地a的长和宽都小于等于另一块土地b的长和宽，那么a是可以去掉不考虑的
我们对土地处理完后按照长升序排序，相同的长再按照宽降序排序，这样选取的土地一定是一段段连续的区间，且最大长和宽分别在区间的开头和结尾，然后列出dp柿子，然后斜率优化

还是那个问题，注意正负号和不等号的方向

Code

---

#include <stdio.h>
#include <string.h>
#include <algorithm>
#define rep(i,st,ed) for (int i=st;i<=ed;++i)
#define drp(i,st,ed) for (int i=st;i>=ed;--i)
#define fi first
#define se second

typedef long long LL;
typedef std:: pair<LL,LL> pair;
typedef long double db;

const int INF=0x7fffffff;
const int N=200005;

int queue[N];
LL f[N];

pair a[N];

db slope(int k,int j) {
    db ret=f[j]-f[k];
    ret=ret/(a[k+1].se-a[j+1].se);
    return -ret;
}

void solve(int n) {
    int head=1,tail=1;
    rep(i,1,n) {
        while (head<tail&&slope(queue[head],queue[head+1])<a[i].fi) head++;
        f[i]=f[queue[head]]-a[queue[head]+1].se*a[i].fi;
        while (head<tail&&slope(queue[tail-1],queue[tail])>slope(queue[tail],i)) tail--;
        queue[++tail]=i;
    }
    printf("%lld\n", f[n]);
}

int main(void) {
    freopen("bzoj1597.in","r",stdin);
    int n; scanf("%d",&n);
    rep(i,1,n) {
        scanf("%lld%lld",&a[i].fi,&a[i].se);
        a[i].se*=-1;
    }
    std:: sort(a+1,a+n+1);
    LL mx=0; int cnt=0;
    drp(i,n,1) if (-a[i].se<=mx) {
        a[i].fi=INF;
        cnt++;
    } else mx=-a[i].se;
    std:: sort(a+1,a+n+1);
    solve(n-cnt);
    return 0;
}