結果
問題 | No.622 点と三角柱の内外判定 |
ユーザー |
|
提出日時 | 2017-07-03 02:34:13 |
言語 | PyPy2 (7.3.15) |
結果 |
AC
|
実行時間 | 87 ms / 1,500 ms |
コード長 | 1,466 bytes |
コンパイル時間 | 2,456 ms |
コンパイル使用メモリ | 76,800 KB |
実行使用メモリ | 76,928 KB |
最終ジャッジ日時 | 2024-11-27 11:39:48 |
合計ジャッジ時間 | 6,243 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 32 |
ソースコード
#!/usr/bin/python#works with Py2/3 and use numpy only if available.import copy,mathtry:import numpydef mul(a,b):return numpy.array(a).dot(b)except ImportError:def mul(a,b):assert len(a[0])==len(b)r=[[0]*len(b[0]) for _ in a]for i in range(len(a)):for j in range(len(b[0])):for k in range(len(a[0])):r[i][j]+=a[i][k]*b[k][j]return r#def dist(a,b):# return math.sqrt(sum((a[i]-b[i])**2 for i in range(len(a))))def rotateTriangle(*_):d=copy.deepcopy(list(_))for i in range(len(d)-1,-1,-1):d[i][0]-=d[0][0]d[i][1]-=d[0][1]d[i][2]-=d[0][2]d=list(zip(*d))thetay=-math.atan2(d[2][1],d[0][1])d=mul([[math.cos(thetay),0,-math.sin(thetay)],[0,1,0],[math.sin(thetay),0,math.cos(thetay)]],d)thetaz=-math.atan2(d[1][1],d[0][1])d=mul([[math.cos(thetaz),-math.sin(thetaz),0],[math.sin(thetaz),math.cos(thetaz),0],[0,0,1]],d)thetax=-math.atan2(d[2][2],d[1][2])d=mul([[1,0,0],[0,math.cos(thetax),-math.sin(thetax)],[0,math.sin(thetax),math.cos(thetax)]],d)d=list(zip(*d))return [e[0:2] for e in d]import sysa=[list(map(float,sys.stdin.readline().split())) for _ in range(4)]b=rotateTriangle(*a)for i in range(3):edgex=b[(i+1)%3][0]-b[i][0]edgey=b[(i+1)%3][1]-b[i][1]vx=b[3][0]-b[i][0]vy=b[3][1]-b[i][1]if edgex*vy-edgey*vx<0:print('NO')exit()print('YES')