POJ 1755 Triathlon【半平面交】

题意：在铁人三项比赛中，给出每个人游泳，骑车，跑步的速度，你可以任意安排游泳，骑车，跑步的路程（路程非0），若某人可以获得第一（严格），则输出Yes，否则No

分析：将三个路程其中一个设为1，剩下两个是x,y.对于每个人i，列出方程Ai*x+Bi*y+Ci<Aj*x+Bj*x+Cj。求半平面交，若解集为空，No。否则Yes。

#include<stdio.h>
#include<string.h>
#include<math.h>
#define maxn 110
#define INF 1e30
#define eps 1e-8
struct node
{
	double u,v,w;
}man[maxn];
int n;
struct point
{
	double x,y;
	point(){}
	point(double a,double b){
		x=a,y=b;
	}
}p[maxn],q[maxn];
int size,s;
int sig(double k)
{
	return (k<-eps)?-1:(k>eps);
}
void init()
{
	p[1]=point(0,0);
	p[2]=point(INF,0);
	p[3]=point(INF,INF);
	p[4]=point(0,INF);
	p[5]=p[1];
	p[0]=p[4];
	size=4;
}
double ab(double x)
{ return x>0?x:-x; }
point intersect(const point &x,const point &y,double a,double b,double c)
{
	double u=ab(a*x.x+b*x.y+c);
	double v=ab(a*y.x+b*y.y+c);
	return point((x.x*v+y.x*u)/(u+v),(x.y*v+y.y*u)/(u+v));
}
void cut(double a,double b,double c)
{
	s=0;
	int i;
	for(i=1;i<=size;i++)
	{
		if(sig(a*p[i].x+b*p[i].y+c)<1)
			q[++s]=p[i];
		else 
		{
			if(sig(a*p[i-1].x+b*p[i-1].y+c)<0)
				q[++s]=intersect(p[i-1],p[i],a,b,c);
			if(sig(a*p[i+1].x+b*p[i+1].y+c)<0)
				q[++s]=intersect(p[i+1],p[i],a,b,c);
		}
	}
	for(i=1;i<=s;i++)
		p[i]=q[i];
	p[s+1]=p[1];
	p[0]=p[s];
	size=s;
}
int win(int x)
{
	int i;
	for(i=0;i<n;i++)
		if(i!=x&&man[i].u>=man[x].u&&man[i].v>=man[x].v&&man[i].w>=man[x].w)
			return false;
	init();
	double a,b,c;
	for(i=0;i<n;i++)
		if(i!=x)
		{
			a=1/man[x].u-1/man[i].u;
			b=1/man[x].v-1/man[i].v;
			c=1/man[x].w-1/man[i].w;
			cut(a,b,c);
			if(size<3)
				return 0;
		}
	return size>=3;
}
int main()
{
	int i;
	while(scanf("%d",&n)!=EOF)
	{
		for(i=0;i<n;i++)
			scanf("%lf%lf%lf",&man[i].u,&man[i].v,&man[i].w);
		for(i=0;i<n;i++)
		{
			if(win(i))
				printf("Yes\n");
			else printf("No\n");
		}
	}
	return 0;
}