#include <bits/stdc++.h>

using namespace std;

#define pb push_back
#define mp make_pair
#define mt make_tuple
#define fi first
#define se second
#define all(a) (a).begin(), (a).end()
#define endl "\n"
#define rep(i, a, b) for (auto i = (a); i < (b); ++i)
#define rrep(i, a, b) for (auto i = (a); i > (b); --i)
#define UNIQUE(v) (v).erase(unique((v).begin(), (v).end()), (v).end())
typedef long long ll;
typedef unsigned long long ull;

#ifdef LOCAL_DEFINE
struct InitInput
{
  InitInput()
  {
    FILE *stream1;
    stream1 = freopen("in.txt", "r", stdin);
    assert(stream1 != nullptr);
    cerr << "This problem is not interactive" << endl;
  }
} LOCAL_INPUT;
struct LOCAL_OUTPUT
{
  LOCAL_OUTPUT()
  {
    FILE *stream2;
    const char *outputfile = "out.txt";
    stream2 = freopen(outputfile, "w", stdout);
    assert(stream2 != nullptr);
    cerr << "output [ " << outputfile << " ]" << endl;
  }
} /*LOCAL_OUTPUT*/;
#define show(x) cerr<<#x<<" = "<<(x)<<" (line "<<__LINE__<< ")"<<endl
#define showA(a, n) \
  do{for(int i=0;i<(n);i++)cerr<<"("<<i<<" = "<<(a)[i]<<") ";cerr<<endl;}while(0)
#define showA2(a, n, m) \
  do{for(int i=0;i<(n);i++){for(int j=0;j<(m);j++){cerr<<"("<<i<<", "<<j<<" = "<<(a)[i][j]<<") ";}cerr<<endl;}cerr<<endl;}while(0)
#else
#define show(x)
#define showA(a, n)
#define showA2(a, n, m)
#endif
struct InitAim
{
  InitAim()
  {
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    cout.precision(12);
    cout << fixed;
#ifdef LOCAL_DEFINE
    cerr << "This problem is not interactive" << endl;
#endif
  }
} aim_cpo;

///////////////////////////////////////////////////////////////////////////////////
// TEMPLATE(data structure)
///////////////////////////////////////////////////////////////////////////////////

template <typename T>
bool chmin(T &a, T b) { return a > b ? (a = b, true) : false; }
template <typename T>
bool chmax(T &a, T b) { return a < b ? (a = b, true) : false; }
template <typename T>
void ln(T i, T n) { cout << (i == n - 1 ? "\n" : " "); }
template <typename T, typename S>
ostream &operator<<(ostream &out, const pair<T, S> &pair1)
{
  out << '(' << pair1.fi << ", " << pair1.se << ')';
  return out;
}
template <typename T>
ostream &operator<<(ostream &out, const vector<T> &vector1)
{
  out << '[';
  for (unsigned int i = 0; i < vector1.size(); i++)
  {
    out << vector1[i];
    if (i == vector1.size() - 1)
      out << "]";
    else
      out << ", ";
  }
  return out;
}

#define GCD(a, b) __gcd(a, b)
template <typename T>
inline T LCM(T a, T b)
{
  return a / GCD(a, b) * b;
}
template <typename T>
inline T EXTGCD(T a, T b, T &x, T &y)
{
  T d = a;
  if (b != 0)
  {
    d = EXTGCD(b, a % b, y, x);
    y -= (a / b) * x;
  }
  else
    x = 1, y = 0;
  return d;
}
template <typename T>
inline bool is_prime(T a)
{
  for (int i = 2; i * i <= a; i++)
    if (a % i == 0)
      return true;
  return false;
}
template <typename T, typename S>
inline T Pow(T a, S b)
{
  T res = 1, now = a;
  while (b)
  {
    if (b & 1)
      res *= now;
    b >>= 1;
    now *= now;
  }
  return res;
}

