pku2242The Circumference of the Circle

题目链接：pku2242

题意：已知三角形的顶点坐标，求三角形的外接圆的周长。

方法：已知三角形三边a,b,c，可以用海伦公式求其面积。S=sqrt(p*(p-a)*(p-b)*(p-c))，其中p=(a+b+c)/2.

知道a,b,c及外接圆半径R,则三角形的面积为 S=a*b*c/(4*R)。本题由此公式即可算出外接圆半径。

代码一：

#include <iostream> #include <cmath> using namespace std; const double pi=3.141592653589793; struct point {double x, y;}; double dist(point a,point b) { return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y)); } int main() { point p[3]; double a,b,c,pp,s,R; while(cin>>p[0].x>>p[0].y>>p[1].x>>p[1].y>>p[2].x>>p[2].y) { a=dist(p[0],p[1]); b=dist(p[1],p[2]); c=dist(p[2],p[0]); pp=(a+b+c)/2; s=sqrt(pp*(pp-a)*(pp-b)*(pp-c)); R=a*b*c/(4*s); printf("%.2lf/n",2*pi*R); } }

代码二：用模板

#include <iostream> #include <cmath> using namespace std; const double pi=3.141592653589793; struct point {double x, y;}; typedef point triangle[3]; typedef struct { point centre; double r; } circle; double dist(point a,point b) { return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y)); } double triangle_area(triangle t) { return fabs(t[0].x*t[1].y+t[1].x*t[2].y+t[2].x*t[0].y-t[1].x*t[0].y-t[2].x*t[1].y-t[0].x*t[2].y)/2; } //已知三角形的坐标，求其外接圆的圆心，半径 circle circle_of_triangle(triangle t) { circle tmp; double a,b,c,c1,c2,ax,bx,cx,ay,by,cy; a=dist(t[0],t[1]); b=dist(t[1],t[2]); c=dist(t[2],t[0]); tmp.r=a*b*c/triangle_area(t)/4; ax=t[0].x; ay=t[0].y; bx=t[1].x; by=t[1].y; cx=t[2].x; cy=t[2].y; c1=(ax*ax+ay*ay-bx*bx-by*by)/2; c2=(ax*ax+ay*ay-cx*cx-cy*cy)/2; tmp.centre.x=(c1*(ay-cy)-c2*(ay-by))/ ((ax-bx)*(ay-cy)-(ax-cx)*(ay-by)); tmp.centre.y=(c1*(ax-cx)-c2*(ax-bx))/ ((ay-by)*(ax-cx)-(ay-cy)*(ax-bx)); return tmp; } int main() { triangle t; circle c; while(cin>>t[0].x>>t[0].y>>t[1].x>>t[1].y>>t[2].x>>t[2].y) { c=circle_of_triangle(t); printf("%.2lf/n",2*pi*c.r); } return 0; }