#41746 (C++11) No.66 輝け☆全国たこやき杯

• 問題: No.66 輝け☆全国たこやき杯
• ユーザー: fiord
• 提出日時: 2015-08-05 23:57:07
• 言語: C++11 (gcc 11.4.0)
• 結果: AC
• 実行時間: 4 ms / 5,000 ms
• コード長: 655 bytes
• コンパイル時間: 1,335 ms
• コンパイル使用メモリ: 158,560 KB
• 実行使用メモリ: 5,376 KB
• 最終ジャッジ日時: 2024-07-18 03:33:24
• 合計ジャッジ時間: 1,839 ms
• ジャッジサーバーID （参考情報）: judge2 / judge4

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	2 ms; 5,248 KB
testcase_01	AC	1 ms; 5,376 KB
testcase_02	AC	2 ms; 5,376 KB
testcase_03	AC	2 ms; 5,376 KB
testcase_04	AC	2 ms; 5,376 KB
testcase_05	AC	2 ms; 5,376 KB
testcase_06	AC	2 ms; 5,376 KB
testcase_07	AC	2 ms; 5,376 KB
testcase_08	AC	3 ms; 5,376 KB
testcase_09	AC	4 ms; 5,376 KB

ソースコード

#include <bits/stdc++.h>
using namespace std;
double p[10][1<<10];
double calc(double a,double b){	return a*a/(a*a+b*b);	}
int main(){
	int m;	cin>>m;
	double s[1<<m];
	for(int i=0;i<(1<<m);i++)	cin>>s[i];
	for(int i=0;i<(1<<m);i+=2){
		p[0][i]=calc(s[i],s[i+1]);	p[0][i+1]=1-p[0][i];
	}
	for(int i=1;i<m;i++){
		for(int j=0;j<(1<<m);j++){
			int num=(j>>i);
			if(num%2==0){
				for(int k=((num+1)<<i);k<((num+2)<<i);k++){
					p[i][j]+=calc(s[j],s[k])*p[i-1][k];
				}
			}
			else{
				for(int k=((num-1)<<i);k<(num<<i);k++){
					p[i][j]+=calc(s[j],s[k])*p[i-1][k];
				}
			}
			p[i][j]*=p[i-1][j];
		}
	}
	printf("%.15lf\n",p[m-1][0]);
	return 0;
}