1298 圆与三角形

1298 圆与三角形 

题目来源：  HackerRank

基准时间限制：1 秒 空间限制：131072 KB 分值: 0  难度：基础题

 收藏

 关注

给出圆的圆心和半径，以及三角形的三个顶点，问圆同三角形是否相交。相交输出"Yes"，否则输出"No"。（三角形的面积大于0）。

Input

第1行：一个数T，表示输入的测试数量(1 <= T <= 10000)，之后每4行用来描述一组测试数据。
4-1：三个数，前两个数为圆心的坐标xc, yc，第3个数为圆的半径R。(-3000 <= xc, yc <= 3000, 1 <= R <= 3000）
4-2：2个数，三角形第1个点的坐标。
4-3：2个数，三角形第2个点的坐标。
4-4：2个数，三角形第3个点的坐标。(-3000 <= xi, yi <= 3000）

Output

共T行，对于每组输入数据，相交输出"Yes"，否则输出"No"。

Input示例

2
0 0 10
10 0
15 0
15 5
0 0 10
0 0
5 0
5 5

Output示例

Yes
No

#include<stdio.h>
#include<math.h>
int t;
double r;
struct point{
double x;
double y;
};
point o,a,b,c;
int length(point m,point n)
{
if((m.x-n.x)*(m.x-n.x)+(m.y-n.y)*(m.y-n.y)>r*r*1.0)
return 1;//在圆外 
return 0; 
}
double fork(point o,point m,point n)//向量om，on叉乘   
{  
    m.x-=o.x;  
    m.y-=o.y;  
    n.x-=o.x;  
    n.y-=o.y;  
    return m.x*n.y-m.y*n.x;  
}  
double lenth(point m,point n)
{
return sqrt((m.x-n.x)*(m.x-n.x)+(m.y-n.y)*(m.y-n.y)); 
}
double h(point o,point m,point p)
{
double len_mp=lenth(m,p);  
    point t_mp={(p.x-m.x),(p.y-m.y)};//向量mp  
     point t_ok={-t_mp.y,t_mp.x};//mp的法向量
      
    //下面利用叉乘判断最近距离的垂点是否在线段上  
    if(fork({0,0},t_ok,{p.x-o.x,p.y-o.y})*fork({0,0},t_ok,{m.x-o.x,m.y-o.y})>0)  
    {  
        return lenth(o,m)<lenth(o,p)?lenth(o,m):lenth(o,p);  
    }  
    double s=fabs(fork(o,m,p));//平行四边形面积   
    return s/len_mp;
}
int check()
{
//都在圆内
if(length(o,a)==0&&length(o,b)==0&&length(o,c)==0)
return 0;
//有在圆外的，也有在圆内的 
if(length(o,a)==0||length(o,b)==0||length(o,c)==0)
return 1;
//都在圆外
if(h(o,a,b)<=r||h(o,a,c)<=r||h(o,b,c)<=r)
return 1;
return 0;

}
int main()
{
scanf("%d",&t);
while(t--)
{
scanf("%lf %lf %lf",&o.x,&o.y,&r);
scanf("%lf %lf %lf %lf %lf %lf",&a.x,&a.y,&b.x,&b.y,&c.x,&c.y);
if(check())
printf("Yes\n");
else
printf("No\n");
}
return 0;
 }