結果

問題 No.132 点と平面との距離
コンテスト
ユーザー IL_msta
提出日時 2015-07-31 15:58:09
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 1,029 ms / 5,000 ms
コード長 1,359 bytes
コンパイル時間 851 ms
コンパイル使用メモリ 92,116 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-07-17 22:55:34
合計ジャッジ時間 2,540 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 3
権限があれば一括ダウンロードができます

ソースコード

diff # 

#define _USE_MATH_DEFINES
 
#include <iostream>
#include <iomanip>
#include <sstream>
 
#include <algorithm>
#include <cmath>
 
#include <string>
//#include <array>
#include <list>
#include <queue>
#include <vector>
#include <complex>
#include <set>
#include <map>
 
/////////
#define REP(i, x, n) for(int i = x; i < n; i++)
#define rep(i,n) REP(i,0,n)
#define P(p) cout<<(p)<<endl;
 
#define PII pair<int,int>
/////////
typedef long long LL;
typedef long double LD;
/////////
using namespace::std;
/////////

LD Len(vector<LD> p,vector<LD> q1,vector<LD> q2,vector<LD> q3){
	LD ans=0;
	vector<LD> L0(3),E1(3),E2(3),G(3);
	LD e3;
	rep(i,3){
		L0[i] = p[i] -q1[i];
		E1[i] = q2[i]-q1[i];
		E2[i] = q3[i]-q1[i];
	}
	
	rep(i,3){
		G[i] = E1[(i+1)%3]*E2[(i+2)%3]-E1[(3+i-1)%3]*E2[(3+i-2)%3];
	}
	e3 = sqrt(G[0]*G[0]+G[1]*G[1]+G[2]*G[2]);
	rep(i,3){
		ans += L0[i]*G[i];
	}
	ans /= e3;
	return ans>=0?ans:-ans;
}

int main(void){
    std::cin.tie(0); 
    std::ios::sync_with_stdio(false);
    std::cout << std::fixed;//
    cout << setprecision(16);//
	
	int N;
	cin>>N;
	vector<LD> p(3);
	vector< vector<LD>> q(300,vector<LD>(3));
	cin>>p[0]>>p[1]>>p[2];
	rep(i,N){
		cin>>q[i][0]>>q[i][1]>>q[i][2];
	}
	LD ans = 0;
	for(int i=0;i<N-2;++i){
		for(int j=i+1;j<N-1;++j){
			for(int k=j+1;k<N;++k){
				ans += Len(p,q[i],q[j],q[k]);
			}
		}
	}
	P(ans);
	return 0;
}
0