UVA 11181 Probability|Given

#include <map>
#include <set>
#include <list>
#include <cmath>
#include <ctime>
#include <deque>
#include <stack>
#include <queue>
#include <cctype>
#include <cstdio>
#include <string>
#include <vector>
#include <climits>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
#define LL long long
#define PI 3.1415926535897932626
using namespace std;
int gcd(int a, int b) {return a % b == 0 ? b : gcd(b, a % b);}
#define MAXN 25
double p[MAXN],ans[MAXN];
int res[MAXN];
int N,R;
double pb,pa;
void add()
{
    bool vis[MAXN];
    double tmp=1.0;
    memset(vis,false,sizeof(vis));
    for (int i=1;i<=R;i++)
    {
        vis[res[i]]=true;
        tmp*=p[res[i]];
    }
    for (int i=1;i<=N;i++)
      if (!vis[i])  tmp*=(1-p[i]);
    for (int i=1;i<=R;i++)
        ans[res[i]]+=tmp;
    pb+=tmp;
}
void dfs(int cur,int pre)
{
    if (cur>R)
    {
        //for (int i=1;i<=R;i++) printf("%d ",res[i]);puts("");
        add();
        return ;
    }
    for (int i=pre+1;i<=N;i++)
    {
        if (R-cur>N-i) break;
        res[cur]=i;
        dfs(cur+1,i);
    }
}
int main()
{
    int kase=1;
    while (scanf("%d%d",&N,&R)!=EOF)
    {
        if (N==0 && R==0) break;
        for (int i=1;i<=N;i++)
            scanf("%lf",&p[i]);
        memset(ans,0,sizeof(ans));
        pb=0;
        dfs(1,0);
        //printf("%lf\n",pb);
        printf("Case %d:\n",kase++);
        for (int i=1;i<=N;i++)
            printf("%.6lf\n",ans[i]/pb);
    }
    return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值