计算几何 线段的各种情况 [好题]

队内比赛又悲剧了，咱们队各种不在状态啊~ 我悲剧的计算几何就这么悲催咯~ 不过是个好题，可以用来做模板了。线段相交分为几种情况，1不相交，2,重叠，3相交端点在线上，4相交端点不在线上。

为什么还是这么弱啊~~ 和我不想编代码有关系！加油啊~狐狸

#include<stdio.h>
#include<cmath>
#include<string.h>
#define eps 1e-7

struct node
{
       int x,y;
};

struct segment
{
       node node1,node2;
}line[101];

int sign( double x )
{
    if( fabs(x)<eps )
    return 0;
    else if( x>0 )
    return 1;
    else return -1;
}

int xmult( node p,node a,node b ){
       return sign( (a.x-p.x)*(b.y-p.y) - (b.x-p.x)*(a.y-p.y) );
}

bool xmultSegment( segment a,segment b )
{
     if( (xmult( a.node1,a.node2,b.node1 )^xmult( a.node1,a.node2,b.node2 ))==-2  )
     {    
          if( (xmult( b.node1,b.node2,a.node1 )^xmult( b.node1,b.node2,a.node2 ))==-2  )
             return true;
     }
     return false;
}

double dist( node p1,node p2 )
{
       return sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y));
}

bool examine( node a,node b,node c )
{
     if( fabs(dist(a,c)+dist(c,b)-dist(a,b))<eps )
         return true;
     return false;
}
bool SamePoint( node a,node b,node c )
{
     if( a.x==c.x && a.y==c.y )
         return true;
     if( b.x==c.x && b.y==c.y )
         return true;
     return false;
}
bool exSame( segment a,segment b )
{
     int index=0;
     if( SamePoint( a.node1,a.node2,b.node1 ) )
         index++;
     if( SamePoint( a.node1,a.node2,b.node2 ) )
         index++;
     if( SamePoint( b.node1,b.node2,a.node1 ) )
         index++;
     if( SamePoint( b.node1,b.node2,a.node2 ) )
         index++;
     if( index==0 )
         return false;
     else
         return true;
}

int inLine( segment a,segment b )
{
    int index=0;
     if( xmult( a.node1,a.node2,b.node1 )==0 && examine( a.node1,a.node2,b.node1 ) )
         index++;
     if( xmult( a.node1,a.node2,b.node2 )==0 && examine( a.node1,a.node2,b.node2 ) )
         index++;
     if( xmult( b.node1,b.node2,a.node1 )==0 && examine( b.node1,b.node2,a.node1 ) )
         index++;
     if( xmult( b.node1,b.node2,a.node2 )==0 && examine( b.node1,b.node2,a.node2 ) )
         index++;
    return index;
}

int main()
{
    int T,n;
    int i,j;
    scanf( "%d",&T );
    while( T-- )
    {
           int ans[5];
           memset( ans,0,sizeof(ans) );
           scanf( "%d",&n );
           for( i=1;i<=n;i++ )
                scanf( "%d %d %d %d",&line[i].node1.x,&line[i].node1.y,&line[i].node2.x,&line[i].node2.y );           
           for( i=1;i<=n;i++ )
           for( j=i+1;j<=n;j++ )
           {
                if( xmultSegment( line[i],line[j] ) || inLine( line[i],line[j] ) )
                {
                    if( inLine( line[i],line[j] )>=2 && !exSame( line[i],line[j] ) )
                        ans[2]++;
                    else if( inLine( line[i],line[j] )>=1 )
                        ans[3]++;
                    else
                        ans[4]++;
                }
                else
                    ans[1]++;
           }
           printf( "%d\n%d\n%d\n%d\n\n",ans[1],ans[2],ans[3],ans[4] );
    }
    return 0;
}