C - One Person Game ZOJ - 期望DP-归零条件

最新推荐文章于 2020-12-09 20:47:37 发布

ID_BePosit

最新推荐文章于 2020-12-09 20:47:37 发布

阅读量180

点赞数

CC 4.0 BY-SA版权

分类专栏： DP

本文链接：https://blog.youkuaiyun.com/BePosit/article/details/83478897

DP 专栏收录该内容

105 篇文章

订阅专栏

C - One Person Game
ZOJ - 3329
转自一篇写的非常好的博客：
原文：https://blog.youkuaiyun.com/sd_stjean/article/details/75175232
只需按照公式去求解A，B即可，我的代码中dp数组只是预处理了一下个各个分数的概率。

#include<iostream>
#include<stdio.h>
#include<cstring>
using namespace std;
#define maxn 666
double dp[maxn],a,b,c,k1,k3;
double A[maxn],B[maxn],k2;
int t,n;
int main()
{
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d%lf%lf%lf%lf%lf%lf",&n,&k1,&k2,&k3,&a,&b,&c);
        memset(dp,0,sizeof(dp));
        memset(A,0,sizeof(A));
        memset(B,0,sizeof(B));
        dp[0]=1.0/k1/k2/k3;
        for(int i=1; i<=k1; i++)
            for(int j=1; j<=k2; j++)
                for(int k=1; k<=k3; k++)
                    if(i!=a||j!=b||k!=c)
                        dp[i+j+k]+=dp[0];
        for(int i=n; i>=0; i--)
        {
            A[i]+=dp[0];
            B[i]+=1;
            for(int j=3; j<=k1+k2+k3; j++)
            {
                A[i]+=A[i+j]*dp[j];
                B[i]+=B[i+j]*dp[j];
            }
        }
        printf("%.12lf\n",B[0]/(1-A[0]));
    }
    return 0;
}