結果
問題 | No.66 輝け☆全国たこやき杯 |
ユーザー |
![]() |
提出日時 | 2014-11-13 23:32:08 |
言語 | C++11 (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 10 ms / 5,000 ms |
コード長 | 1,912 bytes |
コンパイル時間 | 1,544 ms |
コンパイル使用メモリ | 160,064 KB |
実行使用メモリ | 6,820 KB |
最終ジャッジ日時 | 2024-12-31 10:17:57 |
合計ジャッジ時間 | 2,191 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 10 |
ソースコード
#include<bits/stdc++.h>using namespace std;#define REP(i,a,b) for(i=a;i<b;i++)#define rep(i,n) REP(i,0,n)#define mygc(c) (c)=getchar_unlocked()#define mypc(c) putchar_unlocked(c)#define ll long long#define ull unsigned llvoid reader(int *x){int k,m=0;*x=0;for(;;){mygc(k);if(k=='-'){m=1;break;}if('0'<=k&&k<='9'){*x=k-'0';break;}}for(;;){mygc(k);if(k<'0'||k>'9')break;*x=(*x)*10+k-'0';}if(m)(*x)=-(*x);}void reader(int *x, int *y){reader(x);reader(y);}void reader(int *x, int *y, int *z){reader(x);reader(y);reader(z);}void reader(ll *x){int k,m=0;*x=0;for(;;){mygc(k);if(k=='-'){m=1;break;}if('0'<=k&&k<='9'){*x=k-'0';break;}}for(;;){mygc(k);if(k<'0'||k>'9')break;*x=(*x)*10+k-'0';}if(m)(*x)=-(*x);}int reader(char c[]){int i,s=0;for(;;){mygc(i);if(i!=' '&&i!='\n'&&i!='\r'&&i!='\t'&&i!=EOF) break;}c[s++]=i;for(;;){mygc(i);if(i==' '||i=='\n'||i=='\r'||i=='\t'||i==EOF) break;c[s++]=i;}c[s]='\0';return s;}void writer(int x, char c){int s=0,m=0;char f[10];if(x<0)m=1,x=-x;while(x)f[s++]=x%10,x/=10;if(!s)f[s++]=0;if(m)mypc('-');while(s--)mypc(f[s]+'0');mypc(c);}void writer(ll x, char c){int s=0,m=0;char f[20];if(x<0)m=1,x=-x;while(x)f[s++]=x%10,x/=10;if(!s)f[s++]=0;if(m)mypc('-');while(s--)mypc(f[s]+'0');mypc(c);}void writer(const char c[]){int i;for(i=0;c[i]!='\0';i++)mypc(c[i]);}int M;int S[100000];double dp[100000], nx[10000];double winp(int a, int b){return pow(S[a],2) / (pow(S[a],2) + pow(S[b],2));}int main(){int i, j, k, per = 1, st, ed;reader(&M);rep(i,1<<M) reader(S+i);rep(i,1<<M) dp[i] = 1;rep(k,M){rep(i,1<<M) nx[i] = 0;per *= 2;rep(i,1<<M){st = i / per * per;ed = st + per;REP(j,st,ed) if(j/(per/2) != i/(per/2)){nx[i] += dp[i] * dp[j] * winp(i,j);}}rep(i,1<<M) dp[i] = nx[i];}printf("%.15f\n",dp[0]);return 0;}