Bloodsucker ZOJ-3551 期望DP

最新推荐文章于 2021-05-01 22:51:21 发布

西伯利亚松鼠

最新推荐文章于 2021-05-01 22:51:21 发布

阅读量765

点赞数

分类专栏： # 概率/期望DP

本文链接：https://blog.csdn.net/qq_46173805/article/details/115606444

版权

概率/期望DP 专栏收录该内容

14 篇文章 0 订阅

订阅专栏

Bloodsucker

在这里插入图片描述

solution

$dp[i+1]=(dp[i]+1)*\frac{i*(n-i)p}{C_n^2}+(dp[i+1]+1)*(\frac{i(n-i)(1-p)}{C_n^2}+1-\frac{i(n-i)}{C_n^2})$

$dp[i+1]=dp[i]+\frac{C_n^2}{i(n-i)p}$

code

/*SiberianSquirrel*/
/*CuteKiloFish*/
#include <bits/stdc++.h>
using namespace std;

long double dp[101000];

void solve() {
    int o; cin >> o;
    while(o --) {
        int n; cin >> n;
        long double p; cin >> p;
        dp[1] = 0;
        for(int i = 1; i < n; ++ i) {
            dp[i + 1] = dp[i] + (1.0 * n * (n - 1) / 2.0) / (1.0 * i * (n - i) * p);
        }
        cout << fixed << setprecision(3) << dp[n] << endl;
    }
}

signed main() {
    ios_base::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
#ifdef ACM_LOCAL
    freopen("input", "r", stdin);
    freopen("output", "w", stdout);
#endif
    solve();
    return 0;
}