ZOJ 3329 One Person Game (概率DP)

     那么  dp[0]=B[0]/(1-A[0]);

转载自kuangbin博客，链接kuangbin博客

代码如下：

#include <iostream>
#include <string.h>
#include <math.h>
#include <queue>
#include <algorithm>
#include <stdlib.h>
#include <map>
#include <set>
#include <stdio.h>
using namespace std;
#define LL __int64
#define pi acos(-1.0)
const int mod=100000000;
const int INF=0x3f3f3f3f;
const double eqs=1e-8;
int p[1000];
double A[1000], B[1000];
int main()
{
    int t, i, j, k, x, n, a, b, c, k1, k2, k3;
    double p0;
    scanf("%d",&t);
    while(t--){
        scanf("%d%d%d%d%d%d%d",&n,&k1,&k2,&k3,&a,&b,&c);
        memset(p,0,sizeof(p));
        for(i=1;i<=k1;i++){
            for(j=1;j<=k2;j++){
                for(k=1;k<=k3;k++){
                    if(i==a&&j==b&&k==c) continue ;
                    p[i+j+k]++;
                }
            }
        }
        p0=1.0/(k1*k2*k3);
        memset(A,0,sizeof(A));
        memset(B,0,sizeof(B));
        for(i=n;i>=0;i--){
            for(j=1;j<=k1+k2+k3;j++){
                A[i]+=A[i+j]*p[j]*1.0/(k1*k2*k3);
                B[i]+=B[i+j]*p[j]*1.0/(k1*k2*k3);
            }
            A[i]+=p0;
            B[i]+=1;
        }
        printf("%.15f\n",B[0]/(1-A[0]));
    }
    return 0;
}

ZOJ 3329 One Person Game (概率DP)

原文：http://blog.csdn.net/scf0920/article/details/42772315