1487. 麻球繁衍
★ 输入文件:tribbles.in
输出文件:
tribbles.out
评测插件
时间限制:3 s 内存限制:256 MB
【题目描述】
万有引力定律:
“使物体相互靠近的力的大小与物体的质量成正比——而物体的质量又由同一种力决定。这是一个有趣并且有益的例子,说明了科学是如何用A证明B,再用B证明A的。”——安布罗斯·比尔斯(美国讽刺作家——译者注)。
你有一坨K个毛球(<星际迷航>中的种族——译者注)。这种毛球只会存活一天。在死亡之前,一个毛球有P_i的概率生出i个毛球(i=0,1,...,n-1)。m天后所有毛球都死亡的概率是多少?(包含在第m天前全部死亡的情况)
【输入格式】
输入包含多组数据。
输入文件的第1行是一个正整数N,表示数据组数。
每组数据的第1行有3个正整数n(1<=n<=1000),k(0<=k<=1000),m(0<=m<=1000)。
接下来有n行,给出P_0,P_1,...,P_n-1。
【输出格式】
对于第i组数据,输出"Case #i: ",后面是第m天后所有毛球均已死亡的概率。
【样例输入】
4
3 1 1
0.33
0.34
0.33
3 1 2
0.33
0.34
0.33
3 1 2
0.5
0.0
0.5
4 2 2
0.5
0.0
0.0
0.5
【样例输出】
Case #1: 0.3300000
Case #2: 0.4781370
Case #3: 0.6250000
Case #4: 0.3164063
【提示】
如果你的输出与标准答案相差不超过10^-5,那么你的答案就被认为是正确的。
【来源】
刘汝佳,《算法竞赛入门经典训练指南》表2.8
概率论水题,不解释了。。。
<pre name="code" class="cpp"> #include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
double ans;
int N;
int n,k,m;
double a[1001],num[1001];
int main()
{
freopen("tribbles.in","r",stdin);
freopen("tribbles.out","w",stdout);
scanf("%d",&N);
for (int q=1;q<=N;q++)
{
scanf("%d%d%d",&n,&k,&m);
for (int i=0;i<n;i++)
scanf("%lf",&a[i]);
memset(num,0,sizeof(num));
num[m]=a[0];
for (int i=m-1;i>0;i--)
{
num[i]=a[0];
for(int j=1;j<n;j++)
num[i]+=a[j]*pow(num[i+1],j);
}
ans=1;
for(int i=1;i<=k;i++)ans*=num[1];
printf("Case #%d: %.7lf\n",q,ans);
}
}