poj 1755 Triathlon 半平面交求不等式的 是否为空集-------构造有向直线

double V[Max_N] , U[Max_N] , W[Max_N] ;
int main(){
    int t, i , j , ln ,n , k;
    double a,b,c;
    scanf("%d",&n) ;
    for(i = 0 ; i< n ; i++)
        scanf("%lf%lf%lf", &V[i], &U[i], &W[i]) ;
    for(i = 0 ; i< n ; i++){
        k = 0 ;
        double val = 10000 ;
        for(j = 0 ; j < n; j++){ //构造有向直线 ax + by + c > 0 , 有向（b, -a）对应是逆时针
            if( j != i){
                if(V[i] <= V[j] && U[i] <= U[j] && W[i] <= W[j]) { //特判
                    break;
                }
                if( V[i] >= V[j] && U[i] >= U[j] && W[i] >= W[j] ) //特判
                    continue ;
                a = (val / V[j] - val / W[j] ) - (val / V[i] - val / W[i]) ;
                b = (val / U[j] - val / W[j] ) - (val / U[i] - val / W[i]) ;
                c = val / W[j] - val / W[i] ;
                if(fabs(a) < fabs(b)) // 防止b 为0
                    List[k].s = Point( 0 , - c / b) ;
                else
                    List[k].s = Point( - c / a , 0) ;
                List[k].e = List[k].s + Point(b, -a) ;
                k++;
             }
        }
        if(j != n) {printf("No\n") ; continue ;}
        //由题意增加的三个平面 x > 0 y > 0 , x + y < 1(总长度为1)
        List[k++] = pVector( Point(0, 0) , Point( 1, 0) ) ;
        List[k++] = pVector(Point( 0 , 1) , Point(0 , 0 ) ) ;
        List[k++] = pVector(Point( 0 , 1) , Point(-1 , 2 ) ) ;
        if(  halfPanelCross(List , k)  ) // 并列冠军不算，故用的面积计算的， 但是测试不用面积用顶点也ac?
            printf("Yes\n") ;
        else printf("No\n") ;
    }
    return 0;
}