hdu4454 计算几何 三分

传送门

题意：给定一个圆和一个矩形、一个点，求点先到圆再到矩形的最短距离和

思路：点是定点，实际上是求圆上的点到这个定点和矩形的最小距离和，所以遍历圆上的点，三分就可以了，在求与矩形的最小的距离的时候，分成八部分（可以说十分特别麻烦了T T）

完整代码：

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <cmath>
using namespace std;
const double eps=1e-9;
const double PI=acos(-1.0);
struct point
{
    double x,y;
}o,p,px,p1,p2,p3,p4;
double r;
double dis(point p1,point p2)
{
    return (p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y);
}
double dis2(point p1,point p2)
{
    return sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y));
}
double mindis(point pp)
{
    if(pp.x<=p1.x&&pp.y>=p2.y)
        return dis2(pp,p3);
    if(pp.x>=p1.x&&pp.x<=p2.x&&pp.y>=p2.y)
        return pp.y-p2.y;
    if(pp.x>=p2.x&&pp.y>=p2.y)
        return dis2(pp,p2);
    if(pp.x<=p1.x&&pp.y<=p1.y)
        return dis2(pp,p1);
    if(pp.x>=p1.x&&pp.x<=p2.x&&pp.y<=p1.y)
        return p1.y-pp.y;
    if(pp.x>=p2.x&&pp.y<=p1.y)
        return dis2(pp,p4);
    if(pp.y>=p1.y&&pp.y<=p2.y&&pp.x<=p1.x)
        return p1.x-pp.x;
    if(pp.y>=p1.y&&pp.y<=p2.y&&pp.x>=p2.x)
        return pp.x-p2.x;
}
double Function(double alpha)
{
    point pp;
    pp.x=r*cos(alpha)+o.x;
    pp.y=r*sin(alpha)+o.y;
    return mindis(pp)+dis2(pp,p);
}
double trisection(double l,double r)
{
    while(r-l>eps)
    {
        double ll=l+(r-l)/3;
        double rr=l+2*(r-l)/3;
        double xl=Function(ll);
        double xr=Function(rr);
        if(xl<xr)
            r=rr;
        else l=ll;
    }
    return r;
}
int main()
{
    while(scanf("%lf%lf",&p.x,&p.y)!=-1)
    {
        if(p.x==0&&p.y==0) break;
        scanf("%lf%lf%lf",&o.x,&o.y,&r);
        double t1,t2,t3,t4;
        scanf("%lf%lf%lf%lf",&t1,&t2,&t3,&t4);
        p1.x=min(t1,t3);p1.y=min(t2,t4);
        p2.x=max(t1,t3);p2.y=max(t2,t4);
        p3.x=p1.x;p3.y=p2.y;
        p4.x=p2.x;p4.y=p1.y;
        double left=0;
        double right=2*PI;
        double res=trisection(left,right);
        double ans=Function(res);
        printf("%.2f\n",ans);
    }
    return 0;
}