A051717 1, followed by denominators of first differences of Bernoulli numbers (B(i)-B(i-1)).

1, 2, 3, 6, 30, 30, 42, 42, 30, 30, 66, 66, 2730, 2730, 6, 6, 510, 510, 798, 798, 330, 330, 138, 138, 2730, 2730, 6, 6, 870, 870, 14322, 14322, 510, 510, 6, 6, 1919190, 1919190, 6, 6, 13530, 13530, 1806, 1806, 690, 690, 282, 282, 46410, 46410, 66, 66, 1590, 1590

Offset: 0

Author

N. J. A. Sloane

Keywords

• nonn
• easy
• nice
• frac

Comments

Equivalently, denominators of Bernoulli twin numbers C(n) (cf. A051716).

The Bernoulli twin numbers C(n) are defined by C(0) = 1, then C(2n) = B(2n) + B(2n-1), C(2n+1) = -B(2n+1) - B(2n), where B() are the Bernoulli numbers A027641/A027642. The definition is due to Paul Curtz.

Denominators of column 1 of table described in A051714/A051715.

Examples

			Bernoulli numbers: 1, -1/2, 1/6, 0, -1/30, 0, 1/42, 0, -1/30, 0, 5/66, ...
First differences: -3/2, 2/3, -1/6, -1/30, 1/30, 1/42, -1/42, -1/30, ...
Numerators: -3, 2, -1, -1, 1, 1, -1, -1, 1, 5, -5, -691, 691, 7, ...
Denominators: 2, 3, 6, 30, 30, 42, 42, 30, 30, 66, 66, 2730, ...
Sequence of C(n)'s begins: 1, -1/2, -1/3, -1/6, -1/30, 1/30, 1/42, -1/42, -1/30, 1/30, 5/66, -5/66, -691/2730, 691/2730, 7/6, -7/6, ...
		

Links

• G. C. Greubel, Table of n, a(n) for n = 0..5000
• M. Kaneko, The Akiyama-Tanigawa algorithm for Bernoulli numbers, J. Integer Sequences, 3 (2000), #00.2.9.

Crossrefs

Cf. A027641, A027642, A051714, A051715, A051716, A129825, A129826.

Cf. A129724.

For numerators see A172083.

Programs

Magma

f:= func< n | Bernoulli(n) + Bernoulli(n-1) >;
function A051717(n)
  if n eq 0 then return 1;
  elif (n mod 2) eq 0 then return Denominator(f(n));
  else return Denominator(-f(n));
  end if;
end function;
[A051717(n): n in [0..50]]; // G. C. Greubel, Apr 22 2023

Maple

C:=proc(n) if n=0 then RETURN(1); fi; if n mod 2 = 0 then RETURN(bernoulli(n)+bernoulli(n-1)); else RETURN(-bernoulli(n)-bernoulli(n-1)); fi; end;

Mathematica

c[0]= 1; c[n_?EvenQ]:= BernoulliB[n] + BernoulliB[n-1]; c[n_?OddQ]:= -BernoulliB[n] - BernoulliB[n-1]; Table[Denominator[c[n]], {n,0,53}] (* Jean-François Alcover, Dec 19 2011 *)
Join[{1},Denominator[Total/@Partition[BernoulliB[Range[0,60]],2,1]]] (* Harvey P. Dale, Mar 09 2013 *)
Join[{1},Denominator[Differences[BernoulliB[Range[0,60]]]]] (* Harvey P. Dale, Jun 28 2021 *)

PARI

a(n)=if(n<3,n+1,denominator(bernfrac(n)-bernfrac(n-1))) \\ Charles R Greathouse IV, May 18 2015

SageMath

def f(n): return bernoulli(n)+bernoulli(n-1)
def A051717(n):
    if (n==0): return 1
    elif (n%2==0): return denominator(f(n))
    else: return denominator(-f(n))
[A051717(n) for n in range(51)] # G. C. Greubel, Apr 22 2023

Extensions

More terms from James Sellers, Dec 08 1999

Edited by N. J. A. Sloane, May 25 2008

Entry revised by N. J. A. Sloane, Apr 22 2021