XDOJ1035--淼·孟德尔

Description

　　作为一个实验室，仅在里面编程序显得很单调，为此WM也经常在实验室做一些物理、化学、生物等实验。一次他受到“植物大战僵尸”游戏的启发，杂交出了能生产金币的植物！这种植物刚开始每株一天能生产K个金币，这样经过D天，WM就能得到K *D个金币了。（-_-）但是WM想要得到更多的金币。每天的工作结束时，WM可以花费P * X个金币对所有植物进行杂交改良，使之一天能多生产X个金币；也可以花费P * Y个金币培育出Y株新植物（这Y株新植物的生产能力与其他植物一样）。
　　一开始WM只有N株每天能生产K个金币的植物，他想通过这些植物赚到M个金币，最少需要几天呢？
注：可以认为培育新植物或者是进行杂交能瞬间完成，即不占用轮次。

Input

输入数据第一行是一正整数T（0＜T≤30），表示有T组测试数据。
每组测试数据只有一行，包括四个整数N, K, P, M（0＜N, K, P≤100, 0＜M≤1000000000），意义如上文描述。

Output

对于每组测试数据，输出一个整数D，表示通过这些植物赚到M个金币最少要D天。

Sample Input

4
2 1 2 10
2 1 2 9
5 4 15 100
1 1 100 1000000000

Sample Output

4
3
5
272

Hint

对于第一组数据2 1 2 10：第一天赚到2 * 1个金币，花费这2个金币去提高生产能力1，接下来3天每天都赚到2 * 2个金币，第4天结束后就有12个金币，满足要求。

 

解题思路:

这个问题的难点是如何找到其中的贪心性质:

（1）如果当前有P个金币、A株植物、植物的生产能力为B，要用这P个金币来提高一天的生产能力并使之尽可能大的方法是：如果A>B的话，用这个P个金币去提高植物的生产能力（B++）；反之则用这P个金币去培育一株新的植物（A++）。简单地说，就是A+B恒定的时候，尽量缩小A和B的差距，使A*B更大。

（2）如果从开始到结束共要花P×Z个金币在提高一天的生产能力上，那么在前期一有金币就用来提高一天的生产能力、直到花了P×Z个金币才停止，这种策略是最快达到目标的。也就是说，如果已经确定要提高一天的生产能，就应该趁早、一有金币就提高。

基于以上两个假设，我们可以枚举最优解总共需要花P*Z（for Z=0 to 2×sqrt（M）），用上面两个假设去计算出每种情况的最少天数。

其实上面是从网上找到答案，对于其中为什么是2×sqrt（M）不明白，而且程序中为了提速，使用的sum>1000的这的条件也不明白。

 

#include <iostream>
#include<cmath>
using namespace std;
const int INF =  300;

int main()
{
    int T;
    cin>>T;
    for(int m=0;m<T;++m)
    {
        int N,K,P,M;
        cin>>N>>K>>P>>M;
        int Z = 2*sqrt(M);
        int minDay = INF;
        int harvest = N*K;
        int sum = 0;
        //枚举
        for(int input=Z;input>=0;--input)
        {
            int temp = P*input;
            int n=N;
            int k = K;
            int days = 1;
            harvest = N*K;
            while(harvest<M)
            {
                if(temp==0)
                {
                    int w = n*k;
                    days += (M-harvest)/w;
                    if((M-harvest)%w!=0)
                        ++days;
                    break;
                }

                if(temp>=P)
                {
                    int q1 = harvest/P;
                    int q2 = temp/P;
                    int minPN = q1>q2?q2:q1;
                    int diff = n>k?(n-k):(k-n);
                    int piff = minPN>diff?minPN-diff:0;
                    diff = minPN>diff?diff:minPN;
                    temp -= minPN*P;
                    harvest -= minPN*P;
                    if(n>k)
                    {
                        k += diff+piff/2;
                        n += piff/2;
                    }
                    else
                    {
                        n += diff+piff/2;
                        k += piff/2;
                    }
                    if(piff%2==1)
                    {
                        n>k?++k:++n;
                    }

                }
                ++days;
                if(days>minDay)
                    break;
                harvest += n*k;
            }
            if(days<minDay)
                minDay = days;
            if(days>minDay)
                ++sum;
            if(sum>1000)
                break;
        }
        cout<<minDay<<endl;
    }
    return 0;
}

 

最后欢迎大家访问我的个人网站： 1024s