OFFSET
0,3
LINKS
Robert P. P. McKone, Table of n, a(n) for n = 0..199
FORMULA
Let b(n) = 2 if n = 2^m or n = 2^m + 1 for some m, otherwise 1; then a(n) = lcm({1,2,...,n}) / b(n). - Peter Luschny, Jun 26 2009
MAPLE
LCM := proc(n) option remember; if n < 2 then 1 else ilcm(n, LCM(n-1)) fi end;
a := proc(n) local i; add(i, i=convert(2*iquo(n+2, 2), base, 2)); `if`(%=1, LCM(n), LCM(n)/2) end: # Peter Luschny, Jun 26 2009
MATHEMATICA
l[1] = 1; l[n_Integer?NonNegative] := l[n] = LCM[n, l[n - 1]];
a[0] = 1; a[n_Integer?NonNegative] := a[n] = Module[{s}, s = Total[IntegerDigits[2*Quotient[n + 2, 2], 2]]; If[s == 1, l[n], l[n]/2]];
Table[a[n], {n, 0, 29}] (* Robert P. P. McKone, Oct 31 2023 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
Offset corrected by Peter Luschny, Jun 26 2009
STATUS
approved