結果
問題 | 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 = floatp = 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) ** 2def 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 = mlyelse:ry = mryreturn max_norm(x, ly), lydef 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 = mlxelse:rx = mrxreturn lx, bsearch_y(lx)[1]x, y = bsearch_x()print(f'{x:.20f} {y:.20f}')