Miller–Rabin primality test: Difference between revisions

Content added Content deleted
Line 1,215: Line 1,215:
miller_rabin(n) := block(
miller_rabin(n) := block(
[v: [2, 3, 5, 7, 11, 13, 17], s, d, q: true],
[v: [2, 3, 5, 7, 11, 13, 17], s, d, q: true, a],
if n < 19 then member(n, v) else (
if n < 19 then member(n, v) else (
[s, d]: facpow(n - 1, 2),
[s, d]: facpow(n - 1, 2),