A090799 - OEIS

%I #6 Mar 15 2015 21:02:28

%S 71,131,191,251,311,431,971,1031,1091,1511,1571,1811,1931,2111,2351,

%T 2411,2711,3371,3491,3671,4091,4211,4391,4871,5051,5231,5351,5471,

%U 5711,6011,6131,6311,6911,7331,7451,7691,8111,8231,8291,8831,9371,9851,10091

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

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

%F prime(i) / denominator(Bernoulli(prime(n)-1) + 1/prime(n))=66

%t Prime[ Select[ Range[ PrimePi[10192]], Denominator[ BernoulliB[ Prime[ # ] - 1] + 1/Prime[ # ]] == 66 &]] (* _Robert G. Wilson v_, Feb 10 2004 *)

%o (PARI) bouay66(n)=bernfrac(n-1)+1/n for (i=1,1000,if(denominator(bouay42(prime(i)))==66,print1(prime(i), ",")))

%K nonn

%O 1,1

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

%E More terms from _Robert G. Wilson v_, Feb 10 2004