#include <bits/stdc++.h>

using namespace std;
using ll = long long;
const double eps=1e-11;
const int MAXN=5e5+7;                                                         
struct Point{                                                                        
    ll x,y;                                                                      
    Point(ll x=0,ll y=0):x(x),y(y){}
};
typedef Point Vector;

Vector operator + (Vector A,Vector B){return Vector(A.x+B.x,A.y+B.y);}
Vector operator - (Vector A,Vector B){return Vector(A.x-B.x,A.y-B.y);}
Vector operator * (Vector A,ll B){return Vector(A.x*B,A.y*B);}
Vector operator / (Vector A,ll B){return Vector(A.x/B,A.y/B);}

int dcmp(ll x){if(x==0)return 0;return (x>0)?1:-1;}
bool operator == (const Vector A,const Vector B){
    return dcmp(A.x-B.x)==0 && dcmp(A.y-B.y)==0;
}
ll Dot(Vector A,Vector B){return A.x*B.x+A.y*B.y;}                     
ll Cross(Vector A,Vector B){return A.x*B.y-B.x*A.y;}              

ll ComputePolygonArea(vector<Point> &points)
{
    int point_num = points.size();
    if(point_num < 3)return 0.0;
    ll s = points[0].y * (points[point_num-1].x - points[1].x);
    for(int i = 1; i < point_num; ++i)
        s += points[i].y * (points[i-1].x - points[(i+1)%point_num].x);
    return s;
}
Point a[MAXN],be;
vector<Point> c;
int t;
int n;
ll be_x,be_y,r;
ll Distance(Point p,Point a,Point b)
{
    Vector v1=b-a,v2=p-a;
    return fabs(Cross(v1,v2));
}
bool check(Point P,Point A,Point B){
    if(A==B) return 1;
    Vector v1=B-A,v2=P-A,v3=P-B;
    if(dcmp(Dot(v1,v2))<0)
    {
        __int128 k2=Dot(v2,v2);
        if(k2>=r*r) return 1;
        return 0; 
    }          
    else if(dcmp(Dot(v1,v3))>0)
    {
        __int128 k3=Dot(v3,v3);
        if(k3>=r*r) return 1;
        return 0; 
    }    
    else
    {
        __int128 k1=Dot(v1,v1);
        __int128 sss=Cross(v1,v2);
        
        sss=sss*sss;
        __int128 rk=r;
        k1=k1*rk*rk;
        if(sss>=k1) return 1;
        return 0;
    }
}
// bool check(Point P,Point A,Point B){
//     if(A==B)return 1;
//     Vector v1=B-A,v2=P-A,v3=P-B;
//     if(dcmp(Dot(v1,v2))<0)
//     {
//         if(Length(v2)>r||fabs(Length(v2)-r)<=eps) return 1;
//         return 0; 
//     }          
//     else if(dcmp(Dot(v1,v3))>0)
//     {
//         if(Length(v3)>r||fabs(Length(v3)-r)<=eps) return 1;
//         return 0;
//     }    
//     else
//     {
//         if(fabs(Cross(v1,v2))>r*Length(v1)||fabs(fabs(Cross(v1,v2))-r*Length(v1))<=eps) return 1;
//         return 0;
//     }
// }
bool check1(Point a,Point b,Point c)
{ 
    __int128 res1=Cross(a,b),res2=Cross(a,c);
    if(res1*res2<0) return 0;
    return 1;
}
int main()
{
    cin.tie(nullptr) ->sync_with_stdio(false);

    cin>>t;
    while(t--)
    {
        cin>>n;
        cin>>be_x>>be_y>>r;
        be={be_x,be_y};
        for(int i=1;i<=n;i++) cin>>a[i].x>>a[i].y;
        c.clear();
        ll ans=0LL;
        
        for(int i=1;i<=n;i++)
        {
            c.clear();
            c.push_back(a[i]);
            int j=i,numk=1;
            int jj=i+1;
            jj%=n;
            if(jj==0) jj=n;
            Point k1=be-a[i];
            Point k2=a[jj]-a[i];
            while(numk<n)
            {
                j++;numk++;
                j%=n;
                if(j==0) j=n;
                if(!check(be,a[i],a[j])) break;
                Point k3=a[j]-a[i];
                if(!check1(k1,k2,k3)) break;
                c.push_back(a[j]);
                ans=max(ans,(ll) ComputePolygonArea(c));
            }
            // cerr<<i<<" "<<j-1<<" asd"<<endl;
            // cout<<check(be,a[3],a[1])<<"  asd"<<endl;
            c.clear();
            c.push_back(a[i]);
            j=i;numk=1;
            jj=i-1;
            jj=(jj%n+n)%n;
            if(jj==0) jj=n;
            k1=be-a[i];
            k2=a[jj]-a[i];
            while(numk<n)
            {
                j--;numk++;
                j=(j%n+n)%n;
                if(j==0) j=n;
                if(!check(be,a[i],a[j])) break;
                Point k3=a[j]-a[i];
                if(!check1(k1,k2,k3)) break;
                c.push_back(a[j]);
                ans=max(ans,(ll) ComputePolygonArea(c));
            }

        }

        cout<<ans<<'\n';
    }
    
    return 0;
}
/*

*/