第１７講 関数の再帰的呼び出しによる３次・４次魔方陣ソフトの改良
第１話 プログラムの改良

今からプログラムの改良に入りますが、
実は前回作ったプログラムは３次・４次魔方陣に限定されたプログラムではありません。
理論的には１００次魔方陣でも作ることのできるプログラムです。
もちろん冒頭の
#pragma once
int n;
int a[10][10];
int s;
は
#pragma once
int n;
int a[100][100];
int s;
と変更しなければなりませんが。
理論的には int a[n][n];としておけば、３
以上でn以下の整数なら何次の魔方陣でも作ることのできる普遍的なプログラムです。
だから、デザインの方も
としておいてもてよいのです。
１５以下としたのは、DataGridViewがまだ１５列しかないからです。
だから、５を入力して忍耐強く待ち続ければいつかは魔方陣は出てきます。
もちろん、
void f(char g){
　　if(s==10)return;
　　int i,j,k,h,x,y,xx,yy,wa;
　　　　　　　　　　．
　　　　　　　　　　．
　　　　　　　　　　．
　　　　　　　　if(h==1){
　　　　　　　　　　if(g<n*n-1){
　　　　　　　　　　　　f(g+1);
　　　　　　　　　　}
　　　　　　　　　　else{
　　　　　　　　　　　　for(j=0;j<n;j++){
　　　　　　　　　　　　　　for(k=0;k<n;k++){
　　　　　　　　　　　　　　　　　w[k]=(a[j][k]).ToString();
　　　　　　　　　　　　　　}
　　　　　　　　　　　　　　dataGridView1->Rows->Add(w);
　　　　　　　　　　　　}
　　　　　　　　　　　　for(j=0;j<16;j++)w[j]=L"";
　　　　　　　　　　　　dataGridView1->Rows->Add(w);
　　　　　　　　　　　　s++;
　　　　　　　　　　　　if(s==10)return;
　　　　　　　　　　}
のif(s==10)return;入れておくことが前提です。このプログラムでは２１億個の魔方陣を作成させるには、１０００年かけても終わらないからです。
DataGridViewは残念ながらリアルタイムに表示させることができず、すべての計算が終わってからしか結果が表示されないのです。
ただし、if(s==1)return;としても、私のパソコンでは６時間待っても結果は表示されませんでした。
もし６を入力したなら１００日待ったも出てこないかもしれません。

ところが、これからちょっととした改良を加えるだけで、
５次魔方陣なら１個出てくるまでの時間は、０．３７４４秒（私の使用パソコンで）まで縮めることができます。
６時間待っても一つも表示されなかったので、
仮に１個出てくるまでの時間が６時間であるとして計算しもて、６万倍以上の作成速度の差があることになります。
４次の場合は、７０４０個のすべてを作り出すスピードは１７倍程度ですから、
次数が一つ増えるだけで大きく作成速度が違ってくることが分かります。
したがって、６次以降については、作成速度が１０００万倍以上は軽く越えるでしょう。

さて、その小改良とは枠番後gの付け方を変えるだけです。

（赤の番号はx、濃紺の番号はy、ピンクの番号はgに対応）

つまり、対角線→逆対角線→その他の順に変更するだけで、５次魔方陣において３万倍以上も作成速度が大きくなるのです。
さて、問題は番号の付け替えをどのようにするかです。
前回においては、x=g%n;　y=g/n;と座標をいちいち計算させていましたが、これは無駄です。
毎回座標を計算し直しているわけですから。
ですから、座標作成関数g1とg2を用意しておいて、まず最初に座標を計算させておきたいものです。
関数を二つ用意する理由は、偶数次と奇数次では作成ルールが違うからです。
上の４次と５次を比較して頂ければそれは理解して頂けると思います。
もし、番号の付け方が

（赤の番号はx、濃紺の番号はy、ピンクの番号はgに対応）
なら、用意する関数は一つでよく、
#pragma once
int n;
int a[10][10];
int x[100],y[100];
int s;
private: System::Void button1_Click(System::Object^ sender, System::EventArgs^ e) {
　 　　　　　　　DateTime^ hj=DateTime::Now;
　　 　　　　　　n=int::Parse(textBox1->Text);
　　　 　　　　　s=0;
　　　　　　　　 g(n);
　　　　 　　　　f(0);
　　　　　　　　　　　・
　　　　　　　　　　　・
　　　　　　　　　　　・
　 　　　　}
　　 　　　void g(char n){
　　　　　　　　 int i;
　　　　　　　　 for(i=0;i<n*n;i++){
　　　　　　　　　　 x[i]=i%n;
　　　　　　　　　　 y[i]=i/n;
　　　　　　　　 }
　　　　　　}
というプログラム例が考えられます。

（赤の番号はx、濃紺の番号はy、ピンクの番号はgに対応）
のように番号付けを行うにはどうしたらよいでしょうか。

３次・４次に限定するなら、
void g1(char n){
　　x[0]=0;
　　x[1]=1;
　　x[2]=2;
　　x[3]=2;
　　x[4]=0;
　　x[5]=1;
　　x[6]=0;
　　x[7]=2;
　　x[8]=1;
　　y[0]=0;
　　y[1]=1;
　　y[2]=2;
　　y[3]=0;
　　y[4]=2;
　　y[5]=0;
　　y[6]=1;
　　y[7]=1;
　　y[8]=2;
}

void g2(char n){
　　x[0]=0;
　　x[1]=1;
　　x[2]=2;
　　x[3]=3;
　　x[4]=3;
　　x[5]=2;
　　x[6]=1;
　　x[7]=0;
　　x[8]=1;
　　x[9]=2;
　　x[10]=0;
　　x[11]=3;
　　x[12]=0;
　　x[13]=3;
　　x[14]=1;
　　x[15]=2;
　　y[0]=0;
　　y[1]=1;
　　y[2]=2;
　　y[3]=3;
　　y[4]=0;
　　y[5]=1;
　　y[6]=2;
　　y[7]=3;
　　y[8]=0;
　　y[9]=0;
　　y[10]=1;
　　y[11]=1;
　　y[12]=2;
　　y[13]=2;
　　y[14]=3;
　　y[15]=3;
}
でもいいわけですが、これは芸がありません。
プログラムにおいて大事なのは、如何にして一般化するかです。
今回の改良で実際に扱える（時間的に可能なのは）５次までです。
６次以降だと軽く１０００万倍以上は速くなっていますが、
それでも６次以降だと太刀打ちできないのです。
ですが、将来のおいては、26次あたりまで扱えるようになります。
講義が進んでいけば、26次でもコンピュータは1秒に数千個の勢いで魔方陣を作れるようになります。
ですから、将来に向けて座標番号付けは是非とも普遍的(一般的）なものにしてほしいのです。
前に示した
　　 　　　void g(char n){
　　　　　　　　 int i;
　　　　　　　　 for(i=0;i<n*n;i++){
　　　　　　　　　　 x[i]=i%n;
　　　　　　　　　　 y[i]=i/n;
　　　　　　　　 }
　　　　　　}
は普遍的なプログラムです。nには自由に整数を入力することができるからです。

の場合も普遍的なプログラムを考えてほしいわけです。
超難解ですよ。
頭が爆発しそうになるほど。

今回は、30分考えて分からなければ第2話を開きましょう。
今回は超ド級の難問だからです。

第１６講第７話へ　　第１７講第２話へ

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