OFFSET
1,2
FORMULA
All primes, all power of 2 are in the sequence. The composites are in A035163.
PROG
(PARI) a000364(n)=subst(bernpol(2*n+1), 'x, 1/4)*4^(2*n+1)*(-1)^(n+1)/(2*n+1);
lista(nn) = {for (n=1, nn, 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 Weston Markham (weston.markham(AT)gmail.com), May 22 2005
STATUS
approved