第２０講 一般種法による魔方陣ソフトの高速化
第１話 一般種とは？

特殊種法の工夫により、５次魔方陣の作成速度は改良版に比べて１８２２倍にもなりました。
計測していませんが、７次以降では軽く１０万倍（一番最初のソフトと比べれば100億倍）を超えると思われますが、
それでも７次以降では歯が立ちませんでした。
さらに、特殊種法では３次と６次はできません。
そこで、さらに改良したいと思います。
それは一般種法による改良です。
前回の特殊種を、一辺種に変更します。
するとかなり高速化されると同時に、３次と６次もできるようになります。
というより理論的にはあらゆる魔方陣（２次の魔方陣は存在しませんので３次以上）を作成できます。
とっても現実的には、この方法をもってしても１０次あたりが限界です。
１０次あたりでも１個作成させるのに１０分程度は時間が必要です。
１１次あたりでは１時間以上必要でしょう。

では、一般種とは何でしょうか。一般種の例を挙げると、

などです。
いずれの対角線、列、行の合計も１０になっていて、そこは特殊種と同じですが、
特殊種との違いは一目瞭然です。
特殊種の場合は対角線上にも列にも行にも同じ数字が重複してはいけませんでしたが、
上の例においては、行も列も対角線上も同じ数字が重複しています。
白以外の色がついている数字は重複の例です。
色以外にも多数重複していることを確認して下さい。
このように対角線、列、行の合計がすべて１０になっていれば、
数字の重複は問題にしない種が、一般種です。
特殊種は一般種の特別の場合であることが分かります。
だから、特殊な種ということで特殊種と名付けているのです。
ならばプログラミングにおいても、
特別な種から出発しないで、一般種から出発した方がよいと思うかもしれません。
ですが、実は一般種を作らせる普遍的なプログラミングを考えることは、
特殊種の場合より遙かに複雑で大変なのです。
というのは効率的に、一般種を作らせるには番号付け

を改良して

（紺はｙ、赤はｘに、ピンクはｉに対応）
としなければなりません。
４次と５次に限定するなら簡単です。
もちろん偶数の場合と奇数の場合では番号付けのルールが違いますから、
偶数用と奇数用の二つは用意しなければならないわけですが、
どんな偶数でもどんな奇数でも通用するプログラミングを用意しなければなりません。
講義が進めば100次魔方陣あたりも作成可能になるからです。
ところが、普遍的なプログラミングをすることは非常にやっかいなのです。

ソースでいえば、i1とi2を変更しなければならないです。
　　 　　　　　void i1(char n){
　　　 　　　　　　char i,j,m;
　　　　 　　　　　m=n/2;
　　　　　 　　　　for(i=0;i<n*n;i++){
　　　　　　 　　　　　if(i<n){
　　　　　　　 　　　　　　x[i]=i;
　　　　　　　　 　　　　　y[i]=i;
　　　　　　　　　 　　}
　　　　　　　　　　 　else if(i<2*n-1){
　　　　　　　　　　　 　　x[i]=2*n-i-1;
　　　　　　　　　　　　 　y[i]=i-n;
　　　　　　　　　　　　　 if(x[i]<=y[i]){
　 　　　　　　　　　　　　　　x[i]--;
　　 　　　　　　　　　　　　　y[i]++;
　　　 　　　　　　　　　　}
　　　　 　　　　　　　}
　　　　　 　　　　　　else if(i<2*n+(n*n-3*n)/2){
　　　　　　 　　　　　　　x[i]=(i-2*n+1)%(n-2);
　　　　　　　 　　　　　　y[i]=(i-2*n+1)/(n-2);
　　　　　　　　 　　　　　if(x[i]>=y[i])x[i]++;
　　　　　　　　　 　　　　if(x[i]>=n-y[i]-1)x[i]++;
　　　　　　　　　　 　}
　　　　　　　　　　　 else if(i<3*n-1+(n*n-3*n)/2){
　 　　　　　　　　　　　　x[i]=i-2*n-(n*n-3*n)/2;
　　 　　　　　　　　　　　y[i]=m;
　　　 　　　　　　　　　　if(x[i]>=y[i])x[i]++;
　　　　 　　　　　　　}
　　　　　 　　　　　　else{
　　　　　　 　　　　　　　x[i]=(i-2*n)%(n-2);
　　　　　　　 　　　　　　y[i]=(i-2*n)/(n-2);
　　　　　　　　 　　　　　if(x[i]>=n-y[i]-1)x[i]++;
　　　　　　　　　 　　　　if(x[i]>=y[i])x[i]++;
　　　　　　　　　　 　}
　 　　　　　　　　}
　　 　　　　　}
　　　 　　　　void i2(char n){
　　　　 　　　　　char i,m;
　　　　　 　　　　m=n/2;
　　　　　　 　　　for(i=0;i<n*n;i++){
　　　　　　　 　　　　if(i<n){
　　　　　　　　 　　　　　x[i]=i;
　　　　　　　　　 　　　　y[i]=i;
　　　　　　　　　　 　}
　　　　　　　　　　　 else if(i<2*n){
　 　　　　　　　　　　　　x[i]=2*n-i-1;
　　 　　　　　　　　　　　y[i]=i-n;
　　　 　　　　　　　　}
　　　　 　　　　　　　else if(i<(m+1)*n){
　　　　　 　　　　　　　　x[i]=(i-2*n)%(n-2);
　　　　　　 　　　　　　　y[i]=(i-2*n)/(n-2);
　　　　　　　 　　　　　　if(x[i]>=y[i])x[i]++;
　　　　　　　　 　　　　　if(x[i]>=n-y[i]-1)x[i]++;
　　　　　　　　　 　　}
　　　　　　　　　　 　else{
　　　　　　　　　　　 　　x[i]=(i-2*n)%(n-2);
　　　　　　　　　　　　 　y[i]=(i-2*n)/(n-2);
　　　　　　　　　　　　　 if(x[i]>=n-y[i]-1)x[i]++;
　 　　　　　　　　　　　　if(x[i]>=y[i])x[i]++;
　　 　　　　　　　　　}
　　　 　　　　　　}
　　　　 　　　}

