41 lines
648 B
Plaintext
41 lines
648 B
Plaintext
global Prime1
|
|
n = 3
|
|
c = 0
|
|
|
|
print "The first 50 Blum integers:"
|
|
while True
|
|
if isSemiprime(n) then
|
|
if Prime1 % 4 = 3 then
|
|
Prime2 = n / Prime1
|
|
if (Prime2 <> Prime1) and (Prime2 % 4 = 3) then
|
|
c += 1
|
|
if c <= 50 then
|
|
print rjust(string(n), 4);
|
|
if c % 10 = 0 then print
|
|
end if
|
|
if c >= 26828 then
|
|
print : print "The 26828th Blum integer is: "; n
|
|
exit while
|
|
end if
|
|
end if
|
|
end if
|
|
end if
|
|
n += 2
|
|
end while
|
|
end
|
|
|
|
function isSemiprime(n)
|
|
d = 3
|
|
c = 0
|
|
while d*d <= n
|
|
while n % d = 0
|
|
if c = 2 then return false
|
|
n /= d
|
|
c += 1
|
|
end while
|
|
d += 2
|
|
end while
|
|
Prime1 = n
|
|
return c = 1
|
|
end function
|