23/07/08 10:24:35.14 nQGRHO9D.net
出題おじさんが、題材を提供してくれるからプログラムネタに改造して解くのは楽しいぞ。
統計処理ソフトのオマケ機能を使ってだから証明や数式解を出すのは無理だが。
出題おじさんのおかげで作図に必要な関数のストックが溜まったので他に流用できて( ・∀・)ｲｲ!!
例
# 三点の座標から外接円の中心と半径を返す
Outcircle <- function(A,B,C){ # circumscribing circle
a1=Re(A) ; a2=Im(A)
b1=Re(B) ; b2=Im(B)
c1=Re(C) ; c2=Im(C)
o1 = (a1^2*(b2-c2) + a2^2*(b2-c2) + a2*(-b1^2-b2^2+c1^2+c2^2)+
b1^2*c2+b2^2*c2-b2*c1^2-b2*c2^2)/(2*(a1*(b2-c2)+a2*(c1-b1)+b1*c2-b2*c1))
o2 = (a1^2*(c1-b1)+a1*(b1^2+b2^2-c1^2-c2^2)+a2^2*(c1-b1)-b1^2*c1+b1*c1^2+b1*c2^2-
b2^2*c1)/(2*(a1*(b2-c2)+ a2*(c1-b1)+ b1*c2- b2*c1))
center=o1+1i*o2
radius=abs(center-A)
c(center=center,radius=radius)
}