結果

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

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
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 1 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 1 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 1 ms
5,376 KB
testcase_20 WA -
testcase_21 AC 2 ms
5,376 KB
testcase_22 AC 1 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 1 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;
};
bool operator==(const point& pa, const point& pb) {
  return pa.x == pb.x && pa.y == pb.y && pa.z == pb.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;
  }

  if (p[0] == p[1] && p[1] == p[2]) {
    assert(false);
    ok = true;
  } else {
    for (int i = 0; i < 3; i++) {
      if (p[i] == p[(i + 1) % 3]) {
        assert(false);
        point a = get_vec_ab(p[i], p[(i + 2) % 3]);
        point v = get_vec_ab(p[i], d);
        auto va = naiseki(v, a);
        auto aa = naiseki(a, a);
        cerr << "va = " << va << endl;
        cerr << "aa = " << aa << endl;
        if (va <= 0.0  || va >= aa) {
          ok = false;
        } else {
          ok = true;
        }
        break;
      }
    }
  }
  cout << (ok ? "YES" : "NO") << endl;
}
0