| author | Christian Urban <urbanc@in.tum.de> | 
| Fri, 22 Nov 2019 16:41:45 +0000 | |
| changeset 323 | 93b6c16dded8 | 
| parent 321 | f6cb539e0ffb | 
| child 325 | 26058bf089ae | 
| permissions | -rw-r--r-- | 
| 67 | 1 | // Scala Lecture 3 | 
| 2 | //================= | |
| 3 | ||
| 320 | 4 | // - last week | 
| 5 | // | |
| 6 | // option type | |
| 7 | // higher-order function | |
| 8 | ||
| 9 | ||
| 323 | 10 | def add(x: Int, y: Int) : Int = x + y | 
| 11 | ||
| 12 | def plus5(x: Int) : Int = add(5, x) | |
| 13 | ||
| 14 | plus5(6) | |
| 15 | ||
| 16 | def add2(x: Int)(y: Int) : Int = x + y | |
| 17 | ||
| 18 | def plus3(y: Int) : Int => Int = add2(3)(y) | |
| 19 | ||
| 20 | plus3(9) | |
| 21 | ||
| 22 | List(1,2,3,4,5).map(add2(3)) | |
| 23 | List(1,2,3,4,5).map(add(3, _)) | |
| 24 | ||
| 25 | type Pos = (Int, Int) | |
| 26 | ||
| 27 | def test(p: Pos) = {
 | |
| 28 |   if (p._1 < 5 && p._2 < 5) {
 | |
| 29 | Some(p) | |
| 30 | } | |
| 31 | } | |
| 32 | ||
| 33 | val l = List((1,2), (5,3), (2,5), (1,3)) | |
| 34 | ||
| 35 | l.map(test).flatten | |
| 320 | 36 | |
| 37 | // Recursion Again ;o) | |
| 38 | //==================== | |
| 39 | ||
| 217 | 40 | |
| 41 | // A Web Crawler / Email Harvester | |
| 42 | //================================= | |
| 43 | // | |
| 44 | // the idea is to look for links using the | |
| 45 | // regular expression "https?://[^"]*" and for | |
| 218 | 46 | // email addresses using yet another regex. | 
| 217 | 47 | |
| 48 | import io.Source | |
| 49 | import scala.util._ | |
| 155 | 50 | |
| 217 | 51 | // gets the first 10K of a web-page | 
| 52 | def get_page(url: String) : String = {
 | |
| 53 |   Try(Source.fromURL(url)("ISO-8859-1").take(10000).mkString).
 | |
| 320 | 54 |     getOrElse { println(s" Problem with: $url"); ""}
 | 
| 217 | 55 | } | 
| 155 | 56 | |
| 217 | 57 | // regex for URLs and emails | 
| 58 | val http_pattern = """"https?://[^"]*"""".r | |
| 59 | val email_pattern = """([a-z0-9_\.-]+)@([\da-z\.-]+)\.([a-z\.]{2,6})""".r
 | |
