NYOJ 83-迷宫寻宝（二）（计算几何->线段相交）

迷宫寻宝（二）

时间限制： 1000 ms  |  内存限制： 10000 KB

难度： 5

描述

一个叫ACM的寻宝者找到了一个藏宝图，它根据藏宝图找到了一个迷宫，这是一个很特别的迷宫，迷宫是一100*100的个正方形区域，里面有很多墙，这些墙都是由一些直线构成的，如下图。

墙把迷宫分隔成很多藏宝室，任何两个藏宝室之间都没有门。

ACM现在准备用开凿设备在相邻两个藏宝室的墙中间凿开一个门，最终取出迷宫中的宝物。

但是，开凿门是一件很费力费时的工作，ACM想开凿尽量少的门取出宝物，现在请你写一个程序来帮助它判断一下最少需要开几个门吧。

输入

第一行输入一个正数N(N<10)表示测试数据组数
每组测试数据的第一行是一个整数n(0<=n<=30),代表了墙的个数，随后的n行里每行有四个整数x1,x2,y1,y2，这四个数分别是代表一个墙的两个端点的坐标。外围的正方形四个顶点固定在(0,0)(0,100)(100,0)(100,100)这四堵个墙不在上面的n个数里。注意，不能在两个线的交点处开凿门。
数据保证任意两个中间墙的交点不在四周的墙上。
输完所有的墙后，输入两个数，x,y（可能不是整数），表示宝藏的坐标。

输出

输出最少需要开凿的门的个数

样例输入

1
7 
20 0 37 100 
40 0 76 100 
85 0 0 75 
100 90 0 90 
0 71 100 61 
0 14 100 38 
100 47 47 100 
54.5 55.4 

样例输出

2

来源

East Central North America 1999 改编

上传者

张云聪

题意：题目翻译过来就是让你在矩形边界上找一点，使得该点与宝藏的位置组成的线段与其他线段相交最少，求最小值。

题解：我们回顾一下计算几何的基础知识，即线段相交问题：

一般线段的空间位置关系分三种情况：

（1）两线段不相交

（2）两线段完全相交

（3）一条线段的端点在另一条线段上（也包含重合）

我们本题解决的只是第二种情况，即不含有线段端点在另一线段上的情况（代码中的online函数可以忽略掉！！）呢线段的完全相交怎么判断呢，我们回顾一下叉积的知识，叉积的其中一个很重要的作用是判断一点在线段或直线的左侧还是右侧，或者在线段上。这便是我们用来判断的非常重要的条件，也就是已知两端点的情况下，枚举第三个点和第四个点，假若第3和第4这两点分别在已知线段的两侧，则这两条线段必相交。。。OK，问题解决了。

#include<map>          
#include<stack>          
#include<queue>          
#include<vector>          
#include<math.h>          
#include<stdio.h>          
#include<iostream>          
#include<string.h>          
#include<stdlib.h>  
#include<algorithm> 
#include<functional>  
using namespace std;          
typedef long long  ll;          
#define inf  1000000000         
#define MOD 1000000007           
#define  maxn  510005
#define  lowbit(x) (x&-x)          
#define  eps 1e-9
int n;
struct node
{
	double x,y;
}a[105],b[105],c,d;
double Cross(node p1,node p2,node p3)
{//叉积判断一点在这条线段的左侧还是右侧
	return (p2.x-p1.x)*(p3.y-p2.y)-(p3.x-p2.x)*(p2.y-p1.y);
}
bool online(node p1,node p2,node p3)
{//判断线段的一个端点是否在另一条线段上
	if(p3.x>=min(p1.x,p2.x) && p3.x<=max(p1.x,p2.x) && p3.y>=min(p1.y,p2.y) && p3.y<=max(p1.y,p2.y))
		return 1;
	return 0;
}
bool judge(node p1,node p2,node p3,node p4)
{
	double d1=Cross(p1,p2,p3);
	double d2=Cross(p1,p2,p4);
	double d3=Cross(p3,p4,p1);
	double d4=Cross(p3,p4,p2);
	if(d1*d2<0 && d3*d4<0)
		return 1;
	if(fabs(d1)<=eps && online(p1,p2,p3)) return 1;
	if(fabs(d2)<=eps && online(p1,p2,p4)) return 1;
	if(fabs(d3)<=eps && online(p3,p4,p1)) return 1;
	if(fabs(d4)<=eps && online(p3,p4,p2)) return 1;
	return 0;
}
int work(node p)
{
	int i,sum=0;
	for(i=1;i<=n;i++)
		if(judge(c,p,a[i],b[i]))
			sum++;
	return sum;
}
int main(void)
{
	int T,i,j,ans,tmp;
	scanf("%d",&T);
	while(T--)
	{
		ans=100000000;
		scanf("%d",&n);
		for(i=1;i<=n;i++)
			scanf("%lf%lf%lf%lf",&a[i].x,&a[i].y,&b[i].x,&b[i].y);
		scanf("%lf%lf",&c.x,&c.y);
		for(i=1;i<100;i++)
		{
			d.x=0;d.y=i;
			tmp=work(d);
			ans=min(ans,tmp);
			d.x=100;d.y=i;
			tmp=work(d);
			ans=min(ans,tmp);
			d.x=i;d.y=0;
			tmp=work(d);
			ans=min(ans,tmp);
			d.x=i;d.y=100;
			tmp=work(d);
			ans=min(ans,tmp);
			//printf("%d\n",ans);
		}
		printf("%d\n",ans+1);
	}
	return 0;
}