HDU-5045 Contest 状态压缩DP求期望

N个人，M道题，M个小时，每个人做一道题需要1个小时。给出一个N*M的矩阵代表每个人做对每道题的概率。然后要求在任何时刻，任意两个人的敲题时间差不能大于1，也就是说，m道题要分成多段长度为n的最优排列，n为10，2^10=1024 1024*1000  状压即可。

#include <iostream>
#include <cstdio>
#include <cmath>
#include <queue>
#include <cstdlib>
#include <cstring>
#include <algorithm>
using namespace std;
int n,m;
double ans;
double p[1111][1111];
double dp[1111][1111];
int main()
{
	int t;
	int cas=1;
	scanf("%d",&t);
	while(t--)
	{
		for(int i=0;i<=1024;i++)
		{
			for(int j=0;j<=1000;j++)
			{
				dp[i][j]=-1;
			}
		}
		ans=-1;
		scanf("%d%d",&n,&m);
		for(int i=1;i<=n;i++)
		{
			for(int j=1;j<=m;j++)
			{
				scanf("%lf",&p[i][j]);
			}
		}
		int s=(1<<n)-1;//N个全部选过了 
		dp[0][0]=0;
		for(int j=1;j<=m;j++)
		{
			for(int i=0;i<s;i++)
			{
				for(int k=1;k<=n;k++)
				{
					int temp=1<<(k-1);//选到第i个人的时候 
					if(dp[i][j-1]<0)
					{
						continue;
					}
					if(temp&i)//这个人已经被选过 
					{
						continue;
					}
					int s1=i|temp;//把这个人标记为选过 
					if(s1==s)//因为时间只差不能超过1，因此是多段n的排列，然后从0开始 
					{
						s1=0;
					}
					dp[s1][j]=max(dp[s1][j],dp[i][j-1]+p[k][j]);
					if(j==m)
					{
						ans=max(ans,dp[s1][j]);
					}
				}
			}
		}
		printf("Case #%d: %.5f\n",cas++,ans);
	}
	return 0;
}