読者です 読者をやめる 読者になる 読者になる

美しき行列式たち 〜 第七回

行列式

美しき行列式たち, 最終回です. 連載にはタグ『行列式』をつけています. 今日はL字型の行列式を扱います. L字型と言えば, 第一回に紹介したL same (仮名)ですね. あれは簡単に証明できました. しかし, 今回の難易度はマックスですよっ♪

18. L increasing power (仮名)

\begin{align} \left| \begin{array}{cccccc} x^0 & x^1 & x^4 & x^9 & \cdots & x^{n^2-2n+1} \\ x^3 & x^2 & x^5 & x^{10} & \vdots & x^{n^2-2n+2} \\ x^8 & x^7 & x^6 & x^{11} & \vdots & x^{n^2-2n+3} \\ x^{15} & x^{14} & x^{13} & x^{12} & \vdots & x^{n^2-2n+4} \\ \vdots & \cdots & \cdots & \cdots & \ddots & \vdots \\ x^{n^2-1} & x^{n^2-2} & x^{n^2-3} & x^{n^2-4} & \cdots & x^{n^2-n} \\ \end{array} \right| = x^{\frac{1}{3}\left(n^3-n\right)} \prod_{k=1}^{n-1}{\left(1-x^{2k}\right)} \end{align}
L字型の行列式は美しいです. そして指数の部分が順々に増えていく感じの行列がこちら. 答えはx多項式, しかも積和となっています. という事は, この行列式=0の解は複素平面上で考えると云々...
(証明)
\begin{align} \left| \left( a_{ij} \right) \right| &= \left| \begin{array}{cccccc} x^0 & x^1 & x^4 & x^9 & \cdots & x^{n^2-2n+1} \\ x^3 & x^2 & x^5 & x^{10} & \vdots & x^{n^2-2n+2} \\ x^8 & x^7 & x^6 & x^{11} & \vdots & x^{n^2-2n+3} \\ x^{15} & x^{14} & x^{13} & x^{12} & \vdots & x^{n^2-2n+4} \\ \vdots & \cdots & \cdots & \cdots & \ddots & \vdots \\ x^{n^2-1} & x^{n^2-2} & x^{n^2-3} & x^{n^2-4} & \cdots & x^{n^2-n} \\ \end{array} \right| \\ &= \left| \begin{array}{cccccc} x^0-x^2 & x^1 & x^4 & x^9 & \cdots & x^{n^2-2n+1} \\ 0 & x^2 & x^5 & x^{10} & \vdots & x^{n^2-2n+2} \\ 0 & x^7 & x^6 & x^{11} & \vdots & x^{n^2-2n+3} \\ 0 & x^{14} & x^{13} & x^{12} & \vdots & x^{n^2-2n+4} \\ \vdots & \cdots & \cdots & \cdots & \ddots & \vdots \\ 0 & x^{n^2-2} & x^{n^2-3} & x^{n^2-4} & \cdots & x^{n^2-n} \\ \end{array} \right|\nonumber \\ &= \left(x^0-x^2\right) \left| \begin{array}{ccccc} x^2 & x^5 & x^{10} & \cdots & x^{n^2-2n+2} \\ x^7 & x^6 & x^{11} & \vdots & x^{n^2-2n+3} \\ x^{14} & x^{13} & x^{12} & \vdots & x^{n^2-2n+4} \\ \vdots & \cdots & \cdots & \ddots & \vdots \\ x^{n^2-2} & x^{n^2-3} & x^{n^2-4} & \cdots & x^{n^2-n} \\ \end{array} \right|\nonumber \\ &= \left(x^0-x^2\right) \left| \begin{array}{ccccc} x^2-x^6 & x^5 & x^{10} & \cdots & x^{n^2-2n+2} \\ 0 & x^6 & x^{11} & \vdots & x^{n^2-2n+3} \\ 0 & x^{13} & x^{12} & \vdots & x^{n^2-2n+4} \\ \vdots & \cdots & \cdots & \ddots & \vdots \\ 0 & x^{n^2-3} & x^{n^2-4} & \cdots & x^{n^2-n} \\ \end{array} \right|\nonumber \\ &= \left(x^0-x^2\right)\left(x^2-x^6\right) \left| \begin{array}{cccc} x^6 & x^{11} & \cdots & x^{n^2-2n+3} \\ x^{13} & x^{12} & \vdots & x^{n^2-2n+4} \\ \vdots & \cdots & \ddots & \vdots \\ x^{n^2-3} & x^{n^2-4} & \cdots & x^{n^2-n} \\ \end{array} \right|\nonumber \\ &= \left(x^0-x^2\right)\left(x^2-x^6\right) \left| \begin{array}{cccc} x^6-x^{12} & x^{11} & \cdots & x^{n^2-2n+3} \\ 0 & x^{12} & \vdots & x^{n^2-2n+4} \\ \vdots & \cdots & \ddots & \vdots \\ 0 & x^{n^2-4} & \cdots & x^{n^2-n} \\ \end{array} \right|\nonumber \\ &= \left(x^0-x^2\right)\left(x^2-x^6\right)\left(x^6-x^{12}\right) \left| \begin{array}{ccc} x^{12} & \cdots & x^{n^2-2n+4} \\ \vdots & \ddots & \vdots \\ x^{n^2-4} & \cdots & x^{n^2-n} \\ \end{array} \right|\nonumber \\ &\vdots \nonumber \\ &= \left(x^0-x^2\right)\left(x^2-x^6\right)\cdots\left(x^{n^2-3n+2}-x^{n^2-n}\right)x^{n^2-n} \nonumber \\ &= \left[\prod_{k=1}^{n-1}{\left(x^{k^2-k}-x^{k^2+k}\right)}\right]x^{n^2-n} \nonumber \\ &= x^{\frac{1}{3}\left(n^3-n\right)} \prod_{k=1}^{n-1}{\left(1-x^{2k}\right)} \end{align}
(証明終)
証明は左の列から順々に, 二行目のx倍を引いて外に出していきました. ええ, これはまだ簡単なんですよ...

