https://blog.csdn.net/mengxiang000000/article/details/69193806
題目大意:
現(xiàn)在有石頭r人锣光,剪刀s人,布p人.
每次會(huì)有兩個(gè)不同陣營(yíng)的人見(jiàn)面惑畴,然后一個(gè)人可能會(huì)出局辐棒。
問(wèn)最終三個(gè)隊(duì)獲勝的幾率。
一個(gè)隊(duì)獲勝意味著其他隊(duì)的人都死了秕噪。
思路:
觀察到r钳降,s,p均小于100.那么很明顯巢价,設(shè)定dp【i】【j】【k】表示石頭還剩下i人牲阁,剪刀還剩下j人固阁,布還剩下k人的幾率。
那么就有:
①dp【i-1】【j】【k】+=dp【i】【j】【k】*這種情況的概率城菊;
②dp【i】【j-1】【k】+=dp【i】【j】【k】*這種情況的概率备燃;
③dp【i】【j】【k-1】+=dp【i】【j】【k】*這種情況的概率;
即遞推出所有相遇的情況
很明顯凌唬,三種情況的概率分別為:
①ik/(ij+jk+ik)并齐;石頭和布一起出現(xiàn)的時(shí)候,石頭才可能死一個(gè)人客税。
②ij/(ij+jk+ik)况褪;石頭和剪子一起出現(xiàn)的時(shí)候,剪子才可能死一個(gè)人更耻。
③jk/(ij+jk+ik)测垛;剪子和布一起出現(xiàn)的時(shí)候,布才可能死一個(gè)人秧均。
那么過(guò)程維護(hù)一下dp數(shù)組食侮,注意三層for的枚舉方向即可。
代碼:
#include<algorithm>
#include <iostream>
#include <cstdlib>
#include <cstring>
#include <cassert>
#include <cstdio>
#include <vector>
#include <string>
#include <cmath>
#include <queue>
#include <stack>
#include <set>
#include <map>
using namespace std;
#define P(a,b,c) make_pair(a,make_pair(b,c))
#define rep(i,a,n) for (int i=a;i<=n;i++)
#define per(i,a,n) for (int i=n;i>=a;i--)
#define all(x) (x).begin(),(x).end()
#define SZ(x) ((int)(x).size())
#define pb push_back
#define mp make_pair
#define fi first
#define se second
//#define INF 0x3f3f3f3f
typedef pair<int,pair<int,int> >pii;//注意空格呦目胡!
typedef long long ll;
const ll mod = 1000000007;
ll gcd(ll a, ll b) { return b ? gcd(b, a%b) : a; }
const int INF=0x7f7f7f7f;
const int maxn = 1e5+5;
double dp[105][105][105];
double cal(int a,int b,int c)
{
double aa=(double)a;
double bb=(double)b;
double cc=(double)c;
if(aa*bb+bb*cc+cc*aa==0)return 0;
else
{
return (aa*bb)/(aa*bb+bb*cc+cc*aa);
}
}
int main()
{
int r,s,p;
while(~scanf("%d%d%d",&r,&s,&p))
{
memset(dp,0,sizeof(dp));
dp[r][s][p]=1;
per(i,0,r)
{
per(j,0,s)
{
per(k,0,p)
{
if(i-1>=0)dp[i-1][j][k]+=dp[i][j][k]*cal(i,k,j);
if(j-1>=0)dp[i][j-1][k]+=dp[i][j][k]*cal(i,j,k);
if(k-1>=0)dp[i][j][k-1]+=dp[i][j][k]*cal(j,k,i);
}
}
}
double a,b,c;
a=0;b=0;c=0;
for(int i=0;i<=r;i++)a+=dp[i][0][0];
for(int i=0;i<=s;i++)b+=dp[0][i][0];
for(int i=0;i<=p;i++)c+=dp[0][0][i];
printf("%.12lf %.12lf %.12lf\n",a,b,c);
}
return 0;
}