結果

問題 No.66 輝け☆全国たこやき杯
コンテスト
ユーザー firiexp
提出日時 2019-10-24 00:37:23
言語 C++17(gcc12)
(gcc 12.4.0 + boost 1.89.0)
コンパイル:
g++-12 -O2 -lm -std=c++17 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
AC  
実行時間 4 ms / 5,000 ms
コード長 1,351 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 1,065 ms
コンパイル使用メモリ 115,508 KB
実行使用メモリ 7,844 KB
最終ジャッジ日時 2026-03-08 21:45:00
合計ジャッジ時間 1,671 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 10
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:33:29: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   33 |     for (auto &&i : v) scanf("%lg", &i);
      |                        ~~~~~^~~~~~~~~~~

ソースコード

diff #
raw source code 

#include <iostream>
#include <algorithm>
#include <iomanip>
#include <map>
#include <set>
#include <queue>
#include <stack>
#include <numeric>
#include <bitset>
#include <cmath>

static const int MOD = 1000000007;
using ll = long long;
using u32 = uint32_t;
using namespace std;

template<class T> constexpr T INF = ::numeric_limits<T>::max()/32*15+208;

template <class T, class U>
vector<T> make_v(U size, const T& init){ return vector<T>(static_cast<size_t>(size), init); }

template<class... Ts, class U>
auto make_v(U size, Ts... rest) { return vector<decltype(make_v(rest...))>(static_cast<size_t>(size), make_v(rest...)); }

template<class T> void chmin(T &a, const T &b){ a = (a < b ? a : b); }
template<class T> void chmax(T &a, const T &b){ a = (a > b ? a : b); }


int main() {
    int k;
    cin >> k;
    vector<double> v(1<<k);
    for (auto &&i : v) scanf("%lg", &i);
    auto dp = make_v(k+1, 1<<k, 0.0);
    fill(dp[0].begin(), dp[0].end(), 1.0);
    int mask = (1<<k)-1;
    for (int i = 0; i < k; ++i) {
        for (int j = 0; j < (1<<k); ++j) {
            int j2 = (j ^ (1<<i)) & mask;
            for (int l = 0; l < (1<<i); ++l) {
                dp[i+1][j] += dp[i][j]*dp[i][j2+l] *(v[j]*v[j])/(v[j]*v[j]+v[j2+l]*v[j2+l]);
            }
        }
        mask ^= (1<<i);
    }
    printf("%.10f\n", dp[k][0]);
    return 0;
}
0