平面の重心座標系における距離の公式を紹介します。
この記事は重心座標について既に知っている人向けです。
重心座標そのものについての説明は行いません。
重心座標系の基準となる三角形を$\triangle ABC$とし、辺$BC,CA,AB$の長さをそれぞれ$a,b,c$とおきます。
$S_A=\frac{1}{2}(b^2+c^2-a^2),\quad S_B=\frac{1}{2}(c^2+a^2-b^2),\quad \,S_C=\frac{1}{2}(a^2+b^2-c^2)\quad$とします。
重心座標が正規化されている必要がある箇所はその旨を記載します。
記載がなければ正規化された重心座標である必要はありません。
この記事では、座標は原則として$(x_1,x_2,x_3)$や$(p_1,p_2,p_3)$のように添え字$1,2,3$を使用し、3つの数をまとめて$\bm{x}$や$\bm{p}$のように太字で書くことにします。
点$P$の座標のひとつが$\bm{p}$であること(座標$\bm{p}$が指し示す点が$P$であること)を、点$P(\bm{p})$のように書き表します。
座標に対して、ベクトルと同様の演算(和と定数倍)を定義します。
つまり、$k\bm{p}+l\bm{q}=(kp_1+lq_1,\,kp_2+lq_2,\,kp_3+lq_3)$ とします。
2点$P(\bm{p}),Q(\bm{q})$間の距離$PQ$の公式は様々な形があります。(作ろうと思えば無限に作れます。)
まずは私が基本の形だと思っている式から載せます。
$$PQ^2=\frac{a^2(p_2q_3-p_3q_2)^2+b^2(p_3q_1-p_1q_3)^2+c^2(p_1q_2-p_2q_1)^2-2S_A(p_3q_1-p_1q_3)(p_1q_2-p_2q_1)-2S_B(p_1q_2-p_2q_1)(p_2q_3-p_3q_2)-2S_C(p_2q_3-p_3q_2)(p_3q_1-p_1q_3)}{(p_1+p_2+p_3)^2\,(q_1+q_2+q_3)^2}$$
はい、長いですね。
ですが、次のように見ればそれほど難しくはありません。
座標$\bm{p},\bm{q}$をベクトルだと思って、その外積(クロス積,ベクトル積)を $\bm{u}=\bm{p}\times\bm{q}$ とおきます。($\bm{u}$は行ベクトルとします。)
また、行列$\bm{A}$を
$$\bm{A}=\mqty(a^2&-S_C&-S_B\\-S_C&b^2&-S_A\\-S_B&-S_A&c^2)$$
で定義します。
すると、fmlBasicの分子は $\bm{u}\,\bm{A}\,^t\bm{u}$ となります。とても簡単ですね。
なお、$\bm{p},\bm{q}$が正規化された重心座標である場合にはfmlBasicの分母が1になり、よりシンプルになります。
式の簡略化のため、
$T_p=p_1+p_2+p_3,\quad T_q=q_1+q_2+q_3$
および
$u_1=p_2q_3-p_3q_2,\quad u_2=p_3q_1-p_1q_3,\quad u_3=p_1q_2-p_2q_1$
とおいておきます。
すると、fmlBasicの右辺は
$\displaystyle \frac{1}{{T_p}\!^2{T_q}\!^2}\bigl(a^2u_1^2+b^2u_2^2+c^2u_3^2-2S_Au_2u_3-2S_Bu_3u_1-2S_Cu_1u_2\bigr)$
となります。
点$O$をユークリッド平面上の(つまり無限遠点ではない)任意の点とすると、重心座標の性質により、
$\displaystyle \overrightarrow{OP}=\frac{p_1}{p_1+p_2+p_3}\overrightarrow{OA}+\frac{p_2}{p_1+p_2+p_3}\overrightarrow{OB}+\frac{p_3}{p_1+p_2+p_3}\overrightarrow{OC}$ および
$\displaystyle \overrightarrow{OQ}=\frac{q_1}{q_1+q_2+q_3}\overrightarrow{OA}+\frac{q_2}{q_1+q_2+q_3}\overrightarrow{OB}+\frac{q_3}{q_1+q_2+q_3}\overrightarrow{OC}$ が成り立ちます。
$PQ^2$をベクトルで計算すると次のようになります。
\begin{align}
&\sq{PQ}\\
=&\abl\ora{OQ}-\ora{OP}\abr^2\\
=&\Biggl\lvert\biggl(\frac{q_1}{q_1+q_2+q_3}\ora{OA}+\frac{q_2}{q_1+q_2+q_3}\ora{OB}+\frac{q_3}{q_1+q_2+q_3}\ora{OC}\biggr)\\
&\qquad-\biggl(\frac{p_1}{p_1+p_2+p_3}\ora{OA}+\frac{p_2}{p_1+p_2+p_3}\ora{OB}+\frac{p_3}{p_1+p_2+p_3}\ora{OC}\biggr)\Biggr\rvert^2\\
=&\Biggl\lvert\biggl(\frac{q_1}{q_1+q_2+q_3}-\frac{p_1}{p_1+p_2+p_3}\biggr)\ora{OA}\\
&\qquad+\biggl(\frac{q_2}{q_1+q_2+q_3}-\frac{p_2}{p_1+p_2+p_3}\biggr)\ora{OB}\\
&\qquad\qquad+\biggl(\frac{q_3}{q_1+q_2+q_3}-\frac{p_3}{p_1+p_2+p_3}\biggr)\ora{OC}\Biggr\rvert^2\\
=&\Biggl\lvert\frac{q_1(p_1+p_2+p_3)-p_1(q_1+q_2+q_3)}{(p_1+p_2+p_3)(q_1+q_2+q_3)}\ora{OA}\\
&\qquad+\frac{q_2(p_1+p_2+p_3)-p_2(q_1+q_2+q_3)}{(p_1+p_2+p_3)(q_1+q_2+q_3)}\ora{OB}\\
&\qquad\qquad+\frac{q_3(p_1+p_2+p_3)-p_3(q_1+q_2+q_3)}{(p_1+p_2+p_3)(q_1+q_2+q_3)}\ora{OC}\Biggr\rvert^2\\
=&\frac{1}{(p_1+p_2+p_3)^2(q_1+q_2+q_3)^2}\\
&\times\biggl\lvert\bigl((p_3q_1-p_1q_3)-(p_1q_2-p_2q_1)\bigr)\ora{OA}\\
&\ \qquad+\bigl((p_1q_2-p_2q_1)-(p_2q_3-p_3q_2)\bigr)\ora{OB}\\
&\ \qquad\qquad+\bigl((p_2q_3-p_3q_2)-(p_3q_1-p_1q_3)\bigr)\ora{OC}\,\biggr\rvert^2\\
=&\frac{1}{{T_p}\!^2{T_q}\!^2}\biggl\lvert(u_2-u_3)\ora{OA}+(u_3-u_1)\ora{OB}+(u_1-u_2)\ora{OC}\,\biggr\rvert^2\\
=&\frac{1}{{T_p}\!^2{T_q}\!^2}\biggl\{(u_2-u_3)^2\sq{OA}+(u_3-u_1)^2\sq{OB}+(u_1-u_2)^2\sq{OC}\\
&\qquad\quad+2(u_3-u_1)(u_1-u_2)\ip{OB}{OC}+2(u_1-u_2)(u_2-u_3)\ip{OC}{OA}+2(u_2-u_3)(u_3-u_1)\ip{OA}{OB}\ \biggr\}\\
=&\frac{1}{{T_p}\!^2{T_q}\!^2}
\begin{aligned}[t]
\biggl\{\ &\Bigl(\sq{OB}+\sq{OC}-2\ip{OB}{OC}\Bigr)u_1^2\\
+&\Bigl(\sq{OC}+\sq{OA}-2\ip{OC}{OA}\Bigr)u_2^2\\
+&\Bigl(\sq{OA}+\sq{OB}-2\ip{OA}{OB}\Bigr)u_3^2\\
-&\Bigl(2\sq{OA}+2\ip{OB}{OC}-2\ip{OC}{OA}-2\ip{OA}{OB}\Bigr)u_2u_3\\
-&\Bigl(2\sq{OB}-2\ip{OB}{OC}+2\ip{OC}{OA}-2\ip{OA}{OB}\Bigr)u_3u_1\\
-&\Bigl(2\sq{OC}-2\ip{OB}{OC}-2\ip{OC}{OA}+2\ip{OA}{OB}\Bigr)u_1u_2\ \biggr\}
\end{aligned}\\
=&\frac{1}{{T_p}\!^2{T_q}\!^2}
\begin{aligned}[t]
\biggl\{\ &\Bigl(\sq{OB}-2\ip{OB}{OC}+\sq{OC}\Bigr)u_1^2\\
+&\Bigl(\sq{OC}-2\ip{OC}{OA}+\sq{OA}\Bigr)u_2^2\\
+&\Bigl(\sq{OA}-2\ip{OA}{OB}+\sq{OB}\Bigr)u_3^2\\
-&\Bigl(\Bigl(\sq{OC}-2\ip{OC}{OA}+\sq{OA}\Bigr)+\Bigl(\sq{OA}-2\ip{OA}{OB}+\sq{OB}\Bigr)-\Bigl(\sq{OB}-2\ip{OB}{OC}+\sq{OC}\Bigr)\Bigr)u_2u_3\\
-&\Bigl(\Bigl(\sq{OA}-2\ip{OA}{OB}+\sq{OB}\Bigr)+\Bigl(\sq{OB}-2\ip{OB}{OC}+\sq{OC}\Bigr)-\Bigl(\sq{OC}-2\ip{OC}{OA}+\sq{OA}\Bigr)\Bigr)u_3u_1\\
-&\Bigl(\Bigl(\sq{OB}-2\ip{OB}{OC}+\sq{OC}\Bigr)+\Bigl(\sq{OC}-2\ip{OC}{OA}+\sq{OA}\Bigr)-\Bigl(\sq{OA}-2\ip{OA}{OB}+\sq{OB}\Bigr)\Bigr)u_1u_2\ \biggr\}
\end{aligned}\\
=&\frac{1}{{T_p}\!^2{T_q}\!^2}
\begin{aligned}[t]
\biggl\{\ &\abl\ora{OC}-\ora{OB}\abr^2u_1^2+\abl\ora{OA}-\ora{OC}\abr^2u_2^2+\abl\ora{OB}-\ora{OA}\abr^2u_3^2\\
-&\Bigl(\abl\ora{OA}-\ora{OC}\abr^2+\abl\ora{OB}-\ora{OA}\abr^2-\abl\ora{OC}-\ora{OB}\abr^2\Bigr)u_2u_3\\
-&\Bigl(\abl\ora{OB}-\ora{OA}\abr^2+\abl\ora{OC}-\ora{OB}\abr^2-\abl\ora{OA}-\ora{OC}\abr^2\Bigr)u_3u_1\\
-&\Bigl(\abl\ora{OC}-\ora{OB}\abr^2+\abl\ora{OA}-\ora{OC}\abr^2-\abl\ora{OB}-\ora{OA}\abr^2\Bigr)u_1u_2\ \biggr\}
\end{aligned}\\
=&\frac{1}{{T_p}\!^2{T_q}\!^2}
\biggl\{\ \sq{BC}u_1^2+\sq{CA}u_2^2+\sq{AB}u_3^2\\
&\qquad\qquad-\Bigl(\sq{CA}+\sq{AB}-\sq{BC}\Bigr)u_2u_3
-\Bigl(\sq{AB}+\sq{BC}-\sq{CA}\Bigr)u_3u_1
-\Bigl(\sq{BC}+\sq{CA}-\sq{AB}\Bigr)u_1u_2\ \biggr\}\\
=&\frac{1}{{T_p}\!^2{T_q}\!^2}\Bigl(\,a^2u_1^2+b^2u_2^2+c^2u_3^2-\bigl(b^2+c^2-a^2\bigr)u_2u_3-\bigl(c^2+a^2-b^2\bigr)u_3u_1-\bigl(a^2+b^2-c^2\bigr)u_1u_2\,\Bigr)\\
=&\frac{1}{{T_p}\!^2{T_q}\!^2}\bigl(a^2u_1^2+b^2u_2^2+c^2u_3^2-2S_Au_2u_3-2S_Bu_3u_1-2S_Cu_1u_2\bigr)
\end{align}
これでfmlBasicが示されました。
その他の距離の公式は円の式を利用したもので、正規化された重心座標を使います。
まず、正規化された重心座標について説明します。
無限遠点ではない点$P(\bm{p})$は $p_1+p_2+p_3\neq0$ を満たします。
そこで、$\bm{p}'=\frac{1}{p_1+p_2+p_3}\bm{p}$ とおくと $p'_1+p'_2+p'_3=1$ が成り立ちます。
このように、 $p'_1+p'_2+p'_3=1$ を満たす重心座標$\bm{p}'$を正規化された重心座標といいます。
円の式を利用した距離の公式について、まずは一般の円の場合です。
ここで使う円の式は実円または虚円または点円を表すものならば何でもいいので、公式はいくらでも作れることになります。
(実円・虚円・点円などの広義の円に関しては下の方に
付録
を書いたので参考にしてください。)
$d_1,d_2,d_3,e_1,e_2,e_3\in\mathbb{R}$として、$x_1,x_2,x_3$の関数$g$を
$$g(\bm{x})=d_1{x_1}\!^2+d_2{x_2}\!^2+d_3{x_3}\!^2+2e_1x_2x_3+2e_2x_3x_1+2e_3x_1x_2$$
とします。
$\bm{x}$を重心座標として$g(\bm{x})=0$が広義の円を表している場合、
$$\frac{1}{a^2}g(0,1,-1)=\frac{1}{b^2}g(-1,0,1)=\frac{1}{c^2}g(1,-1,0)$$
が成り立つので、この値を$K_g$とおくことにします。
$g(\bm{x})=0$が実円または虚円または点円ならば$K_g\neq0$になります。
$g(\bm{x})=0$は実円または虚円または点円を表す式であるとし、$K_g$を上に書いたように定めます。
$\displaystyle \bm{p}'=\frac{1}{p_1+p_2+p_3}\bm{p}$ および $\displaystyle \bm{q}'=\frac{1}{q_1+q_2+q_3}\bm{q}$ とします。
このとき、距離の公式は次のようになります。
$$PQ^2=\frac{1}{K_g}\,g\bigl(\bm{q}'-\bm{p}'\bigr)$$
式の簡略化のため、
$T_p=p_1+p_2+p_3,\quad T_q=q_1+q_2+q_3$
および
$u_1=p_2q_3-p_3q_2,\quad u_2=p_3q_1-p_1q_3,\quad u_3=p_1q_2-p_2q_1$
とおいておきます。
$g\bigl(\bm{q}'-\bm{p}'\bigr)$を計算すると、
\begin{align}
&g\bigl(\bm{q}'-\bm{p}'\bigr)\\
=&g\biggl(\frac{1}{T_q}\bm{q}-\frac{1}{T_p}\bm{p}\biggr)\\
=&d_1\biggl(\frac{q_1}{T_q}-\frac{p_1}{T_p}\biggr)^2+d_2\biggl(\frac{q_2}{T_q}-\frac{p_2}{T_p}\biggr)^2+d_3\biggl(\frac{q_3}{T_q}-\frac{p_3}{T_p}\biggr)^2\\
&+2e_1\biggl(\frac{q_2}{T_q}-\frac{p_2}{T_p}\biggr)\biggl(\frac{q_3}{T_q}-\frac{p_3}{T_p}\biggr)+2e_2\biggl(\frac{q_3}{T_q}-\frac{p_3}{T_p}\biggr)\biggl(\frac{q_1}{T_q}-\frac{p_1}{T_p}\biggr)+2e_3\biggl(\frac{q_1}{T_q}-\frac{p_1}{T_p}\biggr)\biggl(\frac{q_2}{T_q}-\frac{p_2}{T_p}\biggr)\\
=&\frac{1}{{T_p}\!^2{T_q}\!^2}\Bigl\{d_1\bigl(q_1T_p-p_1T_q\bigr)^2+d_2\bigl(q_2T_p-p_2T_q\bigr)^2+d_3\bigl(q_3T_p-p_3T_q\bigr)^2\\
&\qquad\qquad+2e_1\bigl(q_2T_p-p_2T_q\bigr)\bigl(q_3T_p-p_3T_q\bigr)+2e_2\bigl(q_3T_p-p_3T_q\bigr)\bigl(q_1T_p-p_1T_q\bigr)+2e_3\bigl(q_1T_p-p_1T_q\bigr)\bigl(q_2T_p-p_2T_q\bigr)\Bigr\}\\
=&\frac{1}{{T_p}\!^2{T_q}\!^2}
\begin{aligned}[t]
\Bigl\{\ &d_1\bigl((p_3q_1-p_1q_3)-(p_1q_2-p_2q_1)\bigr)^2\\
+&d_2\bigl((p_1q_2-p_2q_1)-(p_2q_3-p_3q_2)\bigr)^2\\
+&d_3\bigl((p_2q_3-p_3q_2)-(p_3q_1-p_1q_3)\bigr)^2\\
+&2e_1\bigl((p_1q_2-p_2q_1)-(p_2q_3-p_3q_2)\bigr)\bigl((p_2q_3-p_3q_2)-(p_3q_1-p_1q_3)\bigr)\\
+&2e_2\bigl((p_2q_3-p_3q_2)-(p_3q_1-p_1q_3)\bigr)\bigl((p_3q_1-p_1q_3)-(p_1q_2-p_2q_1)\bigr)\\
+&2e_3\bigl((p_3q_1-p_1q_3)-(p_1q_2-p_2q_1)\bigr)\bigl((p_1q_2-p_2q_1)-(p_2q_3-p_3q_2)\bigr)\,\Bigr\}
\end{aligned}\\
=&\frac{1}{{T_p}\!^2{T_q}\!^2}\Bigl\{d_1(u_2-u_3)^2+d_2(u_3-u_1)^2+d_3(u_1-u_2)^2\\
&\qquad\qquad+2e_1(u_3-u_1)(u_1-u_2)+2e_2(u_1-u_2)(u_2-u_3)+2e_3(u_2-u_3)(u_3-u_1)\Bigr\}\\
=&\frac{1}{{T_p}\!^2{T_q}\!^2}\Bigl\{(d_2+d_3-2e_1)u_1^2+(d_3+d_1-2e_2)u_2^2+(d_1+d_2-2e_3)u_3^2\\
&\qquad\qquad+2(-d_1-e_1+e_2+e_3)u_2u_3+2(-d_2-e_2+e_3+e_1)u_3u_1+2(-d_3-e_3+e_1+e_2)u_1u_2\Bigr\}\\
=&\frac{1}{{T_p}\!^2{T_q}\!^2}
\begin{aligned}[t]
\Bigl\{\ &(d_2+d_3-2e_1)u_1^2+(d_3+d_1-2e_2)u_2^2+(d_1+d_2-2e_3)u_3^2\\
-&\bigl((d_3+d_1-2e_2)+(d_1+d_2-2e_3)-(d_2+d_3-2e_1)\bigr)u_2u_3\\
-&\bigl((d_1+d_2-2e_3)+(d_2+d_3-2e_1)-(d_3+d_1-2e_2)\bigr)u_3u_1\\
-&\bigl((d_2+d_3-2e_1)+(d_3+d_1-2e_2)-(d_1+d_2-2e_3)\bigr)u_1u_2\,\Bigr\}
\end{aligned}
\end{align}
となって、$g(\bm{x})=0$が広義の円を表しているときは
$d_2+d_3-2e_1=g(0,1,-1)=a^2K_g$ および
$d_3+d_1-2e_2=g(-1,0,1)=b^2K_g$ および
$d_2+d_3-2e_1=g(1,-1,0)=c^2K_g$ が成り立つので、
\begin{align}
=&\frac{1}{{T_p}\!^2{T_q}\!^2}\Bigl\{a^2K_gu_1^2+b^2K_gu_2^2+c^2K_gu_3^2\\
&\qquad\qquad-\bigl(b^2K_g+c^2K_g-a^2K_g\bigr)u_2u_3-\bigl(c^2K_g+a^2K_g-b^2K_g\bigr)u_3u_1-\bigl(a^2K_g+b^2K_g-c^2K_g\bigr)u_1u_2\Bigr\}\\
=&\frac{K_g}{{T_p}\!^2{T_q}\!^2}\bigl(a^2u_1^2+b^2u_2^2+c^2u_3^2-2S_Au_2u_3-2S_Bu_3u_1-2S_Cu_1u_2\bigr)\\
=&K_g\cdot PQ^2
\end{align}
と式変形されます。(最後の等号はfmlBasicによります。)
$g(\bm{x})=0$が実円または虚円または点円を表す式である場合は$K_g\neq0$なので、$g\bigl(\bm{q}'-\bm{p}'\bigr)=K_g\cdot PQ^2$ の両辺を$K_g$で割ると、fmlGCが得られます。
次に、基準三角形によって決まる円の式を使ってfmlGCの距離の公式の具体例を求めてみます。
$\triangle ABC$の外接円の式は $a^2x_2x_3+b^2x_3x_1+c^2x_1x_2=0$ です。
$K_g$の値は$-1$になります。
$\displaystyle \bm{p}'=\frac{1}{p_1+p_2+p_3}\bm{p}$ および $\displaystyle \bm{q}'=\frac{1}{q_1+q_2+q_3}\bm{q}$ とします。
$$PQ^2=-a^2(q'_2-p'_2)(q'_3-p'_3)-b^2(q'_3-p'_3)(q'_1-p'_1)-c^2(q'_1-p'_1)(q'_2-p'_2)$$
$\triangle ABC$の極円の式は $S_A{x_1}\!^2+S_B{x_2}\!^2+S_C{x_3}\!^2=0$ です。
$K_g$の値は$1$になります。
$\displaystyle \bm{p}'=\frac{1}{p_1+p_2+p_3}\bm{p}$ および $\displaystyle \bm{q}'=\frac{1}{q_1+q_2+q_3}\bm{q}$ とします。
$$PQ^2=S_A(q'_1-p'_1)^2+S_B(q'_2-p'_2)^2+S_C(q'_3-p'_3)^2$$
点$A,B,C$の点円の式はそれぞれ
$c^2{x_2}\!^2+b^2{x_3}\!^2+2S_Ax_2x_3=0$,
$a^2{x_3}\!^2+c^2{x_1}\!^2+2S_Bx_3x_1=0$,
$b^2{x_1}\!^2+a^2{x_2}\!^2+2S_Cx_1x_2=0$
です。
$K_g$の値はいずれも$1$になります。
$\displaystyle \bm{p}'=\frac{1}{p_1+p_2+p_3}\bm{p}$ および $\displaystyle \bm{q}'=\frac{1}{q_1+q_2+q_3}\bm{q}$ とします。
$$PQ^2=c^2(q'_2-p'_2)^2+b^2(q'_3-p'_3)^2+2S_A(q'_2-p'_2)(q'_3-p'_3)$$
$\displaystyle \bm{p}'=\frac{1}{p_1+p_2+p_3}\bm{p}$ および $\displaystyle \bm{q}'=\frac{1}{q_1+q_2+q_3}\bm{q}$ とします。
$$PQ^2=a^2(q'_3-p'_3)^2+c^2(q'_1-p'_1)^2+2S_B(q'_3-p'_3)(q'_1-p'_1)$$
$\displaystyle \bm{p}'=\frac{1}{p_1+p_2+p_3}\bm{p}$ および $\displaystyle \bm{q}'=\frac{1}{q_1+q_2+q_3}\bm{q}$ とします。
$$PQ^2=b^2(q'_1-p'_1)^2+a^2(q'_2-p'_2)^2+2S_C(q'_1-p'_1)(q'_2-p'_2)$$
$\triangle ABC$の1頂点を通り別の1頂点で辺に接する円の式は
$a^2{x_3}\!^2+(a^2-b^2)x_3x_1-c^2x_1x_2=0$,
$b^2{x_1}\!^2+(b^2-c^2)x_1x_2-a^2x_2x_3=0$,
$c^2{x_2}\!^2+(c^2-a^2)x_2x_3-b^2x_3x_1=0$,
$a^2{x_2}\!^2-b^2x_3x_1+(a^2-c^2)x_1x_2=0$,
$b^2{x_3}\!^2-c^2x_1x_2+(b^2-a^2)x_2x_3=0$,
$c^2{x_1}\!^2-a^2x_2x_3+(c^2-b^2)x_3x_1=0$
です。
$K_g$の値はいずれも$1$になります。
$\displaystyle \bm{p}'=\frac{1}{p_1+p_2+p_3}\bm{p}$ および $\displaystyle \bm{q}'=\frac{1}{q_1+q_2+q_3}\bm{q}$ とします。
$$PQ^2=a^2(q'_3-p'_3)^2+(a^2-b^2)(q'_3-p'_3)(q'_1-p'_1)-c^2(q'_1-p'_1)(q'_2-p'_2)$$
$$PQ^2=b^2(q'_1-p'_1)^2+(b^2-c^2)(q'_1-p'_1)(q'_2-p'_2)-a^2(q'_2-p'_2)(q'_3-p'_3)$$
$$PQ^2=c^2(q'_2-p'_2)^2+(c^2-a^2)(q'_2-p'_2)(q'_3-p'_3)-b^2(q'_3-p'_3)(q'_1-p'_1)$$
$$PQ^2=a^2(q'_2-p'_2)^2-b^2(q'_3-p'_3)(q'_1-p'_1)+(a^2-c^2)(q'_1-p'_1)(q'_2-p'_2)$$
$$PQ^2=b^2(q'_3-p'_3)^2-c^2(q'_1-p'_1)(q'_2-p'_2)+(b^2-a^2)(q'_2-p'_2)(q'_3-p'_3)$$
$$PQ^2=c^2(q'_1-p'_1)^2-a^2(q'_2-p'_2)(q'_3-p'_3)+(c^2-b^2)(q'_3-p'_3)(q'_1-p'_1)$$
広義の円という用語は一般には反転幾何でよく用いられます。
しかし、ここでは射影幾何における広義の円を考えるので、混同しないようにしてください。
この記事で説明する広義の円は反転幾何に出てくる広義の円とは少し異なります!
$d_1,d_2,d_3,e_1,e_2,e_3\in\mathbb{R}$として、$x_1,x_2,x_3$の関数$g$を
$$g(\bm{x})=d_1{x_1}\!^2+d_2{x_2}\!^2+d_3{x_3}\!^2+2e_1x_2x_3+2e_2x_3x_1+2e_3x_1x_2$$
とします。つまり$g(\bm{x})$は$x_1,x_2,x_3$の2次斉次多項式です。
座標$\bm{x}$が$g(\bm{x})=0$を満たす点$X(\bm{x})$の集合を広義の円錐曲線(広義の二次曲線)といいます。
行列$\displaystyle \mqty(d_1&e_3&e_2\\e_3&d_2&e_1\\e_2&e_1&d_3)$を用いると広義の円錐曲線$g(\bm{x})=0$を次の5つに分類できます。
i.の実円錐曲線は通常の円錐曲線で楕円・放物線・双曲線があります。
ii.の虚円錐曲線は平面上に図形は現れず空集合になります。ただし、座標が虚数になる虚点というものを考えれば、式を満たす虚点の集合として図形が存在することになります。
iii.の1点は実平面上ではただの1点です。ただし、虚点の範囲まで考えれば複素共役な2本の虚直線になります。
iv.の2直線は、交差する2直線・平行な2直線・通常の直線と無限遠直線を合わせた2直線という3種があります。
v.はiv.の2直線が近づいて同じ直線になったものと思ってください。つまり2本の直線が重なって1直線になったものです。
この記事で説明する広義の円は反転幾何に出てくる広義の円とは少し異なります!
2つの虚円点を通る広義の円錐曲線を(射影幾何における)広義の円といいます。(この名称や定義が一般的に使用されているかどうかはわかりませんが、私はそう呼んでいます。)
その条件を整理すると、$g(\bm{x})$の各項の係数が
$$\frac{d_2+d_3-2e_1}{a^2}=\frac{d_3+d_1-2e_2}{b^2}=\frac{d_1+d_2-2e_3}{c^2}$$
を満たすような広義の円錐曲線$g(\bm{x})=0$が広義の円になります。
広義の円には次の5種類があり、広義の円錐曲線の5分類と対応しています。
i.の実円はいわゆる通常の円です。広義の円に対する言葉としての狭義の円がこの実円になります。
ii.の虚円は平面上に図形は現れず空集合になります。ただし、座標が虚数になる虚点というものを考えれば、式を満たす虚点の集合が存在します。
iii.の点円は平面上の1点です。これは円の半径が0になったものと考えてください。
iv.は円の半径が無限大になったものとみなせます。反転幾何においては直線が広義の円に含まれますが、射影幾何においては無限遠直線を含めた2直線になることに注意してください。
v.はiv.における"別の直線"が"同じ無限遠直線"になったものと思ってください。つまり無限遠直線2つが重なったものです。
$g(\bm{x})=0$が広義の円ならば次の等式も成り立ちます。
$$\frac{1}{a^2}g(a^2,-S_C,-S_B)=\frac{1}{b^2}g(-S_C,b^2,-S_A)=\frac{1}{c^2}g(-S_B,-S_A,c^2)=\frac{S^2}{a^2}g(0,1,-1)=\frac{S^2}{b^2}g(-1,0,1)=\frac{S^2}{c^2}g(1,-1,0)$$
ただし、$S=2\cdot\triangle ABC=\frac{1}{2}\sqrt{2(b^2c^2+c^2a^2+a^2b^2)-(a^4+b^4+c^4)}$ です。
参考文献jkやftでは外接円の式を使用したfmlCircumcircleが距離の公式として載っており、どうやらこれが一般的なようです。
fmlCircumcircle以外にも距離の公式はいろいろと作れるということを紹介したくて今回この記事を書きました。