結果

問題 No.132 点と平面との距離
ユーザー mkawa2mkawa2
提出日時 2020-01-21 14:48:18
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
WA  
実行時間 -
コード長 1,155 bytes
コンパイル時間 700 ms
コンパイル使用メモリ 10,876 KB
実行使用メモリ 8,412 KB
最終ジャッジ日時 2023-09-19 05:07:01
合計ジャッジ時間 5,030 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 WA -
testcase_01 WA -
testcase_02 WA -
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys

sys.setrecursionlimit(10 ** 6)
def II(): return int(sys.stdin.readline())
def LF(): return list(map(float, sys.stdin.readline().split()))

def main():
    # 3点Q1,Q2,Q3からできる三角形を底面とし、Pを頂点とする三角すいを考える
    # 底面積(△Q1Q2Q3)と体積が分かれば高さ(点Pと平面Q1Q2Q3の距離)が分かる
    n = II()
    xp, yp, zp = LF()
    qq = [LF() for _ in range(n)]
    qq = [[x - xp, y - yp, z - zp] for x, y, z in qq]
    ans = 0
    for k, (x3, y3, z3) in enumerate(qq):
        for j, (x2, y2, z2) in enumerate(qq[:k]):
            for (x1, y1, z1) in qq[:j]:
                # 体積(の6倍)vを求める
                v = abs(x1 * y2 * z3 + x2 * y3 * z1 + x3 * y1 * z2 - x3 * y2 * z1 - x2 * y1 * z3 - x1 * y3 * z2)
                # 底面積(の2倍)sを求める
                x12, y12, z12 = x2 - x1, y2 - y1, z2 - z1
                x13, y13, z13 = x3 - x1, y3 - y1, z3 - z1
                s = abs(x12 * y13 + y12 * z13 + z12 * x13 - y12 * x13 - z12 * y13 - x12 * z13)
                # h=v/sより高さhを求める
                ans += v / s
    print(ans)

main()
0