結果
問題 | No.132 点と平面との距離 |
ユーザー | mkawa2 |
提出日時 | 2020-01-21 14:14:37 |
言語 | Python3 (3.12.2 + numpy 1.26.4 + scipy 1.12.0) |
結果 |
TLE
|
実行時間 | - |
コード長 | 902 bytes |
コンパイル時間 | 138 ms |
コンパイル使用メモリ | 12,416 KB |
実行使用メモリ | 48,800 KB |
最終ジャッジ日時 | 2024-07-05 16:41:17 |
合計ジャッジ時間 | 12,566 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | TLE | - |
testcase_01 | -- | - |
testcase_02 | -- | - |
ソースコード
import sys sys.setrecursionlimit(10 ** 6) def II(): return int(sys.stdin.readline()) def LF(): return list(map(float, sys.stdin.readline().split())) import numpy as np def main(): # 3点Q1,Q2,Q3からできる三角形を底面とし、Pを頂点とする三角すいを考える # 底面積(△Q1Q2Q3)と体積が分かれば高さ(点Pと平面Q1Q2Q3の距離)が分かる n = II() p = np.array(LF()) qq = [np.array(LF()) for _ in range(n)] ans = 0 for k, q3 in enumerate(qq): for j, q2 in enumerate(qq[:k]): for q1 in qq[:j]: # 体積vを求める v = abs(np.dot(q1 - p, np.cross(q2 - p, q3 - p)) / 6) # 底面積sを求める s = np.linalg.norm(np.cross(q2 - q1, q3 - q1)) / 2 # h=3v/sより高さhを求める ans += 3 * v / s print(ans) main()