R^2−2Rr=d^2 (オイラーの定理)
2r^2(R^2+d^2)=(R^2−d^2)^2 (フースの定理)
フースは双心五角形,六角形,七角形,八角形に関する同様の公式も見つけているが,今回のコラムでは,双心五角形の場合を扱ってみたい.
===================================
【1】アルゴリズム
外接円:x^2+y^2=R^2
上の2点を(xk,yk),(xk+1,yk+1)とおく.
xk^2+yk^2=R^2
xk+1^2+yk+1^2=R^2
この2点を結ぶ直線:
y−yk=mk(x−xk),mk=(yk+1−yk)/(xk+1−xk)
が,
内接円:(x−d)^2+y^2=r^2
に接することから,判別式=0とおいて,
{d+mk(mkxk−yk)}^2−(1+mk^2){d^2+(mkxk−yk)^2−r^2}=0が得られる.
初期値(x0,y0)=(R,0)から始まって,(x1,y1),(x2,y2),・・・が満たす連立方程式を求める.
xk+1=d−rとなるときの(R,r,d)の関係式を求める(k=1,n=5).d=0とおいたとき,
r/R=cos(π/5)
であれば,正解が得られていることが確認できる.
===================================
【2】グレブナー基底
d^6−2d^4rR+8d^2r^3R−3d^4R^2−4d^2r^2R^2+4d^2rR^3+3d^2R^4+4r^2R^4−2rR^5−R^6=0
d=0とおくと
4r^2R^4−2rR^5−R^6=0
となるが,
r/R=cos(π/5)=(1+√5)/4
はこれを満たす.
===================================