#include <bits/stdc++.h>
using namespace std;

int main() {
    int n;
    scanf("%d", &n);
    vector<double> x(n), y(n);
    for (int i = 0; i < n; i++) scanf("%lf%lf", &x[i], &y[i]);

    vector<double> sx(n + 1), sx2(n + 1), sxy(n + 1), sy(n + 1);
    for (int i = 0; i < n; i++) {
        sx[i + 1] = sx[i] + x[i];
        sx2[i + 1] = sx2[i] + x[i] * x[i];
        sxy[i + 1] = sxy[i] + x[i] * y[i];
        sy[i + 1] = sy[i] + y[i];
    }

    int m;
    scanf("%d", &m);
    for (int q = 0; q < m; q++) {
        int l, r;
        double lambda, x;
        scanf("%d%d%lf%lf", &l, &r, &lambda, &x);
        l--;

        int cn = r - l;
        double sum_x = sx[r] - sx[l];
        double sum_x2 = sx2[r] - sx2[l];
        double sum_xy = sxy[r] - sxy[l];
        double sum_y = sy[r] - sy[l];
        double k = (sum_xy * (cn + lambda) - sum_x * sum_y) / ((sum_x2 + lambda) * (cn + lambda) - sum_x * sum_x);
        double b = (sum_y - k * sum_x) / (cn + lambda);
        double res = k * x + b;
        printf("%.17f\n", res);
    }
}