結果

問題 No.1151 チャレンジゲーム
ユーザー 37zigen37zigen
提出日時 2020-08-08 18:14:56
言語 PyPy3
(7.3.15)
結果
RE  
実行時間 -
コード長 851 bytes
コンパイル時間 166 ms
コンパイル使用メモリ 82,388 KB
実行使用メモリ 67,520 KB
最終ジャッジ日時 2024-04-09 21:44:34
合計ジャッジ時間 4,296 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 RE -
testcase_01 RE -
testcase_02 RE -
testcase_03 RE -
testcase_04 RE -
testcase_05 RE -
testcase_06 RE -
testcase_07 RE -
testcase_08 RE -
testcase_09 RE -
testcase_10 RE -
testcase_11 RE -
testcase_12 RE -
testcase_13 RE -
testcase_14 RE -
testcase_15 RE -
testcase_16 RE -
testcase_17 RE -
testcase_18 RE -
testcase_19 RE -
testcase_20 RE -
testcase_21 RE -
testcase_22 RE -
testcase_23 RE -
testcase_24 RE -
testcase_25 RE -
testcase_26 RE -
testcase_27 RE -
testcase_28 RE -
testcase_29 RE -
testcase_30 RE -
testcase_31 RE -
testcase_32 RE -
testcase_33 RE -
testcase_34 RE -
testcase_35 RE -
testcase_36 RE -
testcase_37 RE -
testcase_38 RE -
testcase_39 RE -
testcase_40 RE -
testcase_41 RE -
testcase_42 RE -
testcase_43 RE -
testcase_44 RE -
testcase_45 RE -
testcase_46 RE -
testcase_47 RE -
testcase_48 RE -
testcase_49 RE -
testcase_50 RE -
testcase_51 RE -
testcase_52 RE -
権限があれば一括ダウンロードができます

ソースコード

diff # 

import numpy as np
N=int(input())
A=[int(x) for x in input().split()]
s=[0]*(1<<N)
for i in range(1<<N):
    for j in range(N):
        if (i&(1<<j))>0:
            s[i]+=A[j]
dp={}
def f(a,b,t):
    if (a|b)==(1<<N)-1:
        return 1 if 2*s[a]+t>2*s[b] else 0
    c=int((a|(b<<N))*t)
    if c in dp:
        return dp[c]
    ret=0
    for i in range(N):
        na=a|(1<<i)
        if (na|b)==(a|b):
            continue
        e=1
        for j in range(N):
            nb=b|(1<<j)
            if (a|nb)==(a|b):
                continue
            p,q=A[i],A[j]
            vec=np.array([[-f(b,na,t*(-1))/p+1],[-f(a,nb,t)/q+1]])
            mat=np.array([[1,-(1-1/q)],[-(1-1/p),1]])
            det=(p+q-1)/(p*q)
            mat/=det
            e=min(e,np.dot(mat,vec)[0][0])
        ret=max(ret,e)
    dp[c]=ret
    return ret
print(f(0,0,-1))
0