#include <bits/stdc++.h>


#define N 3010

using namespace std;

struct Point{
    double x, y; Point(double X = 0, double Y = 0) {x = X, y = Y;}
    void in() {cin>>x>>y;}
};
int n;
double rc;
long long sum = 0, ans = 0;
Point nd[N], c; 

Point operator+(Point a, Point b) {
    return Point(a.x + b.x, a.y + b.y);
}
Point operator-(Point a, Point b) {
    return Point(a.x - b.x, a.y - b.y);
}
bool operator==(Point a, Point b) {
    return (a.x == b.x && a.y == b.y);
}
double Dot(Point a, Point b) {
    return (a.x * b.x) + (a.y * b.y);
}
double Cro(Point a, Point b) {
    return (a.x * b.y) - (a.y * b.x);
}
double len(Point a) {
    return sqrt(Dot(a,a));
}
double dis(Point p, Point a, Point b) {
    if(a==b) return len(p - a);
    Point x = p - a, y = p - b, z = b - a;
    if(Dot(x, z) < 0) return len(x);
    if(Dot(y, z) > 0) return len(y);
    // double k = Cro(x, z), kk = len(z);
    return fabs(Cro(x, z) / len(z));
}
double sjx(Point p, Point a, Point b) {
    if(p == a || p == b) return 0;
    Point x = p - a, y = p - b, z = b - a;
    return fabs(Cro(x, z));
}
bool check(Point a, Point b) {
    // cout<<dis(c, a, b)<<endl;
    if(dis(c, a, b) > rc) return 1;
    else  return 0;
}
bool mid(Point p, Point a, Point b) {
    Point x = a - p, y = b - p, z = c - p;
    if(Cro(x, z) * Cro(y, z) < 0) return 0;
    else return 1;
}
void sl() {
    
    cin>>n;
    c.in();
    cin>>rc;
    ans = sum = 0;
    int l = 1, r = 1;
    for(int i = 1;i <= n; i++) nd[i].in();
    for(int i = 1; i <= n; i++) {
        if(l < i) l = i;
        if(r < i) r = i;
        while (check(nd[i], nd[r]) && mid(nd[i], nd[r], nd[l])) {
            sum += sjx(nd[i], nd[l], nd[r]);
            l = r; r++; if(r > n) r = 1; 
        } 
        ans = max(sum, ans);
        sum -= sjx(nd[r - 1], nd[i], nd[i + 1]);
    }
    cout<<ans<<"\n";
}
int main() {
    int T; cin>>T;
    while (T--) {
        sl();
    }
}