結果

問題 No.622 点と三角柱の内外判定
ユーザー Ryuto
提出日時 2017-12-22 22:07:13
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 1,430 bytes
コンパイル時間 142 ms
コンパイル使用メモリ 82,472 KB
実行使用メモリ 54,844 KB
最終ジャッジ日時 2024-12-17 23:46:23
合計ジャッジ時間 2,286 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 23 WA * 9
権限があれば一括ダウンロードができます

ソースコード

diff # 

#!/usr/bin/env python


def cross(A, B):
    X = A[1]*B[2] - A[2]*B[1]
    Y = A[2]*B[0] - A[0]*B[2]
    Z = A[0]*B[1] - A[1]*B[0]
    return [X, Y, Z]


def inverse(A):
    det = A[0][0]*A[1][1]*A[2][2] + A[0][1]*A[1][2]*A[2][0]
    det += A[0][2]*A[1][0]*A[2][1] - A[0][0]*A[1][2]*A[2][1]
    det += - A[0][2]*A[1][1]*A[2][0] - A[0][1]*A[1][0]*A[2][2]
    A11 = (-1)**2 * (A[1][1]*A[2][2] - A[1][2]*A[2][1]) / det
    A12 = (-1)**3 * (A[1][0]*A[2][2] - A[1][2]*A[2][0]) / det
    A13 = (-1)**4 * (A[1][0]*A[2][1] - A[1][1]*A[2][0]) / det
    A21 = (-1)**3 * (A[0][1]*A[2][2] - A[0][2]*A[2][1]) / det
    A22 = (-1)**4 * (A[0][0]*A[2][2] - A[0][2]*A[2][0]) / det
    A23 = (-1)**5 * (A[0][0]*A[2][1] - A[0][1]*A[2][0]) / det
    A31 = (-1)**4 * (A[0][1]*A[1][2] - A[0][2]*A[1][1]) / det
    A32 = (-1)**5 * (A[0][0]*A[1][2] - A[0][2]*A[1][0]) / det
    A33 = (-1)**6 * (A[0][0]*A[1][1] - A[0][1]*A[1][0]) / det
    return [[A11, A21, A31], [A12, A22, A32], [A13, A23, A33]]


A = [int(x) for x in input().split()]
B = [int(x) for x in input().split()]
C = [int(x) for x in input().split()]
D = [int(x) for x in input().split()]

AB = [b - a for (a, b) in zip(A, B)]
AC = [c - a for (a, c) in zip(A, C)]
AD = [d - a for (a, d) in zip(A, D)]

Z = cross(AB, AC)
M = inverse([AB, AC, Z])

Co = [AD[0]*x + AD[1]*y + AD[2]*z for (x, y, z) in zip(M[0], M[1], M[2])]

if Co[0] <= 1 and Co[1] <= 1:
    print('YES')
else:
    print('NO')
0