51Nod 1298 圆与三角形

给出圆的圆心和半径，以及三角形的三个顶点，问圆同三角形是否相交。相交输出”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

计算几何模板真难打啊！
cin还卡了一组样例的时间

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#define eps 1e-8
#define IO ios::sync_with_stdio(false);\
cin.tie(0);\
cout.tie(0);
using namespace std;
struct point
{
    double x,y;
};
typedef struct point point;
double xmult(point p1,point p2,point p0)
{
    return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);
}
double dis(point p1,point p2)
{
    return sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y));
}
double disptoline(point p,point l1,point l2)
{
    return fabs(xmult(p,l1,l2)/dis(l1,l2));
}
int intersect_seg_cricle(point c,double r,point l1,point l2)
{
    double t1=dis(c,l1)-r,t2=dis(c,l2)-r;
    point t=c;
    if(t1<eps||t2<eps)
        return t1>-eps||t2>-eps;
    t.x+=l1.y-l2.y;
    t.y+=l2.x-l1.x;
    return xmult(l1,c,t)*xmult(l2,c,t)<eps&&disptoline(c,l1,l2)-r<eps;
}
int main()
{
    IO;
    int T;
    cin>>T;
    while(T--)
    {
        double r;
        point p,p1,p2,p3;
        cin>>p.x>>p.y>>r;
        cin>>p1.x>>p1.y;
        cin>>p2.x>>p2.y;
        cin>>p3.x>>p3.y;
        int flag1=intersect_seg_cricle(p,r,p1,p2);
        int flag2=intersect_seg_cricle(p,r,p2,p3);
        int flag3=intersect_seg_cricle(p,r,p1,p3);
        if(flag1||flag2||flag3)
            cout<<"Yes"<<endl;
        else
            cout<<"No"<<endl;
    }
}