[51nod]1298 圆与三角形

1298 圆与三角形 题目来源： HackerRank 给出圆的圆心和半径，以及三角形的三个顶点，问圆同三角形是否相交。相交输出"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 <iostream>
#include <stdio.h>
#include <math.h>
using namespace std;

const double eps = 1e-8;
int cmp(double x)
{
if (fabs(x) < eps) return 0;
if (x > 0) return 1;
return -1;
}

struct point
{
double x, y;

point() {}
point(double a, double b):x(a),y(b) {}

void input() {
scanf("%lf%lf", &x, &y);
}
friend point operator - (const point &a, const point &b) {
return point(a.x-b.x, a.y-b.y);
}
double norm() {
return sqrt(x*x + y*y);
}
};

struct cicle
{
point p;
double r;

void input() {
scanf("%lf%lf%lf", &p.x, &p.y, &r);
}
};

double dot(const point &a, const point &b)
{
return a.x*b.x + a.y*b.y;
}

double det(const point &a, const point &b)
{
return a.x*b.y - a.y*b.x;
}

double dist(const point &a, const point &b)
{
return (a-b).norm();
}

double dis_point_segment(const point p, const point s, const point t)
{
if (cmp(dot(p-s,t-s))<0) return (p-s).norm();
if (cmp(dot(p-t,s-t))<0) return (p-t).norm();
return fabs(det(s-p,t-p)/dist(s, t));
}

bool cross(cicle o, point a, point b, point c)
{
double d1, d2, d3;
d1 = dist(o.p, a);
d2 = dist(o.p, b);
d3 = dist(o.p, c);

if (d1<o.r && d2<o.r && d3<o.r)
return false;
if (dis_point_segment(o.p, a, b)>o.r
&& dis_point_segment(o.p, a, c)>o.r
&& dis_point_segment(o.p, b, c)>o.r)
return false;
return true;
}

int main()
{
//freopen("1.txt", "r", stdin);
int T;
scanf("%d", &T);
while (T--) {
cicle o;
o.input();
point a, b, c;
a.input();
b.input();
c.input();
if (cross(o, a, b, c))
printf("Yes\n");
else
printf("No\n");

}

return 0;
}