LightOJ 1030 Discovering Gold（期望）

Description

You are in a cave, a long cave! The cave can be represented by a 1 x N grid. Each cell of the cave can contain any amount of gold.

Initially you are in position 1. Now each turn you throw a perfect 6 sided dice. If you get X in the dice after throwing, you add X to your position and collect all the gold from the new position. If your new position is outside the cave, then you keep throwing again until you get a suitable result. When you reach the N^th position you stop your journey. Now you are given the information about the cave, you have to find out the expected number of gold you can collect using the given procedure.

Input

Input starts with an integer T (≤ 100), denoting the number of test cases.

Each case contains a blank line and an integer N (1 ≤ N ≤ 100) denoting the dimension of the cave. The next line contains N space separated integers. The i^th integer of this line denotes the amount of gold you will get if you come to the i^th cell. You may safely assume that all the given integers will be non-negative and no integer will be greater than 1000.

Output

For each case, print the case number and the expected number of gold you will collect. Errors less than 10^-6will be ignored.

Sample Input
3

1
101

2
10 3

3
3 6 9
Sample Output
Case 1: 101.0000000000
Case 2: 13.000
Case 3: 15

题目链接：http://lightoj.com/volume_showproblem.php?problem=1030

*********************************************

题意：给出T 组实例，每组实例给出n，接着n个数，分别代表i位置有黄金的数量。你掷骰子前进，掷到几就走几步，如果即将到达的位置存在的话。走到了哪个位置就可以得到那个位置的黄金。让你求你到达最后一个位置时收集黄金的期望。

分析：dp保存没点的概率，注意骰子只有6步，所以注意大于6的情况；

数学期望 = 每一点的概率 * 到该点的值。

ＡＣ代码：

 1 #include<stdio.h>
 2 #include<string.h>
 3 #include<math.h>
 4 #include<queue>
 5 #include<algorithm>
 6 #include<time.h>
 7 #include<stack>
 8 using namespace std;
 9 #define N 12000
10 #define INF 0x3f3f3f3f
11
12 double dp[N];
13 int a[N];
14
15 int main()
16 {
17     int n,T,k=1,i,j,x;
18
19     scanf("%d", &T);
20
21     while(T--)
22     {
23         scanf("%d", &n);
24         memset(dp,0,sizeof(dp));
25         double sum=0.0;
26
27         for(i=1;i<=n;i++)
28             scanf("%d", &a[i]);
29
30         dp[1]=1;
31
32         for(i=1;i<=n;i++)
33         {
34            x=(n-i<6?n-i:6);
35             for(j=1;j<=x;j++)
36                 dp[i+j]+=dp[i]*1.0/x;
37
38             sum+=dp[i]*a[i];
39         }
40
41         printf("Case %d: %.10f\n", k++,sum);
42     }
43
44     return 0;
45 }