結果

問題 No.1151 チャレンジゲーム
コンテスト
ユーザー Chanyuh
提出日時 2020-08-21 17:02:19
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 2,376 bytes
コンパイル時間 1,297 ms
コンパイル使用メモリ 124,360 KB
最終ジャッジ日時 2025-01-13 04:45:57
ジャッジサーバーID
（参考情報） 		judge5 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2 WA * 1
other AC * 8 WA * 42
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<iostream>
#include<string>
#include<cstdio>
#include<vector>
#include<cmath>
#include<algorithm>
#include<functional>
#include<iomanip>
#include<queue>
#include<ciso646>
#include<random>
#include<map>
#include<set>
#include<complex>
#include<bitset>
#include<stack>
#include<unordered_map>
#include<utility>
#include<tuple>
#include<cassert>
using namespace std;
typedef long long ll;
typedef unsigned int ui;
const ll mod = 1000000007;
const ll INF = (ll)1000000007 * 1000000007;
typedef pair<int, int> P;
#define stop char nyaa;cin>>nyaa;
#define rep(i,n) for(int i=0;i<n;i++)
#define per(i,n) for(int i=n-1;i>=0;i--)
#define Rep(i,sta,n) for(int i=sta;i<n;i++)
#define Per(i,sta,n) for(int i=n-1;i>=sta;i--)
#define rep1(i,n) for(int i=1;i<=n;i++)
#define per1(i,n) for(int i=n;i>=1;i--)
#define Rep1(i,sta,n) for(int i=sta;i<=n;i++)
typedef long double ld;
const ld eps = 1e-8;
const ld pi = acos(-1.0);
typedef pair<ll, ll> LP;
int dx[4]={1,-1,0,0};
int dy[4]={0,0,1,-1};

int n,N=1;
ld a[15];
ld dp[2][100050];
int po[15];

int inv(int S){
  int T=0;
  rep(i,n){
    if(S%3!=0){
      T+=po[i]*(3-S%3);
    }
    S/=3;
  }
  return T;
}

ld get_dp(int turn,int S){
  if(dp[turn][S]!=-1) return dp[turn][S];
  vector<int> used(n);
  int S_=S;
  bool flag=true;
  rep(i,n){
    used[i]=S_%3;
    if(S_%3==0) flag=false;
    S_/=3;
  }
  if(flag){
    int point=0;
    rep(i,n){
      if(used[i]==1) point+=a[i];
      else point-=a[i];
    }
    if(point>0) return 1;
    else return 0;
  }
  ld ma=0;
  rep(i,n){
    if(used[i]!=0) continue;
    ld mi=1;
    rep(j,n){
      if(used[j]!=0) continue;
      ld p=a[j]/(a[i]+a[j]-1);
      ld res=p*(get_dp(1-turn,S+(1+turn)*po[i]))+(1-p)*get_dp(turn,S+(2-turn)*po[j]);
      mi=min(mi,res);
    }
    ma=max(ma,mi);
  }
  dp[turn][S]=ma;
  return dp[turn][S];
}

void solve(){
  cin >> n;
  po[0]=1;
  rep(i,n){
    po[i+1]=3*po[i];
  }
  //cout << 5 << " " << inv(5) << endl;
  rep(i,n) cin >> a[i];
  rep(i,n) N*=3;
  rep(k,2) rep(i,N) dp[k][i]=-1;
  // rep(p,2) rep(S,N){
  //   cout << p << endl;
  //   int S_=S;
  //   rep(i,n){
  //     cout << S_%3 << " ";
  //     S_/=3;
  //   }
  //   // cout << "" << endl;
  //   // cout << get_dp(p,S) << endl;
  // }
  cout << get_dp(0,0) << endl;
}

int main(){
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout << fixed << setprecision(50);
    solve();
}
0