login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A216639 A027642(6*n+6)/(sequence of period 2:repeat 42,210). 1
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 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,2

COMMENTS

Is a(n) always an integer?  Is there an a(n) ending with 5?

It appears (tested for n <= 800) that a(n) mod 9 is always one of {1, 2, 4, 5, 7, 8}.

There is a similar sequence of ratios A027642(10n+1)/(66*A010686(n)) which starts 1, 1, 217, 41, 1, 172081, 71, 697, 4123, 101, 23, 7055321, 131, 2059, 32767, 697, 1, 21896102683,...

a(n) is always an integer: 42 = 2*3*7 and 1, 2, and 6 divide 12n+6; 210 = 2*3*5*7 and 1, 2, 4, and 6 divide 12n+12. a(n) never ends in 5 (or 0) since 12n+6 is not divisible by 4 hence the (12n+6)-th Bernoulli denominator is not divisible by 5, and Bernoulli denominators are squarefree and hence the (12n+12)-th Bernoulli denominator, divided by 210, cannot be divisible by 5. - Charles R Greathouse IV, Sep 12 2012

The previous comments argue that 3 or 5 are never prime divisors of a(n). In addition (tested up to n <=900), 7 apparently is also a non-divisor of a(n). In summary, the prime divisors appear all to be in A140461. - Jean-Fran├žois Alcover, Sep 17 2012

LINKS

Charles R Greathouse IV, Table of n, a(n) for n = 0..10000

A. Joyal, Les nombres de Bernoulli (in French), 2003.

FORMULA

a(n) = A027642(6*n+6)/(42*A010686(n)).

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

PROG

(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

CROSSREFS

Cf. A002445, A027762, A165734, A165949.

Sequence in context: A119149 A217658 A236999 * A171098 A118845 A050717

Adjacent sequences:  A216636 A216637 A216638 * A216640 A216641 A216642

KEYWORD

nonn

AUTHOR

Paul Curtz, Sep 12 2012

EXTENSIONS

a(20)-a(40) from Charles R Greathouse IV, Sep 12 2012

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 17 21:53 EDT 2021. Contains 343071 sequences. (Running on oeis4.)