A387027 a(n) = lcm({1, 2, ..., n}) * (n + 1) / n for n > 0, a(0) = 1.
1, 2, 3, 8, 15, 72, 70, 480, 945, 2800, 2772, 30240, 30030, 388080, 386100, 384384, 765765, 12972960, 12932920, 245044800, 244432188, 243877920, 243374040, 5587021440, 5577321750, 27841990176, 27800803800, 83288004800, 83181770100, 2409402996000, 2406725881560
Offset: 0
Keywords
Links
- Paolo Xausa, Table of n, a(n) for n = 0..2000
Crossrefs
Cf. A003418.
Programs
-
Maple
A387027 := n -> local k; ifelse(n = 0, 1, ((n+1) * lcm(seq(k, k = 1..n))) / n): seq(A387027(n), n = 0..30);
-
Mathematica
A387027[n_] := If[n == 0, 1, Quotient[(n + 1) LCM @@ Range[1, n], n]]; Table[A387027[n], {n, 0, 30}]
-
Python
from math import lcm def A387027(n): return (n+1)*lcm(*range(1,n+1))//n if n else 1 # Chai Wah Wu, Aug 17 2025
Formula
a(n) = A003418(n) * (n + 1) / n for n >= 1.
a(n) = (n+1)*A002944(n). - R. J. Mathar, Aug 19 2025
Comments