|      1 // A Version with an explicit n-times regular expression; |         | 
|      2 // this keeps the size of the regular expression in the |         | 
|      3 // EVIL1 test-case quite small |         | 
|      4  |         | 
|      5 abstract class Rexp  |         | 
|      6 case object ZERO extends Rexp |         | 
|      7 case object ONE extends Rexp |         | 
|      8 case class CHAR(c: Char) extends Rexp |         | 
|      9 case class ALT(r1: Rexp, r2: Rexp) extends Rexp  |         | 
|     10 case class SEQ(r1: Rexp, r2: Rexp) extends Rexp  |         | 
|     11 case class STAR(r: Rexp) extends Rexp  |         | 
|     12 case class NTIMES(r: Rexp, n: Int) extends Rexp   //explicit constructor for n-times |         | 
|     13  |         | 
|     14  |         | 
|     15 def nullable (r: Rexp) : Boolean = r match { |         | 
|     16   case ZERO => false |         | 
|     17   case ONE => true |         | 
|     18   case CHAR(_) => false |         | 
|     19   case ALT(r1, r2) => nullable(r1) || nullable(r2) |         | 
|     20   case SEQ(r1, r2) => nullable(r1) && nullable(r2) |         | 
|     21   case STAR(_) => true |         | 
|     22   case NTIMES(r, i) => if (i == 0) true else nullable(r) |         | 
|     23 } |         | 
|     24  |         | 
|     25  |         | 
|     26 def der (c: Char, r: Rexp) : Rexp = r match { |         | 
|     27   case ZERO => ZERO |         | 
|     28   case ONE => ZERO |         | 
|     29   case CHAR(d) => if (c == d) ONE else ZERO |         | 
|     30   case ALT(r1, r2) => ALT(der(c, r1), der(c, r2)) |         | 
|     31   case SEQ(r1, r2) =>  |         | 
|     32     if (nullable(r1)) ALT(SEQ(der(c, r1), r2), der(c, r2)) |         | 
|     33     else SEQ(der(c, r1), r2) |         | 
|     34   case STAR(r1) => SEQ(der(c, r1), STAR(r1)) |         | 
|     35   case NTIMES(r1, i) =>  |         | 
|     36     if (i == 0) ZERO else SEQ(der(c, r1), NTIMES(r1, i - 1)) |         | 
|     37 } |         | 
|     38  |         | 
|     39 def ders (s: List[Char], r: Rexp) : Rexp = s match { |         | 
|     40   case Nil => r |         | 
|     41   case c::s => ders(s, der(c, r)) |         | 
|     42 } |         | 
|     43  |         | 
|     44 def matcher(r: Rexp, s: String) : Boolean = nullable(ders(s.toList, r)) |         | 
|     45  |         | 
|     46  |         | 
|     47 // the optional regular expression: one or zero times |         | 
|     48 // this regular expression is still defined in terms of ALT |         | 
|     49 def OPT(r: Rexp) = ALT(r, ONE) |         | 
|     50  |         | 
|     51  |         | 
|     52 // Test Cases |         | 
|     53  |         | 
|     54 // evil regular expressions |         | 
|     55 def EVIL1(n: Int) = SEQ(NTIMES(OPT(CHAR('a')), n), NTIMES(CHAR('a'), n)) |         | 
|     56 val EVIL2 = SEQ(STAR(STAR(CHAR('a'))), CHAR('b')) |         | 
|     57  |         | 
|     58 def time_needed[T](i: Int, code: => T) = { |         | 
|     59   val start = System.nanoTime() |         | 
|     60   for (j <- 1 to i) code |         | 
|     61   val end = System.nanoTime() |         | 
|     62   (end - start) / (i * 1.0e9) |         | 
|     63 } |         | 
|     64  |         | 
|     65  |         | 
|     66  |         | 
|     67 // test: (a?{n}) (a{n}) |         | 
|     68 for (i <- 0 to 1000 by 100) { |         | 
|     69   println(f"$i: ${time_needed(2, matcher(EVIL1(i), "a" * i))}%.5f") |         | 
|     70 } |         | 
|     71  |         | 
|     72 // test: (a*)* b |         | 
|     73 for (i <- 0 to 20) { |         | 
|     74   println(f"$i: ${time_needed(2, matcher(EVIL2, "a" * i))}%.5f") |         | 
|     75 } |         | 
|     76  |         | 
|     77  |         | 
|     78 // the size of a regular expressions - for testing purposes  |         | 
|     79 def size(r: Rexp) : Int = r match { |         | 
|     80   case ZERO => 1 |         | 
|     81   case ONE => 1 |         | 
|     82   case CHAR(_) => 1 |         | 
|     83   case ALT(r1, r2) => 1 + size(r1) + size(r2) |         | 
|     84   case SEQ(r1, r2) => 1 + size(r1) + size(r2) |         | 
|     85   case STAR(r) => 1 + size(r) |         | 
|     86   case NTIMES(r, _) => 1 + size(r) |         | 
|     87 } |         | 
|     88  |         | 
|     89 // EVIL1(n) has now a constant size, no matter |         | 
|     90 // what n is; also the derivative only grows  |         | 
|     91 // moderately  |         | 
|     92  |         | 
|     93 size(EVIL1(1))  // 7 |         | 
|     94 size(EVIL1(3))  // 7 |         | 
|     95 size(EVIL1(5))  // 7 |         | 
|     96 size(EVIL1(7))  // 7 |         | 
|     97 size(EVIL1(20)) // 7 |         | 
|     98  |         | 
|     99 size(ders("".toList, EVIL1(5)))       // 7 |         | 
|    100 size(ders("a".toList, EVIL1(5)))      // 16 |         | 
|    101 size(ders("aa".toList, EVIL1(5)))     // 35 |         | 
|    102 size(ders("aaa".toList, EVIL1(5)))    // 59 |         | 
|    103 size(ders("aaaa".toList, EVIL1(5)))   // 88 |         | 
|    104 size(ders("aaaaa".toList, EVIL1(5)))  // 122 |         | 
|    105 size(ders("aaaaaa".toList, EVIL1(5))) // 151 |         | 
|    106  |         | 
|    107 // but the size of the derivatives can still grow  |         | 
|    108 // quite dramatically in case of EVIL2 |         | 
|    109  |         | 
|    110 size(ders("".toList, EVIL2))       // 5 |         | 
|    111 size(ders("a".toList, EVIL2))      // 12 |         | 
|    112 size(ders("aa".toList, EVIL2))     // 28 |         | 
|    113 size(ders("aaa".toList, EVIL2))    // 58 |         | 
|    114 size(ders("aaaa".toList, EVIL2))   // 116 |         | 
|    115 size(ders("aaaaa".toList, EVIL2))  // 230 |         | 
|    116 size(ders("aaaaaa".toList, EVIL2)) // 456 |         | 
|    117  |         | 
|    118 size(ders(("a" * 20).toList, EVIL2)) // 7340068 |         |