cogs 1487 麻球繁衍

【题目描述】

万有引力定律：

“使物体相互靠近的力的大小与物体的质量成正比——而物体的质量又由同一种力决定。这是一个有趣并且有益的例子，说明了科学是如何用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，那么你的答案就被认为是正确的。

 1 #include <iostream>
 2 #include <algorithm>
 3 #include <cstdio>
 4 #include <cstring>
 5 #include <cstdlib>
 6 #include <cmath>
 7 
 8 using namespace std;
 9 
10 int N, n, k, m;
11 double p[1100], f[1100];
12 
13 int main() {
14     // freopen("in", "r", stdin);
15     freopen("tribbles.in", "r", stdin);
16     freopen("tribbles.out", "w", stdout);
17     scanf("%d", &N);
18     for(int _ = 1 ; _ <= N ; _ ++) {
19         memset(f, 0, sizeof(f));
20         scanf("%d%d%d", &n, &k, &m);
21         for(int i = 0 ; i < n ; i ++) {
22             scanf("%lf", &p[i]);
23         }
24         f[m] = p[0];
25         for(int i = m - 1 ; i >= 1 ; i --) {
26             f[i] = p[0];
27             for(int j = 1 ; j < n ; j ++) {
28                 f[i] += p[j] * pow(f[i + 1], j);
29             }
30         }
31         double ans = 1;
32         for(int i = 1 ; i <= k ; i ++)
33             ans *= f[1];
34         printf("Case #%d: %lf\n", _, ans);
35     }
36 }

View Code

 

转载于:https://www.cnblogs.com/KingSann/articles/7344279.html