34 lines
590 B
Plaintext
34 lines
590 B
Plaintext
fastfunc semiprim n .
|
|
d = 3
|
|
while d * d <= n
|
|
while n mod d = 0
|
|
if c = 2 : return 0
|
|
n /= d
|
|
c += 1
|
|
.
|
|
d += 2
|
|
.
|
|
if c = 1 : return n
|
|
return 0
|
|
.
|
|
print "The first 50 Blum integers:"
|
|
n = 3
|
|
numfmt 4 0
|
|
repeat
|
|
prim1 = semiprim n
|
|
if prim1 <> 0 and prim1 mod 4 = 3
|
|
prim2 = n div prim1
|
|
if prim2 <> prim1 and prim2 mod 4 = 3
|
|
c += 1
|
|
if c <= 50
|
|
write n
|
|
if c mod 10 = 0 : print ""
|
|
.
|
|
.
|
|
.
|
|
until c >= 26828
|
|
n += 2
|
|
.
|
|
print ""
|
|
print "The 26828th Blum integer is: " & n
|