题意:在铁人三项比赛中,给出每个人游泳,骑车,跑步的速度,你可以任意安排游泳,骑车,跑步的路程(路程非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;
}