#943789 (C++23) No.2602 Real Collider

提出ソース

結果

問題	No.2602 Real Collider
ユーザー	👑 emthrm
提出日時	2024-01-12 23:34:51
言語	C++23 (gcc 12.3.0 + boost 1.83.0)
結果	AC
実行時間	206 ms / 2,000 ms
コード長	5,510 bytes
コンパイル時間	9,399 ms
コンパイル使用メモリ	431,160 KB
実行使用メモリ	6,676 KB
最終ジャッジ日時	2024-01-12 23:35:09
合計ジャッジ時間	17,325 ms
ジャッジサーバーID （参考情報）	judge11 / judge14

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	2 ms 6,676 KB
testcase_01	AC	2 ms 6,676 KB
testcase_02	AC	2 ms 6,676 KB
testcase_03	AC	2 ms 6,676 KB
testcase_04	AC	2 ms 6,676 KB
testcase_05	AC	2 ms 6,676 KB
testcase_06	AC	2 ms 6,676 KB
testcase_07	AC	2 ms 6,676 KB
testcase_08	AC	2 ms 6,676 KB
testcase_09	AC	2 ms 6,676 KB
testcase_10	AC	206 ms 6,676 KB
testcase_11	AC	66 ms 6,676 KB
testcase_12	AC	39 ms 6,676 KB
testcase_13	AC	21 ms 6,676 KB
testcase_14	AC	44 ms 6,676 KB
testcase_15	AC	51 ms 6,676 KB
testcase_16	AC	82 ms 6,676 KB
testcase_17	AC	88 ms 6,676 KB
testcase_18	AC	28 ms 6,676 KB
testcase_19	AC	36 ms 6,676 KB
testcase_20	AC	89 ms 6,676 KB
testcase_21	AC	28 ms 6,676 KB
testcase_22	AC	32 ms 6,676 KB
testcase_23	AC	23 ms 6,676 KB
testcase_24	AC	33 ms 6,676 KB
testcase_25	AC	33 ms 6,676 KB
testcase_26	AC	38 ms 6,676 KB
testcase_27	AC	37 ms 6,676 KB
testcase_28	AC	68 ms 6,676 KB
testcase_29	AC	32 ms 6,676 KB
testcase_30	AC	37 ms 6,676 KB
testcase_31	AC	81 ms 6,676 KB
testcase_32	AC	62 ms 6,676 KB
testcase_33	AC	83 ms 6,676 KB
testcase_34	AC	73 ms 6,676 KB
testcase_35	AC	55 ms 6,676 KB
testcase_36	AC	42 ms 6,676 KB
testcase_37	AC	97 ms 6,676 KB
testcase_38	AC	95 ms 6,676 KB
testcase_39	AC	81 ms 6,676 KB
testcase_40	AC	43 ms 6,676 KB
testcase_41	AC	94 ms 6,676 KB
testcase_42	AC	74 ms 6,676 KB
testcase_43	AC	62 ms 6,676 KB
testcase_44	AC	103 ms 6,676 KB
testcase_45	AC	57 ms 6,676 KB
testcase_46	AC	62 ms 6,676 KB
testcase_47	AC	87 ms 6,676 KB
testcase_48	AC	71 ms 6,676 KB
testcase_49	AC	63 ms 6,676 KB
testcase_50	AC	41 ms 6,676 KB
testcase_51	AC	42 ms 6,676 KB
testcase_52	AC	30 ms 6,676 KB
testcase_53	AC	89 ms 6,676 KB
testcase_54	AC	55 ms 6,676 KB
testcase_55	AC	70 ms 6,676 KB
testcase_56	AC	64 ms 6,676 KB
testcase_57	AC	60 ms 6,676 KB
testcase_58	AC	24 ms 6,676 KB
testcase_59	AC	82 ms 6,676 KB
testcase_60	AC	57 ms 6,676 KB
testcase_61	AC	53 ms 6,676 KB
testcase_62	AC	81 ms 6,676 KB
testcase_63	AC	82 ms 6,676 KB
testcase_64	AC	107 ms 6,676 KB
testcase_65	AC	47 ms 6,676 KB
testcase_66	AC	76 ms 6,676 KB
testcase_67	AC	42 ms 6,676 KB
testcase_68	AC	47 ms 6,676 KB
testcase_69	AC	31 ms 6,676 KB
testcase_70	AC	39 ms 6,676 KB
testcase_71	AC	56 ms 6,676 KB
testcase_72	AC	72 ms 6,676 KB
testcase_73	AC	71 ms 6,676 KB
testcase_74	AC	83 ms 6,676 KB
testcase_75	AC	87 ms 6,676 KB
testcase_76	AC	68 ms 6,676 KB
testcase_77	AC	81 ms 6,676 KB
testcase_78	AC	88 ms 6,676 KB
testcase_79	AC	83 ms 6,676 KB
testcase_80	AC	92 ms 6,676 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

