题意:花费n的时间制作一个烟花,然后每个烟花点然后完美的概率就是p*0.0001,然后每当你完成了一个烟花的制作后,你就可以选择是继续花费n的时间再制作一个烟花还是花费m的时间把之前所有制作的烟花点燃。
然后只要在某次的燃放中,有至少一个烟花完美了,那么就可以结束制作的活动了。
问你最小的期望结束制作的时间是多少?
思路:这道题比赛的时候做的时候,因为看到了期望最小,一开始把队友引到了期望DP的误区,确实是这方面解除的题目不够全面,导致拖延了很多思考和调试正解的情况,锅很大qwq。
根据题目叙述我们可以知道,设我们当前制作了x个烟花,那么这x个烟花至少有一个爆炸的概率就是1 -
(
1
−
p
)
x
(1-p)^x
(1−p)x,然后我们就可以把制作x个烟花看做一个事件,这样每次制作x个烟花假如没爆炸的话,我们就继续制作x个烟花,这样这个事件就符合了独立重复的实验的定义了,它的期望也就是我们的时间去乘上1/P,这是伯努利实验的期望公式,其中P就代表1 -
(
1
−
p
)
x
(1-p)^x
(1−p)x,然后很明显这个事件又是一个关于x的函数,即f(x) = (n*x +m) / (1 -
(
1
−
p
)
x
(1-p)^x
(1−p)x),然后这个函数可以求导证明有零点也可以直接惨单峰然后三分。
注意:如果是用浮点数三分的点,那么最后要取一下整,因为根据实际点燃肯定是整数个烟花,只需要比较距它最近的两个点哪个点函数值大就可以了。
代码:
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const double eps = 1e-8;
double p,q,n,m;
double f(double x){
return (n * x + m) / (1.0 - pow(q,x));
}
int main(){
int T;scanf("%d",&T);
while(T--){
scanf("%lf%lf%lf",&n,&m,&p);
p *= 0.0001;
q = 1.0 - p;
double l = 1,r = 1e9;
while(r - l >= eps) {
double mid = l + (r - l) / 3.0;
double midmid = r - (r - l) / 3.0;
if(f(mid) > f(midmid))
l = mid;
else
r = midmid;
}
int t = l;//如果 是浮点数要注意 烟花肯定是整数个 所以 一定是他两侧的整数里小的那个
//直接整数三分也是可以的 而且那样快速幂还快一点 所以就直接整数比较好
printf("%.10f\n",min(f(t),f(t+1)));
}
return 0;
}