Ormiston pairs: Difference between revisions

m
 
Line 625:
<syntaxhighlight>
maxn = 1000000
len isdivsieve[] maxn
proc sievemksieve . .
max = sqrt len isdivsieve[]
for d = 2 to max
if isdivsieve[d] = 0
for i = d * d step d to len isdivsieve[]
isdivsieve[i] = 1
.
.
.
isdivsieve[] &= 0
.
mksieve
sieve
func nextprim n .
repeat
n += 1
until isdivsieve[n] = 0
.
return n
.
#
func digs n .
while n > 0
Line 658 ⟶ 659:
db = da
a = nextprim a
until a > 1000000maxn
da = digs a
if da = db
2,083

edits