|
%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
|
