OFFSET
1,2
MATHEMATICA
Join[{1}, Select[Range[1, 4140001, 2], PowerMod[2, #, #]==EulerPhi[#]&]] (* The program generates the first 6 terms of the sequence. *) (* Harvey P. Dale, Feb 07 2025 *)
PROG
(PARI) is(n)=Mod(2, n)^n==eulerphi(n)&&n%2==1;
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Altug Alkan, Dec 11 2016
STATUS
approved
