HDU6590 2019杭电多校第一场1013 code

考场上想到是把二维平面的点用一条线分成两个集合，但是没想到用凸包去判断是否相交

跟队友讲了题以后，队友提到了凸包但我没反应过来

看了榜以后以为2题已经结束了，后2个小时都在跟隔壁队伍的同学吹水啊。。。。这道水题就没去开了。。。

走在路上和队友又聊起这题，突然发现凸包相交就行了。。。然后回机房写了30分钟就过了。。。

卧槽，好亏啊，以后比赛一定要每道题都认真思考一下，尽管是过得不多的题，几十个人过的题也是可以做的

UPD:

判断凸包相交要两个条件：

1.判断一个凸包所有点是否会有在另一个凸包内或在另一个凸包边上

2.判断两个凸包是否有线段相交(重合一部分)

今天学弟来问我，我发现我没判断线段相交也过了。。。可能随机数据比较水吧。。。

#include<bits/stdc++.h>
#define maxl 1010
#define eps 1e-9
using namespace std;
 
inline int sgn(double x)
{
	if(x>eps) return 1;
	if(x<-eps) return -1;
	return 0;
}
struct point
{
	double x,y;
	point(double a=0,double b=0)
	{
		x=a;y=b;
	}
	point operator - (const point &b)const
	{
		return point(x-b.x,y-b.y);
	}
	bool operator == (const point &b)const
	{
		return sgn(x-b.x)==0 && sgn(y-b.y)==0;	
	}
	bool operator < (const point &b)const
	{
		return x<b.x-eps;
	}
};
inline double det(const point &a,const point &b)
{
	return a.x*b.y-a.y*b.x;
}
inline double dot(const point &a,const point &b)
{
	return a.x*b.x+a.y*b.y;
}
struct line
{
	point s,e;
	double k;
	line(point a=point(),point b=point())
	{
		s=a;e=b;
		k=atan2(e.y-s.y,e.x-s.x);
	} 
};
struct polygon_convex
{
	vector<point> P;
	polygon_convex(int size=0)
	{
		P.resize(size);
	}
}ap,bp;
inline bool cmp(const point &a,const point &b)
{
	if(sgn(a.x-b.x)==0)
		return sgn(a.y-b.y)<0;
	return sgn(a.x-b.x)<0;
}
 
polygon_convex convex_hull(vector<point> a)
{
	polygon_convex res(2*a.size()+5);
	sort(a.begin(),a.end(),cmp);
	a.erase(unique(a.begin(),a.end()),a.end());
	int m=0,l=a.size();
	for(int i=0;i<l;i++)
	{
		while(m>1 && sgn(det(res.P[m-1]-res.P[m-2],a[i]-res.P[m-2]))<=0)
			--m;
		res.P[m++]=a[i];
	}
	int k=m;
	for(int i=l-2;i>=0;i--)
	{
		while(m>k && sgn(det(res.P[m-1]-res.P[m-2],a[i]-res.P[m-2]))<=0)
			--m;
		res.P[m++]=a[i];
	}
	res.P.resize(m);
	if(a.size()>1)
		res.P.resize(m-1);
	return res;
}
int n;
bool flag;
vector <point> a,b;
 
inline void prework()
{
	scanf("%d",&n);
	a.clear();b.clear();
	int x1,x2,y;
	for(int i=1;i<=n;i++)
	{
		scanf("%d%d%d",&x1,&x2,&y);
		if(y==1)
			a.push_back(point(x1,x2));
		else
			b.push_back(point(x1,x2));
	} 
}
 
inline bool point_on_seg(point p,line a)
{
	return sgn(det(p-a.s,p-a.e))==0 && sgn(dot(p-a.s,p-a.e))<=0;
} 

inline bool parallel(line &l1,line &l2)
{
	return !sgn(det(l1.e-l1.s,l2.e-l2.s));
} 

inline bool point_in(point t,vector <point> a)
{
	int num=0,d1,d2,k,n=((int)a.size())-1;
	for(int i=0;i<n;i++)
	{
		if(point_on_seg(t,line(a[i],a[i+1])))
			return true;
		k=sgn(det(a[i+1]-a[i],t-a[i]));
		d1=sgn(a[i].y-t.y);
		d2=sgn(a[i+1].y-t.y);
		if(k>0 && d1<=0 && d2>0) num++;
		if(k<0 && d2<=0 && d1>0) num++;
	}
	return num%2!=0;
} 

 inline bool line_cross(line l1,line l2)
{
	line vec=line(l1.s,l1.e-l1.s);
	if(parallel(l1,l2))//覆盖的情况 
		return point_on_seg(l1.s,l2) || point_on_seg(l1.e,l2) 
			|| point_on_seg(l2.s,l1) || point_on_seg(l2.e,l1);
	if(det((l2.s-vec.s),vec.e)*det(vec.e,(l2.e-vec.s))<0-eps)
		return false;
	if(det((l2.s-l1.s),(l2.e-l1.s))*det((l2.s-l1.e),(l2.e-l1.e))>0+eps)
		return false;
	
	return true;
}

inline void mainwork()
{
	flag=true;
	ap=convex_hull(a);
	bp=convex_hull(b); 
	ap.P.push_back(ap.P[0]);
	bp.P.push_back(bp.P[0]);
	int la=ap.P.size();
	for(int i=0;i<la;i++)
	if(point_in(ap.P[i],bp.P))
		flag=false;
	int lb=bp.P.size();
	for(int i=0;i<lb;i++)
	if(point_in(bp.P[i],ap.P))
		flag=false;
	for(int i=0;i<la-1;i++)
		for(int j=0;j<lb-1;j++)
		if(line_cross(line(ap.P[i],ap.P[i+1]),line(bp.P[j],bp.P[j+1])))
			flag=false; 
	 
}
 
inline void print()
{
	if(flag)
		puts("Successful!");
	else
		puts("Infinite loop!"); 
}
 
int main()
{
	int t;
	scanf("%d",&t);
	for(int i=1;i<=t;i++)
	{
		prework();
		mainwork();
		print();
	}
	return 0;
}