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[][], b[];
 double c[][];
 int main()
 {
     //freopen("data.in", "r", stdin);
     //freopen("data.out", "w", stdout);
     int T, N, K,x, y, count;
     double  ans, num;

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

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

             memset(b, , sizeof(b));

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

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

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