OFFSET
1,2
COMMENTS
How is this related to A085118? - R. J. Mathar, Jul 02 2025
MAPLE
for n from 1 to 180 do modp(2^(n-1), n) ; if % = numtheory[bigomega](n) then printf("%d, ", n) ; end if; end do: # R. J. Mathar, Dec 07 2010
MATHEMATICA
Select[Range[254], PrimeOmega[#] == PowerMod[2, # - 1, #] &] (* Michael De Vlieger, Jul 02 2025 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Juri-Stepan Gerasimov, Dec 07 2010
STATUS
approved
