#include <bits/stdc++.h>
using namespace std;
#define ll long long 
#define bi __int128

const int N=1e6+5;
int n;
ll ans=0,res=0;

struct vec{
  ll x,y;
  vec(ll a=0,ll b=0){x=a,y=b;}
  vec operator +(const vec &a) const{return vec(x+a.x,y+a.y);}
  vec operator -(const vec &a) const{return vec(x-a.x,y-a.y);}
}a[N];
struct Cir{
  vec a;
  ll r;
  Cir(vec x=vec(0,0),ll R=0){a=x,r=R; }
}C;

ll Dot(vec a,vec b){return a.x*b.x+a.y*b.y;}
ll Cro(vec a,vec b){return a.x*b.y-b.x*a.y;}

ll dcmp(ll x){return x==0?0:(x<0?-1:1);}

bool inter(int l,int r){
  vec nw=a[r]-a[l];
  ll s=Cro(nw,C.a-a[l]),len=Dot(nw,nw);  
  if((bi)s*s>(bi)C.r*C.r*len) return true;
  return false;	
}

bool chk(int l,int r){
  if(dcmp(Cro(C.a-a[l],a[r]-a[l]))!=dcmp(Cro(C.a-a[l],a[(r+1)%n]-a[l]))) return false;
//   cerr<<"chk: "<<l<<' '<<r<<endl;
  return inter(l,(r+1)%n);
}

void SOLVE(){
  cin>>n>>C.a.x>>C.a.y>>C.r;
  for(int i=0;i<n;i++) cin>>a[i].x>>a[i].y;

  ans=res=0;

  for(int l=0,r=1;l<n;l++){
	if(!inter(l,(l+1)%n)) continue;
    if(!inter((l-1+n)%n,l)) r=(l+1)%n;
	while(chk(l,r)) ans+=Cro(a[r]-a[l],a[(r+1)%n]-a[l]),r=(r+1)%n;
    // cerr<<"l = "<<l<<' '<<r<<endl;
	
	res=max(res,ans);
	// cerr<<"ans = "<<ans<<endl<<endl;
	ans-=Cro(a[l]-a[r],a[(l+1)%n]-a[r]);
  }
  cout<<res<<'\n';
}

int main(){
  ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
  int _;cin>>_;
  while(_--) SOLVE();
  return 0;
}
/*
1
6
2 4 1
2 0
4 0
6 3
4 6
2 6
0 3
*/