csu 1009 抛硬币

C - 抛硬币 CSU - 1009

James得到了一堆有趣的硬币，于是决定用这些硬币跟朋友们玩个小游戏。在一个N行M列的表格上，每一个第i行第j列的格子上都放有一枚James的硬币，抛该硬币正面朝上的概率为Pij，所有抛硬币事件两两之间是相互独立的。

现在，玩家在M列硬币中，从每一列里各选择1枚，共M枚，构成一组。如此重复选择N组出来，且保证被选择过的硬币不能再选。选好组之后，每组的M枚硬币各抛一次，如果都是正面朝上，则该组胜利，总分赢得1分；否则该组失败，总分不加也不减。请问，如果让你自行选择硬币的分组，游戏总得分的数学期望的最大值是多少？

Input
输入有多组数据。每组数据第一行为N和M，1≤N≤100，1≤M≤10，以空格分隔。接下来有N行，每行M个小数，表示表格中对应的Pij。

输入以N=M=0结束，这组数据不输出结果。

Output
对于每组数据，输出对应游戏总得分的数学期望的最大值，四舍五入精确至4位小数。每组数据的输出占一行。

Sample Input
2 3
1.0 1.0 1.0
0.5 0.4 0.3
0 0
Sample Output
1.0600

题意：输入一个n行m列的矩阵，每个数代表概率，想得到最大期望，而得分结果为1，即每次优先最大的概率相乘，局部最优，用贪心思想。

　　   在没列上排序，选取最大的概率，在行上相乘，再将每行的结果相加即可。

注意：要将矩阵转置才能正确sort排序，因为二维数组a[][]的第一个下标永远是行标。

#include<stdio.h>
#include<algorithm>
#include<math.h>
#include<string.h>
using namespace std;
const int maxn=1005;

double a[maxn][maxn];
double tmp[maxn][maxn];

int main()
{
    int n,m;
    while(~scanf("%d%d",&n,&m),n|m)
    {
        memset(a,0,sizeof(a));
        memset(tmp,0,sizeof(tmp));
        ///输入 n 行 m 列的矩阵
        for(int i=0;i<n;i++)
        {
            for(int j=0;j<m;j++)
            {
                scanf("%lf",&a[i][j]);
            }
        }

        ///转置为 m行 n 列的矩阵
       for(int j=0;j<m;j++)
       {
           for(int i=0;i<n;i++)
           {
               tmp[j][i]=a[i][j];
           }
       }
        ///排序
        for(int i=0;i<m;i++)
        {
            sort(tmp[i],tmp[i]+n);
        }

        ///输出转置后的矩阵变为 m行 n 列的矩阵
        /*for(int i=0;i<m;i++)
        {
            for(int j=0;j<n;j++)
            {
                printf("%.1lf ",tmp[i][j]);
            }
            printf("\n");
        }*/

        double sum=0,mad;
        for(int i=0;i<n;i++)
        {
            mad=1;
            for(int j=0;j<m;j++)
            {
                mad*=tmp[j][i];///在列上相乘
            }
            sum+=mad;///每列结果相加
        }
        printf("%.4lf\n",sum);
    }

    return 0;
}