結果

問題 No.622 点と三角柱の内外判定
ユーザー te-shte-sh
提出日時 2017-12-22 09:53:01
言語 D
(dmd 2.106.1)
結果
AC  
実行時間 2 ms / 1,500 ms
コード長 1,323 bytes
コンパイル時間 1,067 ms
コンパイル使用メモリ 116,856 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-06-12 23:14:41
合計ジャッジ時間 1,873 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
6,812 KB
testcase_01 AC 1 ms
6,944 KB
testcase_02 AC 1 ms
6,940 KB
testcase_03 AC 1 ms
6,940 KB
testcase_04 AC 1 ms
6,944 KB
testcase_05 AC 1 ms
6,940 KB
testcase_06 AC 1 ms
6,940 KB
testcase_07 AC 1 ms
6,944 KB
testcase_08 AC 1 ms
6,940 KB
testcase_09 AC 1 ms
6,940 KB
testcase_10 AC 1 ms
6,940 KB
testcase_11 AC 1 ms
6,940 KB
testcase_12 AC 1 ms
6,940 KB
testcase_13 AC 1 ms
6,940 KB
testcase_14 AC 1 ms
6,940 KB
testcase_15 AC 1 ms
6,940 KB
testcase_16 AC 1 ms
6,944 KB
testcase_17 AC 1 ms
6,944 KB
testcase_18 AC 1 ms
6,940 KB
testcase_19 AC 1 ms
6,940 KB
testcase_20 AC 2 ms
6,944 KB
testcase_21 AC 1 ms
6,944 KB
testcase_22 AC 1 ms
6,944 KB
testcase_23 AC 1 ms
6,944 KB
testcase_24 AC 1 ms
6,944 KB
testcase_25 AC 1 ms
6,940 KB
testcase_26 AC 1 ms
6,940 KB
testcase_27 AC 2 ms
6,940 KB
testcase_28 AC 1 ms
6,944 KB
testcase_29 AC 1 ms
6,940 KB
testcase_30 AC 1 ms
6,944 KB
testcase_31 AC 1 ms
6,944 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import std.algorithm, std.conv, std.range, std.stdio, std.string;
import std.bigint;    // BigInt

alias point = Point3!BigInt;

void main()
{
  auto p = new point[](4);
  foreach (i; 0..4) {
    auto rd = readln.split.to!(BigInt[]), x = rd[0], y = rd[1], z = rd[2];
    p[i] = point(x, y, z);
  }

  foreach (i; 0..3) p[i] = p[i] - p[3];

  auto v = [outerProd(p[0], p[1]), outerProd(p[1], p[2]), outerProd(p[2], p[0])];
  auto u = outerProd(p[1]-p[0], p[2]-p[0]);

  auto ip = v.map!(vi => vi * u).array;
  writeln(ip.all!"a<0" || ip.all!"a>0" ? "YES" : "NO");
}

struct Point3(T)
{
  T x, y, z;
  pure auto opBinary(string op: "+")(Point3!T rhs) const { return Point3!T(x + rhs.x, y + rhs.y, z + rhs.z); }
  pure auto opBinary(string op: "-")(Point3!T rhs) const { return Point3!T(x - rhs.x, y - rhs.y, z - rhs.z); }
  pure auto opBinary(string op: "*")(Point3!T rhs) const { return x * rhs.x + y * rhs.y + z * rhs.z; }
  pure auto opBinary(string op: "*")(T a) const { return Point3!T(x * a, y * a, z * a); }
  pure auto opBinary(string op: "/")(T a) const { return Point3!T(x / a, y / a, z / a); }
  pure auto hypot2() const { return x ^^ 2 + y ^^ 2 + z ^^ 2; }
}

pure auto outerProd(T)(Point3!T p1, Point3!T p2)
{
  return Point3!T(p1.y * p2.z - p1.z * p2.y, p1.z * p2.x - p1.x * p2.z, p1.x * p2.y - p1.y * p2.x);
}
0