A162173 Absolute difference of numerator and denominator of Bernoulli(2*n), n >= 0.

0, 5, 31, 41, 31, 61, 3421, 1, 4127, 43069, 174941, 854375, 236366821, 8553097, 23749461899, 8615841261683, 7709321041727, 2577687858361, 26315271553055396563, 2929993913841553, 261082718496449135581

Offset

0,2

Links

Muniru A Asiru, Table of n, a(n) for n = 0..310

Formula

a(n) = abs(A000367(n) - A002445(n)).

Example

0(=1-1), 5(=abs(1-6)), 31(=abs(-1-30)), 41(=abs(1-42)), 31(=abs(-1-30)), 61(=abs(5-66), 3421(=abs(-691-2730)), 1(=7-6), 4127(=abs(-3617-510)), etc.

Maple

A000367 := proc(n) numer(bernoulli(2*n)) ; end: A002445 := proc(n) denom(bernoulli(2*n)) ; end: A162173 := proc(n) abs( A000367(n)-A002445(n)) ; end: seq(A162173(n), n=0..40) ; # R. J. Mathar, Sep 13 2009

Mathematica

a[n_] := Abs[Numerator[BernoulliB[2n]] - Denominator[BernoulliB[2n]]]; Array[a, 20, 0] (* Amiram Eldar, Nov 25 2018 *)
Abs[Denominator[#]-Numerator[#]]&/@BernoulliB[2 Range[0, 20]] (* Harvey P. Dale, Jun 18 2022 *)

Prog

(PARI) a(n) = my(b=bernfrac(2*n)); abs(numerator(b)-denominator(b)); \\ Michel Marcus, Nov 25 2018
(GAP) List([0..30], n->AbsInt(NumeratorRat(Bernoulli(2*n))-DenominatorRat(Bernoulli(2*n))));  # Muniru A Asiru, Nov 25 2018

Crossrefs

Cf. A000367, A002445.

Sequence in context: A172030 A042837 A354881 * A027761 A166306 A341632

Adjacent sequences: A162170 A162171 A162172 * A162174 A162175 A162176

Keyword

nonn

Author

Juri-Stepan Gerasimov, Jun 27 2009

Extensions

a(13) corrected by Paul Curtz, Sep 07 2009

a(15) and a(16) corrected by R. J. Mathar, Sep 13 2009

Offset 0 from Amiram Eldar, Nov 25 2018

Status

approved