結果

問題 No.622 点と三角柱の内外判定
ユーザー kotatsugamekotatsugame
提出日時 2019-10-12 04:50:56
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 2 ms / 1,500 ms
コード長 955 bytes
コンパイル時間 642 ms
コンパイル使用メモリ 69,760 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-05-05 04:19:22
合計ジャッジ時間 1,453 ms
ジャッジサーバーID
(参考情報)
judge2 / judge5
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
5,248 KB
testcase_01 AC 2 ms
5,376 KB
testcase_02 AC 2 ms
5,376 KB
testcase_03 AC 2 ms
5,376 KB
testcase_04 AC 2 ms
5,376 KB
testcase_05 AC 1 ms
5,376 KB
testcase_06 AC 2 ms
5,376 KB
testcase_07 AC 1 ms
5,376 KB
testcase_08 AC 1 ms
5,376 KB
testcase_09 AC 2 ms
5,376 KB
testcase_10 AC 1 ms
5,376 KB
testcase_11 AC 1 ms
5,376 KB
testcase_12 AC 2 ms
5,376 KB
testcase_13 AC 1 ms
5,376 KB
testcase_14 AC 2 ms
5,376 KB
testcase_15 AC 2 ms
5,376 KB
testcase_16 AC 1 ms
5,376 KB
testcase_17 AC 2 ms
5,376 KB
testcase_18 AC 1 ms
5,376 KB
testcase_19 AC 1 ms
5,376 KB
testcase_20 AC 2 ms
5,376 KB
testcase_21 AC 2 ms
5,376 KB
testcase_22 AC 2 ms
5,376 KB
testcase_23 AC 2 ms
5,376 KB
testcase_24 AC 1 ms
5,376 KB
testcase_25 AC 1 ms
5,376 KB
testcase_26 AC 2 ms
5,376 KB
testcase_27 AC 2 ms
5,376 KB
testcase_28 AC 1 ms
5,376 KB
testcase_29 AC 1 ms
5,376 KB
testcase_30 AC 1 ms
5,376 KB
testcase_31 AC 1 ms
5,376 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp:18:1: warning: ISO C++ forbids declaration of 'main' with no type [-Wreturn-type]
   18 | main()
      | ^~~~

ソースコード

diff #

#include<iostream>
#include<cmath>
using namespace std;
const double EPS=1e-3;
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;
	}
}
0