結果
| 問題 | 
                            No.66 輝け☆全国たこやき杯
                             | 
                    
| コンテスト | |
| ユーザー | 
                             fiord
                         | 
                    
| 提出日時 | 2015-08-05 23:57:07 | 
| 言語 | C++11(廃止可能性あり)  (gcc 13.3.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 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| other | AC * 10 | 
ソースコード
#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;
}
            
            
            
        
            
fiord