Problem: The commanding officers decided to drop a nuclear bomb on the enemy's forces. You are ordered to determine the power of the warhead that needs to be used.

The enemy has N strategically important objects. Their positions are known due to the intelligence service. The aim of the strike is to deactivate at least K important objects of the enemy. The bombing impact point is already determined and has coordinates of [X0; Y0].

The nuclear warhead is marked by the estimated impact radius R ≥ 0. All the buildings that are located closer than R to the bombing epicentre will be destroyed. All the buildings that are located further than R from the epicentre, can also be deactivated with some degree of probability. Let's assume that D is the distance between a building and the epicentre. This building's deactivation probability P(D, R) is calculated according to the following formula: 

<image> We should regard <image> as ea, where e ≈ 2.7182818284590452353602874713527

If the estimated impact radius of the warhead is equal to zero, then all the buildings located in the impact point will be completely demolished and all the rest of important objects will not be damaged.

The commanding officers want the probability of failing the task to be no more than ε. Nuclear warheads are too expensive a luxury, that's why you have to minimise the estimated impact radius of the warhead. 

Input

The first line contains an integer N which represents the number of the enemy's objects (1 ≤ N ≤ 100). The second line contains two integers: K is the required number of deactivated objects, and ε is the maximally permitted probability of not completing the task, given in per mils (1 ≤ K ≤ N, 1 ≤ ε ≤ 999). The third line contains X0 and Y0 which are the coordinates of the strike impact point. The next N lines contain two numbers Xi and Yi each which are the coordinates of every strategically important object. All the coordinates are integer, their absolute values do not exceed 1000.

Let us remind you that there are a thousand per mils in unity (number one).

There can be several objects in one point.

Output

Print the sought estimated impact radius of the warhead. The absolute or relative measure of the inaccuracy of your answer should not exceed 10 - 6.

Examples

Input

1
1 500
5 5
1 2


Output

3.84257761518762740


Input

5
3 100
0 0
3 4
60 70
100 100
10 10
5 12


Output

13.45126176453737600

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
#pragma warning(error : 4700)
#pragma warning(error : 4715)
#pragma warning(error : 4716)
using namespace std;
const long long MILLER_RABIN = 3215031751;
const long double PI = 3.14159265358979323846;
inline int lg2(long long n) {
  int h = 0;
  while ((1ll << h) < n) ++h;
  return h;
}
struct config_io {
  config_io() {
    cin.tie(nullptr);
    ios_base::sync_with_stdio(false);
  }
} cnf_io;
struct config_rand {
  config_rand() {
    srand(chrono::duration_cast<chrono::nanoseconds>(
              chrono::high_resolution_clock::now().time_since_epoch())
              .count());
  }
} cnf_rand;
namespace std {
template <>
struct hash<pair<int, int> > {
  size_t operator()(const pair<int, int> &x) const {
    return (x.first * 31 + x.second) % ((int)1e9 + 7);
  }
};
}  // namespace std
long double dist(int x1, int y1, int x2, int y2) {
  return sqrt((x1 - x2) * (x1 - x2) + (y1 - y2) * (y1 - y2));
}
long double P(long double R, long double D) {
  if (D <= R) return 1;
  return exp(1 - (D * D) / (R * R));
}
long double calc(long double r, int k, vector<int> &x, vector<int> &y) {
  int n = x.size() - 1;
  long double dp[n + 1][n + 1];
  for (int i = 0; i <= n; ++i) {
    for (int j = 0; j <= n; ++j) {
      dp[i][j] = 0;
    }
  }
  dp[0][0] = 1;
  for (int i = 1; i <= n; ++i) {
    for (int j = 0; j <= i; ++j) {
      dp[i][j] =
          dp[i - 1][j] * (1 - P(r, dist(x[0], y[0], x[i], y[i]))) +
          (j > 0 ? dp[i - 1][j - 1] : 0) * P(r, dist(x[0], y[0], x[i], y[i]));
    }
  }
  long double result = 0;
  for (int i = k; i <= n; ++i) {
    result += dp[n][i];
  }
  return result;
}
int main() {
  int n, k;
  long double eps;
  cin >> n >> k >> eps;
  vector<int> x(n + 1), y(n + 1);
  for (int i = 0; i <= n; ++i) {
    cin >> x[i] >> y[i];
  }
  long double l = 0, r = 1e4;
  const long double EPS = 1e-6;
  while (r - l > EPS) {
    long double mid = (l + r) / 2;
    long double prob = calc(mid, k, x, y);
    if (prob > 1 - eps / 1000) {
      r = mid;
    } else {
      l = mid;
    }
  }
  cout << fixed << setprecision(8);
  cout << r << endl;
  return 0;
}