3点を通る円

平面上の3点のデカルト座標が与えられると、それらすべてを通る円の方程式を見つけます。3つのポイントは直線上にありません。

プログラムへの入力の各行には、3つのポイントxとy座標が順番に含まれA(x),A(y),B(x),B(y),C(x),C(y)ます。これらの座標は、スペースで区切られた1,000,000未満の実数になります。

解はの方程式として印刷され(x-h)^2 + (y-k)^2 = r^2ます。値はh、k、r小数点以下3桁で印刷されることになります。方程式のプラス記号とマイナス記号は、数値の前に複数の記号が付かないように、必要に応じて変更する必要があります。

Sample Inputs

7.0 -5.0 -1.0 1.0 0.0 -6.0
1.0 7.0 8.0 6.0 7.0 -2.0

Sample Outputs

(x - 3.000)^2 + (y + 2.000)^2 = 5.000^2
(x - 3.921)^2 + (y - 2.447)^2 = 5.409^2

Pythonの、176の 189文字

import sys,re
for s in sys.stdin:x,y,z=eval(re.sub(r'(\S+) (\S+)',r'\1+\2j,',s));w=z-x;w/=y-x;c=(x-y)*(w-abs(w)**2)/2j/w.imag-x;print'(x%+.3f)^2+(y%+.3f)^2=%.3f^2'%(c.real,c.imag,abs(c+x))

すべての作業を複雑な平面で行います。このページの下から数学に行きます。 -c円の中心です。

C＃-490

using System;class C{static void Main(){Func<string,double>p=s=>double.Parse(s);Func<double,string>t=s=>(s<0?"+ ":"- ")+Math.Abs(s).ToString("F3");foreach(var l in System.IO.File.ReadAllLines("i")){var v=l.Split();double a=p(v[0]),b=p(v[1]),c=p(v[2]),d=p(v[3]),e=p(v[4]),f=p(v[5]),m=(d-b)/(c-a),n=(f-d)/(e-c),x=(m*n*(b-f)+n*(a+c)-m*(c+e))/(2*(n-m)),y=-(x-(a+c)/2)/m+(b+d)/2,r=Math.Sqrt((x-a)*(x-a)+(y-b)*(y-b));Console.WriteLine("(x "+t(x)+")^2+(y "+t(y)+")^2 = "+r.ToString("F3")+"^2");}}}

これにより、ABとBCの間の2行が見つかります。次に、これらの2本の線の二等分線が交差する場所を見つけます。（気づいたのは、@ PeterTaylorが@PeterOfTheCornへのコメントで言及したことです。）

ルビー、192文字

$<.map{|l|a,b,c,d,e,f=l.split.map &:to_f
n=(f-d)/(e-c)
puts"(x%+.3f)^2+(y%+.3f)^2=%.3f^2"%[x=-(n*(a+c)+(n*(b-f)-(c+e))*m=(d-b)/(c-a))/2/n-=m,y=-(x+(a+c)/2)/m-(b+d)/2,((a+x)**2+(b+y)**2)**0.5]}

使用例：

$ echo "7.0 -5.0 -1.0 1.0 0.0 -6.0
1.0 7.0 8.0 6.0 7.0 -2.0" | ruby circle.rb
(x-3.000)^2+(y+2.000)^2=5.000^2
(x-3.921)^2+(y-2.447)^2=5.409^2

ウルフラムアルファ（27）

私は言う、仕事に適切なツールを使用します。

equation circle ([Input1],[Input2]),([Input3],[Input4]),([Input5],[Input6])

ここに例。

Javascript（299）

これを解決することを考えることができる唯一の方法は、3つの未知数に対して3つの方程式を代数的に解いてh、k、rを見つけることでした。

p=prompt().split(' ');a=p[0],b=p[1],c=p[2],d=p[3],e=p[4],f=p[5];h=((a*a+b*b)*(f-d)+(c*c+d*d)*(b-f)+(e*e+f*f)*(d-b))/(a*(f-d)+c*(b-f)+e*(d-b))/2;k=((a*a+b*b)*(e-c)+(c*c+d*d)*(a-e)+(e*e+f*f)*(c-a))/(b*(e-c)+d*(a-e)+f*(c-a))/2;r=Math.sqrt((a-h)*(a-h)+(b-k)*(b-k));alert("(x-"+h+")²+(y-"+k+")²="+r+"²");

I / Oの例：

7.0 -5.0 -1.0 1.0 0.0 -6.0 -> (x-3)²+(y--2)²=5²

1.0 7.0 8.0 6.0 7.0 -2.0 -> (x-3.9210526315789473)²+(y-2.4473684210526314)² =5.409159155551175²

私が見る唯一のバグは、hまたはkが負の場合、の--代わりに出力することです+。