A216639 A027642(6*n+6)/(sequence of period 2:repeat 42,210).
1, 13, 19, 13, 341, 9139, 43, 221, 19, 270413, 1541, 667147, 79, 16211, 6479, 21437, 103, 996151, 1, 11086933, 103759, 20033, 6533, 11341499, 51491, 8545667, 3097, 16211, 59, 34408161359, 1, 4137341, 5826521, 1339, 219666403, 72719023, 223, 2977, 1501, 45423164501, 83
Offset: 0
Keywords
Links
- Charles R Greathouse IV, Table of n, a(n) for n = 0..10000
- A. Joyal, Les nombres de Bernoulli (in French), 2003.
Programs
-
Maple
A010686 := proc(n) op((n mod 2)+1,[1,5]) ; end proc: A216639 := proc(n) A027642(6*n+6)/42/A010686(n) ; end proc: # R. J. Mathar, Sep 21 2012
-
PARI
a(n)=denominator(bernfrac(6*n+6))/if(n%2,210,42) \\ Charles R Greathouse IV, Sep 12 2012
-
PARI
a(n)=my(t=1);fordiv(3*n+3,d,if(isprime(2*d+1),t*=2*d+1));t/if(n%2,105,21) \\ Charles R Greathouse IV, Sep 12 2012
Extensions
a(20)-a(40) from Charles R Greathouse IV, Sep 12 2012
Comments