【期望dp】LightOJ 1030 Discovering Gold

http://acm.hust.edu.cn/vjudge/contest/view.action?cid=84305#problem/B

期望dp,跟之前一篇文章飞行棋的很像，注意i+j>N情况的处理

/*  
LightOj 1030 期望DP   
  
仍然是倒着递推 
题意： 
      类似于飞行棋,从1-n，每次走的步数掷筛子决定，问到达终点走的总步数的期望 
      注意：如果超过N则重新投掷，也就是最终一定会走到第N格
求解：  
      和之前做的一道飞行棋的很像，需要注意这样一个问题：
      举个例子，假设走到第N-2个位置，那个掷筛子后如果超出N则重掷，N-1和第N个位置给分摊1/2的概率
      全期望公式：E(Y) = E(Y|X1)*p(X1) +　E(Y|X2)*p(X2) + ……
      dp[i]表示已经到位置i最终到达终点的期望 
      dp[i] = sigma ((1/6.0) * dp[i+k] ) + 1,k = 1,2,3,4,5,6 
      如果遇到i+j>N的情况将概率分摊到剩下位置就好了,对应的1/6改下，你会发现这取决于剩下好有多少格N-i
*/    
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<vector>
#include<set>
#include<map>
#include<algorithm>
#include<sstream>
#define eps 1e-9
#define pi acos(-1)
#define long long ll
#define M 10
#define N 1010
using namespace std;
const int _max = 1e2 + 10;

int n,v[_max];

double dp[_max];

int main()
{
   #ifndef ONLINE_JUDGE
   freopen("input.txt","r",stdin);
   #endif // ONLINE_JUDGE
   int T;cin>>T;int cnt =1;
   while(T--){
    scanf("%d",&n);
    for(int i = 1; i <= n; ++ i) scanf("%d",v+i);
    memset(dp,0,sizeof(dp));
    dp[n] = v[n];//已经到达n
    for(int i = n - 1; i > 0; -- i){
      int d = min(6,n-i);
      for(int j = 1,k = 1; j <= d; ++ j){
        dp[i]+=dp[i+j]/d;
      }
      dp[i]+=v[i];
    }
    printf("Case %d: %.6f\n",cnt++,dp[1]);
   }
   return 0;
}