hdu 5045 Contest (状压DP)

状压dp，将m道题分段，每段有n道题，依次对每段进行状压dp即可。

#include<cstdio>
#define max(a,b) (a)>(b)?(a):(b)

bool flag[1005];
double p[13][1003],dp[1030][1003];

int check(int x)
{
    int ret=0,k=0;
    while(x>>k){
        if((x>>k)&1) ret++;
        k++;
    }
    return ret;
}

int main()
{
    int t,ncase=0;
    scanf("%d",&t);
    while(t--){
        int n,m,i,j,k,s,tmp,lun=0;
        scanf("%d%d",&n,&m);
        tmp=m;
        for(i=1;i<=n;i++) for(j=1;j<=m;j++)
            scanf("%lf",&p[i][j]);
        dp[0][0]=0.0;
		int ans=0;
        for(i=0;i<n;i++){
            ans |= (1<<i);
        }
        while(m>=n){
            for(j=1;j<=n;j++){
                for(i=1;i<=ans;i++){
                    if(check(i)==j){
			dp[i][j]=-1;
                        for(k=0;k<n;k++){
                            if(((1<<k)&i)){
                                s=i^(1<<k);
                                dp[i][j]=max(dp[i][j],dp[s][j-1]+p[k+1][j+n*lun]);
                            }
                        }
                    }
                }
            }
            //for(i=1;i<=n;i++) dp[0][i]=dp[1023][i];
			dp[0][0]=dp[ans][n];
            m-=n;
			lun++;
        }
	//	ans=0;
      //  for(i=0;i<m;i++) ans|=(1<<i);
        for(j=1;j<=m;j++){
            for(i=1;i<=ans;i++){
                if(check(i)==j){
					dp[i][j]=-1;
                    for(k=0;k<n;k++){
                        if(((1<<k)&i)){
                            s=i^(1<<k);
                            dp[i][j]=max(dp[i][j],dp[s][j-1]+p[k+1][j+n*lun]);
                        }
                    }
                }
            }
        }
		double aa=-1;
		for(i=0;i<=ans;i++) aa=max(aa,dp[i][m]);
        printf("Case #%d: %.5lf\n",++ncase,aa);
    }
    return 0;
}