tests/Minimum_Enclosing_Circle.test.cpp

View this file on GitHub
Last update: 2025-11-28 02:09:51+07:00
Problem: https://judge.yosupo.jp/problem/minimum_enclosing_circle

Depends on

Code

#define PROBLEM "https://judge.yosupo.jp/problem/minimum_enclosing_circle"

#include "../misc/macros.h"
#include "../geometry/MinimumEnclosingCircle.h"

using P = Point<ld>;
void solve() {
  int n;
  cin >> n;
  vector<P> pts(n);
  for (auto& [x, y] : pts) cin >> x >> y;
  auto [o, r] = mec(pts);
  const ld EPS = 1e-10;
  for (int i = 0; i < n; ++i) {
    if (fabsl((o - pts[i]).dist2() - r * r) < EPS) {
      cout << 1;
    } else {
      cout << 0;
    }
  }
}

int main() {
  cin.tie(0)->sync_with_stdio(0);
  cin.exceptions(cin.failbit);
  int tc = 1;
  // cin >> tc;
  for (int i = 1; i <= tc; ++i) {
    solve();
  }
}

#line 1 "tests/Minimum_Enclosing_Circle.test.cpp"
#define PROBLEM "https://judge.yosupo.jp/problem/minimum_enclosing_circle"

#line 1 "misc/macros.h"
// #pragma GCC optimize("Ofast,unroll-loops")       // unroll long, simple loops
// #pragma GCC target("avx2,fma")                   // vectorizing code
// #pragma GCC target("lzcnt,popcnt,abm,bmi,bmi2")  // for fast bitset operation

#include <bits/extc++.h>
#include <tr2/dynamic_bitset>

using namespace std;
using namespace __gnu_pbds;  // ordered_set, gp_hash_table
// using namespace __gnu_cxx; // rope

// for templates to work
#define all(x) (x).begin(), (x).end()
#define sz(x) (int) (x).size()
#define pb push_back
#define eb emplace_back
using i32 = int32_t;
using u32 = uint32_t;
using i64 = int64_t;
using u64 = uint64_t;
using i128 = __int128_t;
using u128 = __uint128_t;
using ld = long double;
using pii = pair<i32, i32>;
using vi = vector<i32>;

// fast map
const int RANDOM = chrono::high_resolution_clock::now().time_since_epoch().count();
struct chash {  // customize hash function for gp_hash_table
  int operator()(int x) const { return x ^ RANDOM; }
};
gp_hash_table<int, int, chash> table;

/* ordered set
    find_by_order(k): returns an iterator to the k-th element (0-based)
    order_of_key(k): returns the number of elements in the set that are strictly less than k
*/
template <class T>
using ordered_set = tree<T, null_type, less<T>, rb_tree_tag, tree_order_statistics_node_update>;

/*  rope
    rope <int> cur = v.substr(l, r - l + 1);
    v.erase(l, r - l + 1);
    v.insert(v.mutable_begin(), cur);
*/
#line 2 "geometry/Point.h"

template <class T>
int sgn(T x) { return (x > 0) - (x < 0); }
template <class T>
struct Point {
  typedef Point P;
  T x, y;
  explicit Point(T x = 0, T y = 0) : x(x), y(y) {}
  bool operator<(P p) const { return tie(x, y) < tie(p.x, p.y); }
  bool operator==(P p) const { return tie(x, y) == tie(p.x, p.y); }
  P operator+(P p) const { return P(x + p.x, y + p.y); }
  P operator-(P p) const { return P(x - p.x, y - p.y); }
  P operator*(T d) const { return P(x * d, y * d); }
  P operator/(T d) const { return P(x / d, y / d); }
  T dot(P p) const { return x * p.x + y * p.y; }
  T cross(P p) const { return x * p.y - y * p.x; }
  T cross(P a, P b) const { return (a - *this).cross(b - *this); }
  T dist2() const { return x * x + y * y; }
  T dist() const { return sqrt(dist2()); }
  // angle to x-axis in interval [-pi, pi]
  T angle() const { return atan2l(y, x); }
  P unit() const { return *this / dist(); }  // makes dist()=1
  P perp() const { return P(-y, x); }        // rotates +90 degrees
  P normal() const { return perp().unit(); }
  // returns point rotated 'a' radians ccw around the origin
  P rotate(ld a) const {
    return P(x * cos(a) - y * sin(a), x * sin(a) + y * cos(a));
  }
  friend ostream& operator<<(ostream& os, P p) {
    return os << "(" << p.x << "," << p.y << ")";
  }
};
#line 2 "geometry/Circumcircle.h"

typedef Point<ld> P;
ld ccRadius(const P& A, const P& B, const P& C) {
  return (B - A).dist() * (C - B).dist() * (A - C).dist() /
         abs((B - A).cross(C - A)) / 2;
}
P ccCenter(const P& A, const P& B, const P& C) {
  P b = C - A, c = B - A;
  return A + (b * c.dist2() - c * b.dist2()).perp() / b.cross(c) / 2;
}
#line 2 "geometry/MinimumEnclosingCircle.h"

pair<P, ld> mec(vector<P> ps) {
  shuffle(all(ps), mt19937(time(0)));
  P o = ps[0];
  ld r = 0, EPS = 1 + 1e-12;
  for (int i = 0; i < sz(ps); ++i) {
    if ((o - ps[i]).dist() > r * EPS) {
      o = ps[i], r = 0;
      for (int j = 0; j < i; ++j) {
        if ((o - ps[j]).dist() > r * EPS) {
          o = (ps[i] + ps[j]) / 2;
          r = (o - ps[i]).dist();
          for (int k = 0; k < j; ++k) {
            if ((o - ps[k]).dist() > r * EPS) {
              o = ccCenter(ps[i], ps[j], ps[k]);
              r = (o - ps[i]).dist();
            }
          }
        }
      }
    }
  }
  return {o, r};
}
#line 5 "tests/Minimum_Enclosing_Circle.test.cpp"

using P = Point<ld>;
void solve() {
  int n;
  cin >> n;
  vector<P> pts(n);
  for (auto& [x, y] : pts) cin >> x >> y;
  auto [o, r] = mec(pts);
  const ld EPS = 1e-10;
  for (int i = 0; i < n; ++i) {
    if (fabsl((o - pts[i]).dist2() - r * r) < EPS) {
      cout << 1;
    } else {
      cout << 0;
    }
  }
}

int main() {
  cin.tie(0)->sync_with_stdio(0);
  cin.exceptions(cin.failbit);
  int tc = 1;
  // cin >> tc;
  for (int i = 1; i <= tc; ++i) {
    solve();
  }
}