#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define I inline
#define int long long
using namespace std;

const double eps = 1e-10;
const int N = 100010;
struct point {
	int x,y;
	point(int xx=0,int yy=0) {
		x=xx; y=yy;
	}
}a[N],pc;
typedef point Vector;
int n;
double xc,yc,R;

I point operator+(point p1,point p2) {return point(p1.x+p2.x,p1.y+p2.y);};
I point operator-(point p1,point p2) {return point(p1.x-p2.x,p1.y-p2.y);};
I point operator*(int x,point p1) {return point(p1.x*x,p1.y*x);}; 

I int sgn(double x) {
	if(x<-eps) return -1;
	if(x>eps) return 1;
	return 0;
}
I int Cro(point A,point B) {
	return A.x*B.y-A.y*B.x;
}
I int area(point A,point B,point C) {
	return Cro(A-B,A-C);
}
I int Len(point A) {
	return A.x*A.x+A.y*A.y;
}
I double Abs(double x) {
	return x<0?-x:x;
}
I double dis_PL(point A,point B,point c) {
	return Abs(1.0*Cro(A-c,B-c)/sqrt(1.0*Len(A-B)));
}
I void solve() {
	cin>>n; cin>>xc>>yc>>R; pc=point(xc,yc);
	for(int i=1;i<=n;++i) cin>>a[i].x>>a[i].y;
	//reverse(a+1,a+1+n);
	int r=2,ans=0,sum=0;
	for(int i=1;i<=n;++i) {
		if(i==r) r=r%n+1;		
		//cout<<area(a[i],a[r%n+1],a[r])<<' '<<dis_PL(a[i],a[r%n+1],pc)<<endl;
		while(sgn(dis_PL(a[i],a[r%n+1],pc)-R)>0&&area(a[i],a[r],a[r%n+1])>=0) sum+=area(a[i],a[r],a[r%n+1]),r=r%n+1;
		ans=max(ans,sum);
		if(i!=n) {
		//	cout<<"red"<<area(a[i],a[r],a[i+1])<<endl;
			sum+=area(a[i],a[r],a[i+1]);
		}
	}
	cout<<ans<<endl;
}
signed main()
{
	//cout<<area(point(0,0),point(1,1),point(1,0))<<endl;
	int T; cin>>T;
	while(T--) {
		solve();
	}
} 
/*
1
5
1 1 1
0 0
1 0
5 0
3 3
0 5

*/