OFFSET
1,1
COMMENTS
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..16384
FORMULA
a(2*n) = 2*A036283(n).
From Antti Karttunen, Dec 03 2018: (Start)
a(n) = Product_{d|n} [(1+d)^(1+A286561(n,1+d))]^A010051(1+d) - after Peter J. Cameron's Mar 25 2002 comment in A006863.
(End)
MAPLE
MATHEMATICA
max = 52; s = Expand[Normal[Series[(-Log[1-x]/x)^z, {x, 0, max}]]]; a[n_, k_] := Denominator[Coefficient[s, x^n*z^k]]; A053657 = Prepend[LCM @@@ Table[a[n, k], {n, max}, {k, n}], 1]; a[n_] := A053657[[n+1]]/A053657[[n]]; Table[a[n], {n, 1, max}] (* Jean-François Alcover, Dec 20 2012 *)
PROG
(PARI) A185633(n) = if(n%2, 2, denominator(bernfrac(n)/(n))); \\ Antti Karttunen, Dec 03 2018
(PARI) A185633(n) = { my(m=1); fordiv(n, d, if(isprime(1+d), m *= (1+d)^(1+valuation(n, 1+d)))); (m); }; \\ Antti Karttunen, Dec 03 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul Curtz, Dec 18 2012
EXTENSIONS
Name edited by Antti Karttunen, Dec 03 2018
STATUS
approved