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;}