| 60 | ||
| 218 | 61 | // val s = "foo bla christian@kcl.ac.uk 1234567" | 
| 62 | // email_pattern.findAllIn(s).toList | |
| 155 | 63 | |
| 217 | 64 | // drops the first and last character from a string | 
| 65 | def unquote(s: String) = s.drop(1).dropRight(1) | |
| 155 | 66 | |
| 217 | 67 | def get_all_URLs(page: String): Set[String] = | 
| 68 | http_pattern.findAllIn(page).map(unquote).toSet | |
| 155 | 69 | |
| 320 | 70 | // a naive version of crawl - searches until a given depth, | 
| 217 | 71 | // visits pages potentially more than once | 
| 320 | 72 | def crawl(url: String, n: Int) : Unit = {
 | 
| 73 | if (n == 0) () | |
| 217 | 74 |   else {
 | 
| 75 | println(s" Visiting: $n $url") | |
| 321 | 76 | val page = get_page(url) | 
| 77 | for (u <- get_all_URLs(page)) crawl(u, n - 1) | |
| 217 | 78 | } | 
| 155 | 79 | } | 
| 80 | ||
| 217 | 81 | // some starting URLs for the crawler | 
| 82 | val startURL = """https://nms.kcl.ac.uk/christian.urban/""" | |
| 320 | 83 | |
| 217 | 84 | crawl(startURL, 2) | 
| 85 | ||
| 323 | 86 | for (x <- List(1,2,3,4,5,6)) println(x) | 
| 318 | 87 | |
| 320 | 88 | // a primitive email harvester | 
| 89 | def emails(url: String, n: Int) : Set[String] = {
 | |
| 90 | if (n == 0) Set() | |
| 91 |   else {
 | |
| 92 | println(s" Visiting: $n $url") | |
| 93 | val page = get_page(url) | |
| 94 | val new_emails = email_pattern.findAllIn(page).toSet | |
| 323 | 95 | new_emails ++ (for (u <- get_all_URLs(page).par) yield emails(u, n - 1)).flatten | 
| 320 | 96 | } | 
| 218 | 97 | } | 
| 98 | ||
| 323 | 99 | emails(startURL, 3) | 
| 218 | 100 | |
| 101 | ||
| 320 | 102 | // if we want to explore the internet "deeper", then we | 
| 103 | // first have to parallelise the request of webpages: | |
| 104 | // | |
| 105 | // scala -cp scala-parallel-collections_2.13-0.2.0.jar | |
| 106 | // import scala.collection.parallel.CollectionConverters._ | |
| 155 | 107 | |
| 108 | ||
| 109 | ||
| 320 | 110 | // another well-known example | 
| 111 | //============================ | |
| 178 | 112 | |
| 320 | 113 | def move(from: Char, to: Char) = | 
| 114 | println(s"Move disc from $from to $to!") | |
| 67 | 115 | |
| 320 | 116 | def hanoi(n: Int, from: Char, via: Char, to: Char) : Unit = {
 | 
| 117 | if (n == 0) () | |
| 118 |   else {
 | |
| 119 | hanoi(n - 1, from, to, via) | |
| 120 | move(from, to) | |
| 121 | hanoi(n - 1, via, from, to) | |
| 122 | } | |
| 123 | } | |
| 67 | 124 | |
| 320 | 125 | hanoi(4, 'A', 'B', 'C') | 
| 67 | 126 | |
| 155 | 127 | |
| 128 | ||
| 217 | 129 | // Jumping Towers | 
| 130 | //================ | |
| 131 | ||
| 132 | ||
| 133 | // the first n prefixes of xs | |
| 134 | // for 1 => include xs | |
| 135 | ||
| 136 | def moves(xs: List[Int], n: Int) : List[List[Int]] = (xs, n) match {
 | |
| 137 | case (Nil, _) => Nil | |
| 323 | 138 | case (_, 0) => Nil | 
| 139 | case (y::ys, n) => xs :: moves(ys, n - 1) | |
| 217 | 140 | } | 
| 141 | ||
| 142 | ||
| 143 | moves(List(5,1,0), 1) | |
| 144 | moves(List(5,1,0), 2) | |
| 145 | moves(List(5,1,0), 5) | |
| 146 | ||
| 147 | // checks whether a jump tour exists at all | |
| 148 | ||
| 149 | def search(xs: List[Int]) : Boolean = xs match {
 | |
| 150 | case Nil => true | |
| 321 | 151 | case x::xs => | 
| 152 | if (xs.length < x) true | |
| 153 | else moves(xs, x).exists(search(_)) | |
| 217 | 154 | } | 
| 155 | ||
| 156 | ||
| 157 | search(List(5,3,2,5,1,1)) | |
| 158 | search(List(3,5,1,0,0,0,1)) | |
| 159 | search(List(3,5,1,0,0,0,0,1)) | |
| 160 | search(List(3,5,1,0,0,0,1,1)) | |
| 161 | search(List(3,5,1)) | |
| 162 | search(List(5,1,1)) | |
| 163 | search(Nil) | |
| 164 | search(List(1)) | |
| 165 | search(List(5,1,1)) | |
| 166 | search(List(3,5,1,0,0,0,0,0,0,0,0,1)) | |
| 167 | ||
| 168 | // generates *all* jump tours | |
| 321 | 169 | // if we are only interested in the shortest one, we could | 
| 217 | 170 | // shortcircut the calculation and only return List(x) in | 
| 171 | // case where xs.length < x, because no tour can be shorter | |
| 172 | // than 1 | |
| 173 | // | |
| 174 | ||
| 175 | def jumps(xs: List[Int]) : List[List[Int]] = xs match {
 | |
| 176 | case Nil => Nil | |
| 321 | 177 |   case x::xs => {
 | 
| 217 | 178 | val children = moves(xs, x) | 
| 320 | 179 | val results = children.map(cs => jumps(cs).map(x :: _)).flatten | 
| 180 | if (xs.length < x) List(x)::results else results | |
| 217 | 181 | } | 
| 182 | } | |
| 183 | ||
| 320 | 184 | jumps(List(5,3,2,5,1,1)).minBy(_.length) | 
| 217 | 185 | jumps(List(3,5,1,2,1,2,1)) | 
| 186 | jumps(List(3,5,1,2,3,4,1)) | |
| 187 | jumps(List(3,5,1,0,0,0,1)) | |
| 188 | jumps(List(3,5,1)) | |
| 189 | jumps(List(5,1,1)) | |
| 190 | jumps(Nil) | |
| 191 | jumps(List(1)) | |
| 192 | jumps(List(5,1,2)) | |
| 193 | moves(List(1,2), 5) | |
| 194 | jumps(List(1,5,1,2)) | |
| 195 | jumps(List(3,5,1,0,0,0,0,0,0,0,0,1)) | |
| 196 | ||
| 197 | jumps(List(5,3,2,5,1,1)).minBy(_.length) | |
| 198 | jumps(List(1,3,5,8,9,2,6,7,6,8,9)).minBy(_.length) | |
| 199 | jumps(List(1,3,6,1,0,9)).minBy(_.length) | |
| 200 | jumps(List(2,3,1,1,2,4,2,0,1,1)).minBy(_.length) | |
| 201 | ||
| 202 | ||
| 203 | ||
| 318 | 204 | |
| 205 | ||
| 206 | ||
| 320 | 207 | // User-defined Datatypes | 
| 208 | //======================== | |
| 209 | ||
| 323 | 210 | abstract class Tree | 
| 211 | case class Leaf(x: Int) extends Tree | |
| 212 | case class Node(s: String, left: Tree, right: Tree) extends Tree | |
| 213 | ||
| 214 | List(Leaf(20), Node("foo", Leaf(1), Leaf(2)))
 | |
