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

提出ソース

結果

問題	No.2602 Real Collider
ユーザー	👑 emthrm
提出日時	2024-01-12 23:22:31
言語	C++23 (gcc 12.3.0 + boost 1.83.0)
結果	WA
実行時間	-
コード長	5,272 bytes
コンパイル時間	2,802 ms
コンパイル使用メモリ	259,052 KB
実行使用メモリ	5,376 KB
最終ジャッジ日時	2024-09-28 00:06:48
合計ジャッジ時間	23,318 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	AC	1 ms 5,376 KB
testcase_05	AC	1 ms 5,376 KB
testcase_06	AC	2 ms 5,376 KB
testcase_07	AC	1 ms 5,376 KB
testcase_08	AC	2 ms 5,376 KB
testcase_09	AC	1 ms 5,376 KB
testcase_10	AC	881 ms 5,376 KB
testcase_11	AC	263 ms 5,376 KB
testcase_12	AC	70 ms 5,376 KB
testcase_13	AC	41 ms 5,376 KB
testcase_14	AC	79 ms 5,376 KB
testcase_15	AC	211 ms 5,376 KB
testcase_16	AC	331 ms 5,376 KB
testcase_17	AC	359 ms 5,376 KB
testcase_18	AC	48 ms 5,376 KB
testcase_19	AC	69 ms 5,376 KB
testcase_20	AC	341 ms 5,376 KB
testcase_21	AC	50 ms 5,376 KB
testcase_22	AC	55 ms 5,376 KB
testcase_23	AC	38 ms 5,376 KB
testcase_24	AC	59 ms 5,376 KB
testcase_25	AC	59 ms 5,376 KB
testcase_26	AC	137 ms 5,376 KB
testcase_27	AC	68 ms 5,376 KB
testcase_28	AC	248 ms 5,376 KB
testcase_29	AC	60 ms 5,376 KB
testcase_30	AC	69 ms 5,376 KB
testcase_31	AC	330 ms 5,376 KB
testcase_32	AC	254 ms 5,376 KB
testcase_33	WA	-
testcase_34	AC	307 ms 5,376 KB
testcase_35	AC	227 ms 5,376 KB
testcase_36	AC	158 ms 5,376 KB
testcase_37	AC	395 ms 5,376 KB
testcase_38	AC	380 ms 5,376 KB
testcase_39	AC	337 ms 5,376 KB
testcase_40	AC	179 ms 5,376 KB
testcase_41	AC	382 ms 5,376 KB
testcase_42	AC	292 ms 5,376 KB
testcase_43	AC	235 ms 5,376 KB
testcase_44	AC	415 ms 5,376 KB
testcase_45	AC	227 ms 5,376 KB
testcase_46	AC	252 ms 5,376 KB
testcase_47	AC	359 ms 5,376 KB
testcase_48	AC	284 ms 5,376 KB
testcase_49	AC	251 ms 5,376 KB
testcase_50	AC	167 ms 5,376 KB
testcase_51	AC	176 ms 5,376 KB
testcase_52	AC	124 ms 5,376 KB
testcase_53	AC	374 ms 5,376 KB
testcase_54	AC	235 ms 5,376 KB
testcase_55	AC	286 ms 5,376 KB
testcase_56	AC	267 ms 5,376 KB
testcase_57	AC	243 ms 5,376 KB
testcase_58	AC	87 ms 5,376 KB
testcase_59	WA	-
testcase_60	AC	220 ms 5,376 KB
testcase_61	AC	206 ms 5,376 KB
testcase_62	AC	325 ms 5,376 KB
testcase_63	AC	321 ms 5,376 KB
testcase_64	AC	436 ms 5,376 KB
testcase_65	AC	187 ms 5,376 KB
testcase_66	AC	316 ms 5,376 KB
testcase_67	AC	157 ms 5,376 KB
testcase_68	AC	178 ms 5,376 KB
testcase_69	AC	112 ms 5,376 KB
testcase_70	AC	149 ms 5,376 KB
testcase_71	AC	236 ms 5,376 KB
testcase_72	AC	291 ms 5,376 KB
testcase_73	AC	282 ms 5,376 KB
testcase_74	AC	345 ms 5,376 KB
testcase_75	AC	355 ms 5,376 KB
testcase_76	AC	282 ms 5,376 KB
testcase_77	AC	316 ms 5,376 KB
testcase_78	AC	361 ms 5,376 KB
testcase_79	AC	318 ms 5,376 KB
testcase_80	AC	356 ms 5,376 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};
}

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;
    cout << ((cx - x) * (cx - x) + (cy - y) * (cy - y) <= cr ? "Yes\n" : "No\n");
  }
  return 0;
}

yukicoder

結果

テストケース

ソースコード