码迷,mamicode.com
首页 > 其他好文 > 详细

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

时间:2018-02-04 19:34:33      阅读:193      评论:0      收藏:0      [点我收藏+]

标签:double   struct   pos   str   ==   code   tor   就是   https   

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

看到这种最大最小的就先排个序嘛= =以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]);
}

但是这篇好短啊OvO 算了短就短吧~
真实原因: 还有不到15min放学 作为蒟蒻应该A不掉另一道题再把这篇blog搞成二合一于是就这么短放在这儿了. 其实更主要的原因是懒这种事情我怎么会告诉你萌呢~

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

标签:double   struct   pos   str   ==   code   tor   就是   https   

原文地址:https://www.cnblogs.com/enzymii/p/8413707.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!