# Maple code from Robert Israel, April 16 2017 # g:= proc(t) local F; # if isprime(t) then return true fi; # F:= ifactors(t,easy)[2]; # if nops(F) <> 1 then return false fi; # if F[1][1]::integer then return true fi; # nops(ifactors(t)) = 1; # end proc: # N:= 10^75: # Res:= {}: # for i from 1 to ilog2(N) do # for j from 1 to floor(log[3](N/2^i)) do # q:= 2^i*3^j; # if isprime(q-1) and g((q-2)/2^padic:-ordp(q-2,2)) then Res:= Res union {q-1} fi; # if isprime(q+1) and g((q+2)/2^padic:-ordp(q+2,2)) then Res:= Res union {q+1} fi # od od: # Res:= sort(convert(Res,list)): # bfile:= "b284037.txt": # for n from 1 to nops(Res) do fprintf(bfile,"%d %d\n",n,Res[n]) od: # fclose(bfile); 1 11 2 13 3 19 4 23 5 37 6 47 7 53 8 73 9 97 10 107 11 163 12 193 13 383 14 487 15 577 16 863 17 1153 18 2593 19 2917 20 4373 21 8747 22 995327 23 1492993 24 1990657 25 5308417 26 28311553 27 86093443 28 6879707137 29 1761205026817 30 2348273369087 31 5566277615617 32 7421703487487 33 21422803359743 34 79164837199873 35 3470494144278527 36 3799912185593857 37 115422332637413377 38 1332669751402954753 39 4803028329503971873 40 161919374795459002367 41 504202701918008951235073 42 1838129271989302091317247 43 2168345519443636233418208968703 44 213879912621499745085421535625217 45 225322047864625245851308284444673 46 28070062609828769223367060340342783 47 4547350142792260614185463775135531007 48 143556623544770914273611162510277214207 49 65249504295362867872710100461288646829408257 50 120560000905886971974292843165430886709592063 51 1637802228039295265610416163917632858514482560563146653697