Problem: A prisoner wants to escape from a prison. The prison is represented by the interior of the convex polygon with vertices P_1, P_2, P_3, …, P_{n+1}, P_{n+2}, P_{n+3}. It holds P_1=(0,0), P_{n+1}=(0, h), P_{n+2}=(-10^{18}, h) and P_{n+3}=(-10^{18}, 0).

<image>

The prison walls P_{n+1}P_{n+2}, P_{n+2}P_{n+3} and P_{n+3}P_1 are very high and the prisoner is not able to climb them. Hence his only chance is to reach a point on one of the walls P_1P_2, P_2P_3,..., P_{n}P_{n+1} and escape from there. On the perimeter of the prison, there are two guards. The prisoner moves at speed 1 while the guards move, remaining always on the perimeter of the prison, with speed v.

If the prisoner reaches a point of the perimeter where there is a guard, the guard kills the prisoner. If the prisoner reaches a point of the part of the perimeter he is able to climb and there is no guard there, he escapes immediately. Initially the prisoner is at the point (-10^{17}, h/2) and the guards are at P_1. 

Find the minimum speed v such that the guards can guarantee that the prisoner will not escape (assuming that both the prisoner and the guards move optimally).

Notes:

  * At any moment, the guards and the prisoner can see each other. 
  * The "climbing part" of the escape takes no time. 
  * You may assume that both the prisoner and the guards can change direction and velocity instantly and that they both have perfect reflexes (so they can react instantly to whatever the other one is doing). 
  * The two guards can plan ahead how to react to the prisoner movements. 

Input

The first line of the input contains n (1 ≤ n ≤ 50).

The following n+1 lines describe P_1, P_2,..., P_{n+1}. The i-th of such lines contain two integers x_i, y_i (0≤ x_i, y_i≤ 1,000) — the coordinates of P_i=(x_i, y_i).

It is guaranteed that P_1=(0,0) and x_{n+1}=0. The polygon with vertices P_1,P_2,..., P_{n+1}, P_{n+2}, P_{n+3} (where P_{n+2}, P_{n+3} shall be constructed as described in the statement) is guaranteed to be convex and such that there is no line containing three of its vertices.

Output

Print a single real number, the minimum speed v that allows the guards to guarantee that the prisoner will not escape. Your answer will be considered correct if its relative or absolute error does not exceed 10^{-6}.

Examples

Input


2
0 0
223 464
0 749


Output


1


Input


3
0 0
2 2
2 4
0 6


Output


1.0823922


Input


4
0 0
7 3
7 4
5 7
0 8


Output


1.130309669


Input


5
0 0
562 248
460 610
281 702
206 723
0 746


Output


1.148649561


Input


7
0 0
412 36
745 180
747 184
746 268
611 359
213 441
0 450


Output


