三角形外心

double x1,Y1,x2,y2,x3,y3;
struct vec
{
    double a,b,c;
    vec(double x,double y,double z):a(x),b(y),c(z){};
};
vec pro(double x1,double Y1,double x2,double y2)
{
    return (vec){x2-x1,y2-Y1,-((x2-x1)*(x2-x1)+(y2-Y1)*(y2-Y1))/2};
}
int main()
{
    while(cin>>x1>>Y1>>x2>>y2>>x3>>y3)
    {
        vec a=pro(x1,Y1,x2,y2),b=pro(x1,Y1,x3,y3);
        double cx=-(b.c*a.b-a.c*b.b)/(b.a*a.b-b.b*a.a);
        double cy=-(b.c*a.a-a.c*b.a)/(b.b*a.a-a.b*b.a);
        double r=sqrt(cx*cx+cy*cy);//半径
        cx+=x1,cy+=Y1;//外心横坐标,纵坐标
        double c=-2*cx,d=-2*cy,e=cx*cx+cy*cy-r*r;//圆一般方程参数
    }
    return 0;
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值