#include<bits/stdc++.h>
#define int long long
using namespace std;
int n,r;
double eps=1e-8;
struct VECTOR
{
	int x,y;
	double operator^(const VECTOR &v)const
	{
		return x*v.y-y*v.x;
	}
};
struct POINT
{
	int x,y;
	VECTOR operator-(const POINT &p)const
	{
		return {x-p.x,y-p.y};
	}	
}p[2005],p0,o={0,0};
bool eq(double x,double y)
{
	return abs(x-y)<=eps;
}
bool check(POINT p1,POINT p2)
{
	VECTOR v1=p2-p1,v2=p0-p1;
	double dis=(v1^v2)/sqrt(v1*v1);
	return dis>r||eq(dis,r);
}
void solve()
{
	cin>>n;
	cin>>p0.x>>p0.y>>r;
	for(int i=0;i<n;i++)
		cin>>p[i].x>>p[i].y;
	int sum=0;
	for(int i=0;i<n;i++)
		sum+=((p[i]-o)^(p[(i+1)%n]-o));
	int j=0,area=0,ans=0;
	for(int i=0;i<n;i++)
	{
		while(check(p[i],p[(j+1)%n]))
		{
			area-=((p[j]-o)^(p[i]-o));
			j=(j+1)%n;
			area+=((p[j]-o)^(p[i]-o))+((p[(j-1+n)%n]-o)^(p[j]-o));
			if(((p[j]-p[i])^(p0-p[i]))>0)ans=max(ans,area);
			else ans=max(ans,sum-area);
		}
		area-=((p[i]-o)^(p[(i+1)%n]-o))+((p[j]-o)^(p[i]-o));
		area+=((p[j]-o)^(p[(i+1)%n]-o));
	}
	cout<<ans<<endl;
}
signed main()
{
	ios::sync_with_stdio(false),cin.tie(0),cout.tie(0);
	int T=1;
	cin>>T;
	while(T--)solve();
	return 0;
}
/*
4
5
1 1 1
0 0
1 0
5 0
3 3
0 5
6
2 4 1
2 0
4 0
6 3
4 6
2 6
0 3
4
3 3 1
3 0
6 3
3 6
0 3
6
4 2 1
2 0
4 0
6 3
4 6
2 6
0 3
*/

answer.code: In function ‘bool check(POINT, POINT)’:
answer.code:29:35: error: no match for ‘operator*’ (operand types are ‘VECTOR’ and ‘VECTOR’)
   29 |         double dis=(v1^v2)/sqrt(v1*v1);
      |                                 ~~^~~
      |                                 |  |
      |                                 |  VECTOR
      |                                 VECTOR
In file included from /usr/include/c++/13/ccomplex:39,
                 from /usr/include/x86_64-linux-gnu/c++/13/bits/stdc++.h:127,
                 from answer.code:1:
/usr/include/c++/13/complex:395:5: note: candidate: ‘template<class _Tp> constexpr std::complex<_Tp> std::operator*(const complex<_Tp>&, const complex<_Tp>&)’
  395 |     operator*(const complex<_Tp>& __x, const complex<_Tp>& __y)
      |     ^~~~~~~~
/usr/include/c++/13/complex:395:5: note:   template argument deduction/substitution failed:
answer.code:29:36: note:   ‘VECTOR’ is not derived from ‘const std::complex<_Tp>’
   29 |         double dis=(v1^v2)/sqrt(v1*v1);
      |                                    ^~
/usr/include/c++/13/complex:404:5: note: candidate: ‘template<class _Tp> constexpr std::complex<_Tp> std::operator*(const complex<_Tp>&, const _Tp&)’
  404 |     operator*(const complex<_Tp>& __x, const _Tp& __y)
      |     ^~~~~~~~
/usr/include/c++/13/complex:404:5: note:   template argument deduction/substitution failed:
answer.code:29:36: note:   ‘VECTOR’ is not derived from ‘const std::complex<_Tp>’
   29 |         double dis=(v1^v2)/sqrt(v1*v1);
      |                                    ^~