/* MOD */
ll MOD = 1000000000LL + 7LL;
class Mint {
public:
  ll x;
  Mint() : x(0) {}
  Mint(const int x) {
    if (x < MOD && x >= 0) this->x = x;
    else this->x = (MOD + x % MOD) % MOD;
  }
  Mint(const ll x) {
    if (x < MOD && x >= 0) this->x = x;
    else this->x = (MOD + x % MOD) % MOD;
  }
  Mint(const Mint& c) {
    this->x = c.x;
  }
  Mint& operator=(const int c) {
    if (c < MOD && c >= 0) this->x = c;
    else this->x = (MOD + c % MOD) % MOD;
    return *this;
  }
  Mint& operator=(const ll c) {
    if (c < MOD && c >= 0) this->x = c;
    else this->x = (MOD + c % MOD) % MOD;
    return *this;
  }
  Mint& operator=(const Mint c) {
    this->x = c.x;
    return *this;
  }
  bool operator==(const Mint &b) {
    return x == b.x;
  }
  bool operator!=(const Mint &b) {
    return !(x == b.x);
  }
  bool operator<(const Mint &b) {
    return x < b.x;
  }
  bool operator>(const Mint &b) {
    return x > b.x;
  }
};
ostream &operator<<(ostream& out, const Mint &mint) {
  out << mint.x;
  return out;
}
istream &operator>>(istream& in, const Mint &mint) {
  in >> mint.x;
  return in;
}
inline Mint add(const Mint a, const Mint b) {
  const ll x = (a.x + b.x) % MOD;
  return Mint(x);
}
inline Mint mul(const Mint a, const Mint b) {
  ll x = (a.x * b.x) % MOD;
  return Mint(x);
}
template <typename T>
Mint ModPow(Mint a, T b) {
  assert(b >= 0);
  Mint res = 1;
  while (b)
  {
    if (b & 1)
      res = mul(res, a);
    b >>= 1;
    a = mul(a, a);
  }
  return res;
}
inline Mint ModInverse(const Mint x) {
  ll a = x.x, b = MOD, u = 1, v = 0;
  while (b) {
    ll t = a / b;
    swap(a -= t * b, b);
    swap(u -= t * v, v);
  }
  Mint res(u);
  return res;
}
inline Mint operator-(const Mint x) {
  return Mint(-x.x);
}
inline Mint operator+(const Mint a, const Mint b) {
  return add(a, b);
}
inline Mint operator*(const Mint a, const Mint b) {
  return mul(a, b);
}
inline Mint operator-(const Mint a, const Mint b) {
  return add(a, -b);
}
template<typename T>
inline Mint operator+(const Mint a, const T b) {
  return add(a, Mint(b));
}
template<typename T>
inline Mint operator*(const Mint a, const T b) {
  return mul(a, Mint(b));
}
template<typename T>
inline Mint operator-(const Mint a, const T b) {
  return add(a, Mint(-b));
}
template<typename T>
inline Mint operator/(const Mint a, const T b) {
  return mul(a, ModInverse(Mint(b)));
}
template<typename T>
inline Mint operator^(const Mint a, const T b) {
  return ModPow(a, b);
}
template<typename T>
inline Mint operator+=(Mint& a, const T b) {
  a = a + b;
  return a;
}
template<typename T>
inline Mint operator-=(Mint& a, const T b) {
  a = a - b;
  return a;
}
template<typename T>
inline Mint operator*=(Mint& a, const T b) {
  a = a * b;
  return a;
}
template<typename T>
inline Mint operator/=(Mint& a, const T b) {
  a = a / b;
  return a;
}
template <typename T>
T EulerTotient(T a)
{
  vector<pair<int, int>> v;
  for (T i = 2; i * i <= a; i++)
  {
    int cnt = 0;
    while (a % i == 0)
    {
      cnt++;
      a /= i;
    }
    if (cnt != 0)
      v.emplace_back(i, cnt);
  }
  if (a != 1)
    v.emplace_back(a, 1);
  //showV(v, (int) v.size());
  T res = 1;
  for (int i = 0; i < (int)v.size(); i++)
  {
    if (v[i].se == 1)
    {
      res *= v[i].fi - 1;
    }
    else
    {
      res *= Pow(v[i].fi, v[i].se) - Pow(v[i].fi, v[i].se - 1);
    }
  }
  return res;
}
template<typename T, typename S>
Mint comb(T a, S b)
{
  chmin(b, a - b);
  Mint res = 1LL, now = a;
  for (ll i = 1; i <= b; i++)
  {
    res *= now;
    res /= i;
    now -= 1;
  }
  return res;
}

