OFFSET
0,3
LINKS
Paolo Xausa, Table of n, a(n) for n = 0..1000
FORMULA
a(n) = Product_{primes p <= n} (1 + (n*(n+1)/2 - 1 - Sum_{i=2..n} s_p(i))/(p-1)), where s_p(i) is the sum of digits of i in base p. - Amiram Eldar, Jan 15 2025
MATHEMATICA
DivisorSigma[0, FoldList[Times, Range[0, 30]!]] (* Paolo Xausa, Jan 06 2025 *)
PROG
(PARI) a(n)=numdiv(prod(k=1, n, k!));
vector(30, n, a(n-1)) \\ Joerg Arndt, Jan 02 2025
(PARI) a(n) = {my(prd = 1); forprime(p = 2, n, prd *= (1 + (n*(n+1)/2 - 1 - sum(i = 2, n, sumdigits(i, p)))/(p-1))); prd; } \\ Amiram Eldar, Jan 15 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Tsuyoshi Hanatate, Dec 31 2024
STATUS
approved
