|
| |
|
|
A090799
|
|
Denominator(Bernoulli(n-1) + 1/n)=66, where n runs through the primes.
|
|
0
|
|
|
|
71, 131, 191, 251, 311, 431, 971, 1031, 1091, 1511, 1571, 1811, 1931, 2111, 2351, 2411, 2711, 3371, 3491, 3671, 4091, 4211, 4391, 4871, 5051, 5231, 5351, 5471, 5711, 6011, 6131, 6311, 6911, 7331, 7451, 7691, 8111, 8231, 8291, 8831, 9371, 9851, 10091
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
|
OFFSET
|
1,1
|
|
|
COMMENTS
|
60 divides a(i+1)-a(i)
|
|
|
LINKS
|
|
|
|
FORMULA
|
prime(i) / denominator(Bernoulli(prime(n)-1) + 1/prime(n))=66
|
|
|
MATHEMATICA
|
Prime[ Select[ Range[ PrimePi[10192]], Denominator[ BernoulliB[ Prime[ # ] - 1] + 1/Prime[ # ]] == 66 &]] (* Robert G. Wilson v, Feb 10 2004 *)
|
|
|
PROG
|
(PARI) bouay66(n)=bernfrac(n-1)+1/n for (i=1, 1000, if(denominator(bouay42(prime(i)))==66, print1(prime(i), ", ")))
|
|
|
CROSSREFS
|
|
|
|
KEYWORD
|
nonn
|
|
|
AUTHOR
|
mohammed bouayoun (bouyao(AT)wanadoo.fr), Feb 10 2004
|
|
|
EXTENSIONS
|
|
|
|
STATUS
|
approved
|
| |
|
|