template <typename T>
class BIT
{
public:
  BIT(int size)
  {
    BITTable.assign(++size, 0);
  }
  T sum(int k)
  {
    T res = 0;
    for (++k; k > 0; k -= k & -k)
    {
      res += BITTable[k];
    }
    return res;
  }
  T sum(int l, int r) //[,]
  {
    if (l == 0)
      return sum(r);
    return sum(r) - sum(l - 1);
  }
  void update(int k, T x)
  { // b[k] += x;
    for (++k; k < (int)BITTable.size(); k += k & -k)
      BITTable[k] += x;
  }

private:
  vector<T> BITTable;
};

template <typename T>
class IntervalTree
{
  using F = function<T(T, T)>;

public:
  IntervalTree(int n, const F func, const T init) : func(func), init(init)
  {
    size = 1;
    while ((int)size < n)
      size <<= 1;
    table.assign(2 * size, init);
  }
  void set(int k, T &x)
  {
    table[size + k] = x;
  }
  void build()
  {
    for (int i = size - 1; i >= 0; --i)
    {
      table[i] = func(table[i * 2], table[i * 2 + 1]);
    }
  }
  void update(int k, const T &x)
  {
    k += size;
    table[k] = x;
    while (k >>= 1)
    {
      table[k] = func(table[k * 2], table[k * 2 + 1]);
    }
  }
  T query(int a, int b)
  {
    T L = init, R = init;
    for (a += size, b += size; a < b; a >>= 1, b >>= 1)
    {
      if (a & 1)
        L = func(L, table[a++]);
      if (b & 1)
        R = func(table[--b], R);
    }
    return func(L, R);
  }
  T operator[](const int k) const
  {
    return table[k + size];
  }

private:
  unsigned int size;
  vector<T> table;
  const F func;
  const T init;
};

template <typename T, typename S = T>
class LazyIntervalTree {
    using F = function<T(T, T)>;
    using G = function<T(T, S, int)>;
    using H = function<S(S, S)>;
public:
    LazyIntervalTree(int n, const F f, const G g, const H h, const T &init, const S OM)
    : f(f), g(g), h(h), init(init), OM(OM) {
      size = 1;
      while (size < n) size <<= 1;
      data.assign(2 * size, init);
      lazy.assign(2 * size, OM);
    }
    void set(int k, const T &x) {
      data[k + size] = x;
    }
    void build() {
      for (int i = size - 1; i > 0; i--) {
        data[i] = f(data[2 * i], data[2 * i + 1]);
      }
    }
    void propagate(int k, int len) {
      if (lazy[k] != OM) {
        if (k < size) {
          lazy[2 * k] = h(lazy[2 * k], lazy[k]);
          lazy[2 * k + 1] = h(lazy[2 * k + 1], lazy[k]);
        }
        data[k] = g(data[k], lazy[k], len);
        lazy[k] = OM;
      }
    }
    T update(int a, int b, const S &x, int k, int l, int r) {
      propagate(k, r - l);
      if (r <= a || b <= l) {
        return data[k];
      } else if (a <= l && r <= b) {
        lazy[k] = h(lazy[k], x);
        propagate(k, r - l);
        return data[k];
      } else {
        return data[k] = f(update(a, b, x, 2 * k, l, (l + r) >> 1),
          update(a, b, x, 2 * k + 1, (l + r) >> 1, r));
      }
    }
    T update(int a, int b, const S &x) {
      return update(a, b, x, 1, 0, size);
    }
    T query(int a, int b, int k, int l, int r) {
      propagate(k, r - l);
      if (r <= a || b <= l) {
        return init;
      } else if (a <= l && r <= b) {
        return data[k];
      } else {
        return f(query(a, b, 2 * k, l, (l + r) >> 1),
          query(a, b, 2 * k + 1, (l + r) >> 1, r));
      }
    }
    T query(int a, int b) {
      return query(a, b, 1, 0, size);
    }
    T operator[](const int &k) {
      return query(k, k + 1);
    }
    void showLazy() {
      for (int i = 0; i < (int)lazy.size(); i++) {
        cout << lazy[i]; ln(i, (int) lazy.size());
      }
    }
private:
    int size;
    vector <T> data;
    vector <S> lazy;
    const F f;
    const G g;
    const H h;
    const T init;
    const S OM;
};

