The cool breeze blows gently, the flowing water ripples steadily.

"Flowing and passing like this, the water isn't gone ultimately; Waxing and waning like that, the moon doesn't shrink or grow eventually."

"Everything is transient in a way and perennial in another."

Kanno doesn't seem to make much sense out of Mino's isolated words, but maybe it's time that they enjoy the gentle breeze and the night sky — the inexhaustible gifts from nature.

Gazing into the sky of stars, Kanno indulges in a night's tranquil dreams. 

There is a set S of n points on a coordinate plane.

Kanno starts from a point P that can be chosen on the plane. P is not added to S if it doesn't belong to S. Then the following sequence of operations (altogether called a move) is repeated several times, in the given order:

  1. Choose a line l such that it passes through at least two elements in S and passes through Kanno's current position. If there are multiple such lines, one is chosen equiprobably. 
  2. Move to one of the points that belong to S and lie on l. The destination is chosen equiprobably among all possible ones, including Kanno's current position (if it does belong to S). 



There are q queries each consisting of two integers (t_i, m_i). For each query, you're to help Kanno maximize the probability of the stopping position being the t_i-th element in S after m_i moves with a proper selection of P, and output this maximum probability. Note that according to rule 1, P should belong to at least one line that passes through at least two points from S.

Input

The first line contains a positive integer n (2 ≤ n ≤ 200) — the number of points in S.

The i-th of the following n lines contains two space-separated integers x_i and y_i (-10^4 ≤ x_i, y_i ≤ 10^4) — the coordinates of the i-th point in S. The input guarantees that for all 1 ≤ i < j ≤ n, (x_i, y_i) ≠ (x_j, y_j) holds.

The next line contains a positive integer q (1 ≤ q ≤ 200) — the number of queries.

Each of the following q lines contains two space-separated integers t and m (1 ≤ t_i ≤ n, 1 ≤ m_i ≤ 10^4) — the index of the target point and the number of moves, respectively.

Output

Output q lines each containing a decimal number — the i-th among them denotes the maximum probability of staying on the t_i-th point after m_i steps, with a proper choice of starting position P.

Your answer will be considered correct if each number in your output differs from the corresponding one in jury's answer by at most 10^{-6}.

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

Example

Input

5
0 0
1 3
2 2
3 1
4 4
10
1 1
2 1
3 1
4 1
5 1
3 2
3 3
3 4
3 5
3 6


Output

0.50000000000000000000
0.50000000000000000000
0.33333333333333331483
0.50000000000000000000
0.50000000000000000000
0.18518518518518517491
0.15226337448559670862
0.14494741655235482414
0.14332164812274550414
0.14296036624949901017

Note

The points in S and possible candidates for line l are depicted in the following figure.

<image>

For the first query, when P = (-1, -3), l is uniquely determined to be 3x = y, and thus Kanno will move to (0, 0) with a probability of \frac 1 2.

