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;
struct Point {
  long double x, y;
  Point operator~() const { return {y, -x}; }
  bool operator<(const Point &P) const {
    return y > P.y || (y == P.y && x < P.x);
  }
  long double operator*(const Point &P) const { return x * P.x + y * P.y; }
  bool operator==(const Point &P) const { return (x == P.x) && (y == P.y); }
  bool operator!=(const Point &P) const { return (x != P.x) || (y != P.y); }
  long double operator^(const Point &P) const { return x * P.y - y * P.x; }
  Point operator+(const Point &P) const { return {x + P.x, y + P.y}; }
  Point operator-(const Point &P) const { return {x - P.x, y - P.y}; }
  Point operator-() const { return {-x, -y}; }
  Point operator/(const long double d) const { return {x / d, y / d}; }
  Point operator*(const long double d) const { return {x * d, y * d}; }
};
long double CK(const Point &A, const Point &B, const Point &C) {
  return ((C - B) ^ (B - A));
}
int N, s;
vector<Point> P, H;
long double M = 50000000000001, E = 0.0000001, two = 2, nn = -1, zero = 0;
bool check(const long double &R) {
  if (H[0].y > two * R) return false;
  long double x1 = H[0].x - sqrtl(H[0].y * (two * R - H[0].y));
  long double x2 = H[0].x + sqrtl(H[0].y * (two * R - H[0].y));
  for (int i = 1; i < H.size() && x1 <= x2; i++) {
    x1 = max(x1, H[i].x - sqrtl(H[i].y * (two * R - H[i].y)));
    x2 = min(x2, H[i].x + sqrtl(H[i].y * (two * R - H[i].y)));
  }
  return x1 <= x2;
}
long double get() {
  long double Ra = 0, Rb = M;
  int C = 0;
  while (Rb - Ra > E && C++ < 10000) {
    if (check((Ra + Rb) / two))
      Rb = (Ra + Rb) / two;
    else
      Ra = (Ra + Rb) / two;
  }
  return (Ra + Rb) / two;
}
int main() {
  ios::sync_with_stdio(false);
  cin.tie(0);
  cin >> N;
  P.resize(N);
  int positive = 0, negative = 0;
  for (int i = 0; i < N; i++) {
    cin >> P[i].x >> P[i].y;
    if (P[i].y > 0)
      positive++;
    else
      negative++;
  }
  if (positive != 0 && negative != 0) {
    cout << -1 << endl;
    return 0;
  }
  if (negative > 0)
    for (int i = 0; i < N; i++) P[i].y *= nn;
  sort(P.begin(), P.end());
  for (int i = 0; i < N; i++) {
    while (H.size() > 1 && CK(H[H.size() - 2], H[H.size() - 1], P[i]) <= zero)
      H.pop_back();
    H.push_back(P[i]);
  }
  s = H.size();
  for (int i = N - 1; i >= 0; i--) {
    while (H.size() > s && CK(H[H.size() - 2], H[H.size() - 1], P[i]) <= zero)
      H.pop_back();
    H.push_back(P[i]);
  }
  H.pop_back();
  printf("%.20Lf\n", get());
  return 0;
}