#include <bits/stdc++.h>
using namespace std;
#define FOR(i,m,n) for(int i=(m);i<(n);++i)
#define REP(i,n) FOR(i,0,n)
using ll = long long;
constexpr int INF = 0x3f3f3f3f;
constexpr long long LINF = 0x3f3f3f3f3f3f3f3fLL;
constexpr double EPS = 1e-8;
constexpr int MOD = 998244353;
// constexpr int MOD = 1000000007;
constexpr int DY4[]{1, 0, -1, 0}, DX4[]{0, -1, 0, 1};
constexpr int DY8[]{1, 1, 0, -1, -1, -1, 0, 1};
constexpr int DX8[]{0, -1, -1, -1, 0, 1, 1, 1};
template <typename T, typename U>
inline bool chmax(T& a, U b) { return a < b ? (a = b, true) : false; }
template <typename T, typename U>
inline bool chmin(T& a, U b) { return a > b ? (a = b, true) : false; }
struct IOSetup {
  IOSetup() {
    std::cin.tie(nullptr);
    std::ios_base::sync_with_stdio(false);
    std::cout << fixed << setprecision(20);
  }
} iosetup;

namespace std {

__int128 gcd(__int128 x, __int128 y) {
  while (y != 0) std::swap(x %= y, y);
  return x;
}

}

template <typename T = long long>
struct Rational {
  T num, den;

  Rational() : num(0), den(1) {}
  Rational(const T num, const T den = 1) : num(num), den(den) {
    // assert(den != 0);
    reduce();
  }

  template <typename Real = long double>
  Real to_real() const { return static_cast<Real>(num) / den; }

  Rational& operator+=(const Rational& x) {
    const T g = std::gcd(den, x.den);
    num = num * (x.den / g) + x.num * (den / g);
    den *= x.den / g;
    reduce();
    return *this;
  }
  Rational& operator-=(const Rational& x) { return *this += -x; }

  Rational& operator*=(const Rational& x) {
    const T g1 = std::gcd(num, x.den), g2 = std::gcd(den, x.num);
    num = (num / g1) * (x.num / g2);
    den = (den / g2) * (x.den / g1);
    reduce();
    return *this;
  }
  Rational& operator/=(const Rational& x) {
    return *this *= Rational(x.den, x.num);
  }

  auto operator<=>(const Rational& x) const {
    return num * x.den <=> x.num * den;
  }
  bool operator==(const Rational& x) const {
    return num == x.num && den == x.den;
  }

  Rational& operator++() {
    if ((num += den) == 0) den = 1;
    return *this;
  }
  Rational operator++(int) {
    const Rational res = *this;
    ++*this;
    return res;
  }
  Rational& operator--() {
    if ((num -= den) == 0) den = 1;
    return *this;
  }
  Rational operator--(int) {
    const Rational res = *this;
    --*this;
    return res;
  }

