結果

問題 No.2632 Center of Three Points in Lp Norm
ユーザー tatyamtatyam
提出日時 2024-02-12 07:49:52
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 876 bytes
コンパイル時間 141 ms
コンパイル使用メモリ 82,496 KB
実行使用メモリ 80,468 KB
最終ジャッジ日時 2024-09-28 17:52:56
合計ジャッジ時間 10,244 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 172 ms
78,756 KB
testcase_01 AC 187 ms
77,716 KB
testcase_02 AC 150 ms
77,464 KB
testcase_03 WA -
testcase_04 AC 297 ms
79,700 KB
testcase_05 WA -
testcase_06 AC 244 ms
79,132 KB
testcase_07 AC 154 ms
78,028 KB
testcase_08 AC 177 ms
78,356 KB
testcase_09 AC 240 ms
78,056 KB
testcase_10 AC 288 ms
80,276 KB
testcase_11 AC 269 ms
78,616 KB
testcase_12 AC 135 ms
78,396 KB
testcase_13 AC 261 ms
80,268 KB
testcase_14 AC 146 ms
78,436 KB
testcase_15 WA -
testcase_16 WA -
testcase_17 AC 186 ms
79,028 KB
testcase_18 WA -
testcase_19 WA -
testcase_20 AC 151 ms
78,504 KB
testcase_21 AC 145 ms
78,180 KB
testcase_22 WA -
testcase_23 WA -
testcase_24 AC 186 ms
78,824 KB
testcase_25 AC 197 ms
79,400 KB
testcase_26 AC 144 ms
78,152 KB
testcase_27 AC 137 ms
77,408 KB
testcase_28 WA -
testcase_29 AC 163 ms
77,648 KB
testcase_30 AC 144 ms
77,372 KB
testcase_31 AC 143 ms
77,480 KB
testcase_32 WA -
testcase_33 WA -
testcase_34 AC 144 ms
79,044 KB
testcase_35 WA -
testcase_36 AC 179 ms
79,876 KB
testcase_37 AC 234 ms
80,368 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import random
import math

def norm(x: float, y: float) -> float:
    x = abs(x)
    y = abs(y)
    if x > y:
        x, y = y, x
    if y == 0.0:
        return 0.0
    x /= y
    return pow(1.0 + pow(x, p), 1 / p) * y

def dist(a, b):
    return norm(a[0] - b[0], a[1] - b[1])

get = lambda: list(map(float, input().split()))
p = float(input())
a = get()
b = get()
c = get()

def f(o) -> float:
    D = [dist(a, o), dist(b, o), dist(c, o)]
    D.sort()
    return D[2] - x * D[0]

o = [0.0, 0.0]
d = 1.0
R = 1e-15
T = int(1e5)
r = R ** (10 / T)

for _ in range(T):
    x = _ / T
    theta = (random.random() * 2.0 - 1.0) * math.pi
    while True:
        o2 = [o[0] + d * math.cos(theta), o[1] + d * math.sin(theta)]
        if f(o2) >= f(o):
            break
        o = o2
        d /= r
    d *= r

print(*o)
print(dist(a, o), dist(b, o), dist(c, o), file=open(2, "w"))
0