結果

問題 No.1151 チャレンジゲーム
ユーザー t9unkubjt9unkubj
提出日時 2024-10-01 02:57:37
言語 C++23
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 57 ms / 2,000 ms
コード長 6,790 bytes
コンパイル時間 3,899 ms
コンパイル使用メモリ 270,144 KB
実行使用メモリ 7,936 KB
最終ジャッジ日時 2024-10-01 02:57:44
合計ジャッジ時間 7,206 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,816 KB
testcase_01 AC 2 ms
6,820 KB
testcase_02 AC 2 ms
6,816 KB
testcase_03 AC 2 ms
6,816 KB
testcase_04 AC 2 ms
6,820 KB
testcase_05 AC 2 ms
6,816 KB
testcase_06 AC 2 ms
6,820 KB
testcase_07 AC 2 ms
6,820 KB
testcase_08 AC 2 ms
6,816 KB
testcase_09 AC 2 ms
6,820 KB
testcase_10 AC 2 ms
6,816 KB
testcase_11 AC 2 ms
6,820 KB
testcase_12 AC 2 ms
6,820 KB
testcase_13 AC 18 ms
6,816 KB
testcase_14 AC 18 ms
6,816 KB
testcase_15 AC 18 ms
6,816 KB
testcase_16 AC 18 ms
6,816 KB
testcase_17 AC 18 ms
6,820 KB
testcase_18 AC 18 ms
6,820 KB
testcase_19 AC 17 ms
6,820 KB
testcase_20 AC 18 ms
6,816 KB
testcase_21 AC 18 ms
6,820 KB
testcase_22 AC 18 ms
6,820 KB
testcase_23 AC 18 ms
6,820 KB
testcase_24 AC 18 ms
6,816 KB
testcase_25 AC 17 ms
6,816 KB
testcase_26 AC 18 ms
6,816 KB
testcase_27 AC 18 ms
6,816 KB
testcase_28 AC 55 ms
7,808 KB
testcase_29 AC 56 ms
7,808 KB
testcase_30 AC 56 ms
7,808 KB
testcase_31 AC 55 ms
7,680 KB
testcase_32 AC 56 ms
7,808 KB
testcase_33 AC 55 ms
7,808 KB
testcase_34 AC 57 ms
7,808 KB
testcase_35 AC 57 ms
7,808 KB
testcase_36 AC 56 ms
7,808 KB
testcase_37 AC 55 ms
7,808 KB
testcase_38 AC 55 ms
7,808 KB
testcase_39 AC 55 ms
7,936 KB
testcase_40 AC 57 ms
7,808 KB
testcase_41 AC 57 ms
7,808 KB
testcase_42 AC 56 ms
7,808 KB
testcase_43 AC 56 ms
7,808 KB
testcase_44 AC 56 ms
7,808 KB
testcase_45 AC 55 ms
7,808 KB
testcase_46 AC 56 ms
7,680 KB
testcase_47 AC 56 ms
7,808 KB
testcase_48 AC 50 ms
7,808 KB
testcase_49 AC 57 ms
7,808 KB
testcase_50 AC 57 ms
7,808 KB
testcase_51 AC 56 ms
7,680 KB
testcase_52 AC 56 ms
7,808 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#ifdef t9unkubj
#include"template.h"
//#include"template_no_debug.h"
#else
#undef _GLIBCXX_DEBUG
#pragma GCC optimize("O3")
#define dbg(...) 199958
using namespace std;
#include<bits/stdc++.h>
using uint=unsigned;
using ll=long long;
using ull=unsigned long long;
using ld=long double;
using pii=pair<int,int>;
using pll=pair<ll,ll>;
template<class T>using vc=vector<T>;
template<class T>using vvc=vc<vc<T>>;
template<class T>using vvvc=vvc<vc<T>>;
using vi=vc<int>;
using vvi=vc<vi>;
using vvvi=vc<vvi>;
using vl=vc<ll>;
using vvl=vc<vl>;
using vvvl=vc<vvl>;
template<class T>using smpq=priority_queue<T,vector<T>,greater<T>>;
template<class T>using bipq=priority_queue<T>;
#define rep(i,n) for(ll i=0;i<(ll)(n);i++)
#define REP(i,j,n) for(ll i=(j);i<(ll)(n);i++)
#define DREP(i,n,m) for(ll i=(n);i>=(m);i--)
#define drep(i,n) for(ll i=((n)-1);i>=0;i--)
#define all(x) x.begin(),x.end()
#define rall(x) x.rbegin(),x.rend()
#define mp make_pair
#define pb push_back
#define eb emplace_back
#define fi first
#define se second
#define is insert
#define bg begin()
#define ed end()
void scan(int&a) { cin >> a; }
void scan(ll&a) { cin >> a; }
void scan(string&a) { cin >> a; }
void scan(char&a) { cin >> a; }
void scan(uint&a) { cin >> a; }
void scan(ull&a) { cin >> a; }
void scan(bool&a) { cin >> a; }
void scan(ld&a){ cin>> a;}
template<class T> void scan(vector<T>&a) { for(auto&x:a) scan(x); }
void read() {}
template<class Head, class... Tail> void read(Head&head, Tail&... tail) { scan(head); read(tail...); }
#define INT(...) int __VA_ARGS__; read(__VA_ARGS__);
#define LL(...) ll __VA_ARGS__; read(__VA_ARGS__);
#define ULL(...) ull __VA_ARGS__; read(__VA_ARGS__);
#define STR(...) string __VA_ARGS__; read(__VA_ARGS__);
#define CHR(...) char __VA_ARGS__; read(__VA_ARGS__);
#define DBL(...) double __VA_ARGS__; read(__VA_ARGS__);
#define LD(...) ld __VA_ARGS__; read(__VA_ARGS__);
#define VC(type, name, ...) vector<type> name(__VA_ARGS__); read(name);
#define VVC(type, name, size, ...) vector<vector<type>> name(size, vector<type>(__VA_ARGS__)); read(name);
void print(int a) { cout << a; }
void print(ll a) { cout << a; }
void print(string a) { cout << a; }
void print(char a) { cout << a; }
void print(uint a) { cout << a; }
void print(bool a) { cout << a; }
void print(ull a) { cout << a; }
void print(double a) { cout << a; }
void print(ld a){ cout<< a; }
template<class T> void print(vector<T>a) { for(int i=0;i<(int)a.size();i++){if(i)cout<<" ";print(a[i]);}cout<<endl;}
void PRT() { cout <<endl; return ; }
template<class T> void PRT(T a) { print(a); cout <<endl; return; }
template<class Head, class... Tail> void PRT(Head head, Tail ... tail) { print(head); cout << " "; PRT(tail...); return; }
template<class T,class F>
bool chmin(T &x, F y){
    if(x>y){
        x=y;
        return true;
    }
    return false;
}
template<class T, class F>
bool chmax(T &x, F y){
    if(x<y){
        x=y;
        return true;
    }
    return false;
}
void YesNo(bool b){
    cout<<(b?"Yes":"No")<<endl;
}
void Yes(){
    cout<<"Yes"<<endl;
}
void No(){
    cout<<"No"<<endl;
}
template<class T>
int popcount(T n){
    return __builtin_popcountll(n);
}
template<class T>
T sum(vc<T>&a){
    return accumulate(all(a),T(0));
}
template<class T>
T max(vc<T>&a){
    return *max_element(all(a));
}
template<class T>
T min(vc<T>&a){
    return *min_element(all(a));
}
template<class T>
void unique(vc<T>&a){
    a.erase(unique(all(a)),a.end());
}
vvi readgraph(int n,int m,int off = -1){
    vvi g(n);
    rep(i, m){
        int u,v;
        cin>>u>>v;
        u+=off,v+=off;
        g[u].push_back(v);
        g[v].push_back(u);
    }
    return g;
}
vvi readtree(int n,int off=-1){
    return readgraph(n,n-1,off);
}
template<class T>
vc<T> presum(vc<T> &a){
    vc<T> ret(a.size()+1);
    rep(i,a.size())ret[i+1]=ret[i]+a[i];
    return ret;
}
template<class T, class F>
vc<T> &operator+=(vc<T> &a,F b){
    for (auto&v:a)v += b;
    return a;
}
template<class T, class F>
vc<T> &operator-=(vc<T>&a,F b){
    for (auto&v:a)v-=b;
    return a;
}
template<class T, class F>
vc<T> &operator*=(vc<T>&a,F b){
    for (auto&v:a)v*=b;
    return a;
}
#endif
double pass_time=0;
constexpr int three(int n){
    return (n?three(n-1)*3:1);
}
int at(int x,int i){
    return x/three(i)%3;
}
void solve(){
    INT(n);
    VC(ld,a,n);
    //こう取った時に勝つことのできる確率
    vvc<ld>dp(three(n),vc<ld>(2));
    drep(i,three(n)){
        int a1=0,a2=0;
        vi yet;
        rep(j,n){
            if(at(i,j)==0)yet.pb(j);
            if(at(i,j)==1)a1+=a[j];
            if(at(i,j)==2)a2+=a[j];
        }
        if(yet.empty()){
            rep(j,2)dp[i][j]=a1>a2;
        }else{
        rep(j,2){
        ld mx=0;
            for(auto&tar:yet){
                ld win_min=1;
                for(auto&tar2:yet){
            ld now=0;
            if(j==0){
                ld my=(a[tar2])/ld(a[tar]+a[tar2]-1);
                if(a[tar]==1)my=1;
                if(a[tar]!=1&&a[tar2]==1)my=1/a[tar];
                now+=dp[i+three(tar)][1]*my;
                now+=dp[i+three(tar2)*2][0]*(1-my);
            }
            else {
                ld his=(a[tar])/ld(a[tar]+a[tar2]-1);
                if(a[tar2]==1)his=1;
                if(a[tar2]!=1&&a[tar]==1)his=1.0/a[tar2];
                now+=dp[i+three(tar)][1]*(1-his);
                now+=dp[i+three(tar2)*2][0]*his;
            }
            
            chmin(win_min,now);
                }
            chmax(mx,win_min);
            }
            dp[i][j]=mx;
            
        }
            /*ld win=-1e9;
            int tar=-1;
            for(auto x:yet){
                if(chmax(win,dp[i+three(x)][1]/a[x])){
                    tar=x;
                }
            }
            ld win2=-1e9;
            int tar2=-1;
            for(auto x:yet){
                if(chmax(win2,1-(dp[i+three(x)*2][0]/a[x]))){
                    tar2=x;
                }
            }   
            {
                ld my=(a[tar2])/ld(a[tar]+a[tar2]-1);
                if(a[tar]==1)my=1;
                if(a[tar]!=1&&a[tar2]==1)my=1/a[tar];
                dp[i][0]+=dp[i+three(tar)][1]*my;
                dp[i][0]+=dp[i+three(tar2)*2][0]*(1-my);
            }
            {
                ld his=(a[tar])/ld(a[tar]+a[tar2]-1);
                if(a[tar2]==1)his=1;
                if(a[tar2]!=1&&a[tar]==1)his=1.0/a[tar2];
                dp[i][1]+=dp[i+three(tar)][1]*(1-his);
                dp[i][1]+=dp[i+three(tar2)*2][0]*his;
            }*/
        }
    }
    cout<<fixed<<setprecision(20);
    PRT(dp[0][0]);
}
signed main(){
    cin.tie(0)->sync_with_stdio(0);
    pass_time=clock();
    int t=1;
    //cin>>t;
    while(t--)solve();
    pass_time=clock()-pass_time;
    dbg(pass_time/CLOCKS_PER_SEC);
}
0