結果
問題 | No.2632 Center of Three Points in Lp Norm |
ユーザー | suisen |
提出日時 | 2024-02-16 23:44:18 |
言語 | PyPy3 (7.3.15) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,096 bytes |
コンパイル時間 | 375 ms |
コンパイル使用メモリ | 82,524 KB |
実行使用メモリ | 76,152 KB |
最終ジャッジ日時 | 2024-09-28 22:46:54 |
合計ジャッジ時間 | 5,293 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 22 WA * 12 |
ソースコード
Fp = float p = Fp(input()) ax, ay = map(Fp, input().split()) bx, by = map(Fp, input().split()) cx, cy = map(Fp, input().split()) t = Fp(2) / (Fp(3) + Fp(5) ** (Fp(1) / Fp(2))) 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 + t * (ry - ly) mry = ry - t * (ry - ly) 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 + t * (rx - lx) mrx = rx - t * (rx - lx) 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(f'{x:.20f} {y:.20f}')