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