結果

問題 No.622 点と三角柱の内外判定
コンテスト
ユーザー sekiya9311
提出日時 2017-12-22 19:15:04
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,324 bytes
コンパイル時間 1,596 ms
コンパイル使用メモリ 169,676 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2024-12-17 23:15:25
合計ジャッジ時間 2,552 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 24 WA * 8
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;

struct point {
  long double x, y, z;
};
point p[3], d;

inline point get_vec_ab(point pa, point pb) {
  point ab;
  ab.x = pb.x - pa.x;
  ab.y = pb.y - pa.y;
  ab.z = pb.z - pa.z;
  return ab;
}

inline long double naiseki(point pa, point pb) {
  return pa.x * pb.x + pa.y * pb.y + pa.z * pb.z;
}

inline bool is_parallel(point v_1, point v_2) {
  return powl(v_1.y * v_2.z - v_1.z * v_2.y, 2) 
    + powl(v_1.z * v_2.x - v_1.x * v_2.z, 2) 
    + powl(v_1.x * v_2.y - v_1.y * v_2.x, 2) < 1e9;
}

int main() {
  for (int i = 0; i < 3; i++) {
    cin >> p[i].x >> p[i].y >> p[i].z;
  }
  cin >> d.x >> d.y >> d.z;
  bool ok = true;
  for (int i = 0; i < 3; i++) {
    point v_d = get_vec_ab(p[i], d);
    point v_1 = get_vec_ab(p[i], p[(i + 1) % 3]);
    point v_2 = get_vec_ab(p[i], p[(i + 2) % 3]);
    ok &= naiseki(v_d, v_1) >= 0 && naiseki(v_d, v_2) >= 0;
  }
  for (int i = 0; i < 3; i++) {
    for (int j = i + 1; j < 3; j++) {
      if (p[i].x == p[j].x && p[i].y == p[j].y == p[i].z == p[j].z) {
        ok = false;
      }
    }
  }
  for (int i = 0; i < 3; i++) {
    point v_1 = get_vec_ab(p[i], p[(i + 1) % 3]);
    point v_2 = get_vec_ab(p[i], p[(i + 2) % 3]);
    if (is_parallel(v_1, v_2)) {
      ok = false;
    }
  }
  cout << (ok ? "YES" : "NO") << endl;
}
0