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;
const int INFL = (int)1e9;
const long long int INFLL = (long long int)1e18;
const double INFD = numeric_limits<double>::infinity();
const double PI = 3.14159265358979323846;
bool nearlyeq(double x, double y) { return abs(x - y) < 1e-9; }
bool inrange(int x, int t) { return x >= 0 && x < t; }
long long int rndf(double x) {
  return (long long int)(x + (x >= 0 ? 0.5 : -0.5));
}
long long int floorsqrt(double x) {
  long long int m = (long long int)sqrt(x);
  return m + (m * m <= (long long int)(x) ? 0 : -1);
}
long long int ceilsqrt(double x) {
  long long int m = (long long int)sqrt(x);
  return m + ((long long int)x <= m * m ? 0 : 1);
}
long long int rnddiv(long long int a, long long int b) {
  return (a / b + (a % b * 2 >= b ? 1 : 0));
}
long long int ceildiv(long long int a, long long int b) {
  return (a / b + (a % b == 0 ? 0 : 1));
}
long long int gcd(long long int m, long long int n) {
  if (n == 0)
    return m;
  else
    return gcd(n, m % n);
}
template <typename val_t, typename bsargv_t>
val_t upper_binary_search(val_t l, val_t r, const bsargv_t &bsargv,
                          bool (*evalfunc)(val_t, const bsargv_t &)) {
  if (r - l <= 1e-7 * max(1.0, l)) {
    if (evalfunc(l, bsargv))
      return l;
    else
      return r;
  }
  val_t m = (l + r) / 2;
  if (evalfunc(m, bsargv))
    upper_binary_search<val_t, bsargv_t>(l, m, bsargv, evalfunc);
  else
    return upper_binary_search<val_t, bsargv_t>(m, r, bsargv, evalfunc);
}
struct bsargv_t {
  vector<double> xs;
  vector<double> ys;
};
bool evalfunc(double val, const bsargv_t &bsargv) {
  double max_l = -INFD, min_r = INFD;
  for (int i = 0; i < (int)bsargv.xs.size(); i++) {
    double a = val * 2 - abs(bsargv.ys[i]);
    if (a < 0) return false;
    double b = sqrt(abs(bsargv.ys[i])) * sqrt(abs(a));
    double v0 = bsargv.xs[i] - b;
    double v1 = bsargv.xs[i] + b;
    max_l = max(max_l, v0);
    min_r = min(min_r, v1);
  }
  if (max_l <= min_r)
    return true;
  else
    return false;
}
int main() {
  int n;
  cin >> n;
  vector<double> xs(n), ys(n);
  double max_y = -INFD, min_y = INFD;
  for (int i = 0; i < (int)n; i++) {
    cin >> xs[i] >> ys[i];
    max_y = max(max_y, ys[i]);
    min_y = min(min_y, ys[i]);
  }
  if (min_y * max_y < 0) {
    cout << -1 << endl;
  } else {
    bsargv_t bsargv = {xs, ys};
    double ans =
        upper_binary_search<double, bsargv_t>(0.5, 1e15, bsargv, evalfunc);
    cout << setprecision(12) << ans << endl;
  }
}