|
1 // Part 1 about finding and counting Knight's tours |
|
2 //================================================== |
|
3 |
|
4 //object CW8a { // for preparing the jar |
|
5 |
|
6 type Pos = (Int, Int) // a position on a chessboard |
|
7 type Path = List[Pos] // a path...a list of positions |
|
8 |
|
9 |
|
10 // for measuring time in the JAR |
|
11 def time_needed[T](code: => T) : T = { |
|
12 val start = System.nanoTime() |
|
13 val result = code |
|
14 val end = System.nanoTime() |
|
15 println(f"Time needed: ${(end - start) / 1.0e9}%3.3f secs.") |
|
16 result |
|
17 } |
|
18 |
|
19 // for printing a board |
|
20 def print_board(dim: Int, path: Path): Unit = { |
|
21 println |
|
22 for (i <- 0 until dim) { |
|
23 for (j <- 0 until dim) { |
|
24 print(f"${path.reverse.indexOf((j, dim - i - 1))}%3.0f ") |
|
25 } |
|
26 println |
|
27 } |
|
28 } |
|
29 |
|
30 def is_legal(dim: Int, path: Path, x: Pos): Boolean = |
|
31 0 <= x._1 && 0 <= x._2 && x._1 < dim && x._2 < dim && !path.contains(x) |
|
32 |
|
33 // testcases |
|
34 //assert(is_legal(8, Nil, (3, 4)) == true) |
|
35 //assert(is_legal(8, List((4, 1), (1, 0)), (4, 1)) == false) |
|
36 //assert(is_legal(2, Nil, (0, 0)) == true) |
|
37 |
|
38 |
|
39 def add_pair(x: Pos, y: Pos): Pos = |
|
40 (x._1 + y._1, x._2 + y._2) |
|
41 |
|
42 def moves(x: Pos): List[Pos] = |
|
43 List(( 1, 2),( 2, 1),( 2, -1),( 1, -2), |
|
44 (-1, -2),(-2, -1),(-2, 1),(-1, 2)).map(add_pair(x, _)) |
|
45 |
|
46 // 1 mark |
|
47 |
|
48 def legal_moves(dim: Int, path: Path, x: Pos): List[Pos] = |
|
49 moves(x).filter(is_legal(dim, path, _)) |
|
50 |
|
51 // testcases |
|
52 //assert(legal_moves(8, Nil, (2,2)) == |
|
53 // List((3,4), (4,3), (4,1), (3,0), (1,0), (0,1), (0,3), (1,4))) |
|
54 //assert(legal_moves(8, Nil, (7,7)) == List((6,5), (5,6))) |
|
55 //assert(legal_moves(8, List((4,1), (1,0)), (2,2)) == |
|
56 // List((3,4), (4,3), (3,0), (0,1), (0,3), (1,4))) |
|
57 //assert(legal_moves(8, List((6,6)), (7,7)) == List((6,5), (5,6))) |
|
58 //assert(legal_moves(1, Nil, (0,0)) == List()) |
|
59 //assert(legal_moves(2, Nil, (0,0)) == List()) |
|
60 //assert(legal_moves(3, Nil, (0,0)) == List((1,2), (2,1))) |
|
61 |
|
62 // 2 marks |
|
63 |
|
64 def tcount_tours(dim: Int, path: Path): Int = { |
|
65 if (path.length == dim * dim) 1 |
|
66 else |
|
67 (for (x <- legal_moves(dim, path, path.head)) yield tcount_tours(dim, x::path)).sum |
|
68 } |
|
69 |
|
70 def count_tours(dim: Int, path: Path) = |
|
71 time_needed(tcount_tours(dim: Int, path: Path)) |
|
72 |
|
73 |
|
74 def tenum_tours(dim: Int, path: Path): List[Path] = { |
|
75 if (path.length == dim * dim) List(path) |
|
76 else |
|
77 (for (x <- legal_moves(dim, path, path.head)) yield tenum_tours(dim, x::path)).flatten |
|
78 } |
|
79 |
|
80 def enum_tours(dim: Int, path: Path) = |
|
81 time_needed(tenum_tours(dim: Int, path: Path)) |
|
82 |
|
83 // test cases |
|
84 |
|
85 /* |
|
86 def count_all_tours(dim: Int) = { |
|
87 for (i <- (0 until dim).toList; |
|
88 j <- (0 until dim).toList) yield count_tours(dim, List((i, j))) |
|
89 } |
|
90 |
|
91 def enum_all_tours(dim: Int): List[Path] = { |
|
92 (for (i <- (0 until dim).toList; |
|
93 j <- (0 until dim).toList) yield enum_tours(dim, List((i, j)))).flatten |
|
94 } |
|
95 |
|
96 |
|
97 println("Number of tours starting from (0, 0)") |
|
98 |
|
99 for (dim <- 1 to 5) { |
|
100 println(s"${dim} x ${dim} " + time_needed(0, count_tours(dim, List((0, 0))))) |
|
101 } |
|
102 |
|
103 println("Number of tours starting from all fields") |
|
104 |
|
105 for (dim <- 1 to 5) { |
|
106 println(s"${dim} x ${dim} " + time_needed(0, count_all_tours(dim))) |
|
107 } |
|
108 |
|
109 for (dim <- 1 to 5) { |
|
110 val ts = enum_tours(dim, List((0, 0))) |
|
111 println(s"${dim} x ${dim} ") |
|
112 if (ts != Nil) { |
|
113 print_board(dim, ts.head) |
|
114 println(ts.head) |
|
115 } |
|
116 } |
|
117 */ |
|
118 |
|
119 // 1 mark |
|
120 |
|
121 def first(xs: List[Pos], f: Pos => Option[Path]): Option[Path] = xs match { |
|
122 case Nil => None |
|
123 case x::xs => { |
|
124 val result = f(x) |
|
125 if (result.isDefined) result else first(xs, f) |
|
126 } |
|
127 } |
|
128 |
|
129 // test cases |
|
130 //def foo(x: (Int, Int)) = if (x._1 > 3) Some(List(x)) else None |
|
131 // |
|
132 //first(List((1, 0),(2, 0),(3, 0),(4, 0)), foo) |
|
133 //first(List((1, 0),(2, 0),(3, 0)), foo) |
|
134 |
|
135 |
|
136 // 1 mark |
|
137 |
|
138 def tfirst_tour(dim: Int, path: Path): Option[Path] = { |
|
139 if (path.length == dim * dim) Some(path) |
|
140 else |
|
141 first(legal_moves(dim, path, path.head), (x:Pos) => tfirst_tour(dim, x::path)) |
|
142 } |
|
143 |
|
144 def first_tour(dim: Int, path: Path) = |
|
145 time_needed(tfirst_tour(dim: Int, path: Path)) |
|
146 |
|
147 |
|
148 /* |
|
149 for (dim <- 1 to 8) { |
|
150 val t = first_tour(dim, List((0, 0))) |
|
151 println(s"${dim} x ${dim} " + (if (t == None) "" else { print_board(dim, t.get) ; "" })) |
|
152 } |
|
153 */ |
|
154 |
|
155 // 15 secs for 8 x 8 |
|
156 //val ts1 = time_needed(0,first_tour(8, List((0, 0))).get) |
|
157 |
|
158 // no result for 4 x 4 |
|
159 //val ts2 = time_needed(0, first_tour(4, List((0, 0)))) |
|
160 |
|
161 // 0.3 secs for 6 x 6 |
|
162 //val ts3 = time_needed(0, first_tour(6, List((0, 0)))) |
|
163 |
|
164 // 15 secs for 8 x 8 |
|
165 //time_needed(0, print_board(8, first_tour(8, List((0, 0))).get)) |
|
166 |
|
167 |
|
168 //} |
|
169 |
|
170 |