• [BZOJ2667][cqoi2012][kcoj]模拟工厂


    题目描述 Description
    有一个称为“模拟工厂”的游戏是这样的:在时刻0,工厂的生产力等于1。在每个时刻,你可以提高生产力或者生产商品。如果选择提高生产力,在下一个时刻时工厂的生产力加1;如果选择生产商品,则下一个时刻你所拥有的商品数量增加p,其中p是本时刻工厂的生产力。
    n个订单,可以选择接受或者不接受。第i个订单(ti, gi, mi)要求在时刻ti给买家提供gi个商品,事成之后商品数量减少gi,而收入增加mi元。如果接受订单i,则必须恰好在时刻ti交易,不能早也不能晚。同一时刻可以接受多个订单,但每个订单只能被接受一次。要求最后的总收入最大。
    例如,如果一共有两个订单(5,1,8)和(7,15,3),用如下策略是最优的:时刻0, 1, 2提高生产力(时刻3的生产力为4),然后在时刻3,4生产商品,则在时刻5时将拥有8个商品。此时接受第1个订单(还会剩下7个商品),并且在时刻5,6继续生产商品,则在时刻7时拥有7+4+4=15个商品,正好满足订单2。
    输入描述 Input Description
    输入第一行包含一个整数n,即订单数目。以下n行每行三个整数ti, gi, mi
    输出描述 Output Description
     输出仅一行,为最大总收入。输出保证在32位带符号整数范围内。
    样例输入 Sample Input
     2
    5 1 8
    7 15 3
    样例输出 Sample Output
     11
    数据范围及提示 Data Size & Hint
     

    编号

    1-3

    4-6

    7-10

    n

    <=5

    <=10

    <=15

    ti

    100

    100

    <=100,000

    gi

    10,000

    10,000

    <=109

    mi

    10,000

    10,000

    <=109

     这道题让我特别失落。有时候这就是命运,你必须接受。

    首先n<=15,可以2^n枚举到底选哪个物品,然后再进行判定这堆物品到底可不可以选。

    对于一段时间,肯定是先提高生产力,然后再开始干活,这样的策略一定是没有错的。

    于是我先想:我们能不能先玩命加生产力,一边加一边判断,再加生产力会不会没有时间,这样既能保证能取到这个物品,而且此后我们的生产力还很大。但是继续想了想发现,如果前面浪大了一直不造物品,一直蓄力,直到最后才敢时间造完当前物品,结果发现下一个物品时间已经不够了,这种情况是很可能存在的,所以以上策略失败。

    针对于我刚刚出现的问题,我想,一定要在当前这段时间内保证生产的物品尽可能多(不再追求于生产力),于是我就列了一个式子(v0+x)(t-x)其中v0表示这段时间开始的时候的生产力,x表示停止蓄力的时间,t是这段时间的长度,于是我们只需令这个式子最大即可,这是一个简单二次函数求极值问题,我们程序算完之后这段时间内的最大值,也就是说生产了一坨物品,当前这个物品耗掉一些,之后继续做,感觉这种思路每一段时间内都能得到最大的物品数,看样子是很靠谱的。想到这个做法之后,我还有一个质疑,最大的物品数,生产力并一定是最大的,可能最终的答案也会有所不同。但是我无法证明,其实是因为式子过于繁琐,比较难证。

    然后我就写了一波,然后交了一发。很显然是WA的。心情不爽的我看了正解:

    正解是在我最开始的思路上继续的:既然这么做我们只是考虑眼前的利益,并且会影响到我们以后的收益,那么我们为什么不在这时就把后边所有物体全考虑一遍呢?还是那个式子(v0+x)(ti-x)>=gi是一个不等式,解出来一个区间,对于每一个i都有一个区间,也就是说x在这些区间的交集中都可以满足使得后面的全取上,贪心的原则:生产力尽可能大,于是我们就去这些区间交集的最右段作为我们的生产力,继续做即可。

    至于为什么是对的,我也不知道。

    针对于两种贪心做法,好像每一种都有其中的道理,但是正确的做法只有一个。我抱着自己的贪心策略,懒得(不会)证明,就默认是对的了,如果你运气好,直接想到了方法2,算你命好,如果像我一样在走方案1时遇到困难转而跳转为方案2,就卒了。

    针对于这类事情,我也不知道该怎么办。

    看到贪心题的题解之后,觉得这种贪心策略挺显然的,考试的时候我肯定能想到,是么?考试的时候你可能猜到了另外一种策略,抱有自己的策略一定是正确的观点做题,只能靠命运了。

    管大视野要了数据,等数据来了我手测一下,10个点中,我的策略究竟能跑过几个点,事后再把这道题A了的代码交上来。

    明天还有物理化学考试,再见了。

    我终于回来了,真是调了一年。

    新代码改进:发现在计算区间时没有必要计算左端点,(可能是因为怕最后的maxl>minr),但实际上如果出现了maxl>minr,在解后期不等式一定会出现derta<0的情况,所以我们只需要计算L就可以,而不用算一个pair了。2^n的时候如果数组是从1开始记录的话if(i&(1<<(j-1))),从零的话就是if(i&(1<<(j)))

    #include<iostream>
    #include<algorithm>
    #include<cstdio>
    #include<queue>
    #include<cmath>
    #include<cstring>
    using namespace std;
    typedef long long LL;
    #define mem(a,b) memset(a,b,sizeof(a))
    inline int read()
    {
        int x=0,f=1;char c=getchar();
        while(!isdigit(c)){if(c=='-')f=-1;c=getchar();}
        while(isdigit(c)){x=x*10+c-'0';c=getchar();}
        return x*f;
    }
    const int oo=2147483647;
    const double eps=1e-8;
    struct indent
    {
        int t;LL g,m;
        bool operator < (const indent &s)const {return t<s.t;}
        indent() {}
        indent(int _1,LL _2,LL _3):t(_1),g(_2),m(_3) {}
    }bill[20],a[20];
    int n,all,now,cnt,tmp;
    LL ans,nans;
    LL downfloor(double x){return (LL)floor(x);}
    LL f(LL v0,LL t,LL a)
    {
        LL derta=(v0-t)*(v0-t)-4*(a-v0*t);
        if(derta<0)return -1ll;
        double x2=((double)(t-v0)+sqrt(derta))/2.0;
        return downfloor(x2);;
    }
    bool check(int now)
    {
        LL sum=0,v0=1,t,minr,SUM;
        for(int i=1;i<=cnt;i++)
        {
            minr=a[i].t-a[i-1].t;SUM=0;
            for(int j=i;j<=cnt;j++)
            {
                SUM+=a[j].g;
                t=a[j].t-a[i-1].t;
                 if(SUM>sum)minr=min(minr,f(v0,t,SUM-sum));
                if(minr<0)return 0;
            }
            if(minr<0)return false;
            v0+=minr;
             sum+=((a[i].t-a[i-1].t-minr)*v0-a[i].g);
        }
        return true;
    }
    int main()
    {
        n=read();
        for(int i=1;i<=n;i++)bill[i].t=read(),bill[i].g=read(),bill[i].m=read();
        sort(bill+1,bill+1+n);
        all=1<<n;a[0].t=0;
        for(int i=1;i<all;i++)
        {
            cnt=nans=0;
            for(int j=1;j<=n;j++)if(i&(1<<(j-1)))a[++cnt]=indent(bill[j].t,bill[j].g,bill[j].m),nans+=bill[j].m;
            if(check(i))ans=max(ans,nans);
        }
        printf("%lld
    ",ans);
        return 0;
    }
    /*
    
    5
    12 30 49
    20 17 3
    26 7 12    
    71 1109 1113
    86 1578 3136
    
    3200
    */
    View Code

    最后说一句:垃圾贪心毁我青春。

  • 相关阅读:
    curl: (1) Protocol 'http not supported or disabled in libcurl
    线程-分为两类-用户线程和守护线程
    laypage分页插件的使用
    uploadify上传图片插件的使用
    redis安装
    php连接测试memcached
    pageY、clientY、screenY、offsetY的区别
    audio和video样式兼容
    实现剪切和复制功能
    滚动条样式
  • 原文地址:https://www.cnblogs.com/FYH-SSGSS/p/6528124.html
Copyright © 2020-2023  润新知