

第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)
共T行,对于每组输入数据,相交输出"Yes",否则输出"No"。
2 0 0 10 10 0 15 0 15 5 0 0 10 0 0 5 0 5 5
Yes No
转化成三条线段是否与圆相交即可,
当线段两个端点都在圆上时,相交;当一个在圆内,一个在圆外时,相交;当都在圆内,不相交;当都在圆外,需要判断。
都在圆外时,计算出圆心到线段的最短距离与半径比较:
斜率不存在时,比较容易;
斜率存在时,求出线段所在直线的方程,再利用两点间的距离公式列出一个一元二次方程求出圆心到线段每一个点的距离,判断最小距离与半径的关系。
#include<stdio.h>
#include<algorithm>
#include<math.h>
using namespace std;
int t;
double r;
struct point{
double x;
double y;
};
point o,a,b,c;
double dis(point p)
{
return (p.x-o.x)*(p.x-o.x)+(p.y-o.y)*(p.y-o.y);
}
int check(point p,point q)
{
double d1=dis(p);
double d2=dis(q);
//线段两端在圆上
if((d1==r*r)||(d2==r*r))
return 1;
if((d1<r*r&&d2>r*r)||(d1>r*r&&d2<r*r))
return 1;
//线段两端都在圆内
if((d1<r*r)&&(d2<r*r))
return 0;
//斜率不存在
if(p.x==q.x)
{
double maxy=max(p.y,q.y);
double miny=min(p.y,q.y);
if(miny>o.y)
return 0;
if(maxy<o.y)
return 0;
else if(fabs(p.x-o.x)<=r)//点到直线的距离与半径相比较
return 1;
else
return 0;
}
double k,c;
k=(p.y-q.y)/(p.x-q.x);
c=p.y-k*p.x;
//点到直线的距离
double maxx=max(p.x,q.x);
double minx=min(p.x,q.x);
double A,B,C;//化简后的一元二次方程的系数
double K;//对称轴
A=k*k+1;
B=2*k*(c-o.y)-2*o.x;
C=pow(c-o.y,2)+o.x*o.x;
K=-B/2/A;
if(K>maxx)
{
if(A*maxx*maxx+B*maxx + C <= r*r)
return 1;
else
return 0;
}
else if (K < minx)
{
if(A*minx*minx+B*minx+C<= r*r)
return 1;
else
return 0;
}
else
{
if(A*K*K+B*K+C<=r*r)
return 1;
else
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(a,b)==1||check(a,c)==1||check(b,c)==1)
printf("Yes\n");
else
printf("No\n");
}
return 0;
}