Problem: There is a forest that we model as a plane and live n rare animals. Animal number i has its lair in the point (x_{i}, y_{i}). In order to protect them, a decision to build a nature reserve has been made.

The reserve must have a form of a circle containing all lairs. There is also a straight river flowing through the forest. All animals drink from this river, therefore it must have at least one common point with the reserve. On the other hand, ships constantly sail along the river, so the reserve must not have more than one common point with the river.

For convenience, scientists have made a transformation of coordinates so that the river is defined by y = 0. Check whether it is possible to build a reserve, and if possible, find the minimum possible radius of such a reserve.

Input

The first line contains one integer n (1 ≤ n ≤ 10^5) — the number of animals. 

Each of the next n lines contains two integers x_{i}, y_{i} (-10^7 ≤ x_{i}, y_{i} ≤ 10^7) — the coordinates of the i-th animal's lair. It is guaranteed that y_{i} ≠ 0. No two lairs coincide.

Output

If the reserve cannot be built, print -1. Otherwise print the minimum radius. Your answer will be accepted if absolute or relative error does not exceed 10^{-6}.

Formally, let your answer be a, and the jury's answer be b. Your answer is considered correct if \frac{|a - b|}{max{(1, |b|)}} ≤ 10^{-6}.

Examples

Input

1
0 1


Output

0.5

Input

3
0 1
0 2
0 -3


Output

-1


Input

2
0 1
1 1


Output

0.625

Note

In the first sample it is optimal to build the reserve with the radius equal to 0.5 and the center in (0,\ 0.5).

In the second sample it is impossible to build a reserve.

In the third sample it is optimal to build the reserve with the radius equal to 5/8 and the center in (1/2,\ 5/8).

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
template <typename T>
void _dbg(char const* s, const T& t) {
  cout << s << "=" << t << endl;
}
template <typename T, typename... TA>
void _dbg(char const* s, const T& t, const TA&... ta) {
  while (*s != ',') cout << *s++;
  cout << "=" << t << ",";
  _dbg(s + 1, ta...);
}
template <typename T1, typename T2>
ostream& operator<<(ostream& os, const pair<T1, T2>& p) {
  return os << "{" << p.first << "," << p.second << "}";
}
template <typename T>
ostream& operator<<(ostream& os, const vector<T>& v) {
  os << "{";
  for (auto it = v.begin(); it != v.end(); ++it) {
    if (it != v.begin()) os << ",";
    os << *it;
  }
  return os << "}";
}
template <typename T1, typename T2>
ostream& operator<<(ostream& os, const map<T1, T2>& mp) {
  os << "{";
  for (auto it = mp.begin(); it != mp.end(); ++it) {
    if (it != mp.begin()) os << ",";
    os << *it;
  }
  return os << "}";
}
template <typename T>
ostream& operator<<(ostream& os, const set<T>& s) {
  os << "{";
  for (auto it = s.begin(); it != s.end(); ++it) {
    if (it != s.begin()) os << ",";
    os << *it;
  }
  return os << "}";
}
template <typename T>
ostream& operator<<(ostream& os, const multiset<T>& s) {
  os << "{";
  for (auto it = s.begin(); it != s.end(); ++it) {
    if (it != s.begin()) os << ",";
    os << *it;
  }
  return os << "}";
}
template <size_t N>
ostream& operator<<(ostream& os, const bitset<N>& bs) {
  os << "{";
  for (size_t i = 0; i < N; ++i) {
    cout << bs[i];
  }
  return os << "}";
}
using ll = long long;
using ld = long double;
using pii = pair<int, int>;
using pll = pair<ll, ll>;
using pdd = pair<ld, ld>;
using namespace std;
const long double eps = 1e-9;
inline bool gt(const long double& x, const long double& y) { return x > y; }
inline bool lt(const long double& x, const long double& y) {
  return y > x + eps;
}
inline long double safe_sqrt(const long double& D) {
  return D < 0 ? 0 : sqrt(D);
}
struct pt {
  long double x, y;
  pt() {}
  pt(long double a, long double b) : x(a), y(b) {}
};
const int N = 1e5 + 5;
const int STEPS = 150;
int n;
pt p[N];
inline bool can(long double R) {
  long double l = -1e16 - 1, r = 1e16 + 1;
  for (int i = 0; i < n; i++) {
    long double b = -2 * p[i].x;
    long double c = p[i].x * p[i].x + p[i].y * p[i].y - 2 * p[i].y * R;
    long double D = b * b - 4 * c;
    if (D < 0) return false;
    D = safe_sqrt(D);
    long double x1 = p[i].x - D / 2, x2 = p[i].x + D / 2;
    l = max(l, x1);
    r = min(r, x2);
  }
  return !gt(l, r);
}
void solve() {
  cin >> n;
  bool pos = 0, neg = 0;
  for (int i = 0; i < n; ++i) {
    cin >> p[i].x >> p[i].y;
    if (p[i].y > 0)
      pos = 1;
    else
      neg = 1;
  }
  if (pos && neg) {
    cout << -1 << endl;
    return;
  }
  if (neg) {
    for (int i = 0; i < n; ++i) p[i].y = -p[i].y;
  }
  long double L = 0, R = 1e16;
  for (int step = 0; step < STEPS; step++) {
    long double mid = L + (R - L) / 2.0;
    if (can(mid))
      R = mid;
    else
      L = mid;
  }
  cout << setprecision(16) << fixed << L + (R - L) / 2 << endl;
}
int main() {
  solve();
  return 0;
}