OFFSET
1,3
COMMENTS
See A264388 for the numerators and details about the Dedekind sum s(1,n), as well as references.
FORMULA
a(n) = denominator(binomial(n-1, 2)/(6*n)), n >= 1.
a(n) = denominator(s(1,n)), with s(1,n) = Sum_{r=1..(n-1)} (r/n)*(r/n - floor(r/n)- 1/2), n >= 1, where s(h,k) are the Dedekind sums.
MATHEMATICA
Denominator[Table[Binomial[n-1, 2]/(6n), {n, 50}]] (* Harvey P. Dale, Aug 30 2016 *)
PROG
(Julia)
using Nemo
A264389(n) = denominator(dedekind_sum(1, n))
[A264389(n) for n in 1:70] |> println # Peter Luschny, Mar 13 2018
CROSSREFS
KEYWORD
nonn,easy,frac
AUTHOR
Wolfdieter Lang, Jan 11 2016
STATUS
approved