#include<bits/stdc++.h>
#define F(i, a, b) for(int i = a; i <= b; i ++)
#define int long long
using namespace std;
typedef double db;
const db eps = 1e-9;

int sgn(db x) {
    if(fabs(x) < eps) return 0;
    else return x < 0 ? -1 : 1;
}
struct point {
    int x, y;
    point() {}
    point (int x, int y) : x(x), y(y) {}
    point operator + (point B) {return point(x + B.x, y + B.y);}
    point operator - (point B) {return point(x - B.x, y - B.y);}
}; 
struct line {
    point p1, p2;
    line() {}
    line (point p1, point p2) : p1(p1), p2(p2) {}
};
db dis(point a, point b) {return sqrt(1.0 * (a.x - b.x) * (a.x - b.x) + (a.y - b.y) * (a.y - b.y));}
int cross(point a, point b) {return a.x * b.y - a.y * b.x;}
int polygon_area(vector<point> &p, int n) {
    int area = 0;
    F(i, 0, n - 1) area += cross(p[i], p[(i + 1) % n]);
    return area;
}
struct circle {
    point c;
    int r;
    circle() {}
    circle(int x, int y, int _r) {c = point(x, y); r = _r;}  
} C;

db dis_point_line(point p, line v) {
    return fabs(1.0 * cross(p - v.p1, v.p2 - v.p1) / dis(v.p1, v.p2));
}
int line_circle_relation(line v, circle C) {
    db dst = dis_point_line(C.c, v);
    if(sgn(dst - C.r) < 0) return 0;
    else if(sgn(dst - C.r) == 0) return 1;
    else return 2;
}
int point_line_relation(point p, line v) {return sgn(cross(p - v.p1, v.p2 - v.p1));}

signed main() {
    ios::sync_with_stdio(false);
    cin.tie(0);
    int t; cin >> t;
    while(t --) {
        int n; cin >> n;
        int xx, yy, r; cin >> xx >> yy >> r;
        C = circle(xx, yy, r);
        vector<point> a;
        F(i, 1, n) {
            int x, y; cin >> x >> y;
            a.push_back(point(x, y));
        }
        int sum = polygon_area(a, n);
        int res = 0, now = 0;
        for(int i = 0, j = 0; i < n; i ++) {
            if(i == j) {
                j = (j + 1) % n;
                now = 0;
            } else {
                now += cross(a[j], a[i]);
            }
            while((j + 1) % n != i && line_circle_relation(line(a[i], a[j]), C)) {
                if(point_line_relation(C.c, line(a[i], a[j])) == 1) {
                    res = max(res, sum - now);
                } else {
                    res = max(res, now);
                }
                now -= cross(a[j], a[i]), now += cross(a[j], a[(j + 1) % n]), now += cross(a[(j + 1) % n], a[i]); 
                j = (j + 1) % n;
            }
            now -= cross(a[i], a[(i + 1) % n]), now -= cross(a[j], a[i]);
        }
        cout << res << "\n";
    }
    return 0;
}