結果

問題 No.622 点と三角柱の内外判定
ユーザー wajima_wataruwajima_wataru
提出日時 2017-12-29 18:27:27
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
AC  
実行時間 36 ms / 1,500 ms
コード長 1,558 bytes
コンパイル時間 91 ms
コンパイル使用メモリ 12,928 KB
実行使用メモリ 11,392 KB
最終ジャッジ日時 2024-06-01 06:40:57
合計ジャッジ時間 2,325 ms
ジャッジサーバーID
(参考情報)
judge1 / judge4
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 32 ms
11,392 KB
testcase_01 AC 32 ms
11,392 KB
testcase_02 AC 34 ms
11,264 KB
testcase_03 AC 36 ms
11,264 KB
testcase_04 AC 35 ms
11,264 KB
testcase_05 AC 32 ms
11,264 KB
testcase_06 AC 33 ms
11,264 KB
testcase_07 AC 32 ms
11,264 KB
testcase_08 AC 32 ms
11,264 KB
testcase_09 AC 32 ms
11,392 KB
testcase_10 AC 33 ms
11,264 KB
testcase_11 AC 32 ms
11,264 KB
testcase_12 AC 33 ms
11,392 KB
testcase_13 AC 32 ms
11,264 KB
testcase_14 AC 32 ms
11,264 KB
testcase_15 AC 33 ms
11,264 KB
testcase_16 AC 33 ms
11,136 KB
testcase_17 AC 32 ms
11,264 KB
testcase_18 AC 33 ms
11,264 KB
testcase_19 AC 32 ms
11,264 KB
testcase_20 AC 33 ms
11,264 KB
testcase_21 AC 33 ms
11,392 KB
testcase_22 AC 33 ms
11,264 KB
testcase_23 AC 33 ms
11,392 KB
testcase_24 AC 33 ms
11,264 KB
testcase_25 AC 32 ms
11,264 KB
testcase_26 AC 32 ms
11,264 KB
testcase_27 AC 32 ms
11,264 KB
testcase_28 AC 32 ms
11,264 KB
testcase_29 AC 32 ms
11,264 KB
testcase_30 AC 32 ms
11,392 KB
testcase_31 AC 33 ms
11,264 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

from math import sin, cos, atan2
from operator import *
pos = []
for _ in range(4):
	pos.append(list(map(float, input().split())))


def rotate_x(p, theta):
	q = []
	q.append(p[0])
	q.append(p[1] * cos(theta) - p[2] * sin(theta))
	q.append(p[1] * sin(theta) + p[2] * cos(theta))
	return q

def rotate_y(p, theta):
	q = []
	q.append(p[2] * sin(theta) + p[0] * cos(theta))
	q.append(p[1])
	q.append(p[2] * cos(theta) - p[0] * sin(theta))
	return q

def rotate_z(p, theta):
	q = []
	q.append(p[0] * cos(theta) - p[1] * sin(theta))
	q.append(p[0] * sin(theta) + p[1] * cos(theta))
	q.append(p[2])
	return q

def translate(p, q):
	for i in range(3):
		p[i] += q[i]
	return p


d = list(map(neg, pos[0]))
for i in range(4):
	pos[i] = translate(pos[i], d)

th1 = atan2(pos[1][2], pos[1][0])
for i in range(4):
	pos[i] = rotate_y(pos[i], th1)

th2 = atan2(-pos[1][1], pos[1][0])
for i in range(4):
	pos[i] = rotate_z(pos[i], th2)

th3 = atan2(-pos[2][2], pos[2][1])
for i in range(4):
	pos[i] = rotate_x(pos[i], th3)

AD = [pos[3][0] - pos[0][0], pos[3][1] - pos[0][1]]
BD = [pos[3][0] - pos[1][0], pos[3][1] - pos[1][1]]
CD = [pos[3][0] - pos[2][0], pos[3][1] - pos[2][1]]
AB = [pos[1][0] - pos[0][0], pos[1][1] - pos[0][1]]
BC = [pos[2][0] - pos[1][0], pos[2][1] - pos[1][1]]
CA = [pos[0][0] - pos[2][0], pos[0][1] - pos[2][1]]


def cross(a, b):
	return a[0] * b[1] - b[0] * a[1]


if cross(AB, BD) > 0 and cross(BC, CD) > 0 and cross(CA, AD) > 0:
	print('YES')
elif cross(AB, BD) < 0 and cross(BC, CD) < 0 and cross(CA, AD) < 0:
	print('YES')
else:
	print('NO')
0