| 247 |      1 | // Part 1 about finding and counting Knight's tours
 | 
|  |      2 | //==================================================
 | 
| 220 |      3 | 
 | 
| 247 |      4 | //object CW8a {   // for preparing the jar
 | 
| 222 |      5 | 
 | 
| 220 |      6 | type Pos = (Int, Int)    // a position on a chessboard 
 | 
|  |      7 | type Path = List[Pos]    // a path...a list of positions
 | 
|  |      8 | 
 | 
|  |      9 | 
 | 
| 247 |     10 | // for measuring time in the JAR
 | 
| 220 |     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 | 
 | 
| 247 |     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)
 | 
| 220 |     32 | 
 | 
| 247 |     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 | }
 | 
| 220 |    117 | */
 | 
| 247 |    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 | 
 |