骰子【概率dp】

最新推荐文章于 2024-09-07 22:39:09 发布

Wuliwuliii

最新推荐文章于 2024-09-07 22:39:09 发布

阅读量334

点赞数 1

分类专栏： DP动态规划文章标签：概率 dp

本文链接：https://blog.csdn.net/qq_41730082/article/details/106349320

版权

DP动态规划专栏收录该内容

110 篇文章 0 订阅

订阅专栏

题目链接 P1409 骰子

因为会有人被弹出队列，所以我设置的期望dp为 dp[i][j] ，表示当现在队列中有i个人的时候，第j个人获胜的概率。

于是有当只剩一个人的时候，那个人必胜， dp[1][1] = 1 。

再往下，先看它在队首的情况， $dp[i][1] = \frac{1}{6} + \frac{3}{6} *dp[i][i]$ 也就是直接获胜的概率加上它被弹到队尾时候的概率。

其他的情况呢，也就是不在队首的时候呢， $dp[i][j] = \frac{2}{6} * dp[i - 1][j - 1] + \frac{3}{6} * dp[i][j - 1]$ ，表示的是如果第一个被弹出队列，或者第一个被弹到队尾时候，第j-1个获胜，也就是现在的j获胜，因为j的位置就变成了j-1。

于是，连立这两个方程组，可以得到 dp[i][1] 的求解方式，因为这里i比较的不清晰，接下去的式子我用x代替i表示。

$dp[x][1] = \frac{1}{3} * \frac{1}{2^x - 1} * (2^{x - 1} + 2^{x - 1} * dp[x - 1][x - 1] + \cdots + 2^{k} * dp[x - 1][k] + \cdots + 2^1 * dp[x - 1][1])$

所以，就可以利用一个pow来处理double型来解决这个问题了。

#include <iostream>
#include <cstdio>
#include <cmath>
#include <string>
#include <cstring>
#include <algorithm>
#include <limits>
#include <vector>
#include <stack>
#include <queue>
#include <set>
#include <map>
#include <bitset>
//#include <unordered_map>
//#include <unordered_set>
#define lowbit(x) ( x&(-x) )
#define pi 3.141592653589793
#define e 2.718281828459045
#define INF 0x3f3f3f3f
#define HalF (l + r)>>1
#define lsn rt<<1
#define rsn rt<<1|1
#define Lson lsn, l, mid
#define Rson rsn, mid+1, r
#define QL Lson, ql, qr
#define QR Rson, ql, qr
#define myself rt, l, r
using namespace std;
typedef unsigned long long ull;
typedef unsigned int uit;
typedef long long ll;
const int maxN = 1e3 + 7;
int N, M;
double dp[maxN][maxN];
signed main()
{
    scanf("%d%d", &N, &M);
    dp[1][1] = 1.;
    double tmp;
    for(int i=2; i<=N; i++)
    {
        tmp = 1. / 3. / (pow(2., i) - 1.);
        dp[i][1] = tmp * pow(2., i - 1);
        for(int j=i-1; j>=1; j--) dp[i][1] += tmp * pow(2, j) * dp[i - 1][j];
        for(int j=2; j<=i; j++)
        {
            dp[i][j] = dp[i - 1][j - 1] / 3. + dp[i][j - 1] / 2.;
        }
    }
    printf("%.9lf\n", dp[N][M]);
    return 0;
}