2024-11-13 21:50:44 +01:00
|
|
|
open MiniImp
|
|
|
|
|
|
|
|
|
|
let get_result x =
|
2024-11-15 17:23:04 +01:00
|
|
|
Lexing.from_string x |> Parser.prg Lexer.lex |> Semantics.reduce
|
2024-11-13 21:50:44 +01:00
|
|
|
|
|
|
|
|
(* -------------------------------------------------------------------------- *)
|
|
|
|
|
(* Identity program *)
|
|
|
|
|
let program =
|
2024-11-16 15:40:00 +01:00
|
|
|
"def main with input a output b as b := a"
|
2024-11-13 21:50:44 +01:00
|
|
|
;;
|
|
|
|
|
|
2024-11-16 15:40:00 +01:00
|
|
|
Printf.printf "Identity program: ";
|
|
|
|
|
match get_result program 1 with
|
|
|
|
|
Ok d -> Printf.printf "%d\n" d
|
|
|
|
|
| Error `AbsentAssignment msg -> Printf.printf "error -> %s\n" msg
|
|
|
|
|
| Error `DivisionByZero msg -> Printf.printf "error -> %s\n" msg
|
|
|
|
|
;;
|
2024-11-13 21:50:44 +01:00
|
|
|
|
|
|
|
|
(* -------------------------------------------------------------------------- *)
|
|
|
|
|
(* y not defined program *)
|
|
|
|
|
let program =
|
2024-11-16 15:40:00 +01:00
|
|
|
"def main with input a output b as x := 1; b := a + x + y"
|
2024-11-13 21:50:44 +01:00
|
|
|
;;
|
|
|
|
|
|
2024-11-16 15:40:00 +01:00
|
|
|
Printf.printf "y not defined program: ";
|
|
|
|
|
match get_result program 100 with
|
|
|
|
|
Ok d -> Printf.printf "error: %d\n" d
|
|
|
|
|
| Error `AbsentAssignment msg -> Printf.printf "%s\n" msg
|
|
|
|
|
| Error `DivisionByZero msg -> Printf.printf "error -> %s\n" msg
|
2024-11-13 21:50:44 +01:00
|
|
|
;;
|
|
|
|
|
|
|
|
|
|
(* -------------------------------------------------------------------------- *)
|
|
|
|
|
(* Factorial program *)
|
|
|
|
|
let program =
|
2024-11-16 15:40:00 +01:00
|
|
|
"def main with input a output b as
|
2024-11-13 21:50:44 +01:00
|
|
|
b := 1;
|
2024-11-16 15:40:00 +01:00
|
|
|
for (i := 1, i <= a, i := i + 1) do
|
2024-11-13 21:50:44 +01:00
|
|
|
b := b * i;
|
|
|
|
|
"
|
|
|
|
|
;;
|
|
|
|
|
|
2024-11-16 15:40:00 +01:00
|
|
|
Printf.printf "Factorial program: ";
|
|
|
|
|
match get_result program 10 with
|
|
|
|
|
Ok d -> Printf.printf "%d\n" d
|
|
|
|
|
| Error `AbsentAssignment msg -> Printf.printf "error -> %s\n" msg
|
|
|
|
|
| Error `DivisionByZero msg -> Printf.printf "error -> %s\n" msg
|
2024-11-13 21:50:44 +01:00
|
|
|
;;
|
|
|
|
|
|
|
|
|
|
(* -------------------------------------------------------------------------- *)
|
|
|
|
|
(* Hailstone sequence's lenght program *)
|
|
|
|
|
let program =
|
2024-11-16 15:40:00 +01:00
|
|
|
"def main with input a output b as
|
2024-11-13 21:50:44 +01:00
|
|
|
b := 1;
|
2024-11-16 15:40:00 +01:00
|
|
|
while not a == 1 do (
|
2024-11-13 21:50:44 +01:00
|
|
|
b := b + 1;
|
2024-11-16 15:40:00 +01:00
|
|
|
if ((a % 2) == 1) then a := 3 * a + 1 else a := a / 2
|
|
|
|
|
)
|
|
|
|
|
"
|
2024-11-13 21:50:44 +01:00
|
|
|
;;
|
|
|
|
|
|
2024-11-16 15:40:00 +01:00
|
|
|
Printf.printf "Hailstone sequence's lenght program: ";
|
|
|
|
|
match get_result program 77031 with
|
|
|
|
|
Ok d -> Printf.printf "%d\n" d
|
|
|
|
|
| Error `AbsentAssignment msg -> Printf.printf "error -> %s\n" msg
|
|
|
|
|
| Error `DivisionByZero msg -> Printf.printf "error -> %s\n" msg
|
2024-11-13 21:50:44 +01:00
|
|
|
;;
|
|
|
|
|
|
|
|
|
|
(* -------------------------------------------------------------------------- *)
|
|
|
|
|
(* Sum multiples of 3 and 5 program *)
|
|
|
|
|
let program =
|
2024-11-16 15:40:00 +01:00
|
|
|
"def main with input a output b as
|
2024-11-13 21:50:44 +01:00
|
|
|
b := 0;
|
2024-11-16 15:40:00 +01:00
|
|
|
for (i := 0, i <= a, i := i+1) do
|
|
|
|
|
if (i % 3 == 0 || i % 5 == 0) then b := b + i;
|
|
|
|
|
else skip;
|
|
|
|
|
"
|
2024-11-13 21:50:44 +01:00
|
|
|
;;
|
|
|
|
|
|
2024-11-16 15:40:00 +01:00
|
|
|
Printf.printf "Sum multiples of 3 and 5 program: ";
|
|
|
|
|
match get_result program 12345 with
|
|
|
|
|
Ok d -> Printf.printf "%d\n" d
|
|
|
|
|
| Error `AbsentAssignment msg -> Printf.printf "error -> %s\n" msg
|
|
|
|
|
| Error `DivisionByZero msg -> Printf.printf "error -> %s\n" msg
|
2024-11-13 21:50:44 +01:00
|
|
|
;;
|
|
|
|
|
|
|
|
|
|
(* -------------------------------------------------------------------------- *)
|
|
|
|
|
(* Rand program *)
|
|
|
|
|
let program =
|
2024-11-16 15:40:00 +01:00
|
|
|
"def main with input a output b as b := rand(a)"
|
2024-11-13 21:50:44 +01:00
|
|
|
;;
|
|
|
|
|
|
2024-11-16 15:40:00 +01:00
|
|
|
Printf.printf "Rand program: ";
|
|
|
|
|
match get_result program 10 with
|
|
|
|
|
Ok d -> Printf.printf "%b\n" (d < 10)
|
|
|
|
|
| Error `AbsentAssignment msg -> Printf.printf "error -> %s\n" msg
|
|
|
|
|
| Error `DivisionByZero msg -> Printf.printf "error -> %s\n" msg
|
2024-11-13 21:50:44 +01:00
|
|
|
;;
|
|
|
|
|
|
|
|
|
|
(* -------------------------------------------------------------------------- *)
|
|
|
|
|
(* Fibonacci program *)
|
|
|
|
|
let program =
|
2024-11-16 15:40:00 +01:00
|
|
|
"def main with input n output fnext as
|
2024-11-13 21:50:44 +01:00
|
|
|
fnow := 0;
|
|
|
|
|
fnext := 1;
|
2024-11-16 15:40:00 +01:00
|
|
|
while (n > 1) do (
|
2024-11-13 21:50:44 +01:00
|
|
|
tmp := fnow + fnext;
|
|
|
|
|
fnow := fnext;
|
|
|
|
|
fnext := tmp;
|
2024-11-16 15:40:00 +01:00
|
|
|
n := n - 1;
|
|
|
|
|
)
|
|
|
|
|
"
|
2024-11-13 21:50:44 +01:00
|
|
|
;;
|
|
|
|
|
|
2024-11-16 15:40:00 +01:00
|
|
|
Printf.printf "Fibonacci program: ";
|
|
|
|
|
match get_result program 48 with
|
|
|
|
|
Ok d -> Printf.printf "%d\n" d
|
|
|
|
|
| Error `AbsentAssignment msg -> Printf.printf "error -> %s\n" msg
|
|
|
|
|
| Error `DivisionByZero msg -> Printf.printf "error -> %s\n" msg
|
2024-11-13 21:50:44 +01:00
|
|
|
;;
|
|
|
|
|
|
|
|
|
|
(* -------------------------------------------------------------------------- *)
|
|
|
|
|
(* Miller-Rabin primality test program *)
|
|
|
|
|
let program =
|
2024-11-16 15:40:00 +01:00
|
|
|
"def main with input n output result as
|
|
|
|
|
if (n % 2) == 0 then result := 1
|
|
|
|
|
else (
|
|
|
|
|
|
|
|
|
|
result := 0;
|
|
|
|
|
s := 0;
|
|
|
|
|
while (0 == ((n - 1) / (2 ^ s)) % 2) do (
|
|
|
|
|
s := s + 1
|
|
|
|
|
);
|
|
|
|
|
d := ((n - 1) / 2 ^ s);
|
|
|
|
|
for (i := 20, i > 0, i := i - 1) do (
|
|
|
|
|
a := rand(n - 4) + 2;
|
|
|
|
|
x := powmod(a, d, n);
|
2025-01-16 23:48:23 +01:00
|
|
|
y := 0;
|
2024-11-16 15:40:00 +01:00
|
|
|
for (j := 0, j < s, j := j+1) do (
|
|
|
|
|
y := powmod(x, 2, n);
|
|
|
|
|
if (y == 1 && (not x == 1) && (not x == n - 1)) then
|
|
|
|
|
result := 1;
|
|
|
|
|
else
|
|
|
|
|
skip;
|
|
|
|
|
x := y;
|
|
|
|
|
);
|
|
|
|
|
if not y == 1 then result := 1;
|
|
|
|
|
else skip;
|
|
|
|
|
)
|
|
|
|
|
)
|
|
|
|
|
"
|
2024-11-13 21:50:44 +01:00
|
|
|
;;
|
|
|
|
|
|
|
|
|
|
(* should return 0 because prime *)
|
2024-11-16 15:40:00 +01:00
|
|
|
Printf.printf "Miller-Rabin primality test program 1: ";
|
|
|
|
|
match get_result program 179424673 with
|
|
|
|
|
Ok d -> Printf.printf "%d\n" d
|
|
|
|
|
| Error `AbsentAssignment msg -> Printf.printf "error -> %s\n" msg
|
|
|
|
|
| Error `DivisionByZero msg -> Printf.printf "error -> %s\n" msg
|
2024-11-13 21:50:44 +01:00
|
|
|
;;
|
2024-11-16 15:40:00 +01:00
|
|
|
|
2024-11-13 21:50:44 +01:00
|
|
|
(* should return 1 because not prime *)
|
2024-11-16 15:40:00 +01:00
|
|
|
Printf.printf "Miller-Rabin primality test program 2: ";
|
|
|
|
|
match get_result program 179424675 with
|
|
|
|
|
Ok d -> Printf.printf "%d\n" d
|
|
|
|
|
| Error `AbsentAssignment msg -> Printf.printf "error -> %s\n" msg
|
|
|
|
|
| Error `DivisionByZero msg -> Printf.printf "error -> %s\n" msg
|
2024-11-13 21:50:44 +01:00
|
|
|
;;
|