hdoj 5045 （dp 状态压缩）

题意：

给定n个人 m个题目

下面n*m的矩阵表示每个人解出每道题的概率

我们可以得到一个模长为m的集合{1,2,3,1,2}

代表每道题是谁解出的。

有众多集合，且获得这个集合有一个期望，求期望最大的那个集合 （ 的期望值是多少）

一个限制: 对于集合 {1,1,2,3,1} 这样是不合法的（即从[1,n]题必须是1-n的排列，然后[n+1, 2n]题也是一个排列）

然后状压dp

/* ***********************************************
Author        :kuangbin
Created Time  :2014/9/10 20:04:53
File Name     :D.cpp
************************************************ */

#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>
#include <vector>
#include <queue>
#include <set>
#include <map>
#include <string>
#include <math.h>
#include <stdlib.h>
#include <time.h>
using namespace std;
double dp[1010][1<<10];
bool vis[1010][1<<10];
double p[10][1010];
int n,m;
double solve(int id,int s){
	if(s == (1<<n)-1)s = 0;//使得每过N道题形成一周期
	if(vis[id][s])return dp[id][s];
	vis[id][s] = true;
	dp[id][s] = 0;
	if(id == m)return dp[id][s];
	for(int i = 0;i < n;i++){
		if(s&(1<<i))continue;
		dp[id][s] = max(dp[id][s],(1-p[i][id])*solve(id+1,s|(1<<i))+p[i][id]*(solve(id+1,s|(1<<i))+1));
	}
	return dp[id][s];
}
int main()
{
    //freopen("D.in","r",stdin);
    //freopen("D.out","w",stdout);
    int T;
	int iCase = 0;
	scanf("%d",&T);
	while(T--){
		iCase++;
		scanf("%d%d",&n,&m);
		for(int i = 0;i < n;i++)
			for(int j = 0;j < m;j++)
				scanf("%lf",&p[i][j]);
		memset(vis,false,sizeof(vis));
		printf("Case #%d: %.5lf\n",iCase,solve(0,0));
	}
    return 0;
}