  Rational operator+() const { return *this; }
  Rational operator-() const { return Rational(-num, den); }

  Rational operator+(const Rational& x) const { return Rational(*this) += x; }
  Rational operator-(const Rational& x) const { return Rational(*this) -= x; }
  Rational operator*(const Rational& x) const { return Rational(*this) *= x; }
  Rational operator/(const Rational& x) const { return Rational(*this) /= x; }

  friend std::ostream& operator<<(std::ostream& os, const Rational& x) {
    if (x.den == 1) return os << x.num;
    return os << x.num << '/' << x.den;
  }

 private:
  void reduce() {
    const T g = std::gcd(num, den);
    num /= g;
    den /= g;
    if (den < 0) {
      num = -num;
      den = -den;
    }
  }
};

using rational = Rational<__int128>;

using Point = pair<int, int>;
tuple<rational, rational, rational> smallest_enclosing_circle(
    const Point& p1, const Point& p2, const Point& p3) {
  const auto get_circle = [](const Point& p1, const Point& p2) -> tuple<rational, rational, rational> {
    const auto [p1x, p1y] = p1;
    const auto [p2x, p2y] = p2;
    return {rational(p1x + p2x, 2), rational(p1y + p2y, 2),
            rational((p2x - p1x) * (p2x - p1x) + (p2y - p1y) * (p2y - p1y), 4)};
  };
  auto [cx, cy, cr] = get_circle(p1, p2);
  const auto is_in = [&](const Point& p) -> bool {
    const auto [x, y] = p;
    return (cx - x) * (cx - x) + (cy - y) * (cy - y) <= cr;
  };
  if (!is_in(p3)) {
    tie(cx, cy, cr) = get_circle(p1, p3);
    if (!is_in(p2)) {
      tie(cx, cy, cr) = get_circle(p2, p3);
      if (!is_in(p1)) {
        const int a = (p3.first - p2.first) * (p3.first - p2.first) + (p3.second - p2.second) * (p3.second - p2.second);
        const int b = (p1.first - p3.first) * (p1.first - p3.first) + (p1.second - p3.second) * (p1.second - p3.second);
        const int c = (p2.first - p1.first) * (p2.first - p1.first) + (p2.second - p1.second) * (p2.second - p1.second);
        const int idx = p3.first - p1.first, idy = p3.second - p1.second;
        const int jdx = p2.first - p1.first, jdy = p2.second - p1.second;
        const int s = idx * jdy - idy * jdx;
        cx = rational(__int128{p1.first} * a * (b + c - a) + __int128{p2.first} * b * (c + a - b) + __int128{p3.first} * c * (a + b - c), 4LL * s * s);
        cy = rational(__int128{p1.second} * a * (b + c - a) + __int128{p2.second} * b * (c + a - b) + __int128{p3.second} * c * (a + b - c), 4LL * s * s);
        cr = (cx - p1.first) * (cx - p1.first) + (cy - p1.second) * (cy - p1.second);
      }
    }
  }
  return {cx, cy, cr};
}

#include <boost/multiprecision/cpp_int.hpp>
using namespace boost::multiprecision;

int main() {
  int q, xa, ya, xb, yb, xc, yc; cin >> q >> xa >> ya >> xb >> yb >> xc >> yc;
  const auto [cx, cy, cr] = smallest_enclosing_circle({xa, ya}, {xb, yb}, {xc, yc});
  while (q--) {
    int x, y; cin >> x >> y;
    const rational dx = cx - x, dy = cy - y;
    cout << ((int256_t{dx.num} * dx.num * dy.den * dy.den + int256_t{dy.num} * dy.num * dx.den * dx.den) * cr.den <= int256_t{cr.num} * dx.den * dx.den * dy.den * dy.den ? "Yes\n" : "No\n");
  }
  return 0;
}

yukicoder

結果

テストケース

ソースコード