OFFSET
0,4
COMMENTS
a(n) is the sum of all products of three distinct cubes of positive integers up to n, i.e., the sum of all products of three distinct elements from the set of cubes {1^3, ..., n^3}.
LINKS
Roudy El Haddad, Multiple Sums and Partition Identities, arXiv:2102.00821 [math.CO], 2021.
Roudy El Haddad, A generalization of multiple zeta value. Part 2: Multiple sums. Notes on Number Theory and Discrete Mathematics, 28(2) 2022, 200-233, DOI: 10.7546/nntdm.2022.28.2.200-233. See Theorem 5.1 for m = 3 and p = 3.
Index entries for linear recurrences with constant coefficients, signature (13,-78,286,-715,1287,-1716,1716,-1287,715,-286,78,-13,1). [Typo corrected by Georg Fischer, Sep 30 2022]
FORMULA
a(n) = Sum_{k=3..n} Sum_{j=2..k-1} Sum_{i=1..j-1} k^3*j^3*i^3.
a(n) = n^2 * (n + 1)^2 * (n - 1) * (n - 2) * (35*n^6 + 5*n^5 - 237*n^4 - 77*n^3 + 502*n^2 + 148*n -336)/13440.
a(n) = binomial(n+1,4)*binomial(n+1,2)*(35*n^6 + 5*n^5 - 237*n^4 - 77*n^3 + 502*n^2 + 148*n -336)/280.
PROG
(PARI) {a(n) = n^2 * (n + 1)^2 * (n - 1) * (n - 2) * (35*n^6 + 5*n^5 - 237*n^4 - 77*n^3 + 502*n^2 + 148*n -336)/13440};
(Python)
def A352980(n): return n**2*(n*(n*(n*(n*(n*(n*(n*(n*(n*(35*n - 30) - 347) + 180) + 1365) - 350) - 2541) + 240) + 2160) - 40) - 672)//13440 # Chai Wah Wu, May 15 2022
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Roudy El Haddad, Apr 13 2022
STATUS
approved