結果
| 問題 |
No.622 点と三角柱の内外判定
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2019-10-12 04:50:41 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 955 bytes |
| コンパイル時間 | 668 ms |
| コンパイル使用メモリ | 69,504 KB |
| 実行使用メモリ | 5,248 KB |
| 最終ジャッジ日時 | 2024-11-26 22:02:47 |
| 合計ジャッジ時間 | 1,740 ms |
|
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | AC * 31 WA * 1 |
コンパイルメッセージ
main.cpp:18:1: warning: ISO C++ forbids declaration of 'main' with no type [-Wreturn-type]
18 | main()
| ^~~~
ソースコード
#include<iostream>
#include<cmath>
using namespace std;
const double EPS=1e-4;
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);
h=h*(1/abs(h));
double k=dot(h,a-p);
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;
}
}