POJ Monthly--2006.12.31, flymouse King of Fighters

【题目描述】

两个大牛玩拳皇，N v N的比赛，总共M个角色，给出牛A的N个角色对阵牛C的M个角色时赢的概率，求牛C N场比赛全胜的最大概率

【输入输出】

多case

n <= 10 m <= 100000

p11 p12 p13 .... p1m

...

pn1 pn2 pn3 .... pnm

EOF结束

每个用例输出最大概率。

【题目分析】

其实就是说从一个N×M的矩阵里每一行选一个数，这些数不能在同一列，使得它们的乘积最大

鉴于M太大，先减小下规模

对于牛A的每个角色，牛C肯定会从M个角色里挑选N个胜算最大的角色中的一个，所以最后M个角色里最多挑出N × N个角色，里面肯定有使概率最大的N个不同的角色，这下办法就多了，可以直接集合DP，O(N * N * N * 2 ^ N ) <= 1024000，或者使用KM算法，亦或是费用流，因为是求乘积的，用后两种方法的时候可以先log10一下，最后答案再pow(10)一下，题目给出的矩阵是牛A胜牛C的，记得用1减之。我用个第一种方法，代码比较简单

【code】

#include <stdio.h> #include <algorithm> using namespace std; char v[100000], flag; int n, m, now, index[100000], s[100]; double tmp[10][100000], d[2][1024]; bool cmp(int a, int b) { return tmp[now][a] > tmp[now][b]; } int main() { int i, j, q; while (scanf("%d%d", &n, &m) == 2) { flag++; for (i = 0; i < n; i++) { now = i; for (j = 0; j < m; j++) { index[j] = j; scanf("%lf", &tmp[i][j]); tmp[i][j] = 1.0 - tmp[i][j]; } sort(index, index + m, cmp); for (j = 0; j < n; j++) { v[index[j]] = flag; } } for (i = 0, j = 0; i < m; i++) { if (v[i] == flag) { s[j++] = i; } } m = j; int x = 1 << n; memset(d[0], 0, sizeof d[0]); d[0][0] = 1; for (i = 1, q = 1; i <= m; i++, q = !q) { int p = s[i - 1]; for (j = 0; j < x; j++) { d[q][j] = d[!q][j]; for (int k = 0; k < n; k++) { if (j & (1 << k)) { d[q][j] = max(d[q][j], d[!q][j ^ (1 << k)] * tmp[k][p]); } } } } printf("%.15lf/n", d[!q][x - 1]); } return 0; }