結果
問題 | No.2632 Center of Three Points in Lp Norm |
ユーザー | suisen |
提出日時 | 2024-02-16 23:37:39 |
言語 | PyPy3 (7.3.15) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,208 bytes |
コンパイル時間 | 149 ms |
コンパイル使用メモリ | 82,388 KB |
実行使用メモリ | 83,360 KB |
最終ジャッジ日時 | 2024-09-28 22:39:01 |
合計ジャッジ時間 | 8,439 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 21 WA * 13 |
ソースコード
from decimal import Decimal, getcontext, setcontextFp = floatgetcontext().prec = 12p = Fp(input())ax, ay = map(Fp, input().split())bx, by = map(Fp, input().split())cx, cy = map(Fp, input().split())def norm(x1: Fp, y1: Fp, x2: Fp, y2: Fp):return (abs(x1 - x2) ** p + abs(y1 - y2) ** p) ** (Fp(1) / p)def max_norm(x: Fp, y: Fp):a = norm(x, y, ax, ay)b = norm(x, y, bx, by)c = norm(x, y, cx, cy)return abs(a - b) ** 2 + abs(b - c) ** 2 + abs(c - a) ** 2def bsearch_y(x: Fp):ly, ry = Fp(-1000000), Fp(1000000)for _ in range(100):mly = ly + (ry - ly) / 3mry = ry - (ry - ly) / 3if max_norm(x, mly) > max_norm(x, mry):ly = mlyelse:ry = mryreturn max_norm(x, ly), lydef bsearch_x():lx, rx = Fp(-1000000), Fp(1000000)for _ in range(100):mlx = lx + (rx - lx) / 3mrx = rx - (rx - lx) / 3if bsearch_y(mlx)[0] > bsearch_y(mrx)[0]:lx = mlxelse:rx = mrxreturn lx, bsearch_y(lx)[1]x, y = bsearch_x()# print(norm(x, y, ax, ay))# print(norm(x, y, bx, by))# print(norm(x, y, cx, cy))print(f'{x:.10f} {y:.10f}')