【笔记篇】斜率优化dp（五） USACO08MAR土地购(征)买(用)Land Acquisition

最新推荐文章于 2019-09-05 19:55:30 发布

Enzymii

最新推荐文章于 2019-09-05 19:55:30 发布

阅读量389

点赞数

分类专栏：ＤＰ之神虐杀众生文章标签： dp 优化

本文链接：https://blog.csdn.net/Enzymii/article/details/79170319

版权

ＤＰ之神虐杀众生专栏收录该内容

11 篇文章 0 订阅

订阅专栏

好好的题目连个名字都不统一..

看到这种最大最小的就先排个序嘛= =以x为第一关键字, y为第二关键字排序.
然后有一些 $x_i<=x_{i+1},且y_i<=y_{i+1}$ 的土地就完全可以在买 $i+1$ 的时候顺便把 $i$ 买了.
那么现在就剩下了x递增 y递减的一串. 可以证明一次一起买的应该是连续的一段, 因为中间的y一定比左端点的大, x一定比右端点的小, 那么一定可以在买左右端点的同时把中间的买走.

那就是常见套路了, 方程 $f[i]=f[j]+x[i]*y[$ $j+1$ $]$
然后这个出奇地好化… 都不用数项数怕漏项了… $f[j]$ = $-x[i]$ $*y[i+1]+$ $f[i]$ …
但是吧, 我们按照这个式子找单调性的时候发现这个斜率总是与凸包的边符号相反, 就会出问题..比如找不到该转移的时候了… 反正我遇到这种情况一脸懵逼… 但是后来想想我tm可能是傻…其实就是←_←
我们把方程写成 $f[j]$ = $x[i]$ $*(-y[i+1])+$ $f[i]$ 这样不就ok了么= =
然后就变成了一道大水题~~(然而你不还是去除无用土地的时候去错了么(喂你不要说出来啊很丢人的QAQ))~~此人又日常精分了←_←

代码(这次压过行了, 加了个struct和sort反而降到了18行)

#include <cstdio>
#include <algorithm>
const int N=5e4+5;typedef long long LL;
struct lot{int x,y;}l[N];LL f[N]; int q[N],h,t,n,m;
inline bool operator<(const lot&a,const lot&b){return (a.x==b.x)?a.y<b.y:a.x<b.x;}
double slope(int x,int y){return 1.0*(f[x]-f[y])/(l[y+1].y-l[x+1].y);}
inline int gn(int a=0,char c=0){for(;c<'0'||c>'9';c=getchar());
    for(;c>47&&c<58;c=getchar())a=a*10+c-48;return a;}
int main(){m=gn();for(int i=1;i<=m;++i) l[i].x=gn(),l[i].y=gn();
    std::sort(l+1, l+m+1);
    for(int i=1;i<=m;++i){while(n>0&&l[n].y<=l[i].y) --n;l[++n]=l[i];}
    for(int i=1;i<=n;++i){
        while(h<t&&slope(q[h],q[h+1])<=l[i].x) ++h;
        f[i]=f[q[h]]+1LL*l[i].x*l[q[h]+1].y;
        while(h<t&&slope(q[t],q[t-1])>=slope(q[t],i)) --t;
        q[++t]=i;
    } printf("%lld",f[n]);
}