class UnionFind
{
public:
  explicit UnionFind(int _n) : n(_n)
  {
    par.resize(static_cast<unsigned long>(_n));
    rank.resize(static_cast<unsigned long>(_n));
    sizes.resize(static_cast<unsigned long>(_n));
    for (int i = 0; i < _n; i++)
    {
      par[i] = i;
      rank[i] = 0;
      sizes[i] = 1;
    }
  }
  int find(int a)
  {
    if (par[a] == a)
      return a;
    return par[a] = find(par[a]);
  }
  bool same(int a, int b)
  {
    return find(a) == find(b);
  }
  void unite(int a, int b)
  {
    link(find(a), find(b));
  }
  int size(int a)
  {
    return sizes[find(a)];
  }
  void view()
  {
    for (int i = 0; i < n; i++)
    {
      cout << " par"
           << "[" << i << "]=" << par[i] << ((i == n - 1) ? "\n" : ",");
    }
    for (int i = 0; i < n; i++)
    {
      cout << "size"
           << "[" << i << "]=" << sizes[i] << ((i == n - 1) ? "\n" : ",");
    }
    cout << endl;
  }

private:
  void link(int a, int b)
  {
    if (same(a, b))
      return;
    if (rank[a] > rank[b])
    {
      par[b] = a;
      sizes[a] += sizes[b];
      sizes[b] = 0;
    }
    else
    {
      par[a] = b;
      if (rank[a] == rank[b])
        rank[b]++;
      sizes[b] += sizes[a];
      sizes[a] = 0;
    }
  }
  int n;
  vector<int> par;
  vector<int> rank;
  vector<int> sizes;
};

template<typename T>
class Dinic{
public:
  Dinic(int SIZE) :SIZE(SIZE) {
    v = vector<vector<tuple<int, T, int>>>(SIZE + 10);
  }

  void AddEdge(int from, int to, T cap) {
    v[from].push_back(make_tuple(to, cap, v[to].size()));
    v[to].push_back(make_tuple(from, 0, v[from].size() - 1));
  }

  void bfs(int s) {
    level = vector<int>(SIZE + 10, -1);
    queue<int> q;
    level[s] = 0;
    q.push(s);
    while (!q.empty()) {
      int now = q.front(); q.pop();
      for (int i = 0; i < (int)v[now].size(); i++) {
        int next, nextrv; T nextc; tie(next, nextc, nextrv) = v[now][i];
        if (nextc > 0 && level[next] < 0) {
          level[next] = level[now] + 1;
          q.push(next);
        }
      }
    }
  }

  T dfs(int now, int t, T f) {
    if (now == t) return f;
    for (int &i = iter[now]; i < (int)v[now].size(); i++) {
      int next, nextrv; T nextc; tie(next, nextc, nextrv) = v[now][i];
      if (nextc > 0 && level[now] < level[next]) {
        T d = dfs(next, t, min(f, nextc));
        if (d > 0) {
          get<1>(v[now][i]) -= d;
          get<1>(v[next][nextrv]) += d;
          return d;
        }
      }
    }
    return 0;
  }

