結果

問題 No.622 点と三角柱の内外判定
ユーザー sekiya9311sekiya9311
提出日時 2017-12-22 19:06:18
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 1,200 bytes
コンパイル時間 1,598 ms
コンパイル使用メモリ 167,620 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-05-10 02:59:06
合計ジャッジ時間 2,540 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

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

ソースコード

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 ||naiseki(v_d, v_1) > 0 && naiseki(v_d, v_2) >=0 || 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