HDU4643 GSM 计算几何-二分

题目链接：HDU4643

题目大意：给你n个城市，m个基站，从城市a到城市b，路径为直线，手机总是接受离自己位置最近的基站的信号，问从城市a到城市b中会有几次基站改变的情况。

思路：利用二分的思想，求离a b 最近的基站 ，他们连线的中垂线与线段ab的交点，有几个交点就是改变几次。

代码借鉴：大佬的代码

不大好理解的地方加了注释。

AC代码：

/*
HDU 4643
简单计算几何 二分求交点
2017年8月3日10:52:57 
TLE: ~scanf()读到文件末尾 漏写
AC代码 
*/

#include<stdio.h>
#include<algorithm>
#include<string.h>
#include<math.h>
using namespace std;
const int inf=0x3f3f3f3f;
const double eps=1e-7;
const int maxn=55;
struct Point{
	double x,y;
	Point(){}
	Point(double _x,double _y){
		x=_x;
		y=_y;
	}
}city[maxn],base[maxn];

int n,m;

double dis(Point a,Point b){
	return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}

int belong(Point p){
	int k;//记录距离p最近的基站标号 
	double min_dis=inf;
	for(int i=1;i<=m;i++){
		if(dis(p,base[i])<min_dis){
			min_dis=dis(p,base[i]);
			k=i;
		}
	}
	return  k; 
}

int find(Point l,Point r){
	int k1=belong(l),k2=belong(r);
	//距离两个端点 最近的基站相同，（l,r）区间内交换次数为0 
	if(k1==k2) return 0;
	//找到了两个基站在这条线段上的交点，交换一次 
	if(dis(l,r)<eps) return 1; 
	Point m=Point((l.x+r.x)/2,(l.y+r.y)/2) ;
	return find(l,m)+find(m,r);
}

int main(){
	while(~scanf("%d%d",&n,&m)){
		for(int i=1;i<=n;i++) 
			scanf("%lf%lf",&city[i].x,&city[i].y);
		for(int i=1;i<=m;i++){
			scanf("%lf%lf",&base[i].x,&base[i].y);
		}
		int q;
		scanf("%d",&q);
		while(q--){
			int a,b;
			scanf("%d%d",&a,&b);
			printf("%d\n",find(city[a],city[b]));
		}
	}
	return 0;
	
}