CDOJ 1135 邱老师看电影概率dp

概率DP算法详解

最新推荐文章于 2023-10-28 17:52:33 发布

原创最新推荐文章于 2023-10-28 17:52:33 发布 · 611 阅读

0 ·

CC 4.0 BY-SA版权

CDOJ 同时被 2 个专栏收录

64 篇文章

订阅专栏

数据结构

32 篇文章

订阅专栏

本文介绍了一种概率动态规划算法，并通过具体的代码实现了解决特定问题的过程。从基本的边界条件开始，逐步推导到复杂的状态转移方程，展示了如何利用递推公式计算不同状态的概率。

概率dp，推一个概率公式就行。i==2时那里写错了一个+-号，wa了两发

代码：

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <vector>
#include <algorithm>
using namespace std;
#define maxn 1005
double dp[maxn][maxn];
int w, b;
int main()
{
	//freopen("input.txt", "r", stdin);
	scanf("%d%d", &w, &b);
	dp[0][0] = 0;
	for (int i = 1; i < maxn; ++i)
		dp[i][0] = 1, dp[0][i] = 0;
	for (int i = 1; i < maxn; ++i)
		dp[i][1] = i*1.0 / (i + 1);
	dp[1][2] = 1.0 / 3;
	for (int i = 2; i < maxn; ++i)
		dp[i][2] = i*1.0 / (i + 2) + 2.0 / (i + 2) / (i + 1);
	for (int i = 1; i < maxn; ++i)
	{
		for (int j = 3; j < maxn; ++j)
		{
			dp[i][j] = i*1.0 / (i + j) + j*1.0*(j - 1) / (i + j) / (i + j - 1)*((j - 2)*1.0 / (i + j - 2)*dp[i][j - 3] + i*1.0 / (i + j - 2)*dp[i - 1][j - 2]);
		}
	}
	printf("%.9f\n", dp[w][b]);
	//system("pause");
	//while (1);
	return 0;
}