#include <bits/stdc++.h>
using namespace std;
using db = double;
using ld = long double;
constexpr ld eps = 1e-8;
using vdb = complex<db>;
using vld = complex<ld>;
using cv = const vdb&;
template<class T>
istream& operator >>(istream &in, complex<T> &p) {
	T x, y;
	in >> x >> y;
	p = {x, y};
	return in;
}
#define x real()
#define y imag()
#define C(a, b) (conj(a) * (b)).y
#define D(a, b) (conj(a) * (b)).x
bool onseg(cv h, cv p, cv q) {
	return h.x > min(p.x, q.x) - eps && h.x < max(p.x, q.x) + eps
	    && h.y > min(p.y, q.y) - eps && h.y < max(p.y, q.y) + eps;
}
ld dist(cv o, cv p, cv q) {
	cv h = p + (q - p) * (D(q - p, o - p) / norm(q - p));
	if (onseg(h, p, q)) return abs(o - h);
	else return min(abs(o - p), abs(o - q));
}
int main() {
	int n, q;
	cin >> n >> q;
	vector<vdb> a(n);
	vector<vld> la;
	for (vdb &p : a) {
		cin >> p;
		la.push_back(vld(p.x, p.y));
	}
	ld s0 = 0;
	for (int i = 0; i < n; i++) {
		s0 += C(la[i], la[(i + 1) % n]);
	}
	auto inside = [&](const vld &p) {
		ld s1 = 0;
		for (int i = 0; i < n; i++) {
			s1 += fabs(C(la[i] - p, la[(i + 1) % n] - p));
		}
		return fabs(s0 - s1) < eps;
	};
	while (q--) {
		vld A, B, O;
		ld r;
		cin >> A >> B, O = (A + B) / ld(2), r = abs(O - A);
		vdb dO(O.x, O.y);
		ld d = 1e18;
		if (inside(O)) {
			d = 0;
		} else {
			for (int i = 0; i < n; i++) {
				d = min(d, dist(dO, a[i], a[(i + 1) % n]));
			}
		}
		cout << fixed << setprecision(6) << r * r / 2 + d * d << "\n";
	}
	return 0;
}