7-3 求圆周长

求圆周长

在平面直角坐标系内，给定不共线的三个点的坐标，求经过这3点的那个圆的周长。

输入格式:

输入包含多组测试数据。每个测试占一行，包含6个实数x1,y1,x2,y2,x3,y3，代表了平面上3个点的坐标。题目保证由这3个点决定的那个圆的直径不会超过10^6。

输出格式:

对于每组测试数据，在一行上输出计算的结果。圆的周长精确到百分位。圆周率的近似值取为3.141592653589793。

输入样例:

0.0 -0.5 0.5 0.0 0.0 0.5
0.0 0.0 0.0 1.0 1.0 1.0
5.0 5.0 5.0 7.0 4.0 6.0
0.0 0.0 -1.0 7.0 7.0 7.0
50.0 50.0 50.0 70.0 40.0 60.0
0.0 0.0 10.0 0.0 20.0 1.0
0.0 -500000.0 500000.0 0.0 0.0 500000.0

输出样例:

3.14
4.44
6.28
31.42
62.83
632.24
3141592.65

解题思路：

首先需要明确圆的标准方程为(x-a)² + (y-b)² = r²，其中(a , b)为圆心坐标，已知圆上三个点的坐标为(x1 , y1), (x2 , y2), (x3 , y3),分别带入方程中得：

(x₁ - a)² + (y₁ - b)² = r²    ①
(x₂ - a)² + (y₂ - b)² = r²    ②
(x₃ - a)² + (y₃ - b)² = r²	  ③

由① - ②， ① - ③整理得④、⑤：

(x₁ - x₂)a + (y₁ - y₂)b = [(x₁² - x₂²) + (y₁² - y₂²)] / 2	④
(x₁ - x₃)a + (y₁ - y₃)b = [(x₁² - x₃²) + (y₁² - y₃²)] / 2	⑤

该二元二次方程有解的条件是：

(x₁ - x₂) / (y₁ - y₂) ≠ (x₁ - x₃) / (y₁ - y₃)

定义几个数

A = [(x₁² - x₂²) + (y₁² - y₂²)] / 2
B = [(x₁² - x₃²) + (y₁² - y₃²)] / 2
C = (x₁ - x₂)
D = (y₁ - y₂)
E = (x₁ - x₃)
F = (y₁ - y₃)

即得圆心坐标：

a = (AF -BD) / (CF - DE)
b = (AE - BC) / (DE - CF)

半径 r 为圆心到圆上任意一点的距离，用两点间距离公式进行计算

d = √[ (x₁ - x₂)²+(y₁ - y₂)² ]

得到半径后计算周长C = 2Πr

参考代码：

#include<iostream>
#include<cmath>
#include<iomanip>
using namespace std;
double calculate_round(double x1, double y1, double x2, double y2, double x3, double y3)
{
	double A = ((pow(x1, 2) - pow(x2, 2)) + (pow(y1, 2) - pow(y2, 2))) / 2;
	double B = ((pow(x1, 2) - pow(x3, 2)) + (pow(y1, 2) - pow(y3, 2))) / 2;
	double C = (x1 - x2);
	double D = (y1 - y2);
	double E = (x1 - x3);
	double F = (y1 - y3);
	double a = (A*F - B*D) / (C*F - D*E);
	double	b = (A*E - B*C) / (D*E - C*F);
	double rout = fabs( ( sqrt( pow(  a - x1 , 2 ) + pow(  b - y1 , 2 ) ) ) );
	double round = rout * 2 * 3.141592653589793;
	return round;
}
void excute()
{
	double x1[1000] = { 0 }, y1[1000] = { 0 }, x2[1000] = { 0 }, y2[1000] = { 0 }, x3[1000] = { 0 }, y3[1000] = { 0 }, round[1000] = { 0 };
	int group = 0;
	while (cin >> x1[group] >> y1[group] >> x2[group] >> y2[group] >> x3[group] >> y3[group])
	{
		round[group] = calculate_round(x1[group], y1[group], x2[group], y2[group], x3[group], y3[group]);
		group++;
	}
	for (int i=0;i<group;i++)
	{
		if (i!=group-1)
		{
			cout << fixed << setprecision(2) << round[i] << endl;
		}
		else
		{
			cout << round[i];
		}
	}
}
int main()
{
	excute();
	return 0;
}