POJ 1755 Triathlon

http://poj.org/problem?id=1755

题意:铁人三项，每个人有自己在每一段的速度，求有没有一种3条路线长度都不为0的设计使得某个人能严格获胜?

我们枚举每个人获胜，得到不等式组:s1/v1+s2/v2+s3/v3<s1/v1'+s2/v2'+s3/v3' 是三维半平面交？接着看

我们令两边同时除以s3得到

(s1/s3)*(1/v1)+(s2/s3)*(1/v2)+(1/v3)<(s1/s3)*(1/v1')+(s2/s3)*(1/v2')+(1/v3')

这样s1/s3和s2/s3就变成了唯二的变量了，这样不等式组就变成了二元一次不等式组了，用半平面交就可以解决了。

注意有些地方:

eps要设置到1e-16

如果一个人有2个速度和另一个人相等，直接特判掉第三个，不建直线。

#include<cmath>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
const double eps=1e-;
const double inf=;
struct Point{
 double x,y;
 Point(){}
 Point(double x0,double y0):x(x0),y(y0){}
}p[];
struct Line{
 Point s,e;
 double slop;
 Line(){}
 Line(Point s0,Point e0):s(s0),e(e0){}
}l[],c[];
int tot,n,v1[],v2[],v3[],num[];
int sgn(double x){
 if (x>eps) return ;
 if (x<-eps) return -;
 return ;
}
int read(){
 int t=,f=;char ch=getchar();
 while (ch<''||ch>''){if (ch=='-') f=-;ch=getchar();}
 while (''<=ch&&ch<=''){t=t*+ch-'';ch=getchar();}
 return t*f;
}
double operator *(Point p1,Point p2){
 return p1.x*p2.y-p1.y*p2.x;
}
Point operator -(Point p1,Point p2){
 return Point(p1.x-p2.x,p1.y-p2.y);
}
bool pd(int id){
 for (int i=;i<=n;i++)
  if (i!=id)
 {
  int cnt=,ans1=,ans2=;
  if (v1[i]==v1[id]) cnt++;
  else ans1=v1[id],ans2=v1[i];
  if (v2[i]==v2[id]) cnt++;
  else ans1=v2[id],ans2=v2[i];
  if (v3[i]==v3[id]) cnt++;
  else ans1=v3[id],ans2=v3[i];
  if (cnt==) return ;
  else if (cnt==){
   if (ans1<=ans2) return ;
  }
  num[i]=cnt;
 }
 return ;
}
void build(int id){
    tot=;
    l[++tot].s=Point(,inf);l[tot].e=Point(,);
    l[++tot].s=Point(,);l[tot].e=Point(inf,);
    l[++tot].s=Point(inf,);l[tot].e=Point(,inf);
    for (int i=;i<=n;i++)
     if (i!=id){
     if (num[i]==) continue;
     tot++;
     double A=((double)v1[i]-v1[id])/((double)v1[id]*v1[i]);
     double B=((double)v2[i]-v2[id])/((double)v2[id]*v2[i]);
     double C=((double)v3[i]-v3[id])/((double)v3[id]*v3[i]);
     int a=sgn(A),b=sgn(B),c=sgn(C);
     if (a==){
      l[tot].s.x=-;l[tot].s.y=(-C)/B;
      l[tot].e.x=;l[tot].e.y=(-C)/B;
      if (b>) std::swap(l[tot].s,l[tot].e);
     }else if (b==){
      l[tot].s=Point((-C)/A,-);
      l[tot].e=Point((-C)/A,);
      if (a<) std::swap(l[tot].s,l[tot].e);
     }else{
      l[tot].s=Point(-,-((-)*A+C)/B);
      l[tot].e=Point(,-(A+C)/B);
      if (b>) std::swap(l[tot].s,l[tot].e);
     }
 }
 for (int i=;i<=tot;i++)
  l[i].slop=atan2(l[i].e.y-l[i].s.y,l[i].e.x-l[i].s.x);
}
Point inter(Line p1,Line p2){
 double k1=(p2.e-p1.s)*(p1.e-p1.s);
 double k2=(p1.e-p1.s)*(p2.s-p1.s);
 double t=(k2)/(k1+k2);
 double x=p2.s.x+(p2.e.x-p2.s.x)*t;
 double y=p2.s.y+(p2.e.y-p2.s.y)*t;
 return Point(x,y);
}
bool jud(Line p1,Line p2,Line p3){
 Point p=inter(p1,p2);
 return (p-p3.s)*(p3.e-p3.s)>=eps;
}
bool cmp(Line p1,Line p2){
 if (fabs(p1.slop-p2.slop)<=eps) return (p1.e-p1.s)*(p2.e-p1.s)<=;
 else return p1.slop<p2.slop;
}
bool phi(){
   int cnt=;
   std::sort(l+,l++tot,cmp);
   for (int i=;i<=tot;i++)
    if (fabs(l[i].slop-l[i-].slop)>eps) l[++cnt]=l[i];
   int L=,R=;c[]=l[];c[]=l[];
   for (int i=;i<=cnt;i++){
    while (L<R&&jud(c[R],c[R-],l[i])) R--;
    while (L<R&&jud(c[L],c[L+],l[i])) L++;
    c[++R]=l[i];
   }
   while (L<R&&jud(c[R],c[R-],c[L])) R--;
   while (L<R&&jud(c[L],c[L+],c[R])) L++;
   if (R-L+<) return ;
   c[R+]=c[L];
   cnt=;
   for (int i=L;i<=R;i++)
    p[++cnt]=inter(c[i],c[i+]);
   p[cnt+]=p[];
   double res=;
   for (int i=;i<=cnt;i++)
    res+=p[i]*p[i+];
   res/=2.0;
   res=fabs(res);
   if (res<=eps) return ;
   else return ;
}
bool work(int id){
   build(id);
   if (phi()) return ;
   else return ;
}
int main(){
   n=read();
   for (int i=;i<=n;i++){
    v1[i]=read(),v2[i]=read(),v3[i]=read();
   }
   for (int i=;i<=n;i++){
    if (pd(i)&&work(i)) puts("Yes");
    else puts("No");
   }
   return ;
}