OFFSET
1,1
MATHEMATICA
Select[Range[1000], CompositeQ[#] && #/2^IntegerExponent[#, 2] > 1 && Divisible[Abs[EulerE[2*#]] - 1, #] &] (* Amiram Eldar, Nov 26 2020 *)
PROG
(PARI) a000364(n)=subst(bernpol(2*n+1), 'x, 1/4)*4^(2*n+1)*(-1)^(n+1)/(2*n+1);
lista(nn) = {forcomposite(n=1, nn, if ( n != 2^valuation(n, 2), if (Mod(a000364(n), n) == 1, print1(n, ", ")); ); ); } \\ Michel Marcus, Apr 18 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
Benoit Cloitre, Apr 06 2003
EXTENSIONS
More terms from Hans Havermann, Apr 07 2003
a(23)-a(43) from Amiram Eldar, Nov 26 2020
STATUS
approved