HDU 2017 新生赛 正品概率 (全概率公式 贝叶斯公式)

1.全概率公式:P(A)=P(B1)*P(A|B1)+P(B2)*P(A|B2)+P(B3)*P(A|B3)+......+P(Bn)*P(A|Bn);

2.贝叶斯公式:P(A|B)=P(B|A)*PA/P(B);

Problem Description

袋中有m枚正品硬币，n枚次品硬币（次品硬币两面都有国徽），在袋中任取一枚，将它投掷k次，已知每次得到的都是国徽，那么这枚硬币是正品的概率是多少？

Input

输入包含多组数据，每组占一行包含3个正整m,n,k (1<=m,n,k<=50)。

Output

每组输出一行，包含一个最简分数，硬币为正品的概率。

Sample Input

1 1 1

Sample Output

1/3

设：A硬币为正品的概率 B投K次都为国徽

求P(A|B)=P(B|A)*P(A)/P(B)

P(B|A)(硬币为正品 投k次都为正品的概率)=1/pow(2,k);

P(A)=m/m+n;(正品的概率)

P(B)(通过全概率公式求出) =n/(m+n)*1+1/pow(2,k)*m/(m+n);

所以 P(A|B)=(1/pow(2,k)*m/m+n)/n/(m+n)*1+1/pow(2,k)*m/(m+n)=m/m+n*pow(2,k);

#include <stdio.h>
#include<string.h>
#include<cmath>
#include<algorithm>
using namespace std;	
long long gcd(long long da,long long xiao)
{
	if(xiao==0)return da;
	return gcd(xiao,da%xiao);
}
int main(int argc, char *argv[])
{
	
	long long m,n;
	double k;
	while(scanf("%lld%lld%lf",&m,&n,&k)!=EOF)
	{
		long long ans=1;
		while(k--)
		{
			ans*=2;
		}		
		int x=gcd(m,(m+n*ans));
		printf("%lld/%lld\n",m/x,(m+n*ans)/x);
	}
	return 0;
}