miaozk2006

const   double   EPS                 =   1e-5;           //   ¼ÆË㾫¶È 
const   int   CROSS                   =   0;                 //   Á½Ö±ÏßÏཻ 
const   int   COLINE                 =   1;                 //   Á½Ö±Ïß¹²Ïß 
const   int   PARALLEL             =   2;                 //   Á½Ö±ÏßƽÐР
struct   Point   {                             //   ¶þάµã»òʸÁ¿ 
        double   x,   y; 
        Point()   {} 
        Point(double   x0,   double   y0):   x(x0),   y(y0)   {} 
}; 
struct   Line   {                               //   ¶þάµÄÖ±Ïß»òÏ߶Π
        Point   p1,   p2; 
        Line()   {} 
        Line(Point   p10,   Point   p20):   p1(p10),   p2(p20)   {} 
}; 
inline   bool   EQ(double   x,   double   y)             //   eqaul,   x   ==   y 
{ 
        return   (fabs(x   -   y)   <   EPS); 
} 
//   ¼ÆËãÁ½Ìõ¶þάֱÏߵĽ»µã£¬½á¹ûÔÚ²ÎÊýPÖзµ»Ø 
//   ·µ»Øֵ˵Ã÷ÁËÁ½ÌõÖ±ÏßµÄλÖùØϵ 
//   COLINE       --   ¹²Ïß 
//   PARALLEL   --   ƽÐР
//   CROSS         --   Ïཻ 
int   CalCrossPoint(Line   L1,   Line   L2,   Point&   P) 
{ 
        double   A1,   B1,   C1,   A2,   B2,   C2; 
        A1   =   L1.p2.y   -   L1.p1.y; 
        B1   =   L1.p1.x   -   L1.p2.x; 
        C1   =   L1.p2.x   *   L1.p1.y   -   L1.p1.x   *   L1.p2.y; 
        A2   =   L2.p2.y   -   L2.p1.y; 
        B2   =   L2.p1.x   -   L2.p2.x; 
        C2   =   L2.p2.x   *   L2.p1.y   -   L2.p1.x   *   L2.p2.y; 
        if   (EQ(A1   *   B2,   B1   *   A2))         { 
                if   (EQ(   (A1   +   B1)   *   C2,   (A2   +   B2)   *   C1   ))   { 
                        return   COLINE; 
                }   else   { 
                        return   PARALLEL; 
                } 
        }   else   { 
                P.x   =   (B2   *   C1   -   B1   *   C2)   /   (A2   *   B1   -   A1   *   B2); 
                P.y   =   (A1   *   C2   -   A2   *   C1)   /   (A2   *   B1   -   A1   *   B2); 
                return   CROSS; 
        }