比較的簡単な偶数版に限定してもかなり、難しいことを確認しておきましょう。

（紺はｙ、赤はｘ、ピンクはｉに対応）
ｉ＜２＊ｎの場合は、前の番号付けと同じですので、
ｉ≧２＊ｎに限定してみていきましょう。
ｉとｙの対応を４次の場合で見ると、
ｙ[８]＝０，y[９]＝０，ｙ[１０]＝１，y[１１]＝２，ｙ[１２]＝１，y[１３]＝３，ｙ[１４]＝２，y[１５]＝２
つまり、数列０，０，１，２，１，３，２，２です。
このような数列の一般項は見つけられるのでしょうか。
しかも、ｎによらない普遍的（一般的）法則を見いださなければならないのです。
６次の場合は、
ｙ[１２]＝０，y[１３]＝０，ｙ[１４]＝０，y[１５]＝０，ｙ[１６]＝１，y[１７]＝２，ｙ[１８]＝３，y[１９]＝４，
ｙ[２０]＝１，y[２１]＝１，ｙ[２２]＝１，y[２３]＝２，ｙ[２４]＝３，y[２５]＝５，ｙ[２６]＝２，y[２７]＝２，
ｙ[２８]＝４，y[２９]＝５，ｙ[３０]＝３，y[３１]＝３，ｙ[３２]＝４，y[３３]＝５，ｙ[３４]＝４，y[３５]＝５
今回は数列０，０，０，０，１，２，３，４，１，１，１，２，３，５，２，２，４，５，３，３，４，５，４，５
こんな複雑な数列の一般項を求めることなどできるのでしょうか。
さらに１０次の場合を見てみましょう。

ますます混迷を深めているように見ます。
こんな複雑な迷宮に出口などあるのでしょうか。

今回は、課題ということにはしません。
余りにも難しすぎるからです。
ですが、ファイトのある方は是非挑戦して下さい。
次話で、解答例を示しますが皆さんの方がもっとよいプログラミングをされているかもしれません。
（次話まで待てない方は、
VB講義の応用編の第２話　新しいセル位置番号付けはどうするか？
を参照して下さい。
Visual Basicなので文法が少し違いますが、ここまでついて来られた優秀な皆さんなら解読できると思いますよ。）

第１９講第８話へ 第２０講（再講義）第２話へ



VC++講義第１部へ
vb講義へ
ＶＢ講義基礎へ
初心者のための世界で一番わかりやすいVisual C++入門基礎講座
初心者のための世界で一番わかりやすいVisual Basic入門基礎講座