結果

問題 No.2010 Magical Floor
ユーザー rin204
提出日時 2022-07-15 22:30:51
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 1,098 bytes
コンパイル時間 350 ms
コンパイル使用メモリ 82,432 KB
実行使用メモリ 79,172 KB
最終ジャッジ日時 2024-06-27 19:14:42
合計ジャッジ時間 7,811 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 2 WA * 19
権限があれば一括ダウンロードができます

ソースコード

diff # 

def solve():
    n, x = map(int, input().split())
    H = list(map(int, input().split()))
    A = list(map(int, input().split()))
    ans = x * A[0]
    z = x / 2
    h = H[0]
    a = (h * h + 1) ** 0.5 * A[0]
    for h2, b in zip(H[1:], A[1:]):
        f = lambda x:((h * h + x * x) / (h * h + 1)) ** 0.5 * a + (z - x) * b
        l = 0
        r = z
        for _ in range(100):
            llr = (2 * l + r) / 3
            lrr = (l + 2 * r) / 3
            if f(llr) < f(lrr):
                r = lrr
            else:
                l = llr
        tmp = 2 * min(f(l), f(llr), f(lrr), f(r))
        ans = min(ans, tmp)
        f = lambda x:((h * h + x * x) / (h * h + 1)) ** 0.5 * a + (h2 * h2 + (1 - x) ** 2) ** 0.5 * b
        l = 0
        r = 1
        for _ in range(100):
            llr = (2 * l + r) / 3
            lrr = (l + 2 * r) / 3
            if f(llr) < f(lrr):
                r = lrr
            else:
                l = llr
        tmp = min(f(l), f(llr), f(lrr), f(r))
        a = tmp
        h += h2


    print(ans)

        


for _ in range(int(input())):
    solve()
0