/usr/include/c++/13/complex:413:5: note: candidate: ‘template<class _Tp> constexpr std::complex<_Tp> std::operator*(const _Tp&, const complex<_Tp>&)’
  413 |     operator*(const _Tp& __x, const complex<_Tp>& __y)
      |     ^~~~~~~~
/usr/include/c++/13/complex:413:5: note:   template argument deduction/substitution failed:
answer.code:29:36: note:   ‘VECTOR’ is not derived from ‘const std::complex<_Tp>’
   29 |         double dis=(v1^v2)/sqrt(v1*v1);
      |                                    ^~
In file included from /usr/include/c++/13/valarray:605,
                 from /usr/include/x86_64-linux-gnu/c++/13/bits/stdc++.h:166:
/usr/include/c++/13/bits/valarray_after.h:407:5: note: candidate: ‘template<class _Dom1, class _Dom2> std::_Expr<std::__detail::_BinClos<std::__multiplies, std::_Expr, std::_Expr, _Dom1, _Dom2>, typename std::__fun<std::__multiplies, typename _Dom1::value_type>::result_type> std::operator*(const _Expr<_Dom1, typename _Dom1::value_type>&, const _Expr<_Dom2, typename _Dom2::value_type>&)’
  407 |     _DEFINE_EXPR_BINARY_OPERATOR(*, struct std::__multiplies)
      |     ^~~~~~~~~~~~~~~~~~~~~~~~~~~~
/usr/include/c++/13/bits/valarray_after.h:407:5: note:   template argument deduction/substitution failed:
answer.code:29:36: note:   ‘VECTOR’ is not derived from ‘const std::_Expr<_Dom1, typename _Dom1::value_type>’
   29 |         double dis=(v1^v2)/sqrt(v1*v1);
      |                                    ^~
/usr/include/c++/13/bits/valarray_after.h:407:5: note: candidate: ‘template<class _Dom> std::_Expr<std::__detail::_BinClos<std::__multiplies, std::_Expr, std::_Constant, _Dom, typename _Dom::value_type>, typename std::__fun<std::__multiplies, typename _Dom1::value_type>::result_type> std::operator*(const _Expr<_Dom1, typename _Dom1::value_type>&, const typename _Dom::value_type&)’
  407 |     _DEFINE_EXPR_BINARY_OPERATOR(*, struct std::__multiplies)
      |     ^~~~~~~~~~~~~~~~~~~~~~~~~~~~
/usr/include/c++/13/bits/valarray_after.h:407:5: note:   template argument deduction/substitution failed:
answer.code:29:36: note:   ‘VECTOR’ is not derived from ‘const std::_Expr<_Dom1, typename _Dom1::value_type>’
   29 |         double dis=(v1^v2)/sqrt(v1*v1);
      |                                    ^~
/usr/include/c++/13/bits/valarray_after.h:407:5: note: candidate: ‘template<class _Dom> std::_Expr<std::__detail::_BinClos<std::__multiplies, std::_Constant, std::_Expr, typename _Dom::value_type, _Dom>, typename std::__fun<std::__multiplies, typename _Dom1::value_type>::result_type> std::operator*(const typename _Dom::value_type&, const _Expr<_Dom1, typename _Dom1::value_type>&)’
  407 |     _DEFINE_EXPR_BINARY_OPERATOR(*, struct std::__multiplies)
      |     ^~~~~~~~~~~~~~~~~~~~~~~~~~~~
/usr/include/c++/13/bits/valarray_after.h:407:5: note:   template argument deduction/substitution failed:
answer.code:29:36: note:   ‘VECTOR’ is not derived from ‘const std::_Expr<_Dom1, typename _Dom1::value_type>’
   29 |         double dis=(v1^v2)/sqrt(v1*v1);
      |                                    ^~
/usr/include/c++/13/bits/valarray_after.h:407:5: note: candidate: ‘template<class _Dom> std::_Expr<std::__detail::_BinClos<std::__multiplies, std::_Expr, std::_ValArray, _Dom, typename _Dom::value_type>, typename std::__fun<std::__multiplies, typename _Dom1::value_type>::result_type> std::operator*(const _Expr<_Dom1, typename _Dom1::value_type>&, const valarray<typen...