結果

問題 No.2602 Real Collider
ユーザー tnakao0123tnakao0123
提出日時 2024-01-14 02:21:54
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 2,532 bytes
コンパイル時間 551 ms
コンパイル使用メモリ 44,868 KB
実行使用メモリ 6,676 KB
最終ジャッジ日時 2024-01-14 02:22:00
合計ジャッジ時間 5,320 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge14
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,676 KB
testcase_01 AC 1 ms
6,676 KB
testcase_02 AC 2 ms
6,676 KB
testcase_03 AC 2 ms
6,676 KB
testcase_04 WA -
testcase_05 AC 1 ms
6,676 KB
testcase_06 WA -
testcase_07 WA -
testcase_08 WA -
testcase_09 WA -
testcase_10 WA -
testcase_11 AC 17 ms
6,676 KB
testcase_12 WA -
testcase_13 WA -
testcase_14 WA -
testcase_15 AC 13 ms
6,676 KB
testcase_16 AC 20 ms
6,676 KB
testcase_17 AC 22 ms
6,676 KB
testcase_18 WA -
testcase_19 WA -
testcase_20 AC 25 ms
6,676 KB
testcase_21 WA -
testcase_22 WA -
testcase_23 WA -
testcase_24 WA -
testcase_25 WA -
testcase_26 AC 13 ms
6,676 KB
testcase_27 WA -
testcase_28 AC 21 ms
6,676 KB
testcase_29 WA -
testcase_30 WA -
testcase_31 AC 20 ms
6,676 KB
testcase_32 AC 17 ms
6,676 KB
testcase_33 AC 20 ms
6,676 KB
testcase_34 AC 21 ms
6,676 KB
testcase_35 AC 13 ms
6,676 KB
testcase_36 AC 13 ms
6,676 KB
testcase_37 AC 22 ms
6,676 KB
testcase_38 AC 22 ms
6,676 KB
testcase_39 AC 21 ms
6,676 KB
testcase_40 AC 11 ms
6,676 KB
testcase_41 AC 24 ms
6,676 KB
testcase_42 AC 20 ms
6,676 KB
testcase_43 AC 20 ms
6,676 KB
testcase_44 AC 24 ms
6,676 KB
testcase_45 AC 16 ms
6,676 KB
testcase_46 AC 15 ms
6,676 KB
testcase_47 AC 22 ms
6,676 KB
testcase_48 AC 16 ms
6,676 KB
testcase_49 AC 15 ms
6,676 KB
testcase_50 AC 12 ms
6,676 KB
testcase_51 AC 12 ms
6,676 KB
testcase_52 AC 10 ms
6,676 KB
testcase_53 AC 20 ms
6,676 KB
testcase_54 AC 15 ms
6,676 KB
testcase_55 AC 18 ms
6,676 KB
testcase_56 AC 18 ms
6,676 KB
testcase_57 AC 16 ms
6,676 KB
testcase_58 AC 7 ms
6,676 KB
testcase_59 AC 19 ms
6,676 KB
testcase_60 AC 17 ms
6,676 KB
testcase_61 AC 13 ms
6,676 KB
testcase_62 AC 20 ms
6,676 KB
testcase_63 AC 23 ms
6,676 KB
testcase_64 AC 26 ms
6,676 KB
testcase_65 AC 14 ms
6,676 KB
testcase_66 AC 22 ms
6,676 KB
testcase_67 AC 11 ms
6,676 KB
testcase_68 AC 12 ms
6,676 KB
testcase_69 AC 9 ms
6,676 KB
testcase_70 AC 11 ms
6,676 KB
testcase_71 AC 14 ms
6,676 KB
testcase_72 AC 20 ms
6,676 KB
testcase_73 AC 16 ms
6,676 KB
testcase_74 AC 20 ms
6,676 KB
testcase_75 AC 21 ms
6,676 KB
testcase_76 AC 18 ms
6,676 KB
testcase_77 AC 19 ms
6,676 KB
testcase_78 AC 24 ms
6,676 KB
testcase_79 AC 21 ms
6,676 KB
testcase_80 AC 24 ms
6,676 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

