|
| |
|
|
A067775
|
|
Primes n such that n!*B(n+3) is an integer where B(k) are the Bernoulli numbers B(1)=-1/2 B(2)=1/6 B(4)=-1/30... B(2m+1)=0 m>1.
|
|
2
| |
|
|
2, 5, 11, 17, 23, 29, 31, 41, 47, 53, 59, 61, 71, 73, 83, 89, 101, 107, 113, 131, 137, 139, 149, 151, 157, 167, 173, 179, 181, 191, 197, 199, 211, 227, 233, 239, 241, 251, 257, 263, 269, 271, 281, 283, 293, 311, 317, 331, 337, 347, 353, 359, 367, 373, 383, 389
(list; graph; refs; listen; history; internal format)
|
|
|
|
OFFSET
| 1,1
|
|
|
COMMENTS
| If n is prime n!*B(n-1) is always an integer.
Also primes p such that p + 4 is composite. [From Klaus Brockhaus (klaus-brockhaus(AT)t-online.de), Aug 14 2008]
|
|
|
LINKS
| Klaus Brockhaus, Table of n, a(n) for n=1..1026
Eric Weisstein's World of Mathematics, Bernoulli Number [From Klaus Brockhaus (klaus-brockhaus(AT)t-online.de), Aug 14 2008]
|
|
|
MATHEMATICA
| lst = {}; Do[p = Prime@ n; If[ IntegerQ[ p! BernoulliB[p + 3]], AppendTo[lst, p]], {n, 77}]; lst (* Robert G. Wilson v (rgwv(AT)rgwv.com), Aug 19 2008 *)
|
|
|
CROSSREFS
| Cf. A049591.
Sequence in context: A181580 A118754 A140559 * A140552 A138644 A164921
Adjacent sequences: A067772 A067773 A067774 * A067776 A067777 A067778
|
|
|
KEYWORD
| nonn
|
|
|
AUTHOR
| Benoit Cloitre (benoit7848c(AT)orange.fr), Feb 06 2002
|
| |
|
|
