Problem: 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.

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
const int DEBUG = 0;
using namespace std;
long long int gcd(long long int x, long long int y) {
  while (y != 0) {
    long long int r = x % y;
    x = y;
    y = r;
  }
  return x;
}
struct line {
  long long int a, b, c;
  line(long long int a, long long int b, long long int c) : a(a), b(b), c(c) {}
  static line cross(pair<long long int, long long int> x,
                    pair<long long int, long long int> y) {
    long long int a = y.second - x.second;
    long long int b = x.first - y.first;
    long long int c = a * x.first + b * x.second;
    long long int g = a;
    g = gcd(g, b);
    g = gcd(g, c);
    a /= g;
    b /= g;
    c /= g;
    if (b < 0) {
      a = -a;
      b = -b;
      c = -c;
    }
    if (b == 0 && a < 0) {
      a = -a;
      c = -c;
    }
    return line(a, b, c);
  }
  bool operator<(const line &other) const {
    if (a != other.a) return a < other.a;
    return pair<long long int, long long int>(b, c) <
           pair<long long int, long long int>(other.b, other.c);
  }
  bool operator==(const line &other) const {
    return a == other.a && b == other.b && c == other.c;
  }
  bool inhabit(pair<long long int, long long int> p) const {
    return a * p.first + b * p.second == c;
  }
};
ostream &operator<<(ostream &os, const line &l) {
  return os << l.a << "x + " << l.b << "y = " << l.c << endl;
}
vector<vector<double> > mul(const vector<vector<double> > &a,
                            const vector<vector<double> > &b) {
  int n = a.size();
  int m = b.size();
  int l = b[0].size();
  vector<vector<double> > ret(n, vector<double>(l, 0));
  for (int i = (int)(0); i < (int)(n); i++) {
    for (int j = (int)(0); j < (int)(m); j++) {
      for (int k = (int)(0); k < (int)(l); k++) {
        ret[i][k] += a[i][j] * b[j][k];
      }
    }
  }
  return ret;
}
int main(void) {
  int n;
  scanf("%d", &n);
  vector<int> x(n), y(n);
  for (int i = (int)(0); i < (int)(n); i++) scanf("%d%d", &x[i], &y[i]);
  vector<line> sigma;
  for (int i = (int)(0); i < (int)(n); i++) {
    for (int j = (int)(i + 1); j < (int)(n); j++) {
      sigma.push_back(
          line::cross(pair<long long int, long long int>(x[i], y[i]),
                      pair<long long int, long long int>(x[j], y[j])));
    }
  }
  sort(sigma.begin(), sigma.end());
  sigma.erase(unique(sigma.begin(), sigma.end()), sigma.end());
  if (DEBUG) {
    for (int i = (int)(0); i < (int)(sigma.size()); i++) {
      cerr << sigma[i];
    }
  }
  int e = sigma.size();
  vector<vector<double> > a(n, vector<double>(n));
  vector<vector<double> > b(n, vector<double>(e));
  vector<vector<double> > c(e, vector<double>(n));
  for (int i = (int)(0); i < (int)(n); i++) {
    int sum = 0;
    for (int j = (int)(0); j < (int)(e); j++) {
      if (sigma[j].inhabit(pair<long long int, long long int>(x[i], y[i]))) {
        sum += 1;
      }
    }
    for (int j = (int)(0); j < (int)(e); j++) {
      if (sigma[j].inhabit(pair<long long int, long long int>(x[i], y[i]))) {
        b[i][j] = 1.0 / sum;
      }
    }
  }
  for (int j = (int)(0); j < (int)(e); j++) {
    int u = 0;
    for (int k = (int)(0); k < (int)(n); k++) {
      if (sigma[j].inhabit(pair<long long int, long long int>(x[k], y[k]))) {
        u++;
      }
    }
    for (int k = (int)(0); k < (int)(n); k++) {
      if (sigma[j].inhabit(pair<long long int, long long int>(x[k], y[k]))) {
        c[j][k] += 1.0 / u;
      }
    }
  }
  for (int i = (int)(0); i < (int)(n); i++) {
    for (int j = (int)(0); j < (int)(e); j++) {
      for (int k = (int)(0); k < (int)(n); k++) {
        a[i][k] += b[i][j] * c[j][k];
      }
    }
  }
  if (DEBUG) {
    cerr << "a:" << endl;
    for (int i = (int)(0); i < (int)(n); i++) {
      for (int j = (int)(0); j < (int)(n); j++) {
        cerr << " " << a[i][j];
      }
      cerr << endl;
    }
    cerr << "b:" << endl;
    for (int i = (int)(0); i < (int)(n); i++) {
      for (int j = (int)(0); j < (int)(e); j++) {
        cerr << " " << b[i][j];
      }
      cerr << endl;
    }
    cerr << "c:" << endl;
    for (int i = (int)(0); i < (int)(e); i++) {
      for (int j = (int)(0); j < (int)(n); j++) {
        cerr << " " << c[i][j];
      }
      cerr << endl;
    }
    vector<vector<double> > uu(a);
    for (int i = (int)(0); i < (int)(30); i++) {
      cerr << "a^" << (i + 1) << ":" << endl;
      for (int i = (int)(0); i < (int)(n); i++) {
        for (int j = (int)(0); j < (int)(n); j++) {
          cerr << " " << uu[i][j];
        }
        cerr << endl;
      }
      uu = mul(uu, a);
    }
  }
  int THR = 30;
  vector<vector<vector<double> > > chino(
      THR, vector<vector<double> >(n, vector<double>(n, 0)));
  for (int i = (int)(0); i < (int)(n); i++) {
    chino[0][i][i] = 1;
  }
  for (int i = (int)(1); i < (int)(THR); i++) {
    chino[i] = mul(chino[i - 1], a);
  }
  int q;
  scanf("%d", &q);
  for (int i = (int)(0); i < (int)(q); i++) {
    int t, m;
    scanf("%d%d", &t, &m);
    t--;
    m = min(m, THR - 1);
    double ans = 0.0;
    for (int i = (int)(0); i < (int)(n); i++) {
      ans = max(ans, chino[m][i][t]);
    }
    for (int i = (int)(0); i < (int)(e); i++) {
      double sum = 0;
      for (int j = (int)(0); j < (int)(n); j++) {
        sum += chino[m - 1][j][t] * c[i][j];
      }
      ans = max(ans, sum);
    }
    printf("%.15f\n", ans);
  }
}