#include <bits/stdc++.h>

using namespace std;

#define y1 fuckcpp

const int MAX_N = 2000 + 5;
const double Pi = acos(-1);

int N, Q;
struct Node{
    int x, y, id;
    double angle;
    
    Node () {}
    Node (int x_, int y_, int id_) {
        x = x_; y = y_; id = id_;
        angle = atan2(y, x);
    }
    
    inline bool operator < (const Node &comp) const {
        return angle < comp.angle;
    }
}node[MAX_N];
double ans[MAX_N][MAX_N], dis[MAX_N][MAX_N];

bool check_pi(double x, double y) {
    // x in [y, y+Pi]
    if (y < -Pi) y += Pi * 2;
    if (y < 0) {
        return x >= y && x <= y + Pi;
    }
    return (x >= y || x <= y - Pi);
}

double query(int lx, int rx, int ly, int ry) {
    return ans[rx][ry] - ans[rx][ly - 1] - ans[lx - 1][ry] + ans[lx - 1][ly - 1];
}

double calc_dis(int x1, int y1, int x2, int y2) {
    return sqrt((x1 - x2) * (x1 - x2) + (y1 - y2) * (y1 - y2));
}

double calc_dis(int i, int j) {
    return calc_dis(node[i].x, node[i].y, node[j].x, node[j].y);
}

bool leftOf(int i, int j, int k) {
    // k is left of ij
    return (node[k].x - node[i].x) * (node[j].y - node[i].y) >= (node[k].y - node[i].y) * (node[j].x - node[i].x);
}

int main() {
    scanf("%d", &N);
    for (int i = 1; i <= N; i ++) {
        int x, y;
        scanf("%d%d", &x, &y);
        node[i] = Node(x, y, i);
    }
    sort(node + 1, node + N + 1);
    for (int i = 1; i <= N; i ++) {
        // fprintf(stderr, "i = %d, x = %d, y = %d, angle = %.3lf\n", i, node[i].x, node[i].y, node[i].angle);
    }
    for (int i = 1; i <= N; i ++) {
        // counter-clockwise
        fprintf(stderr, "i = %d\n", i);
        for (int j = (i < N ? i + 1 : 1), k = -1; j != i && check_pi(node[j].angle, node[i].angle); j = (j < N ? j + 1 : 1)) {
            if (k == -1 || leftOf(i, j, k)) k = j;
            dis[node[i].id][node[j].id] = calc_dis(i, k) + calc_dis(k, j);
            // fprintf(stderr, "i = %d, j = %d, k = %d, id = (%d, %d), res = %.3lf\n", i, j, k, node[i].id, node[j].id, dis[node[i].id][node[j].id]);
        }
        // clockwise
        for (int j = (i > 1 ? i - 1 : N), k = -1; j != i && check_pi(node[j].angle, node[i].angle - Pi); j = (j > 1 ? j - 1 : N)) {
            if (k == -1 || !leftOf(i, j, k)) k = j;
            dis[node[i].id][node[j].id] = calc_dis(i, k) + calc_dis(k, j);
            // fprintf(stderr, "i = %d, j = %d, k = %d, id = (%d, %d), res = %.3lf\n", i, j, k, node[i].id, node[j].id, dis[node[i].id][node[j].id]);
        }
    }
    for (int i = 1; i <= N; i ++) {
        for (int j = i; j <= N; j ++) {
            // fprintf(stderr, "i = %d, j = %d, dis = %.3lf\n", i, j, dis[i][j]);
            for (int k = 1; k <= N; k ++) {
                ans[i][j] = max(ans[i][j], dis[i][k] + dis[k][j]);
            }
            ans[j][i] = ans[i][j];
        }
    }
    for (int i = 1; i <= N; i ++) {
        for (int j = 1; j <= N; j ++) {
            ans[i][j] += ans[i - 1][j] + ans[i][j - 1] - ans[i - 1][j - 1];
        }
    }
    scanf("%d", &Q);
    while (Q --) {
        int a, b, c, d;
        scanf("%d%d%d%d", &a, &b, &c, &d);
        printf("%.10lf\n", query(a, b, c, d) / (b - a + 1) / (d - c + 1));
    }
    return 0;
}

/*
5
7 0
3 3
0 7
-3 3
-7 0
6
1 1 3 3
3 3 4 4
1 1 5 5
5 5 2 2
2 2 4 4
1 5 1 5
*/