/* -*- coding: utf-8 -*-
 *
 * 2602.cc:  No.2602 Real Collider - yukicoder
 */

#include<cstdio>
#include<algorithm>
 
using namespace std;

/* constant */

const double DELTA = 1e-7;

/* typedef */

template <typename T>
struct Pt {
  T x, y;
  Pt() {}
  Pt(T _x, T _y) : x(_x), y(_y) {}
  Pt(const Pt<T> &p) : x(p.x), y(p.y) {}

  Pt<T> operator+(const Pt<T> p) const { return Pt<T>(x + p.x, y + p.y); }
  Pt<T> operator-() const { return Pt<T>(-x, -y); }
  Pt<T> operator-(const Pt<T> p) const { return Pt<T>(x - p.x, y - p.y); }
  Pt<T> operator*(T t) const { return Pt<T>(x * t, y * t); }
  Pt<T> operator/(T t) const { return Pt<T>(x / t, y / t); }
  T dot(Pt<T> v) const { return x * v.x + y * v.y; }
  T cross(Pt<T> v) const { return x * v.y - y * v.x; }
  Pt<T> mid(const Pt<T> p) { return Pt<T>((x + p.x) / 2, (y + p.y) / 2); }
  T d2() { return x * x + y * y; }

  Pt<T> rot90() { return Pt<T>(-y, x); }

  bool operator==(const Pt<T> pt) const { return x == pt.x && y == pt.y; }
  bool operator<(const Pt<T> &pt) const {
    return x < pt.x || (x == pt.x && y < pt.y);
  }
};

typedef Pt<double> pt;

struct CL {
  pt p;
  double t0, t1;
  CL() {}
  CL(const pt& _p, double _t0, double _t1) : p(_p), t0(_t0), t1(_t1) {}
};

/* global variables */

pt readpt() {
  int x, y;
  scanf("%d%d", &x, &y);
  return pt(x, y);
}

/* subroutines */

bool cross_lines(const pt& ap, const pt av, const pt& bp, const pt bv, CL& cl) {
  double op01 = av.cross(bv);
  //if (op01 == 0.0) return false; /* need to handle parallel?? */

  if (op01 == 0.0) {
    pt v = bp - ap;
    if (v.cross(av) != 0.0) return false;

    pt a1 = ap + av;
    pt b1 = bp + bv;

    return
      ((bp - ap).dot(b1 - ap) <= 0.0 ||
       (bp - a1).dot(b1 - a1) <= 0.0 ||
       (ap - bp).dot(a1 - bp) <= 0.0 ||
       (ap - b1).dot(a1 - b1) <= 0.0);
  }

  pt v = bp - ap;
  double op0 = v.cross(av);
  double op1 = v.cross(bv);

  double t0 = op1 / op01;
  double t1 = op0 / op01;

  cl.p = bv * t1 + bp;
  cl.t0 = t0;
  cl.t1 = t1;

  return true;
  //return (0.0 <= cl.t0 && cl.t0 <= 1.0 && 0.0 <= cl.t1 && cl.t1 <= 1.0);
}

/* main */

int main() {
  int qn;
  scanf("%d", &qn);
  pt pa = readpt(), pb = readpt(), pc = readpt();

  pt p0 = (pa + pb) / 2, p1 = (pa + pc) / 2;
  pt v0 = (pb - pa).rot90(), v1 = (pc - pa).rot90();
  CL cl;

  cross_lines(p0, v0, p1, v1, cl);
  pt cp = cl.p;
  double rr = (pa - cp).d2();

  while (qn--) {
    pt pi = readpt();
    if ((pi - cp).d2() <= rr + DELTA) puts("Yes");
    else puts("No");
  }

  return 0;
}
0