poj1066，线段上的点能否到另一点+规范相交

#include <iostream>
#include <cmath>
//为什么tm是多组 
using namespace std;

struct Point{
	double x,y;
	Point(double x=0,double y=0):x(x),y(y){
	};
	double operator ^(Point b){return x*b.y-y*b.x;
	}
	Point operator -(Point b){return Point(x-b.x,y-b.y);
	}
};

struct Seg{
	Point a,b;
	Seg(int a,int b):a(a),b(b){
	};
	Seg(){
	};
};

const int N=40;
Seg seg[N];
Point point[N<<1];

int  n;
Point target;
int cnt,minn;


int sgn(double x){//判正负，防止丢失精度 
	const double neg=1e-8;
	if(fabs(x)<neg) return 0;
	if(x>0) return 1;
	else return -1;
}


bool isIntersect(Point ax,Point ay,Point bx,Point by){//判断线段相交 ，非严格相交，严格相交去= 
	return
		max(ax.x,ay.x)>=min(bx.x,by.x)&&
		max(bx.x,by.x)>=min(ax.x,ay.x)&&
		max(ax.y,ay.y)>=min(bx.y,by.y)&&
		max(bx.y,by.y)>=min(ax.y,ay.y)&&
		sgn((ax-bx)^(bx-by))*sgn((ay-bx)^(bx-by))<=0&&
		sgn((bx-ax)^(ax-ay))*sgn((by-ax)^(ax-ay))<=0;
}


void solve(){
	minn=1<<10;
	
	for(int i=0;i<n<<1;i++){//枚举所以线段（墙）的端点 
		cnt=0;
		for(int j=0;j<n;j++){//枚举所有墙 
			if(i!=j<<1&&i!=j<<1|1&&isIntersect(target,point[i],seg[j].a,seg[j].b)) cnt++;
		}//自己墙的端点不能匹配自己墙 ，否则就要严格相交 
		minn=min(minn,cnt);
	}
	if(n==1||n==0) minn=0;//如果只有一个墙或者没有墙那么在内部肯定是不要穿墙的，但会跳循环，minn可能为无穷大，所以手动 
	
	cout<<"Number of doors = "<<minn+1<<endl;//+最外面的墙 
}


int main(int argc, char** argv) {
	ios::sync_with_stdio(false);
	
	while(cin>>n){
		for(int i=0;i<n;i++){
			cin>>seg[i].a.x>>seg[i].a.y>>seg[i].b.x>>seg[i].b.y;//输入墙的信息 
			point[i<<1]=seg[i].a;//输入墙端点的信息 
			point[i<<1|1]=seg[i].b;
		}
		cin>>target.x>>target.y;//输入目标的坐标 
		
		solve();	
	}
	
	
	return 0;
}