| 320 | 215 | |
| 321 | 216 | sealed abstract class Colour | 
| 320 | 217 | case object Red extends Colour | 
| 218 | case object Green extends Colour | |
| 219 | case object Blue extends Colour | |
| 323 | 220 | case object Yellow extends Colour | 
| 320 | 221 | |
| 222 | ||
| 223 | def fav_colour(c: Colour) : Boolean = c match {
 | |
| 224 | case Green => true | |
| 323 | 225 | case _ => false | 
| 320 | 226 | } | 
| 227 | ||
| 228 | fav_colour(Green) | |
| 229 | ||
| 230 | // ... a tiny bit more useful: Roman Numerals | |
| 231 | ||
| 321 | 232 | sealed abstract class RomanDigit | 
| 320 | 233 | case object I extends RomanDigit | 
| 234 | case object V extends RomanDigit | |
| 235 | case object X extends RomanDigit | |
| 236 | case object L extends RomanDigit | |
| 237 | case object C extends RomanDigit | |
| 238 | case object D extends RomanDigit | |
| 239 | case object M extends RomanDigit | |
| 240 | ||
| 241 | type RomanNumeral = List[RomanDigit] | |
| 242 | ||
| 323 | 243 | List(X,I,M,D) | 
| 320 | 244 | |
| 245 | /* | |
| 246 | I -> 1 | |
| 247 | II -> 2 | |
| 248 | III -> 3 | |
| 249 | IV -> 4 | |
| 250 | V -> 5 | |
| 251 | VI -> 6 | |
| 252 | VII -> 7 | |
| 253 | VIII -> 8 | |
| 254 | IX -> 9 | |
| 255 | X -> 10 | |
| 256 | */ | |
| 257 | ||
| 258 | def RomanNumeral2Int(rs: RomanNumeral): Int = rs match { 
 | |
| 259 | case Nil => 0 | |
| 260 | case M::r => 1000 + RomanNumeral2Int(r) | |
| 261 | case C::M::r => 900 + RomanNumeral2Int(r) | |
| 262 | case D::r => 500 + RomanNumeral2Int(r) | |
| 263 | case C::D::r => 400 + RomanNumeral2Int(r) | |
| 264 | case C::r => 100 + RomanNumeral2Int(r) | |
| 265 | case X::C::r => 90 + RomanNumeral2Int(r) | |
| 266 | case L::r => 50 + RomanNumeral2Int(r) | |
| 267 | case X::L::r => 40 + RomanNumeral2Int(r) | |
| 268 | case X::r => 10 + RomanNumeral2Int(r) | |
| 269 | case I::X::r => 9 + RomanNumeral2Int(r) | |
| 270 | case V::r => 5 + RomanNumeral2Int(r) | |
| 271 | case I::V::r => 4 + RomanNumeral2Int(r) | |
| 272 | case I::r => 1 + RomanNumeral2Int(r) | |
| 273 | } | |
| 274 | ||
| 275 | RomanNumeral2Int(List(I,V)) // 4 | |
| 276 | RomanNumeral2Int(List(I,I,I,I)) // 4 (invalid Roman number) | |
| 277 | RomanNumeral2Int(List(V,I)) // 6 | |
| 278 | RomanNumeral2Int(List(I,X)) // 9 | |
| 279 | RomanNumeral2Int(List(M,C,M,L,X,X,I,X)) // 1979 | |
| 280 | RomanNumeral2Int(List(M,M,X,V,I,I)) // 2017 | |
| 281 | ||
| 282 | ||
| 283 | // String interpolations as patterns | |
| 284 | ||
| 285 | val date = "2019-11-26" | |
| 286 | val s"$year-$month-$day" = date | |
| 287 | ||
| 288 | def parse_date(date: String) : Option[(Int, Int, Int)]= date match {
 | |
| 289 | case s"$year-$month-$day" => Some((day.toInt, month.toInt, year.toInt)) | |
| 290 | case s"$day/$month/$year" => Some((day.toInt, month.toInt, year.toInt)) | |
| 291 | case s"$day.$month.$year" => Some((day.toInt, month.toInt, year.toInt)) | |
| 292 | case _ => None | |
| 293 | } | |
| 318 | 294 | |
| 320 | 295 | parse_date("2019-11-26")
 | 
