BZOJ3963 WF2011 Machine Works

最新推荐文章于 2018-03-25 17:01:21 发布

Fuxey

最新推荐文章于 2018-03-25 17:01:21 发布

阅读量1.1k

点赞数

分类专栏： CDQ分治文章标签： CDQ分治

本文链接：https://blog.csdn.net/Fuxey/article/details/50779371

版权

CDQ分治专栏收录该内容

3 篇文章 0 订阅

订阅专栏

此题有一些坑点，如果你WAWA不止，可以直接看代码后的注意事项

提示：
1. 做过NOI货币兑换那题么？没做过可以先做做那题， BZOJ1492
2. 其实真正有意义的天数并不多，因为不是每天都可以买入卖出。

代码后详细解释：

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <string>
#include <vector>
#include <deque>
#include <stack>
#include <queue>
#include <set>
#include <map>
#include <algorithm>
#include <cassert>

using namespace std;
typedef long long ll;
const double eps = 1e-10;
const int maxn = 1e5+1e2;

int dcmp(double a) { return fabs(a)<eps?0:a<0?-1:1; } 

struct state
{
    ll d , p , r , g;
    void read() { scanf("%lld%lld%lld%lld" , &d , &p , &r , &g); }
    bool operator <(const state& b) const { return d<b.d; }  
    void out(){ cout<<d<<" "<<p<<" "<<r<<" "<<g<<endl; } 
};

struct points
{
    ll x , y;
    points(ll x=0 , ll y=0):x(x),y(y){}
    bool operator <(const points& b)const { return x<b.x || (x==b.x && y<b.y); }
}p[maxn] , ch[maxn];

typedef points Vector;
Vector operator -(Vector a , Vector b) { return Vector(a.x-b.x , a.y-b.y); }
double Cross(Vector a , Vector b) { return 1.0*a.x*b.y-1.0*a.y*b.x; }  


ll n , D , m , t;
ll f[maxn];
state d[maxn];

ll offer(int i , points p) { return p.y+p.x*d[i].d; }

void solve(int l , int r)
{
    if(l==r)
    {
        f[l] = max(f[l-1] , f[l]);
        return ;
    }

    int mid = (l+r)/2;
    solve(l, mid);

    m = t = 0;
    for(int i=l;i<=mid;i++) if(f[i] >= d[i].p) p[++m] = points(d[i].g , f[i]-d[i].p+d[i].r-(d[i].d+1)*d[i].g);

    sort(p+1, p+1+m);
    for(int i=1;i<=m;i++)
    {
        while(t>1 && dcmp(Cross(ch[t]-ch[t-1], p[i]-ch[t-1]))>=0) t--;
        ch[++t] = p[i];
    }

    for(int i=mid+1 , j = 1;i<=r;i++)
    {
        while(j<t && offer(i, ch[j+1])>=offer(i, ch[j])) j++;
        if(j<=t) f[i] = max(f[i] , offer(i, ch[j]));
    }
    solve(mid+1, r);
}

int main(int argc, char *argv[]) {

    int Case=0;
    while(cin>>n>>f[0]>>D && n+f[0]+D)
    {
        for(int i=1;i<=n;i++) d[i].read();
        for(int i=1;i<=n+1;i++) f[i] = 0 ;
        d[++n].d = D+1;
        sort(d+1, d+1+n);

        solve(1, n);
        printf("Case %d: %lld\n" , ++Case , f[n]);
    }

    return 0;
}