（湖南省大学生计算机程序设计竞赛）点到圆弧的距离

点到圆弧的距离

Time Limit: 1 Sec   Memory Limit: 128 MB   Special Judge
Submit: 594   Solved: 101
[ Submit][ Status][ Web Board]

Description

输入一个点P和一条圆弧（圆周的一部分），你的任务是计算P到圆弧的最短距离。换句话说，你需要在圆弧上找一个点，到P点的距离最小。
提示：请尽量使用精确算法。相比之下，近似算法更难通过本题的数据。

Input

输入包含最多10000组数据。每组数据包含8个整数x1, y1, x2, y2, x3, y3, xp, yp。圆弧的起点是A(x1,y1)，经过点B(x2,y2)，结束位置是C(x3,y3)。点P的位置是 (xp,yp)。输入保证A, B, C各不相同且不会共线。上述所有点的坐标绝对值不超过20。

Output

对于每组数据，输出测试点编号和P到圆弧的距离，保留三位小数。你的输出和标准输出之间最多能有0.001的误差。

Sample Input

0 0 1 1 2 0 1 -1
3 4 0 5 -3 4 0 1

Sample Output

Case 1: 1.414
Case 2: 4.000

HINT

Source

湖南省第十届大学生计算机程序设计竞赛

关键在于判定射线属于的区域，尝试过用叉积，但是超过180度的角就不能用叉积判定顺时针逆时针关系；

所以采用数量积得出一种判定方法。

#include <iostream>
#include<stdio.h>
#include<string.h>
#include<math.h>
#define dd double
using namespace std;
struct o
{
    double x,y;
} p[4];
int cmp(dd a,dd b)
{
    if(fabs(a-b)<0.000000001)return 1;
    return 0;
}
int ff(double x1,double y1,double x2,double y2)
{
    double xx=x1*y2-x2*y1;
    if(xx>0)return 1;
    if(xx<0)return -1;
    return 0;
}
double di(double x1,double y1,double x2,double y2)
{
    return sqrt((x2-x1)*(x2-x1)+(y2-y1)*(y2-y1));
}
int main()
{
    double a,b;
    int t=1,uu=1;
    freopen("a.in","r",stdin);
    freopen("out.txt","w",stdout);
    while(~scanf("%lf%lf",&p[1].x,&p[1].y))
    {
        scanf("%lf%lf",&p[2].x,&p[2].y);
        scanf("%lf%lf",&p[3].x,&p[3].y);
        scanf("%lf%lf",&a,&b);
        p[1].x+=0.00000001;
        p[1].y+=0.00000001;
        dd xx1=(p[1].x+p[2].x)/2,
           yy1=(p[1].y+p[2].y)/2,
           xx2=(p[1].x+p[3].x)/2,
           yy2=(p[1].y+p[3].y)/2,
           k1=(p[2].y-p[1].y)/(p[2].x-p[1].x),
           k2=(p[3].y-p[1].y)/(p[3].x-p[1].x),
           m=-1.0/k1,n=xx1/k1+yy1,
           p1=-1.0/k2,q=xx2/k2+yy2,
           x0=(q-n)/(m-p1),y0=m*x0+n,
           r=sqrt((p[2].x-x0)*(p[2].x-x0)+(p[2].y-y0)*(p[2].y-y0));
        double vv=di(a,b,p[1].x,p[1].y),
               bb=di(a,b,p[3].x,p[3].y),
               nn=di(x0,y0,a,b)-r,
               mm=di(x0,y0,a,b)+r;
        double dx1=p[1].x-x0,dy1=p[1].y-y0,
               dx2=p[3].x-x0,dy2=p[3].y-y0,
               dx=a-x0,dy=b-y0,
               dxz=p[2].x-x0,dyz=p[2].y-y0,
               miu,lmt,lon=sqrt(dx*dx+dy*dy),X,Y,S1,S2,S3;
        dx=dx*r/lon;
        dy=dy*r/lon;
        X=dx+dxz;
        Y=dy+dyz;
        S2=X*dx+Y*dy;
        S1=X*dx1+Y*dy1;
        S3=X*dx2+Y*dy2;
        int lis=0;
        if(S2>=S1&&S2>=S3)
            lis=1;
        X=dx1+dx2;
        Y=dy1+dy2;
        S2=X*dx1+Y*dy1;
        S1=X*dx+Y*dy;
        S3=X*dxz+Y*dyz;
        if(S2>=S1&&S2>=S3)
            lis=1;
        printf("Case %d: ",uu++);
        if(lis)
        {
            if(nn>=0.0)printf("%.3lf\n",nn);
            else printf("%.3lf\n",-nn);
        }
        else if(vv<bb)printf("%.3lf\n",vv);
        else printf("%.3lf\n",bb);
    }
    return 0;
}