ZOJ 3502 Contest <状态压缩 概率 DP>

链接：http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=3502

#include <iostream>
#include <cstdio>
#include <string>
#include <cmath>
#include <cstring>
using namespace std;
int T, N;
double p[15][15], dp[1<<15];
string s[1<<15];
const double eps=1e-8;
int Sign( double x )
{
    if( fabs(x)<eps )return 0;
    return x>0?1:-1;
}
int main( )
{
    scanf("%d", &T);
    while(T--){
        scanf("%d",&N );
        for(int i=0; i<N; ++ i){
            for( int j=0; j<N; ++j ){
                scanf("%lf", &p[i][j]);
            }
        }
       int t=1<<N;double q;

        memset(dp, 0, sizeof dp);
        if(N==1)s[1]="A"; // N=1 p=0 不能被更新到~
        for(int i=1; i<t; ++ i){
            for( int j=0; j<N; ++ j ){// 更新到当前位
                if( i&(1<<j)) {
                    q=0;
                    for( int k=0; k<N; ++ k ){ // 已经被更新过的
                        if(i&(1<<k))
                        q=max(q, p[k][j]);
                    }
                }
                q=dp[i^(1<<j)]+q/100.0;
                if( Sign(q-dp[i])>0 || Sign(q-dp[i])==0 && s[i^(1<<j)]<s[i]){
                    dp[i]=q;
                    s[i]=s[i^(1<<j)];
                    s[i]+=('A'+j);
                }
            }
        }

        printf("%.2f
%s
", dp[t-1], s[t-1].c_str());
    }
    return 0;
}