ACM选修hust 1075 组合+数学+期望值

Description

Input

Output

Sample Input

2
2 1
0 1
1 0
3 1
0 1 1
1 0 1
1 1 0

Sample Output

0.500
1.125

题意：给出N个node，当节点为女性，且认识K个男性时，该节点就满足条件，求出所有情况满足条件的节点的期望值。
想法：由于输入矩阵i X j来表示i是否认识j。故可以记录每个i认识的个数。枚举每个i，再进行排列组合计算，满足条件的个数计入count，计算期望值时除以所有情况(2^n)。
    （最开始愚蠢的我写了2的幂指函数，及非打表方式求组合函数，结果没有确认精度损失的问题WA无数次……捂脸）

#include <stdio.h>
#include <iostream>
#include <algorithm>
#include <string.h>
using namespace std;



int a[35][35], b[35];
double c[33][33];
int main()
{
    //freopen("data.in", "r", stdin);
    //freopen("data.out", "w", stdout);
    int T, N, K,x, y, count;
    double  ans, num;

    c[0][0]=1;
    for(int i = 1; i <= 30; i++)
    {
        c[i][0] = c[i][i]=1;
        for(int j = 1; j < i; j++)
        {
                c[i][j] = c[i-1][j-1] + c[i-1][j];
        }
    }
    scanf("%d", &T);

        while(T--)
        {
             scanf("%d%d",&N, &K);

            memset(b, 0, sizeof(b));


            for(int i=0; i<N; i++)
            {
                for(int j=0; j<N; j++)
                {
                    scanf("%d", &a[i][j]);
                        b[i]+=a[i][j];
                }
            }
            //printf("%d",b[4]);
            ans=0;
            for(int i = 0; i < N; i++)
            {
                if(b[i]>=K)
                {
                    num=0;
                    for(int j = K; j <= b[i]; j++)
                    {
                        num+=c[b[i]][j];
                    }
                    for(int j = 1; j <= b[i] + 1; j++)
                        num/=2.00;

                    ans+=num;
                }
            }
            printf("%.3lf
", ans);
        }
        return 0;


    //fclose(stdin);
    //fclose(stdout);
}