OFFSET
0,1
COMMENTS
Difference table of B(n+2):
1/6, 0, -1/30, 0, 1/42, 0, -1/30, ...
-1/6, -1/30, 1/30, 1/42, -1/42, -1/30, ...
2/15, 1/15, -1/105, -1/21, -1/105, ...
-1/15, -8/105, -4/105, 4/105, ...
-1/105, 4/105, 8/105, ...
1/21, 4/105, ...
-1/105, ...
...
a(n) is the denominator of the n-th term of the first column.
a(n+2) is the denominator of the n-th term of the third row.
See A239315(n), which is the table without the first two rows.
Inverse binomial transform: 1/6, -1/6, 2/15, -1/15, -1/105, 1/21, -1/105, -1/15, 7/165, 5/33, -2663/15015, ... .
FORMULA
MATHEMATICA
max = 42; bb = Table[BernoulliB[n+2], {n, 0, max}]; dd = Table[Differences[bb, n], {n, 0, max}]; dd[[All, 1]] // Denominator (* Jean-François Alcover, Apr 09 2015 *)
PROG
(PARI) lista(nn) = {A = vector(nn, n, bernfrac(n+1)); for (i=1, #A-1, for(j=0, i-1, A[i+1]-=binomial(i, j)*A[j+1])); for (i=1, #A, print1(denominator(A[i]), ", ")); } \\ Michel Marcus, Apr 08 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul Curtz, Apr 07 2015
STATUS
approved