結果
問題 | No.132 点と平面との距離 |
ユーザー |
|
提出日時 | 2016-09-07 10:09:05 |
言語 | D (dmd 2.109.1) |
結果 |
AC
|
実行時間 | 404 ms / 5,000 ms |
コード長 | 1,037 bytes |
コンパイル時間 | 2,923 ms |
コンパイル使用メモリ | 167,524 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-06-12 04:12:01 |
合計ジャッジ時間 | 3,546 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 3 |
ソースコード
import std.algorithm, std.array, std.container, std.range, std.bitmanip;import std.numeric, std.math, std.bigint, std.random, core.bitop;import std.string, std.regex, std.conv, std.stdio, std.typecons;struct point {real x, y, z;}void main(){auto n = readln.chomp.to!size_t;auto ai = iota(n + 1).map!(_ => readln.split.map!(to!real)).map!(rd => point(rd[0], rd[1], rd[2])).array;auto p = ai.front;auto qi = ai[1..$];auto bi = new bool[](n);bi[3..$] = true;real s = 0;foreach (i; 0..n - 2)foreach (j; i + 1..n - 1)foreach (k; j + 1..n)s += calc(p, qi[i], qi[j], qi[k]);writefln("%.9f", s);}real calc(point p, point a, point b, point c){auto d = (b.y - a.y) * (c.z - a.z) - (c.y - a.y) * (b.z - a.z);auto e = (b.z - a.z) * (c.x - a.x) - (c.z - a.z) * (b.x - a.x);auto f = (b.x - a.x) * (c.y - a.y) - (c.x - a.x) * (b.y - a.y);auto g = -(d * a.x + e * a.y + f * a.z);return (d * p.x + e * p.y + f * p.z + g).abs / (d ^^ 2 + e ^^ 2 + f ^^ 2).sqrt;}