  T max_flow(int s, int t) {
    T flow = 0;
    for (;;) {
      bfs(s);
      if (level[t] < 0) return flow;
      iter = vector<int>(SIZE + 10, 0);
      int f;
      while ((f = dfs(s, t, INT_MAX)) > 0) {
        flow += f;
      }
    }
  }
private:
  int SIZE;
  vector<vector<tuple<int, T, int>>> v;
  vector<int> level, iter;
};

template< typename flow_t, typename cost_t >
struct PrimalDual {
  const cost_t INF;

  struct edge {
    int to;
    flow_t cap;
    cost_t cost;
    int rev;
    bool isrev;
  };
  vector< vector< edge > > graph;
  vector< cost_t > potential, min_cost;
  vector< int > prevv, preve;

  PrimalDual(int V) : graph(V), INF(numeric_limits< cost_t >::max()) {}

  void add_edge(int from, int to, flow_t cap, cost_t cost) {
    graph[from].emplace_back((edge) {to, cap, cost, (int) graph[to].size(), false});
    graph[to].emplace_back((edge) {from, 0, -cost, (int) graph[from].size() - 1, true});
  }

  cost_t min_cost_flow(int s, int t, flow_t f) {
    int V = (int) graph.size();
    cost_t ret = 0;
    using Pi = pair< cost_t, int >;
    priority_queue< Pi, vector< Pi >, greater< Pi > > que;
    potential.assign(V, 0);
    preve.assign(V, -1);
    prevv.assign(V, -1);

    while(f > 0) {
      min_cost.assign(V, INF);
      que.emplace(0, s);
      min_cost[s] = 0;
      while(!que.empty()) {
        Pi p = que.top();
        que.pop();
        if(min_cost[p.second] < p.first) continue;
        for(int i = 0; i < graph[p.second].size(); i++) {
          edge &e = graph[p.second][i];
          cost_t nextCost = min_cost[p.second] + e.cost + potential[p.second] - potential[e.to];
          if(e.cap > 0 && min_cost[e.to] > nextCost) {
            min_cost[e.to] = nextCost;
            prevv[e.to] = p.second, preve[e.to] = i;
            que.emplace(min_cost[e.to], e.to);
          }
        }
      }
      if(min_cost[t] == INF) return -1;
      for(int v = 0; v < V; v++) potential[v] += min_cost[v];
      flow_t addflow = f;
      for(int v = t; v != s; v = prevv[v]) {
        addflow = min(addflow, graph[prevv[v]][preve[v]].cap);
      }
      f -= addflow;
      ret += addflow * potential[t];
      for(int v = t; v != s; v = prevv[v]) {
        edge &e = graph[prevv[v]][preve[v]];
        e.cap -= addflow;
        graph[v][e.rev].cap += addflow;
      }
    }
    return ret;
  }

  void output() {
    for(int i = 0; i < graph.size(); i++) {
      for(auto &e : graph[i]) {
        if(e.isrev) continue;
        auto &rev_e = graph[e.to][e.rev];
        cout << i << "->" << e.to << " (flow: " << rev_e.cap << "/" << rev_e.cap + e.cap << ")" << endl;
      }
    }
  }
};

template<typename T>
using Matrix = vector<vector<T>>;

template<typename T>
void init(int h, int w, Matrix<T> &matrix) {
    matrix = vector<vector<T>>(h, vector<T>(w));
    /*for (int i = 0; i < n; i++)
    {
        for (int j = 0; j < n; j++)
        {
            cin >> matrix[i][j];
        }
    }*/
}
template<typename T>
Matrix<T> &operator*(Matrix<T> &matrix, Matrix<T> &matrix1) {
    int height = (int)matrix.size();
    int width = (int)matrix[0].size();
    int height2 = (int)matrix1.size();
    int width2 = (int)matrix1[0].size();
    assert(width == height2);
    Matrix<T> res;
    res = vector<vector<T>>(height, vector<T>(width2, 0));
    for (int i = 0; i < height; i++) {
        for (int j = 0; j < width2; j++) {
            for (int k = 0; k < width; k++) {
                res[i][j] += matrix[i][k] * matrix1[k][j];
            }
        }
    }
    matrix.swap(res);
    return matrix;
}
template<typename T>
Matrix<T> Pow(Matrix<T> matrix, ll x) {
    Matrix<T> res = vector<vector<T>>(matrix.size(), vector<T>(matrix.size(), 0));
    for (int i = 0; i < (int)matrix.size(); i++) {
        res[i][i] = 1;
    }
    while (x) {
        if (x & 1) {
            res = res * matrix;
        }
        matrix = matrix * matrix;
        x >>= 1;
    }
    return res;
}


