poj3335---计算几何

//半平面交，就是求有没有内核，有输出“YES”，否则“NO”
//这里都是逆时针方向
#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<math.h>
#include<queue>
#define dist(a,b) sqrt(1.0*(a.x-b.x)*(a.x-b.x)+1.0*(a.y-b.y)*(a.y-b.y))
#define cross(a,b,c) 1.0*(b.x-a.x)*(c.y-a.y)-1.0*(b.y-a.y)*(c.x-a.x)
#define dot(a,b,c) (b.x-a.x)*(c.x-a.x)+(b.y-a.y)*(c.y-a.y)
#define delt(a) fabs(a)<eps?0:a>0?1:-1
#define pi acos(-1.0)
#define eps 1e-8
#define inf 1e20
#define N 105
#define M 1005
using namespace std;
typedef __int64 i64;
struct TPoint
{
	double x,y;
}pt[N],st,ans[N];
int n,m,t;
void scan()
{
	scanf("%d",&n);
	for(int i=0;i<n;i++)
		scanf("%lf%lf",&pt[i].x,&pt[i].y);
	st.x=st.y=0;
}
double getarea()
{
	double area=0;
	pt[n]=pt[0];
	for(int i=0;i<n;i++) area+=cross(st,pt[i],pt[i+1]);
	return area;
}
void cut(TPoint pl[],TPoint a,TPoint b)
{
	int i,k;
	for(i=k=0;i<m;i++)
	{
		double s1=cross(pl[i],a,b),s2=cross(pl[i+1],a,b);
		int d1=delt(s1),d2=delt(s2);
		if(d1>=0) ans[k++]=pl[i];
		else if(d1*d2<0)
		{
			ans[k].x=(s2*pl[i].x-s1*pl[i+1].x)/(s2-s1);
			ans[k++].y=(s2*pl[i].y-s1*pl[i+1].y)/(s2-s1);
		}
	}
	m=k;
}
void solve()
{
	int i,j,k;
	bool tmp1=false,tmp2=false;
	if(delt(getarea())<0) reverse(pt,pt+n);
	pt[n]=pt[0],m=n;
	for(i=0;i<=n;i++) ans[i]=pt[i];
	for(i=0;i<n;i++)
		cut(ans,pt[i],pt[i+1]);
	if(m==0) puts("NO");
	else puts("YES");
}
int main()
{
    scanf("%d",&t);
	while(t--)
	{
		scan();
		solve();
	}
    return 0;
}