ZOJ 1648 Circuit Board(判断线段相交)

#include<stdio.h>
#include<math.h>
#include<string.h>
#define max(a,b) a>b?a:b
#define min(a,b) a<b?a:b
#define N 2001
#define esp 0.000000001
struct point
{
double x,y;
};
struct Line
{
point a,b;
}line
;
double xmult(point a1,point b1,point a2,point b2)
{
return (a1.x-b1.x)*(a2.y-b2.y)-(a2.x-b2.x)*(a1.y-b1.y);
}
int quick_judge(int i,int j)
{
point tmp;
int flag_x,flag_y;
//保证每条 line 中 a 为较小的点(在x轴上较左,在y轴上较下)
if(line[i].a.x>line[i].b.x) tmp=line[i].a,line[i].a=line[i].b,line[i].b=tmp;
if(line[j].a.x>line[j].b.x) tmp=line[j].a,line[j].a=line[j].b,line[j].b=tmp;
//两条线段的 左端点中的较大值 <= 右端点中的较小值，就可以 x 轴上判断相交
double max_left=max(line[i].a.x,line[j].a.x),min_right=min(line[i].b.x,line[j].b.x);
if(max_left > min_right) flag_x=0;
else flag_x=1;
if(line[i].a.y>line[i].b.y) tmp=line[i].a,line[i].a=line[i].b,line[i].b=tmp;
if(line[j].a.y>line[j].b.y) tmp=line[j].a,line[j].a=line[j].b,line[j].b=tmp;
double max_down=max(line[i].a.y,line[j].a.y),min_up=min(line[i].b.y,line[j].b.y);
if(max_down > min_up) flag_y=0;
else flag_y=1;
if(flag_x && flag_y)return 1;  //矩阵相交,不排斥,返回1
else return 0;
}
int judge_cross(int i,int j)
{
if(!quick_judge(i,j))return 0;//排斥，即矩阵不相交，那么线段一定不相交，返回0
if(xmult(line[i].a,line[j].a,line[j].b,line[j].a)*xmult(line[i].b,line[j].a,line[j].b,line[j].a)>0)return 0;
return 1;
}
int main()
{
int i,j,t;
while(~scanf("%d",&t))
{
memset(line,0,sizeof(line));
for(i=0;i<t;i++)
scanf("%lf%lf%lf%lf",&line[i].a.x,&line[i].a.y,&line[i].b.x,&line[i].b.y);
if(t<=1)
{
printf("ok!/n");continue;
}
int cross=0;
for(i=0;i<t;i++)
{
for(j=0;j<t;j++)
{
if(j==i)continue;
if(judge_cross(i,j))
{
cross=1;break;
}
}
if(cross)break;
}
if(!cross)printf("ok!/n");
else printf("burned!/n");
}
return 0;
}