1.134745994

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using ld = long double;
using ul = std::size_t;
template <typename T, ul C>
using ar = std::array<T, C>;
const int MN = 55;
const ld INF = 1e10;
const ld EPS = 1e-9;
template <ul C>
ar<ld, C>& operator*=(ar<ld, C>& a, const ld& v) {
  for (ul i = 0; i < C; ++i) a[i] *= v;
  return a;
}
template <ul C>
ar<ld, C> operator*(ar<ld, C> a, const ld& v) {
  return a *= v;
}
template <ul C>
ar<ld, C>& operator+=(ar<ld, C>& a, const ar<ld, C>& b) {
  for (ul i = 0; i < C; ++i) a[i] += b[i];
  return a;
}
template <ul C>
ar<ld, C> operator+(ar<ld, C> a, const ar<ld, C>& b) {
  return a += b;
}
inline bool equ(const ld& a, const ld& b) { return std::abs(a - b) < EPS; }
struct vec {
 public:
  ld x, y;
  vec(ld _x = 0, ld _y = 0) : x(_x), y(_y) {}
  vec& operator+=(const vec& o) { return x += o.x, y += o.y, *this; }
  vec& operator-=(const vec& o) { return x -= o.x, y -= o.y, *this; }
  vec& operator*=(const ld& o) { return x *= o, y *= o, *this; }
  friend vec operator+(vec a, const vec& b) { return a += b; }
  friend vec operator-(vec a, const vec& b) { return a -= b; }
  friend vec operator*(vec a, const ld& b) { return a *= b; }
  ld mag() const { return hypot(x, y); }
};
struct range {
 public:
  ld l, r;
};
using ranges = std::vector<range>;
ranges unite(const ranges& a, const ranges& b) {
  ranges f(0);
  bool rdy = 0;
  ld l = -INF - 1, r = -INF - 1;
  auto add = [&]() {
    if (rdy)
      f.push_back({l, r});
    else
      rdy = 1;
  };
  range t;
  for (int i = 0, j = 0; i < a.size() || j < b.size();) {
    if (j < b.size() && (i == a.size() || b[j].l < a[i].l))
      t = b[j++];
    else
      t = a[i++];
    if (r + EPS < t.l) add(), l = t.l;
    r = std::max(r, t.r);
  }
  add();
  return f;
}
ranges isect(const ranges& a, const ranges& b) {
  ranges f;
  for (int i = 0, j = 0; i < a.size() && j < b.size();) {
    ld l = std::max(a[i].l, b[j].l);
    ld r = std::min(a[i].r, b[j].r);
    if (l + EPS < r) f.push_back({l, r});
    if (a[i].r < b[j].r)
      ++i;
    else
      ++j;
  }
  return f;
}
int N;
vec p[MN];
ld pf[MN];
ar<ld, 6> sqr(const ar<ld, 3>& a) {
  return {a[0] * a[0], 2 * a[0] * a[1], 2 * a[0] * a[2],
          a[1] * a[1], a[2] * a[2],     2 * a[1] * a[2]};
}
ranges neg(const ar<ld, 3>& a) {
  if (equ(a[2], 0)) {
    if (equ(a[1], 0))
      return a[0] < 0 ? ranges{{-INF, INF}} : ranges();
    else
      return a[1] < 0 ? ranges{{-a[0] / a[1], INF}}
                      : ranges{{-INF, -a[0] / a[1]}};
  } else if (a[2] < 0) {
    ld d = a[1] * a[1] - 4 * a[0] * a[2];
    if (d > EPS) {
      d = sqrt(d) / (2 * -a[2]);
      ld c = -a[1] / (2 * a[2]);
      return ranges{{-INF, c - d}, {c + d, INF}};
    }
    return ranges{{-INF, INF}};
  } else {
    ld d = a[1] * a[1] - 4 * a[0] * a[2];
    if (d > EPS) {
      d = sqrt(d) / (2 * a[2]);
      ld c = -a[1] / (2 * a[2]);
      return ranges{{c - d, c + d}};
    }
    return ranges();
  }
}
ranges eval(const ar<ld, 6>& a) {
  ranges f;
  f = neg({a[0], a[1], a[3]});
  f = unite(f, neg({a[0] + a[2] + a[4], a[1] + a[5], a[3]}));
  if (a[4] > EPS) {
    if (equ(a[5], 0)) {
      ld yv = a[2] / (-2 * a[4]);
      if (0 < yv && yv < 1)
        f = unite(f, neg({a[0] + a[2] * yv + a[4] * yv * yv, a[1], a[3]}));
    } else {
      range xr = {(-2 * a[4] - a[2]) / a[5], -a[2] / a[5]};
      if (a[5] < 0) xr = {xr.r, xr.l};
      f = unite(f, isect(ranges{xr}, neg({a[0] - a[2] * a[2] / (4 * a[4]),
                                          a[1] - 2 * a[2] * a[5] / (4 * a[4]),
                                          a[3] - a[5] * a[5] / (4 * a[4])})));
    }
  }
  return isect(f, ranges{{0, 1}});
}
bool ok(const ld& v) {
  ranges bot, top;
  for (int i = 0; i < N; ++i)
    for (int j = i + 1; j < N; ++j) {
      ld mi = (p[i + 1] - p[i]).mag();
      ld mj = (p[j + 1] - p[j]).mag();
      ar<ld, 6> L = sqr({pf[j] - pf[i], -mi, mj});
      ar<ld, 6> Dx =
          sqr({p[j].x - p[i].x, -p[i + 1].x + p[i].x, p[j + 1].x - p[j].x});
      ar<ld, 6> Dy =
          sqr({p[j].y - p[i].y, -p[i + 1].y + p[i].y, p[j + 1].y - p[j].y});
      L *= -1;
      Dx *= v * v;
      Dy *= v * v;
      ar<ld, 6> fin = L + Dx + Dy;
      ranges bot_ = eval(fin);
      std::swap(fin[1], fin[2]);
      std::swap(fin[3], fin[4]);
      ranges top_ = eval(fin);
      for (range& x : bot_) x.l = pf[i] + mi * x.l, x.r = pf[i] + mi * x.r;
      for (range& x : top_) x.l = pf[j] + mj * x.l, x.r = pf[j] + mj * x.r;
      bot = unite(bot, bot_);
      top = unite(top, top_);
    }
  return isect(bot, top).empty();
}
int main() {
  scanf("%d", &N);
  for (int i = 0; i <= N; ++i) {
    int u, v;
    scanf("%d%d", &u, &v);
    p[i] = {(ld)u, (ld)v};
  }
  for (int i = 0; i < N; ++i) pf[i + 1] = pf[i] + (p[i + 1] - p[i]).mag();
  ld l = 1, r = 5000;
  for (int i = 0; i < 100; ++i) {
    ld m = l + (r - l) / 2;
    if (ok(m))
      r = m;
    else
      l = m;
  }
  printf("%.9lf\n", (double)r);
  return 0;
}