给定三角形三个顶点的坐标,如何求三角形的外心的坐标呢?
例如 :给定a(x1,y1) b(x2,y2) c(x3,y3)求外接圆心坐标O(x,y)
1. 首先,外接圆的圆心是三角形三条边的垂直平分线的交点,我们根据圆心到顶点的距离相等,可以列出以下方程:
(x1-x)(x1-x)-(y1-y)(y1-y)=(x2-x)(x2-x)+(y2-y)(y2-y);
(x2-x)(x2-x)+(y2-y)(y2-y)=(x3-x)(x3-x)+(y3-y)(y3-y);
2.化简得到:
2*(x2-x1)x+2(y2-y1)y=x2^2+y2^2-x1^2-y1^2;
2*(x3-x2)x+2(y3-y2)y=x3^2+y3^2-x2^2-y2^2;
令A1=2*(x2-x1);
B1=2*(y2-y1);
C1=x2^2+y2^2-x1^2-y1^2;
A2=2*(x3-x2);
B2=2*(y3-y2);
C2=x3^2+y3^2-x2^2-y2^2;
即
A1*x+B1y=C1;
A2*x+B2y=C2;
3.最后根据克拉默法则:
x=((C1*B2)-(C2*B1))/((A1*B2)-(A2*B1));
y=((A1*C2)-(A2*C1))/((A1*B2)-(A2*B1));
因此,x,y为最终结果;
对于空间中的三角形,只不过最后解方程组的时候是三元方程组
Ps:克拉默法则可以用向量积和数量积的方法证明,也可以用高等代数的知识证明
可用此知识点解决的问题
http://acm.hdu.edu.cn/showproblem.php?pid=4720点击打开链接
给定三角形三个顶点的坐标,如何求三角形的外心的坐标呢?
例如 :给定a(x1,y1) b(x2,y2) c(x3,y3)求外接圆心坐标O(x,y)
首先,外接圆的圆心是三角形三条边的垂直平分线的交点,我们根据圆心到顶点的距离相等,可以列出以下方程:
(x1-x)(x1-x)-(y1-y)(y1-y)=(x2-x)(x2-x)+(y2-y)(y2-y);
(x2-x)(x2-x)+(y2-y)(y2-y)=(x3-x)(x3-x)+(y3-y)(y3-y);
2.化简得到:
2*(x2-x1)*x+2*(y2-y1)y=x2^2+y2^2-x1^2-y1^2;
2*(x3-x2)*x+2*(y3-y2)y=x3^2+y3^2-x2^2-y2^2;
令A1=2*(x2-x1);
B1=2*(y2-y1);
C1=x2^2+y2^2-x1^2-y1^2;
A2=2*(x3-x2);
B2=2*(y3-y2);
C2=x3^2+y3^2-x2^2-y2^2;
即
A1*x+B1y=C1;
A2*x+B2y=C2;
3.最后根据克拉默法则:
x=((C1*B2)-(C2*B1))/((A1*B2)-(A2*B1));
y=((A1*C2)-(A2*C1))/((A1*B2)-(A2*B1));
因此,x,y为最终结果;
对于空间中的三角形,只不过最后解方程组的时候是三元方程组
时间限制:1.0s 内存限制:256.0MB
提交此题
问题描述
为二维空间中的点设计一个结构体,在此基础上为三角形设计一个结构体。分别设计独立的函数计算三角形的周长、面积、中心和重心。输入三个点,输出这三个点构成的三角形的周长、面积、外心和重心。结果保留小数点后2位数字。
样例输出
与上面的样例输入对应的输出。
例:
数据规模和约定
输入数据中每一个数的范围。
例:doule型表示数据。
#include <iostream>
#include <stdio.h>
#include <iomanip>
#include <cmath>
#include <cstring>
#include <string>
#include<iostream>
#include<algorithm>
#include<cstring>
#include <cstdio>
using namespace std;
int x[4],y[4];
double zc()
{
double a=(x[1]-x[2])*(x[1]-x[2])+(y[1]-y[2])*(y[1]-y[2]);
double b=(x[2]-x[3])*(x[2]-x[3])+(y[2]-y[3])*(y[2]-y[3]);
double c=(x[3]-x[1])*(x[3]-x[1])+(y[3]-y[1])*(y[3]-y[1]);
return sqrt(a)+sqrt(b)+sqrt(c);
}
double mj()
{
double a=(x[1]-x[2])*(x[1]-x[2])+(y[1]-y[2])*(y[1]-y[2]);
double b=(x[2]-x[3])*(x[2]-x[3])+(y[2]-y[3])*(y[2]-y[3]);
double c=(x[3]-x[1])*(x[3]-x[1])+(y[3]-y[1])*(y[3]-y[1]);
a=sqrt(a);
b=sqrt(b);
c=sqrt(c);
double s=(a+b+c)/2;
s=s*(s-a)*(s-b)*(s-c);
return sqrt(s);
}
void wx(double &a,double &b)
{
double A1=2*(x[2]-x[1]);
double B1=2*(y[2]-y[1]);
double C1=(x[2]*x[2]+y[2]*y[2]-x[1]*x[1]-y[1]*y[1]);
double A2=2*(x[3]-x[2]);
double B2=2*(y[3]-y[2]);
double C2=x[3]*x[3]+y[3]*y[3]-x[2]*x[2]-y[2]*y[2];
a=((C1*B2)-(C2*B1))/((A1*B2)-(A2*B1));
b=((A1*C2)-(A2*C1))/((A1*B2)-(A2*B1));
}
void zx(double &c,double &d)
{
c=(x[1]+x[2]+x[3])/3.0;
d=(y[1]+y[2]+y[3])/3.0;
}
int main()
{
for(int i=1;i<=3;i++)
cin>>x[i]>>y[i];
cout<< setiosflags(ios::fixed)<< setiosflags(ios::right)<<setprecision(2);
cout<<zc()<<endl;
cout<<mj()<<endl;
double a,b;
wx(a,b);
cout<<a<<' '<<b<<endl;
double c,d;
zx(c,d);
cout<<c<<' '<<d<<endl;
}