#include <bits/stdc++.h>

using namespace std;

typedef long long ll;

template<typename T> inline void read(T &s) {
	char c = getchar(); int f = 1; s = 0;
	while (!isdigit(c)) { if (c == '-') f = -1; c = getchar(); }
	while (isdigit(c)) s = (s << 3) + (s << 1) + (c ^ 48), c = getchar();
	if (f < 0) s = -s;
}

const int N = 1e5 + 5;
const double eps = 1e-9;

bool equal(double x, double y) { return fabs(x - y) < eps; }

int n;
struct point {
	double x, y;
} p[N];
struct cir {
	point o;
	double r;
} c;
struct line {
	double k, b;
	double gety(double x) {
		return k * x + b;
	}
};

double area(point u, point v, point w) {
	double X1 = v.x - u.x, Y1 = v.y - u.y;
	double X2 = w.x - u.x, Y2 = w.y - u.y;
	double S = fabs(X1 * Y2 - Y1 * X2);
	return S;
}
double dist(point u, point v) {
	return sqrt((u.x - v.x) * (u.x - v.x) + (u.y - v.y) * (u.y - v.y));
}
point intersect(line u, line v) {
	point w;
	w.x = (v.b - u.b) / (u.k - v.k);
	w.y = u.gety(w.x);
	return w;
}
line vertical(line l, point u) {
	line v;
	v.k = -(double)1.0 / l.k;
	v.b = u.y - u.x * v.k;
	return v;
}
line getline(point u, point v) {
	line l;
	l.k = (v.y - u.y) / (v.x - u.x);
	l.b = u.y - l.k * u.x;
	return l;
}
bool intri(point u, point v, point w, point o) {
	double S = area(u, v, w);
	double S1 = area(u, v, o);
	double S2 = area(u, w, o);
	double S3 = area(v, w, o);
	return equal(S, S1 + S2 + S3);
}

ll ans;

void init() {
	ans = 0;
}
int nxt(int id) {
	++id;
	if (id > n) id -= n;
	return id;
}
void solve() {
	init();
	scanf("%d", &n);
	scanf("%lf%lf%lf", &c.o.x, &c.o.y, &c.r);
	for (int i = 1; i <= n; ++i) scanf("%lf%lf", &p[i].x, &p[i].y);
	if (n == 3) { puts("0"); return ; }
	int L = 1, R = 1;
	double now = 0;
	for (L = 1; L <= n; ++L) {
		if (L == R) R = nxt(R);
		while (true) {
			int G = nxt(R);
			if (nxt(G) == L) break;
			if (intri(p[L], p[R], p[G], c.o)) break;
			line l1 = getline(p[L], p[G]);
			line l2 = vertical(l1, c.o);
			point o = intersect(l1, l2);
			if (c.r - dist(c.o, o) > eps) break;
			now += area(p[L], p[R], p[G]);
			// printf("add (%d, %d, %d)\n", L, R, G);
			// printf("area(%d, %d, %d) = %.0lf\n", L, R, G, area(p[L], p[R], p[G]));
			ans = max(ans, (ll)(now + 0.5));
			R = G;
			// printf("L = %d, R = %d, now = %.5lf\n", L, R, now);
		}
		if (R != nxt(L))
			now -= area(p[L], p[nxt(L)], p[R])/*, printf("del (%d, %d, %d)\n", L, nxt(L), R)*/;
	}
	printf("%lld\n", ans);
}

int main() {
	int T; read(T);
	while (T--) solve();
	return 0;
}