///////////////////////////////////////////////////////////////////////////////////
// MAIN
///////////////////////////////////////////////////////////////////////////////////
// 735134400 約数が1344個ある高度合成数(<= 1e9)
// 897612484786617600 約数が103680個ある高度合成数(<= 1e18)
// 苦手分野 重複組み合わせ
// LLの数値をつかう時は最後にLLをつける癖をつけよう
// Mintの数値にしたいときは、Mint(i)みたいにキャストしてね

namespace geo
{
#define IL inline
typedef double D;
typedef bool B;
typedef complex<D> P;
#define X real()
#define Y imag()
ostream &operator<<(ostream &out, const P &point1) // Pの出力
{
  out << '(' << point1.X << ", " << point1.Y << ')';
  return out;
}
B comp(const P &l, const P &r)
{
  return (l.X == r.X) ? l.Y < r.Y : l.X < r.X;
}
typedef pair<P, P> L;  //line
typedef pair<P, P> LS; //line segment
typedef pair<P, D> C;  //circle
typedef vector<P> Poly;
const D EPS = 1e-8;
//Decompotision Macro
#define DCl(a, b, l) \
  P a, b;            \
  tie(a, b) = l
#define DCc(a, b, c) \
  P a;               \
  D b;               \
  tie(a, b) = c

//A dot B
IL D dot(P a, P b) { return a.X * b.X + a.Y * b.Y; }
//A cross B
IL D cross(P a, P b) { return a.X * b.Y - a.Y * b.X; }
IL D ccw(P a, P b, P c)
{
  return cross(b - a, c - a);
}

// 未検証!!
IL P rotation(P a, double angle) {
  a *= exp(complex<double>(.0, angle * M_PI / 180.0));
  return a;
}

IL int sgn(P a, P b, P c)
{
  if (cross(b - a, c - a) > EPS)
    return 1;
  if (cross(b - a, c - a) < -EPS)
    return -1;
  if (dot(b - a, c - a) < -EPS)
    return 2; // b-a-c
  if (dot(a - b, c - b) < -EPS)
    return -2; // a-b-c
  return 0;
}
IL D signed_distance_P_L(P p, L l)
{
  DCl(a, b, l);
  return ccw(a, b, p) / abs(b - a);
}

IL D distance_P_L(P p, L l) // 点と直線の距離
{
  return abs(signed_distance_P_L(p, l));
}

IL D distance_P_P(P p1, P p2) // 二点間距離
{
  return abs(p1 - p2);
}

IL D distance_P_LS(P p, LS ls) // 点と線分の距離
{
  DCl(a, b, ls);
  if (dot(b - a, p - a) < EPS)
    return abs(p - a);
  if (dot(a - b, p - b) < EPS)
    return abs(p - b);
  return abs(cross(b - a, p - a)) / abs(b - a);
}

IL B is_L_L(L p, L q)
{
  DCl(a, b, p);
  DCl(c, d, q);
  return abs(cross(a - b, c - d)) > EPS ||
         abs(cross(a - b, d - b)) < EPS;
}
IL B is_L_LS(L l, LS ls)
{
  DCl(a, b, l);
  DCl(c, d, ls);
  return ccw(a, b, c) * ccw(a, b, d) < EPS;
}
IL B is_LS_LS(LS p, LS q)
{
  DCl(a, b, p);
  DCl(c, d, q);
  return sgn(a, b, c) * sgn(a, b, d) <= 0 && sgn(c, d, a) * sgn(c, d, b) <= 0;
}

IL D distance_LS_LS(LS p, LS q)
{
  if (is_LS_LS(p, q))
    return 0;
  return min({distance_P_LS(p.fi, q),
              distance_P_LS(p.se, q),
              distance_P_LS(q.fi, p),
              distance_P_LS(q.se, p)});
}

//caution: have to use isintersected_L_L before this function
IL P cross_L_L(L p, L q)
{
  DCl(a, b, p);
  D d1 = signed_distance_P_L(a, q);
  D d2 = signed_distance_P_L(b, q);
  return (a * d2 - b * d1) / (d2 - d1);
}

IL B is_C_C(C p, C q)
{
  DCc(cp, rp, p);
  DCc(cq, rq, q);
  D d = abs(cp - cq);
  return abs(rp - rq) - d < EPS &&
         rp + rq - d > -EPS;
}
IL LS cross_C_C(C p, C q)
{
  DCc(cp, rp, p);
  DCc(cq, rq, q);
  D d = abs(cp - cq);
  D rc = (d * d + rp * rp - rq * rq) / (2 * d);
  D rs = sqrt(rp * rp - rc * rc);
  P diff = (cq - cp) / d;
  return {cp + P(rc, rs) * diff, cp + P(rc, -rs) * diff};
}

int in_C(P q, C c) // 円での点の内外判定(円周上も判定)
{
  DCc(p, r, c);
  if (abs(r - abs(p - q)) < EPS)
    return 2; // 円周
  return r - abs(p - q) > -EPS ? 1 : 0;
}

}