| 296 | parse_date("26/11/2019")
 | |
| 297 | parse_date("26.11.2019")
 | |
| 298 | ||
| 299 | ||
| 300 | // User-defined Datatypes and Pattern Matching | |
| 301 | //============================================= | |
| 302 | ||
| 303 | // trees | |
| 304 | ||
| 323 | 305 | |
| 306 | ||
| 307 | // expressions | |
| 308 | ||
| 321 | 309 | sealed abstract class Exp | 
| 320 | 310 | case class N(n: Int) extends Exp // for numbers | 
| 311 | case class Plus(e1: Exp, e2: Exp) extends Exp | |
| 312 | case class Times(e1: Exp, e2: Exp) extends Exp | |
| 313 | ||
| 314 | def string(e: Exp) : String = e match {
 | |
| 315 | case N(n) => s"$n" | |
| 316 |   case Plus(e1, e2) => s"(${string(e1)} + ${string(e2)})" 
 | |
| 317 |   case Times(e1, e2) => s"(${string(e1)} * ${string(e2)})"
 | |
| 318 | } | |
| 319 | ||
| 320 | val e = Plus(N(9), Times(N(3), N(4))) | |
| 321 | println(string(e)) | |
| 322 | ||
| 323 | def eval(e: Exp) : Int = e match {
 | |
| 324 | case N(n) => n | |
| 325 | case Plus(e1, e2) => eval(e1) + eval(e2) | |
| 326 | case Times(e1, e2) => eval(e1) * eval(e2) | |
| 327 | } | |
| 328 | ||
| 329 | println(eval(e)) | |
| 330 | ||
| 331 | def simp(e: Exp) : Exp = e match {
 | |
| 332 | case N(n) => N(n) | |
| 333 |   case Plus(e1, e2) => (simp(e1), simp(e2)) match {
 | |
| 334 | case (N(0), e2s) => e2s | |
| 335 | case (e1s, N(0)) => e1s | |
| 336 | case (e1s, e2s) => Plus(e1s, e2s) | |
| 337 | } | |
| 338 |   case Times(e1, e2) => (simp(e1), simp(e2)) match {
 | |
| 339 | case (N(0), _) => N(0) | |
| 340 | case (_, N(0)) => N(0) | |
| 341 | case (N(1), e2s) => e2s | |
| 342 | case (e1s, N(1)) => e1s | |
| 343 | case (e1s, e2s) => Times(e1s, e2s) | |
| 344 | } | |
| 345 | } | |
| 346 | ||
| 347 | ||
| 348 | val e2 = Times(Plus(N(0), N(1)), Plus(N(0), N(9))) | |
| 349 | println(string(e2)) | |
| 350 | println(string(simp(e2))) | |
| 351 | ||
| 352 | ||
| 353 | // Tokens and Reverse Polish Notation | |
| 321 | 354 | sealed abstract class Token | 
| 320 | 355 | case class T(n: Int) extends Token | 
| 356 | case object PL extends Token | |
| 357 | case object TI extends Token | |
| 358 | ||
| 359 | def rp(e: Exp) : List[Token] = e match {
 | |
| 360 | case N(n) => List(T(n)) | |
| 361 | case Plus(e1, e2) => rp(e1) ::: rp(e2) ::: List(PL) | |
| 362 | case Times(e1, e2) => rp(e1) ::: rp(e2) ::: List(TI) | |
| 363 | } | |
| 364 | println(string(e2)) | |
| 365 | println(rp(e2)) | |
| 366 | ||
| 367 | def comp(ls: List[Token], st: List[Int]) : Int = (ls, st) match {
 | |
| 368 | case (Nil, st) => st.head | |
| 369 | case (T(n)::rest, st) => comp(rest, n::st) | |
| 370 | case (PL::rest, n1::n2::st) => comp(rest, n1 + n2::st) | |
| 371 | case (TI::rest, n1::n2::st) => comp(rest, n1 * n2::st) | |
| 372 | } | |
| 373 | ||
| 374 | comp(rp(e), Nil) | |
| 375 | ||
| 376 | def proc(s: String) : Token = s match {
 | |
| 377 | case "+" => PL | |
| 378 | case "*" => TI | |
| 379 | case _ => T(s.toInt) | |
| 380 | } | |
| 381 | ||
| 382 | comp("1 2 + 4 * 5 + 3 +".split(" ").toList.map(proc), Nil)
 | |
