結果

問題 No.66 輝け☆全国たこやき杯
ユーザー EmKjpEmKjp
提出日時 2015-03-22 10:46:45
言語 C++11
(gcc 13.3.0)
結果
AC  
実行時間 3 ms / 5,000 ms
コード長 1,271 bytes
コンパイル時間 649 ms
コンパイル使用メモリ 81,632 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-06-29 00:00:07
合計ジャッジ時間 1,304 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 10
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:54:22: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   54 |     FOR(i,1<<M) scanf("%lld",&S[i]);
      |                 ~~~~~^~~~~~~~~~~~~~

ソースコード

diff # 

#include<iostream>
#include<sstream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<string>
#include<vector>
#include<cmath>
#include<set>
#include<map>
#include<stack>
#include<queue>
#include<numeric>
#include<functional>
#include<complex>

using namespace std;
#define BET(a,b,c) ((a)<=(b)&&(b)<(c))
#define FOR(i,n) for(int i=0,i##_end=(int(n));i<i##_end;i++)
#define SZ(x) (int)(x.size())
#define ALL(x) (x).begin(),(x).end()
#define MP make_pair
#define FOR_EACH(it,v) for(__typeof(v.begin()) it=v.begin(),it_end=v.end() ; it != it_end ; it++)
typedef vector<int> VI;
typedef vector<VI> VVI;

int M;
long long S[1024];

double memo[1024][11];

double win(int p, int q){
    return 1.0 * S[p] * S[p] / (1.0 * S[p] * S[p] + 1.0 * S[q] * S[q]);
}

double dfs(int pos, int M){
    if(M == 0) return 1.0;
    double& ret = memo[pos][M];
    if(!std::isnan(ret)) return ret;
    double val = 0;
    int r = 1<<M;
    int s = pos / r * r;
    if(pos % r < r / 2) s += r / 2;
    FOR(i,r / 2) {
        val += dfs(s + i, M - 1) * win(pos, s + i);
    }
    val *= dfs(pos, M - 1);
    return ret = val;
}

int main()
{
    cin>>M;
    FOR(i,1<<M) scanf("%lld",&S[i]);
    memset(memo , -1 , sizeof(memo));
    printf("%.10f\n", dfs(0, M));
    return 0;
}
0