int n;
int x[51], y[51];

int main() {
  using namespace geo;
  cin >> n;
  rep(i, 0, n) cin >> x[i] >> y[i];
  double res = 2000.;
  rep(i, 0, n) {
    rep(j, i + 1, n) {
      P p1(x[i], y[i]), p2(x[j], y[j]);
      C tmp1 = mp(p1, distance_P_P(p1, p2));
      C tmp2 = mp(p2, distance_P_P(p1, p2));
      L l1 = cross_C_C(tmp1, tmp2);
      L l2 = mp(p1, p2);
      P center = cross_L_L(l1, l2);
      C now = mp(center, distance_P_P(center, p1));
      bool ok = true;
      rep(k, 0, n) {
        if (in_C(P(x[k], y[k]), now) == 0) {
          ok = false;
        }
      }
      if (ok) {
        if (res - now.se > EPS) res = now.se;
      }
    }
  }
  double ans = res;
  show(res / 2.);
  rep(i, 0, n) {
    rep(j, i + 1, n) {
      rep(k, j + 1, n) {
        P p1(x[i], y[i]), p2(x[j], y[j]), p3(x[k], y[k]);
        double dd = distance_P_L(p3, L(p1, p2));
        if (dd <= EPS && -EPS <= dd) continue;
        C tmp1 = mp(p1, distance_P_P(p1, p2));
        C tmp2 = mp(p2, distance_P_P(p1, p2));
        L l1 = cross_C_C(tmp1, tmp2);
        C tmp3 = mp(p1, distance_P_P(p1, p3));
        C tmp4 = mp(p3, distance_P_P(p1, p3));
        L l2 = cross_C_C(tmp3, tmp4);
        P center = cross_L_L(l1, l2);
        C now = mp(center, distance_P_P(center, p1));
        bool ok = true;
        rep(m, 0, n) {
          if (in_C(P(x[m], y[m]), now) == 0) {
            ok = false;
            break;
          }
        }
        if (ok) {
          if (ans - now.se > EPS) ans = now.se;
        }
      }
    }
  }
  cout << ans << endl;
#ifdef LOCAL_DEFINE
  cerr << "Time elapsed: " << 1.0 * clock() / CLOCKS_PER_SEC << "s." << endl;
  show(MOD);
#endif
  return 0;
}

///////////////////////////////////////////////////////////////////////////////////
// NOTE
///////////////////////////////////////////////////////////////////////////////////
/*
*/

Measured memory_bytes: 256