OFFSET
1,2
COMMENTS
Also, numbers n such that B(n)*n == 6 (mod n), where B(n) is the n-th Bernoulli number. Equivalently, SUM[prime p, (p-1) divides n] n/p == -6 (mod n). There are no other terms below 10^15. - Max Alekseyev, Aug 26 2013
MATHEMATICA
Select[Range[10000], Mod[Sum[PowerMod[i, #, #], {i, #}], #] == 6 &]
PROG
(PARI) is(n)=Mod(sumdiv(n, d, if(isprime(d+1), n/(d+1))), n)==-6 \\ Charles R Greathouse IV, Nov 13 2013
CROSSREFS
KEYWORD
nonn,more
AUTHOR
José María Grau Ribas, Jun 24 2013
EXTENSIONS
1,3,4 prepended by Max Alekseyev, Aug 26 2013
STATUS
approved