|
composite_small := proc (n::integer)
description "procedure to determine if n has a prime factor less than 100";
if igcd(2305567963945518424753102147331756070, n) = 1 then
return false else return true end if;
end proc;
# begin initialization section
p := [0, 2, 6, 8, 12, 18, 20, 26, 30, 32]; o := [1271, 1691]; m := 2310;
# end initialization section
with(ArrayTools); os := Size(o, 2); ps := Size(p, 2);
loopstop := 10^11; loopstart := 0;
print(11);
for n from loopstart to loopstop do
for a to os do
counter := 0; wc := 0; wd := 0;
while `and`(wd > -10, wd < ps) do
wd := wd+1;
if composite_small(m*n+o[a]+p[wd]) = false then wd := wd+1 else wd := -10 end if;
end do;
if wd >= 9 then while `and`(counter >= 0, wc < ps) do
wc := wc+1;
if isprime(m*n+o[a]+p[wc]) then counter := counter+1 else counter := -1 end if end do;
end if;
if counter = ps then print(m*n+o[a]); end if;
end do;
end do;
# Matt C. Anderson, Apr 30 2015
| 