結果

問題 No.622 点と三角柱の内外判定
ユーザー rickythetarickytheta
提出日時 2017-12-22 01:20:37
言語 C++11
(gcc 13.3.0)
結果
AC  
実行時間 2 ms / 1,500 ms
コード長 1,894 bytes
コンパイル時間 1,302 ms
コンパイル使用メモリ 159,852 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-12-16 07:52:41
合計ジャッジ時間 2,168 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 32
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;

typedef long long ll;
typedef vector<int> vi;
typedef vector<ll> vl;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;

typedef int _loop_int;
#define REP(i,n) for(_loop_int i=0;i<(_loop_int)(n);++i)
#define FOR(i,a,b) for(_loop_int i=(_loop_int)(a);i<(_loop_int)(b);++i)
#define FORR(i,a,b) for(_loop_int i=(_loop_int)(b)-1;i>=(_loop_int)(a);--i)

#define DEBUG(x) cout<<#x<<": "<<x<<endl
#define DEBUG_VEC(v) cout<<#v<<":";REP(i,v.size())cout<<" "<<v[i];cout<<endl
#define ALL(a) (a).begin(),(a).end()

#define CHMIN(a,b) a=min((a),(b))
#define CHMAX(a,b) a=max((a),(b))

// mod
const ll MOD = 1000000007ll;
#define FIX(a) ((a)%MOD+MOD)%MOD

// floating
typedef double Real;
const Real EPS = 1e-11;
#define EQ0(x) (abs(x)<EPS)
#define EQ(a,b) (abs(a-b)<EPS)
typedef complex<Real> P;

Real x[4], y[4], z[4];

P ps[4];

int ccw(P a, P b, P c){
  b -= a;
  c -= a;
  return imag(conj(b)*c) < 0.0;
}

int main(){
  REP(i,4)cin>>x[i]>>y[i]>>z[i];
  // origin point(2)
  REP(i,4)if(i!=2){
    x[i]-=x[2];
    y[i]-=y[2];
    z[i]-=z[2];
  }
  x[2] = y[2] = z[2] = 0.0;
  // (xy) rotate point(1) - origin
  {
    P rot = P(x[1], y[1]);
    rot /= abs(rot);
    REP(i,4){
      P p = P(x[i], y[i]) / rot;
      x[i] = real(p);
      y[i] = imag(p);
    }
  }
  // y[1] == 0
  // (xz) rotate point(1) - origin
  {
    P rot = P(x[1], z[1]);
    rot /= abs(rot);
    REP(i,4){
      P p = P(x[i], z[i]) / rot;
      x[i] = real(p);
      z[i] = imag(p);
    }
  }
  // z[1] == 0
  // (yz) rotate point(0) - origin
  {
    P rot = P(y[0], z[0]);
    rot /= abs(rot);
    REP(i,4){
      P p = P(y[i], z[i]) / rot;
      y[i] = real(p);
      z[i] = imag(p);
    }
  }
  // ok
  REP(i,4)ps[i] = P(x[i], y[i]);
  int xs = 0;
  REP(i,3)xs += ccw(ps[i], ps[(i+1)%3], ps[3]);
  if(xs==0 || xs==3){
    puts("YES");
  }else{
    puts("NO");
  }
  return 0;
}
0