poj1329解题报告

赤裸裸的求三角形外接圆，即根据三点坐标确定圆心坐标和圆的半径。

圆的半径=a*b*c/(4*s)(a,b,c分别为三角形的三边长，s是三角形的面积)。

s=a*b*sinA/2

2*r*sinA=a   (画图，构造直角三角形，同弧对同圆周角可证)

由上两式则可推出半径的公式。

 

圆心坐标。设三角形的三个顶点坐标为（ax,ay,(bx,by),(cx,cy);圆心坐标为(rx,ry).

由圆心到三顶点距离相等，可以得到：

(rx-ax)^2+(ry-ay)^2=(rx-bx)^2+(ry-by)^2

(rx-ax)^2+(rx-ay)^2=(rx-cx)^2+(ry-cy)^2

 

 

以上两式消去ry,化简得：

C1(ay-cy)-C2(ay-by)+(ay-by)(cx-ax)rx-(ay-cy)(bx-ax)rx=0

其中C1=ax^2+ay^2-(bx^2+by^2);C2=ax^2+ay^2-(cx^2+cy^2).

同理消去rx，可以解出ry.

 

下面的代码中，Distance原为distance,但是G++中已经定义了distance，导致重名，编译报错。

感谢此文http://hi.baidu.com/%BF%C6%B9%A4ers/blog/item/9b4f13111fd9527acb80c4e9.html解惑。

#include<iostream> #include<stdio.h> #include<string.h> #include<math.h> using namespace std; const double eps=1e-6; struct POINT { double x; double y; }; struct TRIANGLE { POINT a; POINT b; POINT c; }; struct CIRCLE { POINT center; double r; }; double Distance(POINT p1,POINT p2) { return sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y)); } double areaOftriangle(TRIANGLE t) { return fabs(t.a.x*t.b.y+t.b.x*t.c.y+t.c.x*t.a.y-t.b.x*t.a.y-t.c.x*t.b.y-t.a.x*t.c.y)/2; } CIRCLE circleOftriangle(TRIANGLE t) { CIRCLE tmp; double a,b,c; a=Distance(t.a,t.b); b=Distance(t.b,t.c); c=Distance(t.a,t.c); tmp.r=a*b*c/areaOftriangle(t)/4; double ax,ay,bx,by,cx,cy; ax=t.a.x;ay=t.a.y; bx=t.b.x;by=t.b.y; cx=t.c.x;cy=t.c.y; double c1,c2; c1=(ax*ax+ay*ay-bx*bx-by*by)/2; c2=(ax*ax+ay*ay-cx*cx-cy*cy)/2; tmp.center.x=-(c1*(ay-cy)-c2*(ay-by))/ ((ax-bx)*(ay-cy)-(ax-cx)*(ay-by)); tmp.center.y=-(c1*(ax-cx)-c2*(ax-bx))/ ((ay-by)*(ax-cx)-(ay-cy)*(ax-bx)); return tmp; } int main() { TRIANGLE t; CIRCLE circle; while(scanf("%lf%lf%lf%lf%lf%lf",&t.a.x,&t.a.y,&t.b.x,&t.b.y,&t.c.x,&t.c.y)!=EOF) { circle=circleOftriangle(t); if(fabs(circle.center.x)<eps)printf("x^2"); else if(circle.center.x<0)printf("(x - %.3lf)^2 + ",-circle.center.x); else if(circle.center.x>0)printf("(x + %.3lf)^2 + ",circle.center.x); if(fabs(circle.center.y)<eps)printf("y^2"); else if(circle.center.y<0)printf("(y - %.3lf)^2 = ",-circle.center.y); else if(circle.center.y>0)printf("(y + %.3lf)^2 = ",circle.center.y); printf("%.3f^2/n",circle.r); double c=2*circle.center.x; double d=2*circle.center.y; double e=circle.center.x*circle.center.x+circle.center.y*circle.center.y-circle.r*circle.r; printf("x^2 + y^2"); if(c<0)printf(" - %.3lfx",-c); else printf(" + %.3lfx",c); if(d<0)printf(" - %.3lfy",-d); else printf(" + %.3lfy",d); if(e<0)printf(" - %.3lf = 0",-e); else printf(" + %.3lf = 0",e); printf("/n/n"); } return 0; }