結果
問題 | 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, setcontext Fp = float getcontext().prec = 12 p = 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) ** 2 def bsearch_y(x: Fp): ly, ry = Fp(-1000000), Fp(1000000) for _ in range(100): mly = ly + (ry - ly) / 3 mry = ry - (ry - ly) / 3 if max_norm(x, mly) > max_norm(x, mry): ly = mly else: ry = mry return max_norm(x, ly), ly def bsearch_x(): lx, rx = Fp(-1000000), Fp(1000000) for _ in range(100): mlx = lx + (rx - lx) / 3 mrx = rx - (rx - lx) / 3 if bsearch_y(mlx)[0] > bsearch_y(mrx)[0]: lx = mlx else: rx = mrx return 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}')