codeforces 678E Another Sith Tournament 概率dp

奉上官方题解

然后直接写的记忆化搜索

#include <cstdio>
#include <iostream>
#include <ctime>
#include <vector>
#include <cmath>
#include <map>
#include <queue>
#include <algorithm>
#include <cstring>
using namespace std;
typedef long long LL;
const int N=5e2+;
const int INF=0x3f3f3f3f;
const int mod=1e9+;
double p[][],dp[<<][];
int n;
double gao(int s,int cur){
  if(dp[s][cur]!=-1.0) return dp[s][cur];
  for(int i=;i<n;++i){
     if(i!=cur&&(s&(<<i))){
       dp[s][cur]=max(dp[s][cur],p[cur][i]*gao(s^(<<i),cur)+p[i][cur]*gao(s^(<<cur),i));
     }
  }
  return dp[s][cur];
}
int main(){
   scanf("%d",&n);
   for(int i=;i<n;++i){
    for(int j=;j<n;++j)
      scanf("%lf",&p[i][j]);
   }
   for(int i=;i<(<<n);++i)
    for(int j=;j<;++j)
       dp[i][j]=-;
   for(int i=;i<n;++i)
    dp[(<<i)][i]=;
   dp[][]=;
   double ans=;
   for(int i=;i<n;++i){
     ans=max(ans,gao((<<n)-,i));
   }
   printf("%.10f\n",ans);
   return ;
}