結果

問題 No.132 点と平面との距離
ユーザー masa
提出日時 2015-01-21 01:08:34
言語 C++11
(gcc 13.3.0)
結果
AC  
実行時間 47 ms / 5,000 ms
コード長 1,333 bytes
コンパイル時間 683 ms
コンパイル使用メモリ 67,308 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-06-22 23:30:21
合計ジャッジ時間 1,143 ms
ジャッジサーバーID
(参考情報)
judge4 / judge1
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
other AC * 3
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#include <iostream>
#include <cstdio>
#include <vector>
#include <algorithm>
#include <utility>
#include <cmath>
using namespace std;
typedef struct {
double x;
double y;
double z;
} Point;
Point subtract(Point, Point);;
Point cross_product(Point, Point);
double calc(Point, Point, Point);
Point p;
int main() {
int n;
cin >> n;
int x, y, z;
vector<Point> q(n);
cin >> p.x >> p.y >> p.z;
for (int i = 0; i < n; i++) {
cin >> q[i].x >> q[i].y >> q[i].z;
}
Point ab, ac;
double ans = 0.0;
for (int i = 0; i < n; i++) {
for (int j = i + 1; j < n; j++) {
for (int k = j + 1; k < n; k++) {
ans += calc(q[i], q[j], q[k]);
}
}
}
printf("%.10f\n", ans);
return 0;
}
Point subtract(Point a, Point b) {
Point c;
c.x = a.x - b.x;
c.y = a.y - b.y;
c.z = a.z - b.z;
return c;
}
Point cross_product(Point a, Point b) {
Point c;
c.x = a.y * b.z - a.z * b.y;
c.y = a.z * b.x - a.x * b.z;
c.z = a.x * b.y - a.y * b.x;
return c;
}
double calc(Point p1, Point p2, Point p3) {
Point p12, p13, cross;
p12 = subtract(p2, p1);
p13 = subtract(p3, p1);
cross = cross_product(p12, p13);
double d;
d = -1.0 *(cross.x * p1.x + cross.y * p1.y + cross.z * p1.z);
return abs(cross.x * p.x + cross.y * p.y + cross.z * p.z + d) / sqrt(cross.x * cross.x + cross.y * cross.y + cross.z * cross.z);
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0