ログインしてさらにmixiを楽しもう

コメントを投稿して情報交換!
更新通知を受け取って、最新情報をゲット!

魔方陣についてコミュの?4次元方陣〜M次元方陣

  • mixiチェック
  • このエントリーをはてなブックマークに追加
?で↓示した3次元方陣(立方陣)を表した★式
http://mixi.jp/view_bbs.pl?id=2824935

その式を拡張すれば、例えばNxNxNxNの4次元方陣と
 いうのも N≧17以上の素数と、その素数同士の
 積の場合について次のように作れます。

★
T[i][j][k][L]=

N*N*N* mod(  i + 2*j +4*k +8*L +C1 , N)
+ N*N* mod( 2*i + 4*j +8*k + L +C2 , N)
+  N* mod( 4*i + 8*j + k +2*L +C3 , N)
+   mod( 8*i +  j +2*k +4*L +C4 , N)
+1

ここで
 Nは17以上の素数、mod(X,N)はXをNで割った余り。
 0≦i,j,k,L<N
 0≦C1,C2,C3,C4<N C1〜C4は任意の定数。

これで4次元空間をこの超立方体で充填した場合
任意の場所(i,j,k,L)から、それに隣接する方向
(i+α,j+β,K+γ,L+δ) ここで、-1≦α,β,γ,δ≦1
ただしα,β,γ,δすべてが0でない組み合わせとして、
3の4乗−1= 80通りの4次元的な全方向について、
直線上にならぶ、連続したN個の数字の和が定和になります。

この証明は?での記事と同様に
(1)各直線上に連続して並ぶN個の数字の和が等しくなること
(2)重複なく1からNの4乗までの数字から成り立っていること。
の2点から次のように証明できます。

(1)は例えば、α+2β+4γ+8δの取り得る値が
±1から±15までであり、これが具合のいいことに
 ● 0を含まないこと および Nに対して素であること。
 これで、modの項は、各項について、(α,β,γ,δ)方向に
 連続したN個の並びについては、0からN−1を含んで、
 その合計は N*(N−1)/2で一定になります。


(2)については
   0≦i,j,k,L<Nの範囲内では
   あるT[i][j][k][L]と一致するものは、他にないから
   全体では、1からNの4乗までの数字からなることを

   もしT[i][j][k][L]=T[i'][j'][k'][L']ならば
     (i,j,k,L)=(i',j',k',L')が成立することを

Nの剰余系のもとで、
行列
 1248
 2481
 4812
 8124
が逆行列を持つことから?と同様に証明できるからです。

〜〜〜〜〜〜〜〜〜〜〜〜〜〜〜〜〜〜〜〜

同様にして5次元以上:M次元での方陣というものも定義できます。式の形としては同じことですが、


M次元ということで、表記の仕方を変更します。

 i,j,k,L....に替えて X0,X1,X2,...XM-1
 N*N*N といった累乗を^を使って 例えばN^3という具合に
表記すると次式のようになります。


M次元方陣の(X0,X1,X2,...XM-1)座標の値:

H[X0][X1][X2]...[XM-1]=

N^(M-1)* mod (X0+2*X1+4*X2+...+2^(M-1)*XM-1+C1,N)

+N^(M-2)* mod(2*X0+4*X1+8*X2+...+2^(M-1)*XM-2+ 2^0*XM-1+C2,N)

:
;
+ mod(2^(M-1)*X0 +X1 +2*X2+...+ 2^(M-2)*Xm-1 +CM,N)
+1



以上です。
 どうも途中の項がテキストでは表現しずらくなりました(^^ゞ

・・・・後で、Σを使った式に書きなおしてみたいと思います。
(どなたか、先にやっていただいて、ココにコメントくださるとありがたいです。式を記した写真を貼るのも手かなと思っているところです)





・・・・・以上で、私がこのコミュを設置した時からの
目論見であった魔方陣の製法記事は全て載せました。



========================
製法についての記事を?から連番で11/13現在?まで
連番で示しています。 目次として載せておきます
 現在のページを◆で示しておきます。

初めての方は?から順にどうぞ。

・4方陣
? http://mixi.jp/view_bbs.pl?id=2013306
? http://mixi.jp/view_bbs.pl?id=2017172
・5方陣について
? http://mixi.jp/view_bbs.pl?id=2102712
? http://mixi.jp/view_bbs.pl?id=2166770
? http://mixi.jp/view_bbs.pl?id=2370334
・合成数の大きさの方陣
? http://mixi.jp/view_bbs.pl?id=1948761
・立方陣について
 4立方陣
? http://mixi.jp/view_bbs.pl?id=2393919
? http://mixi.jp/view_bbs.pl?id=2552298
? http://mixi.jp/view_bbs.pl?id=2730817
 8立方陣
? http://mixi.jp/view_bbs.pl?id=2730921
・7以上の素数サイズの立方陣 
? http://mixi.jp/view_bbs.pl?id=2750487
・11以の素数サイズの完全立方陣
? http://mixi.jp/view_bbs.pl?id=2824935
・5x5x5,6x6x6の立方陣も存在します。
? http://mixi.jp/view_bbs.pl?id=2892234
・4次元方陣〜M次元方陣
?◆http://mixi.jp/view_bbs.pl?id=2894212

コメント(3)

Xと×が混在しているようです。
>+N^(M-2)* mod(2*X0+4*X1+8*X2+...+2^(M-1)*×M-2+ 2^0*×M-1+C2,N)
+N^(M-2)* mod(2*X0+4*X1+8*X2+...+2^(M-1)*XM-2+ 2^0*XM-1+C2,N)
 御指摘ありがとうございました。訂正しておきます。
ちなみに 行列
 1248
 2481
 4812
 8124

の逆行列は、実数では 1/15*

 -1002
 002-1
 02-10
 2-100

という具合です。剰余系だと
 a00b
 00ba
 0ba0
 ba00

15a≡-1 (mod N)
15b≡2 (mod N)

という具合です。 例えばN=17の時、a=9 b=16です。

ログインすると、みんなのコメントがもっと見れるよ

mixiユーザー
ログインしてコメントしよう!

魔方陣について 更新情報

魔方陣についてのメンバーはこんなコミュニティにも参加しています

星印の数は、共通して参加しているメンバーが多いほど増えます。

人気コミュニティランキング