19. L increasing number (仮名)

\begin{align} \left| \begin{array}{cccccc} 0 & 1 & 4 & 9 & \cdots & n^2-2n+1 \\ 3 & 2 & 5 & 10 & \vdots & n^2-2n+2 \\ 8 & 7 & 6 & 11 & \vdots & n^2-2n+3 \\ 15 & 14 & 13 & 12 & \vdots & n^2-2n+4 \\ \vdots & \cdots & \cdots & \cdots & \ddots & \vdots \\ n^2-1 & n^2-2 & n^2-3 & n^2-4 & \cdots & n^2-n \\ \end{array} \right| = (-2)^{n-2}(n-1)!\left[\sum_{k=1}^{n-1}\frac{1}{k}-2n(n-1)\right] \end{align}
L字型の最終形. 先程の行列の指数部がそのまんま行列の要素になった形です. とても難しいです. もし, 腕に自信があるならば, 今直ぐブラウザーのタブを閉じて計算を始めてください. 計算の基本は『いっぱい0を作る』『\sim 5\times 5くらい小さな行列で様子を見る』です.
(証明)
{\small\begin{eqnarray}  \left|\left(a_{ij}\right)\right|  &=& \left| \begin{array}{cccccccc} 0           & 1           & 4           & 9           & 16          & \cdots & (n-2)^2     & (n-1)^2 \\ 3           & 2           & 5           & 10          & 17          & \vdots & (n-2)^2+1   & (n-1)^2+1 \\ 8           & 7           & 6           & 11          & 18          & \vdots & (n-2)^2+2   & (n-1)^2+2 \\ 15          & 14          & 13          & 12          & 19          & \vdots & (n-2)^2+3   & (n-1)^2+3 \\ 24          & 23          & 22          & 21          & 20          & \vdots & (n-2)^2+4   & (n-1)^2+4 \\ \vdots      & \cdots      & \cdots      & \cdots      & \cdots      & \ddots & \vdots      & \vdots \\ (n-1)^2-1   & (n-1)^2-2   & (n-1)^2-3   & (n-1)^2-4   & (n-1)^2-5   & \cdots & (n-1)(n-2)  & n^2-n-1 \\ n^2-1       & n^2-2       & n^2-3       & n^2-4       & n^2-5       & \cdots & n(n-1)+1    & n(n-1) \\ \end{array} \right| \\ &=& \left| \begin{array}{cccccccc} 0           & 1           & 4           & 9           & 16          & \cdots & (n-2)^2     & (n-1)^2 \\ 3           & 1           & 1           & 1           & 1           & \vdots & 1           & 1       \\ 8           & 6           & 2           & 2           & 2           & \vdots & 2           & 2       \\ 15          & 13          & 9           & 3           & 3           & \vdots & 3           & 3       \\ 24          & 22          & 18          & 12          & 4           & \vdots & 4           & 4       \\ \vdots      & \cdots      & \cdots      & \cdots      & \cdots      & \ddots & \vdots      & \vdots  \\ (n-1)^2-1   & (n-1)^2-3   & (n-1)^2-7   & (n-1)^2-13  & (n-1)^2-21  & \cdots & n-2       & n-2     \\ n^2-1       & n^2-3       & n^2-7       & n^2-13      & n^2-21      & \cdots & 3(n-1)      & n-1     \\ \end{array} \right| \\ &=& \left| \begin{array}{ccccccccc} -(n-1)^2    & -n(n-2)     & -(n+1)(n-3) & -(n+2)(n-4) & -(n+3)(n-5) & \cdots & -(2n-4)2    & -(2n-3)     & (n-1)^2 \\ 2           & 0           & 0           & 0           & 0           & \vdots & 0           & 0           & 1       \\ 6           & 4           & 0           & 0           & 0           & \vdots & 0           & 0           & 2       \\ 12          & 10          & 6           & 0           & 0           & \vdots & 0           & 0           & 3       \\ 20          & 18          & 14          & 8           & 0           & \vdots & 0           & 0           & 4       \\ \vdots      & \cdots      & \cdots      & \cdots      & \cdots      & \ddots & \vdots      & \vdots      & \vdots  \\ (n-1)(n-2)  & n(n-3)      & (n+1)(n-4)  & (n+2)(n-5)  & (n+3)(n-6)  & \cdots & 2n-4        & 0           & n-2     \\ n(n-1)      & (n+1)(n-2)  & (n+2)(n-3)  & (n+3)(n-4)  & (n+4)(n-5)  & \cdots & (2n-3)2     & 2n-2        & n-1     \\ \end{array} \right| \\ &=& \left| \begin{array}{ccccccccc} n-1         & n-2         & n-3         & n-4         & n-5         & \cdots & 2           & 1           & n(n-1)  \\ 2           & 0           & 0           & 0           & 0           & \vdots & 0           & 0           & 1       \\ 6           & 4           & 0           & 0           & 0           & \vdots & 0           & 0           & 2       \\ 12          & 10          & 6           & 0           & 0           & \vdots & 0           & 0           & 3       \\ 20          & 18          & 14          & 8           & 0           & \vdots & 0           & 0           & 4       \\ \vdots      & \cdots      & \cdots      & \cdots      & \cdots      & \ddots & \vdots      & \vdots      & \vdots  \\ (n-1)(n-2)  & n(n-3)      & (n+1)(n-4)  & (n+2)(n-5)  & (n+3)(n-6)  & \cdots & 2n-4        & 0           & n-2     \\ n(n-1)      & (n+1)(n-2)  & (n+2)(n-3)  & (n+3)(n-4)  & (n+4)(n-5)  & \cdots & (2n-3)2     & 2n-2        & n-1     \\ \end{array} \right| \\ &=& \left| \begin{array}{ccccccccc} 0           & 0           & 0           & 0           & 0           & \cdots & 0           & 1           & n(n-1)  \\ 2           & 0           & 0           & 0           & 0           & \vdots & 0           & 0           & 1       \\ 6           & 4           & 0           & 0           & 0           & \vdots & 0           & 0           & 2       \\ 12          & 10          & 6           & 0           & 0           & \vdots & 0           & 0           & 3       \\ 20          & 18          & 14          & 8           & 0           & \vdots & 0           & 0           & 4       \\ \vdots      & \cdots      & \cdots      & \cdots      & \cdots      & \ddots & \vdots      & \vdots      & \vdots  \\ (n-1)(n-2)  & n(n-3)      & (n+1)(n-4)  & (n+2)(n-5)  & (n+3)(n-6)  & \cdots & 2n-4        & 0           & n-2     \\ -(n-2)(n-1) & -(n-3)(n-2) & -(n-4)(n-3) & -(n-5)(n-4) & -(n-6)(n-5) & \cdots & -2          & 2n-2        & n-1     \\ \end{array} \right| \\ &\vdots \\ &=& \left| \begin{array}{cccccccc} 0      & 0      & 0      & 0      &  \cdots  & 0      & 1      & n(n-1)  \\ 2      & 0      & 0      & 0      &  \vdots  & 0      & 0      & 1       \\ -6     & 4      & 0      & 0      &  \vdots  & 0      & 0      & 2       \\ 0      & -8     & 6      & 0      &  \vdots  & 0      & 0      & 3       \\ 0      & 0      & -10    & 8      &  \ddots  & 0      & 0      & 4       \\ \vdots & \cdots & \cdots & \ddots &  \ddots  & \ddots & \vdots & \vdots  \\ 0      & 0      & 0      & \cdots &  -2(n-1) & 2n-4   & 0      & n-2     \\ 0      & 0      & 0      & \cdots &  0       & -2     & 2n-2   & n-1     \\ \end{array} \right| \\ &=& (-1)^{n-1} n(n-1)  \prod_{k=1}^{n-1} 2k + (-1)^n \left| \begin{array}{ccccccc} 2      & 0      & 0      & 0      &  \vdots  & 0      & 1       \\ -6     & 4      & 0      & 0      &  \vdots  & 0      & 2       \\ 0      & -8     & 6      & 0      &  \vdots  & 0      & 3       \\ 0      & 0      & -10    & 8      &  \ddots  & \vdots & 4       \\ \vdots & \cdots & \cdots & \ddots &  \ddots  & 0      & \vdots  \\ 0      & 0      & 0      & \cdots &  -2(n-1) & 2n-4   & n-2     \\ 0      & 0      & 0      & \cdots &  0       & -2     & n-1     \\ \end{array} \right| \end{eqnarray} }
最後の\;\vdots\;の所では, n-2列目を1 \sim n-3列目からうまく引き, 次にn-3列目を1 \sim n-4列目からうまく引き, ... 2列目を1列目から引いていくことで, 下に0が並ぶようにした.
さて, ここで新たに出てきたn-1 \times n-1行列式を計算する.
\begin{eqnarray} &&\left| \begin{array}{ccccccc} 2      & 0      & 0      & 0      &  \vdots  & 0      & 1       \\ -6     & 4      & 0      & 0      &  \vdots  & 0      & 2       \\ 0      & -8     & 6      & 0      &  \vdots  & 0      & 3       \\ 0      & 0      & -10    & 8      &  \ddots  & \vdots & 4       \\ \vdots & \cdots & \cdots & \ddots &  \ddots  & 0      & \vdots  \\ 0      & 0      & 0      & \cdots &  -2(n-1) & 2n-4   & n-2     \\ 0      & 0      & 0      & \cdots &  0       & -2     & n-1     \\ \end{array} \right| \\ &=& (n-1) \prod_{k=1}^{n-2}2k \\ &-& (n-2) \left(\prod_{k=1}^{n-3}2k\right) (-2) \\ &+& (n-3) \left(\prod_{k=1}^{n-4}2k\right) (-2(n-1))(-2) \\ &-& (n-4) \left(\prod_{k=1}^{n-5}2k\right) (-2(n-2))(-2(n-1))(-2) \\ &\vdots \\ &+& (-1)^{n-4} 3 (2 \cdot 4) (-10)(-12)\cdots (-2(n-2))(-2(n-1))(-2) \\ &+& (-1)^{n-3} 2 (2) (-8)(-10)(-12)\cdots (-2(n-2))(-2(n-1))(-2) \\ &+& (-1)^{n-2} 1 (-6)(-8)(-10)(-12)\cdots (-2(n-2))(-2(n-1))(-2) \\ &=& (n-1) \prod_{k=1}^{n-2}2k \\ &-& \left(\prod_{k=1}^{n-1}2k\right) (-1) / (2(n-1)) \\ &+& \left(\prod_{k=1}^{n-1}2k\right) (-1)(-1) / (2(n-2)) \\ &-& \left(\prod_{k=1}^{n-1}2k\right) (-1)(-1)(-1) / (2(n-3)) \\ &\vdots& \\ &+& \left(\prod_{k=1}^{n-1}2k\right) (-1)^{n-4} (-1)^{n-5} (-1) / 8 \\ &+& \left(\prod_{k=1}^{n-1}2k\right) (-1)^{n-3} (-1)^{n-4} (-1) / 6 \\ &+& \left(\prod_{k=1}^{n-1}2k\right) (-1)^{n-2} (-1)^{n-3} (-1) / 4 \\ &=& \left(\prod_{k=1}^{n-1}2k\right) / 2     + \left(\prod_{k=1}^{n-1}2k\right)  / (2(n-1))     + \left(\prod_{k=1}^{n-1}2k\right)  / (2(n-2))     + \left(\prod_{k=1}^{n-1}2k\right)  / (2(n-3))     + \cdots     \\ &+& \left(\prod_{k=1}^{n-1}2k\right)   / 8     + \left(\prod_{k=1}^{n-1}2k\right)   / 6     + \left(\prod_{k=1}^{n-1}2k\right)   / 4 \\ &=& \left(\prod_{k=1}^{n-1}2k\right) \left(\sum_{k=1}^{n-1} \frac 1 {2k}\right) \end{eqnarray}
これを最初の式に代入して
\begin{align} \left|\left(a_{ij}\right)\right|  &= (-1)^{n-1} n(n-1)  \prod_{k=1}^{n-1} 2k + (-1)^n \left(\prod_{k=1}^{n-1}2k\right) \left(\sum_{k=1}^{n-1} \frac 1 {2k} \right) \\ &=  (-1)^n \left(\prod_{k=1}^{n-1}2k\right) \left[\sum_{k=1}^{n-1} \frac 1 {2k} - n (n-1) \right] \\ &=  (-2)^{n-2} (n-1)! \left[\sum_{k=1}^{n-1} \frac 1 k - 2 n (n-1) \right]  \end{align}
より解を得ることができた.
(証明終)
ふぅ... 滅茶苦茶大変です. 私は, この行列式のために一日がかりで付き合いました. 最初にこの行列式を計算したのは, 遥か二年以上も前になります. その時は, 一週間ずっとこの行列式に付きっ切りで考えてました. 当時は漸化式を立てて解き, 上の結果にたどり着いたんです. ところが, あまり綺麗な漸化式ではなく, 一般式を予想→数学的帰納法っていうものだったので, 何とかならないものかーと今回は別の方法を取りました. ともあれ, 今回の連載の中で最強クラスで難しいです.

