結果

問題 No.66 輝け☆全国たこやき杯
ユーザー kjnh10kjnh10
提出日時 2020-08-29 16:52:15
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 203 ms / 5,000 ms
コード長 3,319 bytes
コンパイル時間 2,244 ms
コンパイル使用メモリ 204,064 KB
実行使用メモリ 4,384 KB
最終ジャッジ日時 2023-08-09 07:52:22
合計ジャッジ時間 3,074 ms
ジャッジサーバーID
(参考情報)
judge14 / judge13
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,380 KB
testcase_01 AC 2 ms
4,380 KB
testcase_02 AC 1 ms
4,380 KB
testcase_03 AC 1 ms
4,380 KB
testcase_04 AC 2 ms
4,380 KB
testcase_05 AC 2 ms
4,376 KB
testcase_06 AC 4 ms
4,376 KB
testcase_07 AC 12 ms
4,376 KB
testcase_08 AC 47 ms
4,376 KB
testcase_09 AC 203 ms
4,384 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#line 2 "header.hpp"

//%snippet.set('header')%
//%snippet.fold()%
#ifndef HEADER_H
#define HEADER_H

// template version 2.0
using namespace std;
#include <bits/stdc++.h>

// varibable settings
const long long INF = 1e18;
template <class T> constexpr T inf = numeric_limits<T>::max() / 2.1;

#define _overload3(_1, _2, _3, name, ...) name
#define _rep(i, n) repi(i, 0, n)
#define repi(i, a, b) for (ll i = (ll)(a); i < (ll)(b); ++i)
#define rep(...) _overload3(__VA_ARGS__, repi, _rep, )(__VA_ARGS__)
#define _rrep(i, n) rrepi(i, 0, n)
#define rrepi(i, a, b) for (ll i = (ll)((b)-1); i >= (ll)(a); --i)
#define r_rep(...) _overload3(__VA_ARGS__, rrepi, _rrep, )(__VA_ARGS__)
#define each(i, a) for (auto &&i : a)
#define all(x) (x).begin(), (x).end()
#define sz(x) ((int)(x).size())
#define pb(a) push_back(a)
#define mp(a, b) make_pair(a, b)
#define mt(...) make_tuple(__VA_ARGS__)
#define ub upper_bound
#define lb lower_bound
#define lpos(A, x) (lower_bound(all(A), x) - A.begin())
#define upos(A, x) (upper_bound(all(A), x) - A.begin())
template <class T> inline void chmax(T &a, const T &b) { if ((a) < (b)) (a) = (b); }
template <class T> inline void chmin(T &a, const T &b) { if ((a) > (b)) (a) = (b); }
template <typename X, typename T> auto make_table(X x, T a) { return vector<T>(x, a); }
template <typename X, typename Y, typename Z, typename... Zs> auto make_table(X x, Y y, Z z, Zs... zs) { auto cont = make_table(y, z, zs...); return vector<decltype(cont)>(x, cont); }

#define cdiv(a, b) (((a) + (b)-1) / (b))
#define is_in(x, a, b) ((a) <= (x) && (x) < (b))
#define uni(x) sort(all(x)); x.erase(unique(all(x)), x.end())
#define slice(l, r) substr(l, r - l)

typedef long long ll;
typedef long double ld;
using vl = vector<ll>;
using vvl = vector<vl>;
using pll = pair<ll, ll>;

template <typename T>
using PQ = priority_queue<T, vector<T>, greater<T>>;
void check_input() { assert(cin.eof() == 0); int tmp; cin >> tmp; assert(cin.eof() == 1); }

#if defined(PCM) || defined(LOCAL)
#else
#define dump(...) ;
#define dump_1d(...) ;
#define dump_2d(...) ;
#define cerrendl ;
#endif

#endif /* HEADER_H */
//%snippet.end()%
#line 2 "solve.cpp"
template<class T=ll> using vec = vector<T>;
struct Fast { Fast() { std::cin.tie(0); ios::sync_with_stdio(false); } } fast;


int solve() {
    ll m;cin>>m;

    vec<double> s(1<<m);
    rep(i, 1<<m){
        cin>>s[i];
    }

    vec<vec<double>> dp = make_table(m+1, 1<<m, double(0));
    rep(i, 1<<m) dp[0][i] = 1;

    auto p = [&](ll i, ll j) -> double{
        // prob i win to j;
        return s[i]*s[i] / (s[i]*s[i] + s[j]*s[j]);
    };
    rep(x, 1, m+1){
        rep(y, 1<<m){
            // update dp[x][y]
            cerrendl;
            dump(x, y);
            ll b1 = y / (1 << (x-1));
            ll b2 = y / (1 << (x));
            dump(b1, b2);
            rep(rival, 1<<m){
                ll r1 = rival / (1 << (x-1));
                ll r2 = rival / (1 << (x));
                dump(rival, r1, r2);
                if (r1 != b1 && r2 == b2) {
                    dump(rival);
                    dp[x][y] += dp[x-1][y] * dp[x-1][rival] * p(y, rival);
                }
            }
        }
    }
    dump(dp);
    printf("%.12f\n", dp[m][0]);

    return 0; 
}


int main(){/*{{{*/
    solve();
    check_input();
    return 0;
}/*}}}*/
0