progs/fun/defs.fun
author Christian Urban <christian.urban@kcl.ac.uk>
Sun, 19 Oct 2025 09:51:35 +0200
changeset 1012 c01dfa3ff177
parent 1010 ae9ffbf979ff
permissions -rw-r--r--
updated
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
224
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
     1
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
     2
def zero(x) = 0;
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
     3
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
     4
def suc(x) = x + 1;
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
     5
656
cfc0e730bcda updated
Christian Urban <urbanc@in.tum.de>
parents: 626
diff changeset
     6
def pred(x) =
cfc0e730bcda updated
Christian Urban <urbanc@in.tum.de>
parents: 626
diff changeset
     7
  if x == 0 then x else x - 1;
224
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
     8
656
cfc0e730bcda updated
Christian Urban <urbanc@in.tum.de>
parents: 626
diff changeset
     9
def add(x, y) =
cfc0e730bcda updated
Christian Urban <urbanc@in.tum.de>
parents: 626
diff changeset
    10
  if x == 0 then y else suc(add(x - 1, y));
224
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    11
1010
ae9ffbf979ff updated
Christian Urban <christian.urban@kcl.ac.uk>
parents: 735
diff changeset
    12
def add2(x, y) =
ae9ffbf979ff updated
Christian Urban <christian.urban@kcl.ac.uk>
parents: 735
diff changeset
    13
  if x == 0 then y else add2(x - 1, y + 1); 
ae9ffbf979ff updated
Christian Urban <christian.urban@kcl.ac.uk>
parents: 735
diff changeset
    14
  
656
cfc0e730bcda updated
Christian Urban <urbanc@in.tum.de>
parents: 626
diff changeset
    15
def mult(x, y) =
cfc0e730bcda updated
Christian Urban <urbanc@in.tum.de>
parents: 626
diff changeset
    16
  if x == 0 then 0 else add(y, mult(x - 1, y));
224
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    17
656
cfc0e730bcda updated
Christian Urban <urbanc@in.tum.de>
parents: 626
diff changeset
    18
def pow(x, y) =
cfc0e730bcda updated
Christian Urban <urbanc@in.tum.de>
parents: 626
diff changeset
    19
  if y == 0 then 1 else mult(x, pow(x, y - 1));
224
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    20
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    21
def fib(n) = if n == 0 then 0 
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    22
             else if n == 1 then 1 
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    23
             else fib(n - 1) + fib(n - 2);
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    24
656
cfc0e730bcda updated
Christian Urban <urbanc@in.tum.de>
parents: 626
diff changeset
    25
def fact(n) =
cfc0e730bcda updated
Christian Urban <urbanc@in.tum.de>
parents: 626
diff changeset
    26
  if n == 0 then 1 else n * fact(n - 1);
224
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    27
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    28
def ack(m, n) = if m == 0 then n + 1
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    29
                else if n == 0 then ack(m - 1, 1)
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    30
                else ack(m - 1, ack(m, n - 1));
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    31
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    32
def stack_test(x) = x + 1 + 2 + 3 + 4 + 5;
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    33
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    34
def div(x, y) = x / y;   //integer division
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    35
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    36
def rem(x, y) = x % y;   //remainder
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    37
656
cfc0e730bcda updated
Christian Urban <urbanc@in.tum.de>
parents: 626
diff changeset
    38
def gcd(a, b) =
cfc0e730bcda updated
Christian Urban <urbanc@in.tum.de>
parents: 626
diff changeset
    39
  if b == 0 then a else gcd(b, a % b);
224
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    40
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    41
def is_prime_aux(n, i) = 
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    42
  if n % i == 0 then 0
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    43
  else if (i * i) <= n then is_prime_aux(n, i + 1)  
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    44
  else 1;
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    45
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    46
def is_prime(n) = if n == 2 then 1 else is_prime_aux(n, 2);
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    47
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    48
def primes(n) = 
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    49
  if n == 0 then 0
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    50
  else if is_prime(n) == 1 then (write n; primes(n - 1)) 
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    51
  else primes(n - 1);
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    52
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    53
def is_collatz(n) =
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    54
  if n == 1 then 1
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    55
  else if n % 2 == 0 then is_collatz(n / 2)
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    56
  else is_collatz(3 * n + 1);  
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    57
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    58
def collatz_aux(n, i) = 
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    59
  if i > n then 0
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    60
  else if is_collatz(i) == 1 then (write i; collatz_aux(n, i + 1)) 
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    61
  else collatz_aux(n, i + 1);
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    62
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    63
def collatz(n) = collatz_aux(n, 1);
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    64
656
cfc0e730bcda updated
Christian Urban <urbanc@in.tum.de>
parents: 626
diff changeset
    65
def facT(n, acc) =
cfc0e730bcda updated
Christian Urban <urbanc@in.tum.de>
parents: 626
diff changeset
    66
  if n == 0 then acc else facT(n - 1, n * acc);
224
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    67
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    68
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    69
//zero(3)
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    70
//suc(8)
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    71
//pred(7)
656
cfc0e730bcda updated
Christian Urban <urbanc@in.tum.de>
parents: 626
diff changeset
    72
//write(add(3, 4))
224
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    73
//mult(4,5)
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    74
//pow(2, 3)
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    75
//fib(20)
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    76
//(write(fact(5)) ; fact(6))
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    77
//(write(1) ; 2)
656
cfc0e730bcda updated
Christian Urban <urbanc@in.tum.de>
parents: 626
diff changeset
    78
//write(ack(3, 12))   // for tail-rec test
224
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    79
//stack_test(0)
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    80
//(write (div(11, 3)); rem(11, 3))
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    81
//gcd(54, 24)
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    82
//is_prime(2)
695
484b74bc057e deleted init function from boilerplate code
Christian Urban <urbanc@in.tum.de>
parents: 656
diff changeset
    83
primes(1000)
484b74bc057e deleted init function from boilerplate code
Christian Urban <urbanc@in.tum.de>
parents: 656
diff changeset
    84
//primes(1000000)
484b74bc057e deleted init function from boilerplate code
Christian Urban <urbanc@in.tum.de>
parents: 656
diff changeset
    85
//collatz(4000)
224
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    86
//collatz(5000)
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    87
//facT(6, 1)
Christian Urban <christian dot urban at kcl dot ac dot uk>
parents:
diff changeset
    88
1010
ae9ffbf979ff updated
Christian Urban <christian.urban@kcl.ac.uk>
parents: 735
diff changeset
    89