結果
問題 |
No.622 点と三角柱の内外判定
|
ユーザー |
|
提出日時 | 2019-10-12 04:56:26 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 2 ms / 1,500 ms |
コード長 | 973 bytes |
コンパイル時間 | 713 ms |
コンパイル使用メモリ | 68,992 KB |
実行使用メモリ | 5,248 KB |
最終ジャッジ日時 | 2024-11-26 22:12:10 |
合計ジャッジ時間 | 1,503 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 32 |
コンパイルメッセージ
main.cpp:19:1: warning: ISO C++ forbids declaration of 'main' with no type [-Wreturn-type] 19 | main() | ^~~~
ソースコード
#include<iostream> #include<cmath> using namespace std; #define double long double const double EPS=1e-6; struct Point{ double x,y,z; Point(double x_=0,double y_=0,double z_=0):x(x_),y(y_),z(z_){} void input(){cin>>x>>y>>z;} Point operator+(Point a){return Point{x+a.x,y+a.y,z+a.z};} Point operator-(Point a){return Point{x-a.x,y-a.y,z-a.z};} Point operator*(double a){return Point{x*a,y*a,z*a};} }; double dot(Point a,Point b){return a.x*b.x+a.y*b.y+a.z*b.z;} Point cross(Point a,Point b){return Point{a.y*b.z-a.z*b.y,a.z*b.x-a.x*b.z,a.x*b.y-a.y*b.x};} double norm(Point a){return dot(a,a);} double abs(Point a){return sqrt(dot(a,a));} double area(Point a,Point b,Point c){return abs(cross(b-a,c-a))/2;} main() { Point p,a,b,c; a.input(); b.input(); c.input(); p.input(); Point h=cross(b-a,c-a); double k=dot(h,a-p)/norm(h); Point x=p+h*k; if(area(x,a,b)+area(x,b,c)+area(x,c,a)>area(a,b,c)+EPS) { cout<<"NO"<<endl; } else { cout<<"YES"<<endl; } }