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 | 
         |