Wieferich primes: Difference between revisions
added RPL
(Added Easylang) |
(added RPL) |
||
(3 intermediate revisions by 3 users not shown) | |||
Line 559:
return 1
.
func
if isprim p = 0
return 0
Line 575:
print "Wieferich primes less than 5000: "
for i = 2 to 5000
if
print i
.
Line 996:
=={{header|Python}}==
<syntaxhighlight lang="python">
# Wieferich-Primzahlen
MAX: int = 5_000
# Berechnet a^n mod m
def pow_mod(a: int, n: int, m: int) -> int:
assert n >= 0 and m != 0, "pow_mod(a, n, m), n >= 0, m <> 0"
n -= 1
else:
a = (a*a)%m
n //= 2
return res%m
def is_prime(n: int) -> bool:
for i in range(2, int(n**0.5) + 1):
if n
return False
return True
def
if
return False
▲ p2 = p ** 2
▲ while p > 1:
▲ q = (2 * q) % p2
▲ p -= 1
▲ if q == 1:
return True
else:
return False
if __name__ == '__main__':
print(f"Wieferich primes less than
for i in range(2,
if
print(i)
</syntaxhighlight>
{{out}}
<pre>Wieferich primes less than 5000:
Line 1,120 ⟶ 1,131:
Found 2 Wieferich primes that are < 5,000
</pre>
=={{header|RPL}}==
{{works with|RPL|HP49-C}}
« { } 2
'''WHILE''' DUP 5000 < '''REPEAT'''
'''IF''' 2 OVER 1 - ^ 1 - OVER SQ MOD NOT '''THEN''' SWAP OVER + SWAP '''END'''
NEXTPRIME
'''END''' DROP
» '<span style="color:blue">TASK</span>' STO
{{out}}
<pre>
1: { 1093 3511 }
</pre>
Line 1,248 ⟶ 1,272:
{{libheader|Wren-math}}
{{libheader|Wren-big}}
<syntaxhighlight lang="
import "./big" for BigInt
var primes = Int.primeSieve(5000)
|