Rectangle and Circle（HDU-1221）

---

1 题意

  在直角坐标系下给定一个圆和一个边与坐标轴平行的矩形，判断两个图形的边是否相交。
  链接：link。

2 思路

  命题1：两个图形边相交的充要条件是，圆心与矩形边的最短距离小于等于圆的半径，且最大距离大于等于圆的半径。
  这里暂不证明上述命题，只做简单解释，第一个条件保证与圆最近的一条边与圆相交或相切或者在园内，第二个条件保证园没有包含整个矩形。

  下面介绍如何求点到与坐标轴平行的直线的最短距离和最大距离。
  最短距离：

1. 对于一个点和与x轴平行的线段的最短距离，如果点在线段的左侧就是点与线段左端点的距离；如果在右侧那么最短距离就是点与右端点的距离；如果在中间那么就是点与线段的纵坐标之差。
2. 相反对于y轴平行的线段也做类似的分析。

  最大距离：即点到线段端点距离的最大值。

2.1 时间复杂度分析

  需要计算点到直线的距离八次，所以时间复杂度为$\mathcal{O}(1)$。

2.2 实现

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
double dis(double x1,double y1,double x2,double y2){
    return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
}
double dismin(double x,double y,double x1,double y1,double x2,double y2){
    if(x1==x2){
        if(y1<y2) swap(y1,y2);
        if(y>y1) return dis(x1,y1,x,y);
        else if(y<y2) return dis(x1,y2,x,y);
        else return fabs(x-x1);
    }
    else if(y1==y2){
        if(x1<x2) swap(x1,x2);
        if(x>x1) return dis(x1,y1,x,y);
        else if(x<x2) return dis(x2,y1,x,y);
        else return fabs(y-y1);
    }
}
double dismax(double x,double y,double x1,double y1,double x2,double y2){
    return max(dis(x1,y1,x,y),dis(x2,y2,x,y));
}
int main(){
    int T;cin>>T;
    while(T--){
        double x,y,r,x1,y1,x2,y2,dmin[4],dmax[4];
        cin>>x>>y>>r>>x1>>y1>>x2>>y2;
        dmin[0]=dismin(x,y,x1,y1,x1,y2);
        dmin[1]=dismin(x,y,x1,y1,x2,y1);
        dmin[2]=dismin(x,y,x2,y1,x2,y2);
        dmin[3]=dismin(x,y,x1,y2,x2,y2);
        dmax[0]=dismax(x,y,x1,y1,x1,y2);
        dmax[1]=dismax(x,y,x1,y1,x2,y1);
        dmax[2]=dismax(x,y,x2,y1,x2,y2);
        dmax[3]=dismax(x,y,x1,y2,x2,y2);
        if(*max_element(dmax,dmax+4)>=r&&*min_element(dmin,dmin+4)<=r) cout<<"YES"<<endl;
        else cout<<"NO"<<endl;
    }
    return 0;
}