結果
問題 | No.66 輝け☆全国たこやき杯 |
ユーザー |
|
提出日時 | 2019-04-15 16:37:17 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 4 ms / 5,000 ms |
コード長 | 1,773 bytes |
コンパイル時間 | 1,570 ms |
コンパイル使用メモリ | 168,656 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-09-22 07:57:20 |
合計ジャッジ時間 | 2,076 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 10 |
ソースコード
// includes#include <bits/stdc++.h>// macros#define ll long long int#define pb emplace_back#define mk make_pair#define pq priority_queue#define FOR(i, a, b) for(int i=(a);i<(b);++i)#define rep(i, n) FOR(i, 0, n)#define rrep(i, n) for(int i=((int)(n)-1);i>=0;i--)#define all(x) (x).begin(),(x).end()#define sz(x) ((int)(x).size())#define UNIQUE(v) v.erase(unique(v.begin(), v.end()), v.end())#define FI first#define SE secondusing namespace std;// typestypedef pair<int, int> P;typedef pair<ll, int> Pl;typedef pair<ll, ll> Pll;typedef pair<double, double> Pd;// constantsconst int inf = 1e9;const ll linf = 1LL << 50;const double EPS = 1e-10;const int mod = 1e9 + 7;// solvetemplate <class T>bool chmax(T &a, const T &b){if(a < b){a = b; return 1;} return 0;}template <class T>bool chmin(T &a, const T &b){if(a > b){a = b; return 1;} return 0;}template <typename T> istream &operator>>(istream &is, vector<T> &vec){for(auto &v: vec)is >> v; return is;}double s[1026];double dp[1026][11];int main(int argc, char const* argv[]){ios_base::sync_with_stdio(false);cin.tie(0);int m;cin >> m;rep(i, (1 << m)){cin >> s[i];}rep(i, (1 << m))dp[i][0] = 1.;for(int j = 1; j <= m; j++){for(int i = 0; i < (1 << m); i++){int l = (i / (1 << j)) * (1 << j);int r = l + (1 << j);if(i < (l + r) / 2){for(int k = (l + r) / 2; k < r; k++){dp[i][j] += s[i] * s[i] / (s[i] * s[i] + s[k] * s[k]) * dp[k][j-1] * dp[i][j-1];}}else{for(int k = l; k < (l + r) / 2; k++){dp[i][j] += s[i] * s[i] / (s[i] * s[i] + s[k] * s[k]) * dp[k][j-1] * dp[i][j-1];}}}}cout << setprecision(20);cout << dp[0][m] << endl;return 0;}