結果

問題 No.132 点と平面との距離
ユーザー Lay_ecLay_ec
提出日時 2015-01-20 23:54:07
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 48 ms / 5,000 ms
コード長 1,098 bytes
コンパイル時間 795 ms
コンパイル使用メモリ 83,272 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-06-22 23:24:51
合計ジャッジ時間 1,269 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 3 ms
5,248 KB
testcase_01 AC 16 ms
5,376 KB
testcase_02 AC 48 ms
5,376 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <string>
#include <vector>
#include <cmath>
#include <algorithm>
#include <cstdlib>
#include <ctime>
#include <cstdio>
#include <functional>
#include <set>
#include <sstream>
#include <map>
#include <queue>

using namespace std;
 
const double eps=1e-10;

vector<double> X,Y,Z;
double px,py,pz;


double dist(int i,int j,int k)
{
	double vijx,vijy,vijz;
	double vikx,viky,vikz;

	vijx=X[j]-X[i];
	vijy=Y[j]-Y[i];
	vijz=Z[j]-Z[i];

	vikx=X[k]-X[i];
	viky=Y[k]-Y[i];
	vikz=Z[k]-Z[i];

	//a(x-X[i])+b(y-Y[i])+c(z-Z[i])=0
	//ax+by+cz+(-aX[i]-bY[i]-cZ[i])=0

	double a=vijy*vikz - vijz*viky;
	double b=vijz*vikx - vijx*vikz;
	double c=vijx*viky - vijy*vikx;
	double d=-a*X[i]-b*Y[i]-c*Z[i];

	double res=abs(a*px+b*py+c*pz+d)/sqrt(a*a+b*b+c*c);

	return res;

}

int main()
{

	int n;
	cin>>n;

	cin>>px>>py>>pz;

	double res=0.0;

	X.resize(n);
	Y.resize(n);
	Z.resize(n);
	for(int i=0;i<n;i++){
		cin>>X[i]>>Y[i]>>Z[i];
	}

	for(int i=0;i<n;i++){
		for(int j=i+1;j<n;j++){
			for(int k=j+1;k<n;k++){
				res+=dist(i,j,k);
			}
		}
	}

	printf("%.12f\n",res);

	return 0;
}
0