結果

問題 No.132 点と平面との距離
コンテスト
ユーザー 🍮かんプリン
提出日時 2019-09-13 19:20:09
言語 C++14
(gcc 13.3.0 + boost 1.89.0)
結果
AC  
実行時間 41 ms / 5,000 ms
コード長 1,219 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 1,423 ms
コンパイル使用メモリ 166,712 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-07-04 03:53:26
合計ジャッジ時間 1,916 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 3
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include "bits/stdc++.h"
#define ALL(obj) (obj).begin(),(obj).end()
#define RALL(obj) (obj).rbegin(),(obj).rend()
#define REP(i, n) for(int i = 0; i < (int)(n); i++)
#define REPR(i, n) for(int i = (int)(n); i >= 0; i--)
#define FOR(i,n,m) for(int i = (int)(n); i < int(m); i++)
using namespace std;
typedef long long ll;
const int MOD = 1e9 + 7;
const int INF = MOD - 1;
const ll LLINF = 4e18;

double dist(double p[], double a[], double b[], double c[]) {
    double A = (b[1] - a[1])*(c[2] - a[2]) - (c[1] - a[1])*(b[2] - a[2]); 
    double B = (b[2] - a[2])*(c[0] - a[0]) - (c[2] - a[2])*(b[0] - a[0]); 
    double C = (b[0] - a[0])*(c[1] - a[1]) - (c[0] - a[0])*(b[1] - a[1]);
    double D = -(A*a[0] + B * a[1] + C * a[2]);
    return abs(A*p[0] + B * p[1] + C * p[2] + D) / sqrt(A*A + B * B + C * C);
}

int main() {
    int n; cin >> n;
    double p[3];
    cin >> p[0] >> p[1] >> p[2];
    double v[300][3];
    REP(i, n) {
        cin >> v[i][0] >> v[i][1] >> v[i][2];
    }
    double ans = 0;
    REP(i, n - 2) {
        FOR(j, i + 1, n - 1) {
            FOR(k, j + 1, n) {
                ans += dist(p, v[i], v[j], v[k]);
            }
        }
    }
    printf("%.10f", ans);
    getchar(); getchar();
}
0