A090799 Denominator(Bernoulli(n-1) + 1/n)=66, where n runs through the primes.

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

Offset

1,1

Comments

60 divides a(i+1)-a(i)

Links

Table of n, a(n) for n=1..43.

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

Sequence in context: A115395 A142647 A152951 * A044194 A044575 A039543

Adjacent sequences: A090796 A090797 A090798 * A090800 A090801 A090802

Keyword

nonn

Author

mohammed bouayoun (bouyao(AT)wanadoo.fr), Feb 10 2004

Extensions

More terms from Robert G. Wilson v, Feb 10 2004

Status

approved