题目大意:
铁人三项分连续三段:游泳 自行车 赛跑
已知各选手在每个单项中的速度v[i],u[i],w[i]
设计每个单项的长度 可以让某个特定的选手获胜
判断哪些选手有可能获得冠军
输出n行 有可能获得冠军为Yes 不可能为No
设赛程总长为1,游泳x 自行车y,则赛跑为1-x-y
若选手 i 可以打败选手 j 则
x / v[ i ] + y / u[ i ] + ( 1-x-y ) / w[ i ] < x / v[ j ] + y / u[ j ] + ( 1-x-y ) / w[ j ]
整理成 a*x+b*y+c >0 的形式 那么得到
a = ( 1 / v[ j ] - 1 / w[ j ] ) - ( 1 / v[ i ] - 1 / w[ i ] )
b = ( 1 / u[ j ] - 1 / w[ j ] ) - ( 1 / u[ i ] - 1 / w[ i ] )
c = 1 / w[ j ] - 1 / w[ i ]
最后加上三个固定约束
x > 0 , y > 0 , 1 - x - y > 0
只要 i 与其他所有 j 存在满足的解 就为Yes
![](https://i-blog.csdnimg.cn/blog_migrate/8f900a89c6347c561fdf2122f13be562.gif)
![](https://i-blog.csdnimg.cn/blog_migrate/961ddebeb323a10fe0623af514929fc1.gif)
#include <bits/stdc++.h> using namespace std; const double eps=1e-10; double add(double a,double b) { if(abs(a+b)<eps*(abs(a)+abs(b))) return 0; return a+b; } struct P { double x,y; P(){} P(double _x,double _y):x(_x),y(_y){} P operator - (P p) { return P(add(x,-p.x),add(y,-p.y)); } P operator + (P p) { return P(add(x,p.x),add(y,p.y)); } P operator / (double d) { return P(x/d,y/d); } P operator * (double d) { return P(x*d,y*d); } double dot (P p) { return add(x*p.x,y*p.y); } double det (P p) { return add(x*p.y,-y*p.x); } void read(){ scanf("%lf%lf",&x,&y); } }; struct L { P p,v; double ang; L(){} L(P _p,P _v):p(_p),v(_v){ ang=atan2(v.y,v.x); } bool operator < (const L& b)const { return ang<b.ang; } }l[105]; int v[105],u[105],w[105]; int n, cnt; bool onLeft(L l,P p) { return (l.v).det(p-l.p)>0; } /// p在l的左边 P ins(L a,L b) { return a.p+a.v*((b.v).det(a.p-b.p)/(a.v).det(b.v)); } /// a与b的交点 int insHp() { sort(l,l+cnt); vector <P> pi(2*cnt); vector <L> li(2*cnt); int head,tail; li[head=tail=0]=l[0]; for(int i=1;i<cnt;i++) { while(head<tail && !onLeft(l[i],pi[tail-1])) tail--; while(head<tail && !onLeft(l[i],pi[head])) head++; li[++tail]=l[i]; if(abs((li[tail].v).det(li[tail-1].v))<eps) { tail--; if(onLeft(li[tail],l[i].p)) li[tail]=l[i]; } if(head<tail) pi[tail-1]=ins(li[tail],li[tail-1]); } while(head<tail && !onLeft(li[head],pi[tail-1])) tail--; if(tail-head<=1) return 0; pi[tail]=ins(li[tail],li[head]); return tail-head+1; } /// 半平面交 返回最后得到的多边形的顶点数 void solve() { for(int i=0;i<n;i++) { cnt=0; bool ok=1; double k=10000; /// 数据范围来说 设k=10000 for(int j=0;j<n;j++) { if(i==j) continue; if(v[i]<=v[j] && u[i]<=u[j] && w[i]<=w[j]) { ok=0; break; /// 必败 } if(v[i]>=v[j] && u[i]>=u[j] && w[i]>=w[j]) continue; /// 必胜 double a=(k/v[j]-k/w[j])-(k/v[i]-k/w[i]); double b=(k/u[j]-k/w[j])-(k/u[i]-k/w[i]); double c=k/w[j]-k/w[i]; /// 数值过小会产生精度误差 P v=P(b,-a); if(abs(a)>abs(b)) l[cnt]=L(P(-c/a,0),v); else l[cnt]=L(P(0,-c/b),v); cnt++; } if(ok) { l[cnt++]=L(P(0,0),P(0,-1)); l[cnt++]=L(P(0,0),P(1,0)); l[cnt++]=L(P(0,1),P(-1,1)); // 三个固定约束 if(!insHp()) ok=0; // 半平面交无解 说明必败 } if(ok) printf("Yes\n"); else printf("No\n"); } } int main() { while(~scanf("%d",&n)) { for(int i=0;i<n;i++) scanf("%d%d%d",&v[i],&u[i],&w[i]); solve(); } return 0; } //22648