| 217 | 383 | |
| 384 | ||
| 385 | ||
| 386 | ||
| 387 | // Sudoku | |
| 388 | //======== | |
| 389 | ||
| 390 | // THE POINT OF THIS CODE IS NOT TO BE SUPER | |
| 391 | // EFFICIENT AND FAST, just explaining exhaustive | |
| 392 | // depth-first search | |
| 393 | ||
| 155 | 394 | |
| 395 | val game0 = """.14.6.3.. | |
| 396 | |62...4..9 | |
| 397 | |.8..5.6.. | |
| 398 | |.6.2....3 | |
| 399 | |.7..1..5. | |
| 400 | |5....9.6. | |
| 401 | |..6.2..3. | |
| 402 | |1..5...92 | |
| 403 |               |..7.9.41.""".stripMargin.replaceAll("\\n", "")
 | |
| 53 | 404 | |
| 155 | 405 | type Pos = (Int, Int) | 
| 406 | val EmptyValue = '.' | |
| 407 | val MaxValue = 9 | |
| 408 | ||
| 409 | val allValues = "123456789".toList | |
| 410 | val indexes = (0 to 8).toList | |
| 411 | ||
| 412 | ||
| 413 | def empty(game: String) = game.indexOf(EmptyValue) | |
| 414 | def isDone(game: String) = empty(game) == -1 | |
| 415 | def emptyPosition(game: String) = | |
| 416 | (empty(game) % MaxValue, empty(game) / MaxValue) | |
| 417 | ||
| 67 | 418 | |
| 155 | 419 | def get_row(game: String, y: Int) = | 
| 420 | indexes.map(col => game(y * MaxValue + col)) | |
| 421 | def get_col(game: String, x: Int) = | |
| 422 | indexes.map(row => game(x + row * MaxValue)) | |
| 423 | ||
| 424 | def get_box(game: String, pos: Pos): List[Char] = {
 | |
| 425 | def base(p: Int): Int = (p / 3) * 3 | |
| 426 | val x0 = base(pos._1) | |
| 427 | val y0 = base(pos._2) | |
| 428 | val ys = (y0 until y0 + 3).toList | |
| 429 | (x0 until x0 + 3).toList.flatMap(x => ys.map(y => game(x + y * MaxValue))) | |
| 430 | } | |
| 431 | ||
| 217 | 432 | //get_row(game0, 0) | 
| 433 | //get_row(game0, 1) | |
| 218 | 434 | //get_col(game0, 0) | 
| 435 | //get_box(game0, (3, 1)) | |
| 217 | 436 | |
| 437 | ||
| 155 | 438 | // this is not mutable!! | 
| 439 | def update(game: String, pos: Int, value: Char): String = | |
| 440 | game.updated(pos, value) | |
| 441 | ||
| 442 | def toAvoid(game: String, pos: Pos): List[Char] = | |
| 443 | (get_col(game, pos._1) ++ get_row(game, pos._2) ++ get_box(game, pos)) | |
| 444 | ||
| 445 | def candidates(game: String, pos: Pos): List[Char] = | |
| 218 | 446 | allValues.diff(toAvoid(game, pos)) | 
| 155 | 447 | |
| 448 | //candidates(game0, (0,0)) | |
| 449 | ||
| 450 | def pretty(game: String): String = | |
| 218 | 451 |   "\n" + (game.sliding(MaxValue, MaxValue).mkString("\n"))
 | 
