結果

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

ソースコード

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;
}

long double gaiseki(point a,  point b){
  return powl(a.y *b.z - a.z * b.y, 2 )+ powl(a.z * b.x - a.x * b.z, 2) + powl(a.x * b.y - a.y * b.x, 2);
}

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++) {
    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 (gaiseki(v_1, v_2) < -1e9) {
      ok = false;
    }
  }
  cout << (ok ? "YES" : "NO") << endl;
}
0