OFFSET
0,2
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..300 (terms 0..25 from R. J. Mathar)
FORMULA
a(n) = n!*(1 + Sum_{i=0..n-1} Bernoulli(i)). - R. J. Mathar, Feb 20 2008
EXAMPLE
MAPLE
A129716 := proc(n) n!*(1+add(bernoulli(i), i=0..n-1)); end: seq(A129716(n), n=0..40) ; # R. J. Mathar, Feb 20 2008
MATHEMATICA
max = 21; Accumulate[ Table[ If[n == 0, 1, BernoulliB[n-1]], {n, 0, max}]]*Range[0, max]! (* Jean-François Alcover, Mar 04 2013 *)
PROG
(PARI) vector(26, n, (n-1)!*(1 + sum(j=0, n-2, bernfrac(j))) ) \\ G. C. Greubel, Dec 03 2019
(Magma) [1] cat [Factorial(n)*(1 + (&+[Bernoulli(k): k in [0..n-1]]) ): n in [1..25]]; // G. C. Greubel, Dec 03 2019
(Sage) [factorial(n)*(1 + sum(bernoulli(k) for k in (0..n-1)) ) for n in (0..25)] # G. C. Greubel, Dec 03 2019
(GAP) List([0..25], n-> Factorial(n)*(1 + Sum([0..n-1], j-> Bernoulli(j)) ) ); # G. C. Greubel, Dec 03 2019
CROSSREFS
KEYWORD
sign
AUTHOR
Paul Curtz, Jun 02 2007
EXTENSIONS
More terms from R. J. Mathar, Feb 20 2008
STATUS
approved