| 155 | 452 | |
| 218 | 453 | |
| 155 | 454 | def search(game: String): List[String] = {
 | 
| 455 | if (isDone(game)) List(game) | |
| 456 |   else {
 | |
| 457 | val cs = candidates(game, emptyPosition(game)) | |
| 320 | 458 | cs.map(c => search(update(game, empty(game), c))).toList.flatten | 
| 67 | 459 | } | 
| 460 | } | |
| 461 | ||
| 217 | 462 | search(game0).map(pretty) | 
| 463 | ||
| 464 | val game1 = """23.915... | |
| 465 | |...2..54. | |
| 466 | |6.7...... | |
| 467 | |..1.....9 | |
| 468 | |89.5.3.17 | |
| 469 | |5.....6.. | |
| 470 | |......9.5 | |
| 471 | |.16..7... | |
| 472 |               |...329..1""".stripMargin.replaceAll("\\n", "")
 | |
| 473 | ||
| 474 | ||
| 475 | // game that is in the hard category | |
| 476 | val game2 = """8........ | |
| 477 | |..36..... | |
| 478 | |.7..9.2.. | |
| 479 | |.5...7... | |
| 480 | |....457.. | |
| 481 | |...1...3. | |
| 482 | |..1....68 | |
| 483 | |..85...1. | |
| 484 |               |.9....4..""".stripMargin.replaceAll("\\n", "")
 | |
