結果
問題 | No.132 点と平面との距離 |
ユーザー |
![]() |
提出日時 | 2016-05-05 16:09:08 |
言語 | C++11 (gcc 13.3.0) |
結果 |
CE
(最新)
AC
(最初)
|
実行時間 | - |
コード長 | 1,385 bytes |
コンパイル時間 | 492 ms |
コンパイル使用メモリ | 61,136 KB |
最終ジャッジ日時 | 2024-11-14 19:43:10 |
合計ジャッジ時間 | 969 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
コンパイルエラー時のメッセージ・ソースコードは、提出者また管理者しか表示できないようにしております。(リジャッジ後のコンパイルエラーは公開されます)
ただし、clay言語の場合は開発者のデバッグのため、公開されます。
ただし、clay言語の場合は開発者のデバッグのため、公開されます。
コンパイルメッセージ
main.cpp: In function ‘double calc_area(double*, double*)’: main.cpp:29:12: error: ‘sqrt’ was not declared in this scope 29 | return sqrt(ret); | ^~~~
ソースコード
#include <iostream>#include <vector>#include <algorithm>#include <iomanip>using namespace std;#define RREP(i,s,e) for (i = s; i >= e; i--)#define rrep(i,n) RREP(i,(int)(n)-1,0)#define REP(i,s,e) for (i = s; i <= e; i++)#define rep(i,n) REP(i,0,(int)(n)-1)#define INF 100000000typedef long long ll;double det(double* p, double* q, double* r) {int i;double ret = 0;rep (i,3) ret += p[i]*q[(i+1)%3]*r[(i+2)%3];rep (i,3) ret -= p[i]*q[(i+2)%3]*r[(i+1)%3];return ret;}double calc_area(double* p, double* q) {int i;double ret = 0;ret += (p[1]*q[2]-p[2]*q[1]) * (p[1]*q[2]-p[2]*q[1]);ret += (p[2]*q[0]-p[0]*q[2]) * (p[2]*q[0]-p[0]*q[2]);ret += (p[0]*q[1]-p[1]*q[0]) * (p[0]*q[1]-p[1]*q[0]);return sqrt(ret);}int main() {int i, j, k, l, n;double p[3], a[3], b[3];double xyz[300][3];double ans = 0;cin >> n;rep (i,3) cin >> p[i];rep (i,n) {cin >> xyz[i][0] >> xyz[i][1] >> xyz[i][2];rep (j,3) xyz[i][j] -= p[j];}rep (i,n) REP (j,i+1,n-1) REP (k,j+1,n-1) {double vol = abs(det(xyz[i],xyz[j],xyz[k])/6);rep (l,3) a[l] = xyz[j][l] - xyz[i][l];rep (l,3) b[l] = xyz[k][l] - xyz[i][l];double area = calc_area(a,b) / 2;ans += vol * 3 / area;}cout << fixed << setprecision(10) << ans << endl;return 0;}