For the third query, when P = (2, 2), l is chosen equiprobably between x + y = 4 and x = y. Kanno will then move to the other four points with a probability of \frac 1 2 ⋅ \frac 1 3 = \frac 1 6 each, or stay at (2, 2) with a probability of \frac 1 3.
#include <bits/stdc++.h>
using namespace std;
void readi(int &x) {
  int v = 0, f = 1;
  char c = getchar();
  while (!isdigit(c) && c != '-') c = getchar();
  if (c == '-')
    f = -1;
  else
    v = v * 10 + c - '0';
  while (isdigit(c = getchar())) v = v * 10 + c - '0';
  x = v * f;
}
void readll(long long &x) {
  long long v = 0ll, f = 1ll;
  char c = getchar();
  while (!isdigit(c) && c != '-') c = getchar();
  if (c == '-')
    f = -1;
  else
    v = v * 10 + c - '0';
  while (isdigit(c = getchar())) v = v * 10 + c - '0';
  x = v * f;
}
void readc(char &x) {
  char c;
  while ((c = getchar()) == ' ')
    ;
  x = c;
}
void writes(string s) { puts(s.c_str()); }
void writeln() { writes(""); }
void writei(int x) {
  if (x < 0) {
    putchar('-');
    x = abs(x);
  }
  if (!x) putchar('0');
  char a[25];
  int top = 0;
  while (x) {
    a[++top] = (x % 10) + '0';
    x /= 10;
  }
  while (top) {
    putchar(a[top]);
    top--;
  }
}
void writell(long long x) {
  if (x < 0) {
    putchar('-');
    x = abs(x);
  }
  if (!x) putchar('0');
  char a[25];
  int top = 0;
  while (x) {
    a[++top] = (x % 10) + '0';
    x /= 10;
  }
  while (top) {
    putchar(a[top]);
    top--;
  }
}
inline long long inc(int &x) { return ++x; }
inline long long inc(long long &x) { return ++x; }
inline long long inc(int &x, long long y) { return x += y; }
inline long long inc(long long &x, long long y) { return x += y; }
inline double inc(double &x, double y) { return x += y; }
inline long long dec(int &x) { return --x; }
inline long long dec(long long &x) { return --x; }
inline long long dec(int &x, long long y) { return x -= y; }
inline long long dec(long long &x, long long y) { return x -= y; }
inline double dec(double &x, double y) { return x -= y; }
inline long long mul(int &x) { return x = ((long long)x) * x; }
inline long long mul(long long &x) { return x = x * x; }
inline long long mul(int &x, long long y) { return x *= y; }
inline long long mul(long long &x, long long y) { return x *= y; }
inline double mul(double &x, double y) { return x *= y; }
inline long long divi(const int &x) {
  long long ans, l, r, mid;
  ans = 0;
  l = 0;
  r = 0x3fffffff;
  while (l < r) {
    mid = (l + r) / 2;
    if (mid * mid <= x) {
      ans = mid;
      l = mid + 1;
    } else
      r = mid;
  }
  return ans;
}
inline long long divi(const long long &x) {
  long long ans, l, r, mid;
  ans = 0;
  l = 0;
  r = 0x3fffffff;
  while (l < r) {
    mid = (l + r) / 2;
    if (mid * mid <= x) {
      ans = mid;
      l = mid + 1;
    } else
      r = mid;
  }
  return ans;
}
inline long long divi(int &x, long long y) { return x /= y; }
inline long long divi(long long &x, long long y) { return x /= y; }
inline double divi(double &x, double y) { return x /= y; }
inline long long mod(int &x, long long y) { return x %= y; }
inline long long mod(long long &x, long long y) { return x %= y; }
struct matrix {
  double a[205][205];
} f[17];
long long n, m, i, j, k, px[205], py[205], x, y;
matrix mul(matrix x, matrix y) {
  int i, j, k;
  matrix ans;
  memset((ans.a), (0), (sizeof((ans.a))));
  if ((1) <= ((n)))
    for (((i)) = (1); ((i)) <= ((n)); ((i))++)
      if ((1) <= ((n)))
        for (((j)) = (1); ((j)) <= ((n)); ((j))++)
          if ((1) <= ((n)))
            for (((k)) = (1); ((k)) <= ((n)); ((k))++)
              ans.a[i][j] += x.a[i][k] * y.a[k][j];
  return ans;
}
double cnt[205], g[2][205];
vector<vector<long long> > lne;
void ind() {
  cin >> n;
  if ((1) <= ((n)))
    for (((i)) = (1); ((i)) <= ((n)); ((i))++) cin >> px[i] >> py[i];
  if ((1) <= ((n)))
    for (((i)) = (1); ((i)) <= ((n)); ((i))++)
      if ((1) <= ((i - 1)))
        for (((j)) = (1); ((j)) <= ((i - 1)); ((j))++) {
          long long dx = px[i] - px[j], dy = py[i] - py[j];
          vector<long long> tis;
          if ((1) <= ((n)))
            for (((k)) = (1); ((k)) <= ((n)); ((k))++) {
              if ((px[i] - px[k]) * dy == (py[i] - py[k]) * dx) {
                tis.push_back(k);
              }
            }
          if (tis.size() >= 2) lne.push_back(tis);
        }
  stable_sort(lne.begin(), lne.end());
  lne.resize(unique(lne.begin(), lne.end()) - lne.begin());
  for (__typeof((lne).begin()) it = (lne).begin(); it != (lne).end(); it++) {
    for (__typeof((*it).begin()) it2 = (*it).begin(); it2 != (*it).end(); it2++)
      cnt[*it2] += 1.0;
  }
  for (__typeof((lne).begin()) it = (lne).begin(); it != (lne).end(); it++) {
    for (__typeof((*it).begin()) it2 = (*it).begin(); it2 != (*it).end();
         it2++) {
      for (__typeof((*it).begin()) it3 = (*it).begin(); it3 != (*it).end();
           it3++) {
        f[0].a[*it2][*it3] += 1.0 / cnt[*it2] / it->size();
      }
    }
  }
}
int main() {
  ios_base::sync_with_stdio(false);
  ;
  ind();
  if ((1) <= ((15)))
    for (((i)) = (1); ((i)) <= ((15)); ((i))++) f[i] = mul(f[i - 1], f[i - 1]);
  cin >> m;
  while (m--) {
    cin >> y >> x;
    memset((g), (0), (sizeof((g))));
    g[0][y] = 1;
    x--;
    int z = 0;
    if ((0) <= (15))
      for ((i) = (0); (i) <= (15); (i)++) {
        if ((x >> i) & 1) {
          z ^= 1;
          memset((g[z]), (0), (sizeof((g[z]))));
          if ((1) <= ((n)))
            for (((j)) = (1); ((j)) <= ((n)); ((j))++)
              if ((1) <= ((n)))
                for (((k)) = (1); ((k)) <= ((n)); ((k))++)
                  g[z][j] += g[!z][k] * f[i].a[j][k];
        }
      }
    double res = 0;
    for (__typeof((lne).begin()) it = (lne).begin(); it != (lne).end(); it++) {
      double sum = 0;
      for (__typeof((*it).begin()) it2 = (*it).begin(); it2 != (*it).end();
           it2++) {
        sum += g[z][*it2];
      }
      sum /= it->size();
      res = max(res, sum);
    }
    printf("%.50f\n", res);
  }
  return 0;
}
