链接
https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=2122
题解
条件概率
P
(
A
∣
B
)
=
P
(
A
B
)
P
(
B
)
P(A|B)=\frac{P(AB)}{P(B)}
P(A∣B)=P(B)P(AB),其中
P
(
A
∣
B
)
P(A|B)
P(A∣B)表示
A
A
A在
B
B
B的条件下发生的概率,那也就是说在所有满足条件
B
B
B的事件中,也满足条件
A
A
A的事件所占的比例,那也就是
P
(
A
B
)
P
(
B
)
\frac{P(AB)}{P(B)}
P(B)P(AB)
这题就直接爆搜出所有的情况,算一下
P
(
A
B
)
P(AB)
P(AB)和
P
(
B
)
P(B)
P(B)就好了
代码
#include <bits/stdc++.h>
#define eps 1e-8
#define maxn 100
#define cl(x) memset(x,0,sizeof(x))
using namespace std;
typedef long long ll;
ll n, r, s[maxn];
double p[maxn], ans[maxn];
ll read(ll x=0)
{
ll c, f=1;
for(c=getchar();!isdigit(c);c=getchar())if(c=='-')f=-f;
for(;isdigit(c);c=getchar())x=x*10+c-48;
return f*x;
}
void dfs(ll now, ll cnt, double P)
{
ll i;
if(now>n)
{
if(cnt!=r)return;
for(i=1;i<=n;i++)if(s[i]==1)ans[i]+=P;
ans[0]+=P;
return;
}
if(n-now<r-cnt-1 or cnt>r)return;
dfs(now+1,cnt,P*(1.0-p[now]));
s[now]=1;
dfs(now+1,cnt+1,P*p[now]);
s[now]=0;
}
int main()
{
ll i, kase=0;
while(n=read(), r=read(), n or r)
{
for(i=1;i<=n;i++)scanf("%lf",p+i);
cl(ans);
dfs(1,0,1.0);
printf("Case %lld:\n",++kase);
for(i=1;i<=n;i++)printf("%.6lf\n",ans[i]/ans[0]);
}
return 0;
}