［PDF］Nクイーン問題

Ｎクイーン問題
Ｎクイーン問題
０．目次
１．Ｎクイーン問題
１．１
１．２
１．３
バックトラックに基づく解法
バックトラックに基づく解法の改良
順列生成に基づく解法
- 1 -
Ｎクイーン問題
１．Ｎクイーン問題
Ｎ個のクイーンをＮ × Ｎのチェス盤上に、互いにきき筋 (縦、横、斜め方向 )に
入らないように置く方法を考察する。
●クイーンのきき筋
Ｑ
● Nクイーン問題の解の数。
N
3
4
5
6
7
8
9
10
解の数
0
2
10
4
40
92
352
724
N
11
12
13
14
15
16
17
18
解の数
2680
14200
73712
365596
2279184
14772512
95815104
666090624
- 2 -
N
19
20
21
22
23
24
解の数
4968057848
39029188884
314666222712
2691008701644
24233937684440
227514171973736
Ｎクイーン問題
１．１
バックトラックに基づく解法
チェス盤のマス目を指定するために、 n× n の行列とみなし、 i行目、 j列目の
マス目を (i,j) で指定する。
→
j列目
↓
i行目
●バックトラックの考え方
（１）左の列から右の列へ順にひとつずつクイーンを置いていく。
（２）これからクイーンを置こうとする列のすべてのマス目のうち、
・すでに置かれたクイーンのきき筋に入らないマス目があれば、
そのマス目にクイーンを置き、右隣の列にクイーンを置くことを
続ける。
・そのようなマス目が複数個あるときはその中で一番下 (行番号の
大きい方を下とする )のマス目を選ぶ。
・すべてのマス目がすでに置かれたクイーンのきき筋に入っている
ときは、前の列に戻り、クイーンの置かれたマス目のひとつ上から
一番上のマス目までに対して（２）と同様のことを行う。
（３）こうして、最後の列にクイーンが置ければ解になる。
もし、前の列に戻れなくなったら終了となる。
- 3 -
Ｎクイーン問題
1
2
3
4
状態 (1)
1 2 3 4
□□□□
□□□□
□□□□
●□□□
1
2
3
4
状態 (6)
1 2 3 4
□●□×
□□□×
□□●×
●□□×
1
2
3
4
状態 (11)解
1 2 3 4
□●□□
□□□●
●□□×
□□●×
1
2
3
4
状態 (16)
1 2 3 4
□□●□
●□×□
□□×□
□●×□
1
2
3
4
状態 (21)
1 2 3 4
●□□□
□□□□
□□□□
□●□□
1
2
3
4
状態 (26)
1 2 3 4
●□×□
□□×□
□●×□
□□×□
1
2
3
4
状態 (2)
1 2 3 4
□□□□
□●□□
□×□□
●×□□
1
2
3
4
状態 (7)
1 2 3 4
□●×□
□□×□
□□□□
●□□□
1
2
3
4
状態 (12)
1 2 3 4
□●□×
□□□□
●□□□
□□●□
1
2
3
4
状態 (17)解
1 2 3 4
□□●□
●□□□
□□□●
□●□×
1
2
3
4
状態 (22)
1 2 3 4
●□□□
□□●□
□□×□
□●×□
1
2
3
4
状態 (27)
1 2 3 4
●×□□
□×□□
□□□□
□□□□
1
2
3
4
状態 (3)
1 2 3 4
□□×□
□●×□
□□×□
●□×□
1
2
3
4
状態 (8)
1 2 3 4
□□□□
□□□□
●□□□
□□□□
1
2
3
4
状態 (13)
1 2 3 4
□●×□
□□×□
●□×□
□□□□
1
2
3
4
状態 (18)
1 2 3 4
□□●×
●□□×
□□□□
□●□□
1
2
3
4
状態 (23)
1 2 3 4
●□□×
□□●×
□□□×
□●□×
1
2
3
4
状態 (28)終了
1 2 3 4
□□□□
□□□□
□□□□
□□□□
解の候補は、 n**n個になる。
- 4 -
1
2
3
4
状態 (4)
1 2 3 4
□●□□
□□□□
□□□□
●□□□
1
2
3
4
状態 (9)
1 2 3 4
□●□□
□×□□
●×□□
□×□□
1
2
3
4
状態 (14)
1 2 3 4
□□□□
●□□□
□□□□
□□□□
1
2
3
4
状態 (19)
1 2 3 4
□×□□
●×□□
□×□□
□□□□
1
2
3
4
状態 (24)
1 2 3 4
●□×□
□□□□
□□□□
□●□□
1
2
3
4
状態 (5)
1 2 3 4
□●□□
□□□□
□□●□
●□×□
1
2
3
4
状態 (10)
1 2 3 4
□●□□
□□□□
●□□□
□□●□
1
2
3
4
状態 (15)
1 2 3 4
□□□□
●□□□
□□□□
□●□□
1
2
3
4
状態 (20)
1 2 3 4
●□□□
□□□□
□□□□
□□□□
1
2
3
4
状態 (25)
1 2 3 4
●□□□
□□□□
□●□□
□□□□
Ｎクイーン問題
●４×４のチェス盤における探索過程
(1,4)
(2,4)
(2,3)
(2,2)
(2,1)
(1,3)
(2,4)
(2,3)
(2,2)
(2,1)
(3,4)
(3,3)
(3,2)
(3,1)
(3,4)
(3,3)
(3,4)
(4,4)
(4,3)
(4,2)
(4,1)
(4,4)
(4,3)
(4,2)
(4,1)
解
(3,3)
(3,2)
(3,1)
(1,2)
(2,4)
(3,4)
(3,3)
(3,2)
(3,1)
(4,4)
(4,3)
(4,2)
(4,1)
(2,3)
(2,2)
(2,1)
(1,1)
(2,4)
(3,4)
(3,3)
(3,2)
(3,1)
(2,3)
(3,4)
(3,3)
(3,2)
(3,1)
(2,2)
(2,1)
- 5 -
(4,4)
(4,3)
(4,2)
(4,1)
解
Ｎクイーン問題
●バックトラックに基づくＮクイーン問題プログラム
(i,j)にクイーン Qを置く場合、チェックするマス目を青で示す。
→
j
↓
i
Q
● Ｎクイーン問題プログラム (nq111.c)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
/* << nq111.c >> */
/* Ｎクイーン問題 */
#include <stdio.h>
#define N 30
/* チェス盤サイズの最大値。 */
int b[N+1][N+1], /* チェス盤を表す配列。
b[i][j]=0：マス目 (i,j)にクイーンがある場合。
b[i][j]=1：マス目 (i,j)にクイーンがない場合。 */
count,
/* 配置の個数。 */
n;
/* チェス盤のサイズ。 */
void queen(int i, int j);
int main() {
int i,j;
while( 1 ) {
/* nの読み込み。 */
scanf("%d",&n);
if( (n <= 0)||(n > N) ) { break; }
/* チェス盤等の初期設定。 */
for( i=1; i<=n; i++ ) {
for( j=1; j<=n; j++ ) { b[i][j] = 0; }
}
count = 0;
/* バックトラック。 */
queen(n,0);
printf("%d× %dのチェス盤
解の数： %d\n",n,n,count);
}
}
/* マス目 (i,j)にクイーンを置く操作。 */
void queen(int i, int j) {
int p,q;
- 6 -
Ｎクイーン問題
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
/* 横方向のマス目を調べる。 */
for( q=j-1; q>=1; q-- ) {
/* すでにある場合、戻る。 */
if( b[i][q] == 1 ) { return; }
}
/* 斜め上方向のマス目を調べる。 */
p = i-1; q = j-1;
while( (p >= 1)&&(q >= 1) ) {
/* すでにある場合、戻る。 */
if( b[p][q] == 1 ) { return; }
p--; q--;
}
/* 斜め下方向のマス目を調べる。 */
p = i+1; q = j-1;
while( (p <= n)&&(q >= 1) ) {
/* すでにある場合、戻る。 */
if( b[p][q] == 1 ) { return; }
p++; q--;
}
/* きき筋に他のクイーンがないので、 j列目のマス目 (i,j)にクイーンを置く。 */
b[i][j] = 1;
if( j == n ) {
/* 解の表示。 */
count++; printf("[%d] \n",count);
for( p=1; p<=n; p++ ) {
for( q=1; q<=n; q++ ) { printf("%3d",b[p][q]); }
printf("\n");
}
} else {
/* 右隣の列 j+1の n個のマス目を調べる。 */
for( p=n; p>=1; p-- ) {
queen(p,j+1);
}
}
/* j列目のマス目 (i,j)からクイーンを取り除く。 */
b[i][j] = 0;
}
- 7 -
Ｎクイーン問題
実行結果
% cc nq111.c
% a.out
4
[1]
0 1 0 0
0 0 0 1
1 0 0 0
0 0 1 0
[2]
0 0 1 0
1 0 0 0
0 0 0 1
0 1 0 0
4× 4のチェス盤解の数： 2
8
＜＜途中省略＞＞
8× 8のチェス盤解の数： 92
0
実行時間
実行時間測定
% time ./a.out
13 0
13× 13のチェス盤解の数： 73712
2.579u 0.001s 0:04.73 54.3%
0+0k 0+8io 0pf+0w
% time ./a.out
14 0
14× 14のチェス盤解の数： 365596
16.881u 0.000s 0:18.67 90.4%
0+0k 0+8io 0pf+0w
% time ./a.out
15 0
15× 15のチェス盤解の数： 2279184
116.365u 0.000s 1:58.95 97.8%
0+0k 0+8io 0pf+0w
- 8 -
Ｎクイーン問題
１．２
バックトラックに基づく解法の改良
横方向、斜め方向を配列で表現し、すでにクイーンが置かれているかいないか
の判断を早める。
横方向
左斜め上方向
(i-j+nが一定 )
→j
↓1
i 2
3
4
1
2
3
4
1
4
5
6
7
2
3
4
5
6
3
2
3
4
5
4
1
2
3
4
左斜め下方向
(i+jが一定）
1
2
3
4
1
2
3
4
5
2
3
4
5
6
3
4
5
6
7
4
5
6
7
8
●バックトラックに基づくＮクイーン問題プログラムの改良
(i,j)にクイーンを置く場合、チェックするマス目。
→
j列目
↓
i行目
Ｑ
青線
ue[i-j+n]
赤線
gyo[i]
緑線
sita[i+j]
● Ｎクイーン問題プログラム (nq121.c)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
/* << nq121.c >> */
/* Ｎクイーン問題 */
#include <stdio.h>
#define N 30
/* チェス盤サイズの最大値。 */
int gyo[N+1],
/* gyo[i]=0：行 iにクイーンがない。
gyo[i]=1：行 iにクイーンがある。 */
ue[2*N],
/* ue[i-j+n]=0:左斜め上方向 i-j+nにクイーンがない。
ue[i-j+n]=1:左斜め上方向 i-j+nにクイーンがある。 */
sita[2*N+1], /* sita[i+j]=0：左斜め下方向 i+jにクイーンがない。
sita[i+j]=1：左斜め下方向 i+jにクイーンがある。 */
kai[N+1],
/* 解：列 jの行 kai[j]にクイーンがある。 */
count,
/* 得られた配置の個数。 */
n;
/* チェス盤のサイズ。 */
void queen(int i, int j);
int main() {
int k;
while( 1 ) {
/* nの読み込み。 */
- 9 -
Ｎクイーン問題
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
scanf("%d",&n);
if( (n <= 0)||(n > N) ) { break; }
/* チェス盤等の初期設定。 */
for( k=1; k<=n; k++ ) { gyo[k] = 0; }
for( k=1; k<=2*n-1; k++ ) { ue[k] = 0; }
for( k=2; k<=2*n; k++ ) { sita[k] = 0; }
count = 0;
/* バックトラック。 */
for( k=n; k>=1; k-- ) { queen(k,1); }
printf("%d× %dのチェス盤解の数： %d\n",n,n,count);
}
}
/* マス目 (i,j)にクイーンを置く操作。 */
void queen(int i, int j) {
int k;
/* 横方向のマス目を調べ、すでにある場合、戻る。 */
if( gyo[i] == 1 ) { return; }
/* 左斜め上方向のマス目を調べ、すでにある場合、戻る。 */
if( ue[i-j+n] == 1 ) { return; }
/* 左斜め下方向のマス目を調べ、すでにある場合、戻る。 */
if( sita[i+j] == 1 ) { return; }
/* きき筋に他のクイーンがないので、マス目 (i,j)にクイーンを置く。 */
kai[j] = i;
gyo[i] = 1;
ue[i-j+n] = 1;
sita[i+j] = 1;
if( j == n ) {
/* 解の表示。 */
count++; printf("[%d]:",count);
for( k=1; k<=n; k++ ) {
printf("%3d",kai[k]);
}
printf("\n");
} else {
/* 右隣の列 j+1の n個のマス目を調べる。 */
for( k=n; k>=1; k-- ) {
queen(k,j+1);
}
}
/* マス目 (i,j)のクイーンを取り除く。 */
kai[j] = 0;
- 10 -
Ｎクイーン問題
71
72
73
74
gyo[i] = 0;
ue[i-j+n] = 0;
sita[i+j] = 0;
}
実行結果
% cc nq121.c
% ./a.out
5
[1]: 5 3 1 4 2
[2]: 5 2 4 1 3
＜＜途中省略＞＞
[10]: 1 3 5 2 4
5× 5のチェス盤解の数： 10
11
＜＜途中省略＞＞
11× 11のチェス盤解の数： 2680
0
実行時間
実行時間測定
% time ./a.out
14 0
14× 14のチェス盤解の数： 365596
4.077u 0.000s 0:06.58 61.8%
0+0k 0+0io 0pf+0w
% time ./a.out
15 0
15× 15のチェス盤解の数： 2279184
26.523u 0.000s 0:28.89 91.7%
0+0k 0+8io 0pf+0w
% time ./a.out
16 0
16× 16のチェス盤解の数： 14772512
183.400u 0.000s 3:08.56 97.2%
0+0k 0+8io 0pf+0w
- 11 -
Ｎクイーン問題
１．３
部分順列生成に基づく解法
解の候補を n**nから n!に縮小する。
斜め方向を配列で表現し、効率をあげる。
●考え方：部分順列の生成
集合 {x(1),...,x(k)}上のすべての順列を [x(1),...,x(k)]で表す。
k個の要素から 1個ずつ抽出していく。
[x(1),x(2),...,x(k)] = x(1)[x(2),x(3),...,x(k)]
∪ x(2)[x(1),x(3),...,x(k)]
∪ x(3)[x(2),x(1),...,x(k)]
・・・・・・・・・・・・・
∪ x(k-1)[x(2),x(3),...,x(1),x(k)]
∪ x(k)[x(2),x(3),...,x(k-1),x(1)]
抽出する順序に任意性があるので種々の生成法が考えられる。
- 12 -
Ｎクイーン問題
●考え方１
抽出を繰り返し行う。未決定の先頭要素を残りの要素の先頭から
順に交換していく。
・樹形図１
[1234]
1[234]
2[134]
3[214]
4[231]
12[34]
123[4]
124[3]
13[24]
132[4]
134[2]
14[32]
143[2]
142[3]
21[34]
213[4]
214[3]
23[14]
231[4]
234[1]
24[31]
243[1]
241[3]
32[14]
321[4]
324[1]
31[24]
312[4]
314[2]
34[12]
341[2]
342[1]
42[31]
423[1]
421[3]
43[21]
432[1]
431[2]
41[32]
413[2]
412[3]
- 13 -
Ｎクイーン問題
・樹形図２：要素の交換・移動状況
要素を交換した後、すぐ元に戻すことを繰り返す。
このようにすると、順列に関するプログラムを記述しやすくなる。
赤い×は要素の交換、青い×は元に戻すことを意味する。
[1234]
1[234]
12[34]
123[4]
1234
124[3]
1243
123[4]
13[24]
132[4]
1324
134[2]
143[2]
1342
1432
142[3]
1423
21[34]
213[4]
2134
23[14]
214[3]
231[4]
234[1]
2143
2314
2341
243[1]
241[3]
2431
2413
321[4]
3214
324[1]
312[4]
314[2]
341[2]
342[1]
3241
3124
3142
3412
3421
423[1]
4231
421[3]
432[1]
431[2]
413[2]
412[3]
4213
4321
4312
4132
4123
12[34]
14[32]
12[34]
2[134]
1[234]
21[34]
24[31]
21[34]
3[214]
32[14]
1[234]
31[24]
34[12]
4[231]
42[31]
1[234]
43[21]
41[32]
- 14 -
Ｎクイーン問題
● 部分順列生成プログラム (nq131.c)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
/* << nq131.c >> */
/* 集合 {1,2,… ,n}から r個取り出し ,１列に並べる部分順列を生成する。 */
#include <stdio.h>
#define N 30 /* 集合の最大要素数。 */
int p[N+1], /* 部分順列を保存する配列。 */
n,
/* 集合の要素数。 */
r;
/* 取り出す個数。 */
double count; /* 部分順列の個数。 */
void perm(int k);
int main() {
int i;
while( 1 ) {
/* n,rの読み込み。 */
scanf("%d%d",&n,&r);
if( (n <= 0)||(n > N)||(r > n) ) { break; }
/* 初期設定。 */
for( i=1; i<=n; i++ ) { p[i] = i; }
count = 0;
/* 部分順列の生成。 */
perm(1);
printf("n=%d r=%d count=%16.0f\n",n,r,count);
}
}
/* 部分順列の生成。 */
void perm(int k) {
int i,w;
if( k > r ) {
/* 部分順列の表示。 */
count++; printf("%6.0f: ",count);
for( i=1; i<=r; i++ ) { printf("%2d ",p[i]); }
printf("\n");
return;
}
for( i=k; i<=n; i++ ) {
/* 交換。 */
w = p[k]; p[k] = p[i]; p[i] = w;
perm(k+1);
/* 交換を一度戻す。 */
w = p[k]; p[k] = p[i]; p[i] = w;
}
}
- 15 -
Ｎクイーン問題
実行結果
% cc nq131.c
% ./a.out
3 1
1: 1
2: 2
3: 3
n=3 r=1 count=
3 2
1: 1 2
2: 1 3
3: 2 1
4: 2 3
5: 3 2
6: 3 1
n=3 r=2 count=
3 3
1: 1 2
2: 1 3
3: 2 1
4: 2 3
5: 3 2
6: 3 1
n=3 r=3 count=
0 0
3
6
3
2
3
1
1
2
6
実行時間
実行時間測定
% time ./a.out
11 11 0 0
n=11 r=11 count=
39916800
1.202u 0.001s 0:05.08 23.6%
0+0k 0+8io 0pf+0w
% time ./a.out
12 12 0 0
n=12 r=12 count=
479001600
14.415u 0.000s 0:17.93 80.3%
0+0k 0+8io 0pf+0w
% time ./a.out
13 13 0 0
n=13 r=13 count=
6227020800
187.563u 0.011s 3:10.93 98.2%
0+0k 0+8io 0pf+0w
- 16 -
Ｎクイーン問題
データの移動に着目すると、効率のよいプログラムが得られる。
・樹形図３
[1234]
1[234]
12[34]
123[4]
124[3]
13[24]
132[4]
134[2]
14[32]
143[2]
142[3]
2[134]
3[214]
4[231]
21[34]
213[4]
214[3]
23[14]
231[4]
234[1]
24[31]
243[1]
241[3]
32[14]
321[4]
324[1]
31[24]
312[4]
314[2]
34[12]
341[2]
342[1]
42[31]
423[1]
421[3]
43[21]
432[1]
431[2]
41[32]
413[2]
412[3]
- 17 -
Ｎクイーン問題
● 部分順列生成プログラム (nq132.c)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
/* << nq132.c >> */
/* 集合 {1,2,… ,n}から r個取り出し ,１列に並べる部分順列を生成する。 */
#include <stdio.h>
#define N 30 /* 集合の最大要素数。 */
int p[N+1], /* 部分順列を保存する配列。 */
n,
/* 集合の要素数。 */
r;
/* 取り出す個数。 */
double count; /* 部分順列の個数。 */
void perm(int k);
int main() {
int i;
while( 1 ) {
/* n,rの読み込み。 */
scanf("%d%d",&n,&r);
if( (n <= 0)||(n > N)||(r > n) ) { break; }
/* 初期設定。 */
for( i=1; i<=n; i++ ) { p[i] = i; }
count = 0;
/* 部分順列の生成。 */
perm(1);
printf("n=%d r=%d count=%16.0f\n",n,r,count);
}
}
/* 部分順列の生成。 */
void perm(int k) {
int i,w;
if( k > r ) {
/* 部分順列の表示。 */
count++; printf("%6.0f: ",count);
for( i=1; i<=r; i++ ) { printf("%2d ",p[i]); }
printf("\n");
return;
}
w = p[k];
for( i=k; i<=n; i++ ) {
p[k] = p[i]; p[i] = w;
perm(k+1);
p[i] = p[k];
}
p[k] = w;
}
- 18 -
Ｎクイーン問題
実行結果
% cc nq132.c
% ./a.out
3 1
1: 1
2: 2
3: 3
n=3 r=1 count=
3 2
1: 1 2
2: 1 3
3: 2 1
4: 2 3
5: 3 2
6: 3 1
n=3 r=2 count=
3 3
1: 1 2
2: 1 3
3: 2 1
4: 2 3
5: 3 2
6: 3 1
n=3 r=3 count=
0 0
3
6
3
2
3
1
1
2
6
実行時間
実行時間測定
% time ./a.out
11 11 0 0
n=11 r=11 count=
39916800
1.018u 0.000s 0:04.65 21.7%
0+0k 0+8io 0pf+0w
% time ./a.out
12 12 0 0
n=12 r=12 count=
479001600
12.219u 0.000s 0:15.54 78.5%
0+0k 0+8io 0pf+0w
% time ./a.out
13 13 0 0
n=13 r=13 count=
6227020800
158.758u 0.000s 2:42.25 97.8%
0+0k 0+8io 0pf+0w
- 19 -
Ｎクイーン問題
● Ｎクイーン問題プログラム (nq133.c)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
/* << nq133.c >> */
/* Ｎクイーン問題 */
#include <stdio.h>
#define N 30
/* チェス盤サイズの最大値。 */
int p[N+1],
/* 部分順列を保存する配列。 */
count,
/* 生成される解の個数。 */
n,
/* 要素数。 */
ue[2*N],
/* ue[i-j+n]=0:左斜め上方向 i-j+nにクイーンがない。
ue[i-j+n]=1:左斜め上方向 i-j+nにクイーンがある。 */
sita[2*N+1]; /* sita[i+j]=0:左斜め下方向 i+jにクイーンがない。
sita[i+j]=1:左斜め下方向 i+jにクイーンがある。 */
void perm(int k);
int main() {
int i;
while( 1 ) {
/* nの読み込み。 */
scanf("%d",&n);
if( (n <= 0)||(n > N) ) { break; }
/* 部分順列生成の初期設定。 */
for( i=1; i<=n; i++ ) { p[i] = i; }
/* チェス盤等の初期設定。 */
for( i=1; i<=2*n-1; i++ ) { ue[i] = 0; }
for( i=2; i<=2*n; i++ ) { sita[i] = 0; }
count = 0;
/* バックトラック。 */
perm(1);
printf("チェス盤 %d× %dの解の数： %d \n",n,n,count);
}
}
/* 配置の生成。 */
void perm(int k) {
int i,w;
if( k > n ) {
/* 解の表示。 */
count++; printf("[%d]:",count);
for( i=1; i<=n; i++ ) {
printf("%3d",p[i]);
}
printf("\n");
return;
}
- 20 -
Ｎクイーン問題
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
w = p[k];
for( i=k; i<=n; i++ ) {
p[k] = p[i]; p[i] = w;
/* マス目 (p[k],k)にクイーンが置けることがわかったので置く。 */
if( (ue[p[k]-k+n] == ①
ue[p[k]-k+n] = ③
;
sita[p[k]+k] = ④
;
)&&(sita[p[k]+k] == ②
perm(k+1);
/* マス目 (p[k],k)からクイーンを取り除く。 */
ue[p[k]-k+n] = ⑤
;
sita[p[k]+k] = ⑥
;
}
p[i] = p[k];
}
p[k] = w;
}
実行結果
% cc nq133.c
% ./a.out
5
[1]: 1 3 5 2 4
[2]: 1 4 2 5 3
[3]: 2 4 1 3 5
[4]: 2 5 3 1 4
[5]: 3 1 4 2 5
[6]: 3 5 2 4 1
[7]: 4 2 5 3 1
[8]: 4 1 3 5 2
[9]: 5 2 4 1 3
[10]: 5 3 1 4 2
チェス盤 5× 5の解の数： 10
6
[1]: 2 4 6 1 3 5
[2]: 3 6 2 5 1 4
[3]: 4 1 5 2 6 3
[4]: 5 3 1 6 4 2
チェス盤 6× 6の解の数： 4
0
- 21 -
) ) {
Ｎクイーン問題
実行時間
実行時間測定
% time ./a.out
14 0
チェス盤 14× 14の解の数： 365596
1.874u 0.000s 0:03.50 53.4%
0+0k 0+8io 0pf+0w
% time ./a.out
15 0
チェス盤 15× 15の解の数： 2279184
11.896u 0.000s 0:14.08 84.4%
0+0k 0+8io 0pf+0w
% time ./a.out
16 0
チェス盤 16× 16の解の数： 14772512
80.297u 0.000s 1:21.92 98.0%
0+0k 0+8io 0pf+0w
- 22 -

Download Report