結果

問題 No.66 輝け☆全国たこやき杯
コンテスト
ユーザー pekempey
提出日時 2015-08-20 17:15:33
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 6 ms / 5,000 ms
コード長 877 bytes
コンパイル時間 1,280 ms
コンパイル使用メモリ 159,344 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-07-18 10:56:02
合計ジャッジ時間 1,637 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 10
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#define rep(i, a) for (int i = 0; i < (a); i++)
#define rep2(i, a, b) for (int i = (a); i < (b); i++)
#define repr(i, a) for (int i = (a) - 1; i >= 0; i--)
#define repr2(i, a, b) for (int i = (b) - 1; i >= (a); i--)
using namespace std;
typedef long long ll;
const ll inf = 1e9;
const ll mod = 1e9 + 7;

int M;
double S[101010];

double dp[20][101010];

double sq(double x) { 
	return x * x;
}

int main() {
	cin >> M;
	rep (i, 1 << M) cin >> S[i];

	rep (i, 1 << M) dp[0][i] = 1;

	rep (i, M) {
		rep (j, 1 << M) {
			int l = j / (1 << i + 1) * (1 << i + 1);
			int r = l + (1 << i + 1);
			int l2 = j / (1 << i) * (1 << i);
			int r2 = l2 + (1 << i);
			
			rep2 (k, l, r) if (k < l2 || r2 <= k) {
				double p = sq(S[j]) / (sq(S[j]) + sq(S[k]));
				dp[i + 1][j] += dp[i][j] * dp[i][k] * p; 
			}
		}	
	}

	printf("%.20f\n", dp[M][0]);

	return 0;
}
0