| 485 | ||
| 486 | // game with multiple solutions | |
| 487 | val game3 = """.8...9743 | |
| 488 | |.5...8.1. | |
| 489 | |.1....... | |
| 490 | |8....5... | |
| 491 | |...8.4... | |
| 492 | |...3....6 | |
| 493 | |.......7. | |
| 494 | |.3.5...8. | |
| 495 |               |9724...5.""".stripMargin.replaceAll("\\n", "")
 | |
| 496 | ||
| 497 | ||
| 498 | search(game1).map(pretty) | |
| 499 | search(game3).map(pretty) | |
| 500 | search(game2).map(pretty) | |
| 501 | ||
| 502 | // for measuring time | |
| 503 | def time_needed[T](i: Int, code: => T) = {
 | |
| 504 | val start = System.nanoTime() | |
| 505 | for (j <- 1 to i) code | |
| 506 | val end = System.nanoTime() | |
| 507 | ((end - start) / 1.0e9) + " secs" | |
| 508 | } | |
| 509 | ||
| 510 | time_needed(1, search(game2)) | |
| 511 | ||
| 320 | 512 | |
| 513 | ||
| 514 | ||
| 515 | // Tail recursion | |
| 516 | //================ | |
| 517 | ||
| 518 | ||
| 519 | def fact(n: Long): Long = | |
| 520 | if (n == 0) 1 else n * fact(n - 1) | |
| 521 | ||
| 522 | def factB(n: BigInt): BigInt = | |
| 523 | if (n == 0) 1 else n * factB(n - 1) | |
| 524 | ||
| 525 | factB(100000) | |
| 526 | ||
| 527 | fact(10) //ok | |
| 528 | fact(10000) // produces a stackoverflow | |
| 529 | ||
| 530 | def factT(n: BigInt, acc: BigInt): BigInt = | |
| 531 | if (n == 0) acc else factT(n - 1, n * acc) | |
| 532 | ||
| 533 | factT(10, 1) | |
| 534 | println(factT(100000, 1)) | |
| 535 | ||
| 536 | // there is a flag for ensuring a function is tail recursive | |
| 537 | import scala.annotation.tailrec | |
| 538 | ||
| 539 | @tailrec | |
| 540 | def factT(n: BigInt, acc: BigInt): BigInt = | |
| 541 | if (n == 0) acc else factT(n - 1, n * acc) | |
| 542 | ||
| 543 | ||
| 544 | ||
| 545 | // for tail-recursive functions the Scala compiler | |
| 546 | // generates loop-like code, which does not need | |
| 547 | // to allocate stack-space in each recursive | |
| 548 | // call; Scala can do this only for tail-recursive | |
| 549 | // functions | |
| 550 | ||
| 155 | 551 | // tail recursive version that searches | 
| 158 | 552 | // for all solutions | 
| 553 | ||
| 155 | 554 | def searchT(games: List[String], sols: List[String]): List[String] = games match {
 | 
| 555 | case Nil => sols | |
| 556 |   case game::rest => {
 | |
| 557 | if (isDone(game)) searchT(rest, game::sols) | |
| 558 |     else {
 | |
| 559 | val cs = candidates(game, emptyPosition(game)) | |
| 560 | searchT(cs.map(c => update(game, empty(game), c)) ::: rest, sols) | |
| 561 | } | |
| 562 | } | |
| 67 | 563 | } | 
| 564 | ||
| 158 | 565 | searchT(List(game3), List()).map(pretty) | 
| 566 | ||
| 567 | ||
| 155 | 568 | // tail recursive version that searches | 
| 569 | // for a single solution | |
| 158 | 570 | |
| 155 | 571 | def search1T(games: List[String]): Option[String] = games match {
 | 
| 67 | 572 | case Nil => None | 
| 155 | 573 |   case game::rest => {
 | 
| 574 | if (isDone(game)) Some(game) | |
| 575 |     else {
 | |
| 576 | val cs = candidates(game, emptyPosition(game)) | |
| 577 | search1T(cs.map(c => update(game, empty(game), c)) ::: rest) | |
| 578 | } | |
| 579 | } | |
| 67 | 580 | } | 
| 581 | ||
| 158 | 582 | search1T(List(game3)).map(pretty) | 
| 217 | 583 | time_needed(10, search1T(List(game3))) | 
| 584 | ||
| 158 | 585 | |
| 155 | 586 | // game with multiple solutions | 
| 587 | val game3 = """.8...9743 | |
| 588 | |.5...8.1. | |
| 589 | |.1....... | |
| 590 | |8....5... | |
| 591 | |...8.4... | |
| 592 | |...3....6 | |
| 593 | |.......7. | |
| 594 | |.3.5...8. | |
| 595 |               |9724...5.""".stripMargin.replaceAll("\\n", "")
 | |
