hdu15558(线段相交判定）

#include <iostream>
#include <cstdio>
#include <cstring>

using namespace std;
#define SIZE 1001

struct point{

       double x,y;

}line[SIZE][2];

int num[SIZE];
int p[SIZE];


inline bool across(int i,int j)
{
       double x1,x2,x3,y1,y2,y3;
       x1=line[i][0].x-line[j][0].x;
       y1=line[i][0].y-line[j][0].y;
       x2=line[j][1].x-line[j][0].x;
       y2=line[j][1].y-line[j][0].y;
       x3=line[i][1].x-line[j][0].x;
       y3=line[i][1].y-line[j][0].y;
       return (x1*y2-x2*y1)*(x2*y3-x3*y2)>=0;
}

inline bool isRecIntersect(int i,int j)
{
       return max(line[i][0].x,line[i][1].x)>=min(line[j][0].x,line[j][1].x) &&
           max(line[i][0].y,line[i][1].y)>=min(line[j][0].y,line[j][1].y);
       }

inline bool isIntersect(int i,int j)
{
    return isRecIntersect(i,j)&&isRecIntersect(j,i)&&across(i,j)&&across(j,i);

}

int find(int x)
{
    return x == p[x] ? x:p[x] = find(p[x]);
}


void merge(int x,int y)
{
     int fx=find(x),fy=find(y);

     if(fx!=fy&&isIntersect(x,y))
     {
        p[fx] = fy;
        num[fy] += num[fx];
     }
}

void init(int n) {

    for(int i = 1; i  <= n; i++) {

        num[i] = 1;
        p[i] = i;
    }
}

int main()
{
    int cas;
    cin>>cas;
    while(cas--)
    {

        int c;
        int i;
        int n=1;
        cin>>c;
        init(c);
        while(c--)
        {
            char cmd[2];
            cin>>cmd;
            if(cmd[0]=='P')
            {
                cin>>line[n][0].x>>line[n][0].y>>line[n][1].x>>line[n][1].y;
                for(i=1;i<n;i++)
                    merge(i,n);
                n++;
            }
            else
            {
                cin>>i;
                cout<<num[find(i)]<<endl;
            }
        }
        if(cas)//if not the last cas
            puts("");
    }

    return 0;
}