1. 程式人生 > >CF540 D 概率 DP

CF540 D 概率 DP

names || main long long eof pair using 出現 continue

石頭剪刀布三種生物(?)

隨機戰鬥,最後三方各自只有自己方存活下來的概率是多少。

局面的轉移明顯,註意任選兩方決定戰鬥時有可能出現選了同個種類的,因此註意排除掉同種的組合,也就是條件概率什麽的。

/** @Date    : 2017-08-27 18:47:20
  * @FileName: D.cpp
  * @Platform: Windows
  * @Author  : Lweleth ([email protected])
  * @Link    : https://github.com/
  * @Version : $Id$
  */
#include <bits/stdc++.h>
#define LL long long
#define PII pair<int ,int>
#define MP(x, y) make_pair((x),(y))
#define fi first
#define se second
#define PB(x) push_back((x))
#define MMG(x) memset((x), -1,sizeof(x))
#define MMF(x) memset((x),0,sizeof(x))
#define MMI(x) memset((x), INF, sizeof(x))
using namespace std;

const int INF = 0x3f3f3f3f;
const int N = 1e5+20;
const double eps = 1e-8;


double f[110][110][110];


int main()
{
	int x, y, z;
	while(cin >> x >> y >> z)
	{
		MMF(f);
		f[x][y][z] = 1.00000000000;
		for(int i = x; i >= 0; i--)
		{
			for(int j = y; j >= 0; j--)
			{
				for(int k = z; k >= 0; k--)
				{
					if(!(i || j)|| !(i || k) || !(j || k))
						continue;
					if(i > 0)
						f[i - 1][j][k] += f[i][j][k] * (1.000000 * i * k)/(1.000000*(double)(i*k + j*i + k*j));
					if(j > 0)
						f[i][j - 1][k] += f[i][j][k] * (1.000000 * j * i)/(1.000000*(double)(i*k + j*i + k*j));
					if(k > 0)
						f[i][j][k - 1] += f[i][j][k] * (1.000000 * k * j)/(1.000000*(double)(i*k + j*i + k*j));
					//cout << f[i - 1][j][k] << "~" << f[i][j - 1][k] <<"~"<< f[i][j][k - 1] << endl;
				}
			}
		}
		double ans1, ans2, ans3;
		ans1 = ans2 = ans3 = 0;
		for(int i = 1; i <= x; i++)
			ans1 += f[i][0][0];
		for(int i = 1; i <= y; i++)
			ans2 += f[0][i][0];
		for(int i = 1; i <= z; i++)
			ans3 += f[0][0][i];
		printf("%.9lf %.9lf %.9lf\n", ans1, ans2, ans3);
	}
    return 0;
}

CF540 D 概率 DP