結果

問題 No.132 点と平面との距離
コンテスト
ユーザー nebukuro09
提出日時 2016-11-10 17:37:06
言語 D
(dmd 2.109.1)
結果
AC  
実行時間 1,645 ms / 5,000 ms
コード長 986 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 2,747 ms
コンパイル使用メモリ 163,120 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-06-12 05:03:15
合計ジャッジ時間 5,530 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 3
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

import std.stdio;
import std.array;
import std.string;
import std.conv;
import std.algorithm;
import std.typecons;
import std.range;
import std.random;
import std.math;
import std.container;
import std.numeric;

void main() {
  int N = readln.chomp.to!int;
  auto P = readln.split.map!(to!real).array;
  auto Q = iota(N).map!(_ => readln.split.map!(to!real).array).array;
  real ans = 0.0;

  foreach (i; iota(N))
    foreach (j; iota(i+1, N))
      foreach (k; iota(j+1, N)) {
        auto ij = [Q[j][0]-Q[i][0], Q[j][1]-Q[i][1], Q[j][2]-Q[i][2]];
        auto ik = [Q[k][0]-Q[i][0], Q[k][1]-Q[i][1], Q[k][2]-Q[i][2]];
        auto cross = [ij[1]*ik[2]-ij[2]*ik[1],
                      ij[2]*ik[0]-ij[0]*ik[2],
                      ij[0]*ik[1]-ij[1]*ik[0]];
        cross[] /= sqrt(cross[0]^^2+cross[1]^^2+cross[2]^^2);
        auto ip = [P[0]-Q[i][0], P[1]-Q[i][1], P[2]-Q[i][2]];
        ans += abs(cross[0]*ip[0]+cross[1]*ip[1]+cross[2]*ip[2]);
      }
  writef("%.09f", ans);
}
0