結果

問題 No.622 点と三角柱の内外判定
コンテスト
ユーザー merom686
提出日時 2017-12-22 01:14:13
言語 C++14
(gcc 13.3.0 + boost 1.89.0)
結果
WA  
実行時間 -
コード長 1,290 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 735 ms
コンパイル使用メモリ 75,504 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-12-16 07:52:02
合計ジャッジ時間 1,787 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 26 WA * 6
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <iostream>
#include <vector>
#include <string>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <cmath>
using namespace std;

struct P {
    double inner(const P& p) const {
        double t = 0;
        for (int i = 0; i < 3; i++) {
            t += x[i] * p.x[i];
        }
        return t;
    }
    P& operator-=(const P& p) {
        for (int i = 0; i < 3; i++) {
            x[i] -= p.x[i];
        }
        return *this;
    }
    P& operator*=(const double t) {
        for (int i = 0; i < 3; i++) {
            x[i] *= t;
        }
        return *this;
    }
    double x[3];
};

int main() {
    P p[4];
    for (int i = 0; i < 4; i++) {
        cin >> p[i].x[0] >> p[i].x[1] >> p[i].x[2];
    }

    for (int i = 0; i < 4; i++) {
        if (i != 2) p[i] -= p[2];
    }
    p[2] -= p[2];

    double r[2], s[2];
    for (int i = 0; i < 2; i++) {
    hoge:
        double t0 = p[i].inner(p[i]);
        double t1 = p[i].inner(p[3]);
        r[i] = sqrt(t0);
        s[i] = t1 / r[i];
        if (i == 0 && s[i] == 0) {
            swap(p[0], p[1]);
            goto hoge;
        }
        if (i == 0) p[3] -= (p[0] *= s[0] / r[0]);
    }

    cout << ((s[0] / r[0] + s[1] / r[1] < 1 && s[0] > 0 && s[1] > 0) ? "YES" : "NO") << endl;

    return 0;
}
0