| 596 | ||
| 158 | 597 | searchT(List(game3), Nil).map(pretty) | 
| 155 | 598 | search1T(List(game3)).map(pretty) | 
| 67 | 599 | |
| 77 
3cbe3d90b77f
updated
 Christian Urban <christian dot urban at kcl dot ac dot uk> parents: 
73diff
changeset | 600 | // Moral: Whenever a recursive function is resource-critical | 
| 158 | 601 | // (i.e. works with large recursion depth), then you need to | 
| 77 
3cbe3d90b77f
updated
 Christian Urban <christian dot urban at kcl dot ac dot uk> parents: 
73diff
changeset | 602 | // write it in tail-recursive fashion. | 
| 
3cbe3d90b77f
updated
 Christian Urban <christian dot urban at kcl dot ac dot uk> parents: 
73diff
changeset | 603 | // | 
| 155 | 604 | // Unfortuantely, Scala because of current limitations in | 
| 605 | // the JVM is not as clever as other functional languages. It can | |
| 77 
3cbe3d90b77f
updated
 Christian Urban <christian dot urban at kcl dot ac dot uk> parents: 
73diff
changeset | 606 | // only optimise "self-tail calls". This excludes the cases of | 
| 
3cbe3d90b77f
updated
 Christian Urban <christian dot urban at kcl dot ac dot uk> parents: 
73diff
changeset | 607 | // multiple functions making tail calls to each other. Well, | 
| 
3cbe3d90b77f
updated
 Christian Urban <christian dot urban at kcl dot ac dot uk> parents: 
73diff
changeset | 608 | // nothing is perfect. | 
| 
3cbe3d90b77f
updated
 Christian Urban <christian dot urban at kcl dot ac dot uk> parents: 
73diff
changeset | 609 | |
| 
3cbe3d90b77f
updated
 Christian Urban <christian dot urban at kcl dot ac dot uk> parents: 
73diff
changeset | 610 | |
| 67 | 611 | |
| 612 | ||
| 71 | 613 | |
| 67 | 614 |