総括

七回に渡り美しき行列式たちをお送りしました. いかがでしたか. ここに紹介したものはほとんど教科書に載っているもので, もっと綺麗な証明があるかもしれません. ぜひ教えて下さい. 感想でもなんでもいいですよ. コメントお待ちしてます.
余談ですが, 今回のこの連載は2週間に渡り書き溜め, それを1週間で吐き出しました. (まぁぶっちゃけ昔書いた数式とかも少しはあるんですが... 最近平日は忙しいし... 吐き出しているときは殆どその日は書きくわえてません.) 行の中でないような数式は全て外部のtexファイルに書き, それを読んで改行を除去したりしてインポートするような手段を取りました. その仕組みのお陰で, 七回分全てを1つのファイルにしても300行ほど. とても作業がしやすかったです. これから普段LaTeX使うときもこうしようかなと思います. ありがとうございましたー & おつかれさまでした(自分)
今はあまり書く気が起こりませんが, 連載の中で紹介されていない美しい行列式が出てきたら延長戦に入るかも? 何かあったら教えてください. n\times nの行列で(つまり, 3 \times 3とかのトリビアルなものではなくて), 式が閉じていて綺麗な... まぁなんでもいいですわー. ではでは〜♪


美しき行列式たち 〜 第一回 - プログラムモグモグ

  • 1. Vandermonde determinants
  • 2. L same (仮名)


美しき行列式たち 〜 第二回 - プログラムモグモグ

  • 3. Circulant matrix
  • 4. Circulant matrix - arithmetic progression
  • 5. Circulant matrix - one two three


美しき行列式たち 〜 第三回 - プログラムモグモグ

  • 6. Right ones (仮名)
  • 7. Diagonal shift (仮名)
  • 8. Diagonal no-shift (仮名)
  • 9. Diagonal addition (仮名)


美しき行列式たち 〜 第四回 - プログラムモグモグ

  • 10. Just polynomial (仮名)
  • 11. Diagonal ones with wing (仮名)


美しき行列式たち 〜 第五回 - プログラムモグモグ

  • 12. Almost one except for diagonal (仮名)
  • 13. Tridiagonal matrix
  • 14. Diagonal a others b (仮名)
  • 15. X (仮名)


美しき行列式たち 〜 第六回 - プログラムモグモグ

  • 16. Cauchy determinants
  • 17. Hilbert matrix


美しき行列式たち 〜 第七回 - プログラムモグモグ

  • 18. L increasing power (仮名)
  • 19. L increasing number (仮名)