OFFSET
1,1
COMMENTS
Partial sums of A005476.
a(n) is the dot product of the vectors of the first n positive integers and the next n integers. - Michel Marcus, Sep 02 2020
REFERENCES
A. H. Beiler, Recreations in the Theory of Numbers, Dover, N.Y., 1964, pp. 194-196.
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..1000
Index entries for linear recurrences with constant coefficients, signature (4,-6,4,-1).
FORMULA
G.f.: x*(2+3*x)/(1-x)^4.
a(n) = A132121(n,1). - Reinhard Zumkeller, Aug 12 2007
a(n) = Sum_{i=1..n} Sum_{j=1..n} i+min(i,j). - Enrique Pérez Herrero, Jan 15 2013
a(n) = Sum_{i=1..n} i*(n+i). - Charlie Marion, Apr 10 2013
Sum_{n>=1} 1/a(n) = 36 - 3*Pi*5^(3/4)*phi^(3/2)/4 - 15*sqrt(5)*log(phi)/4 - 75*log(5)/8 = 0.66131826232008423794478..., where phi = A001622 = (1 + sqrt(5))/2 is the golden ratio. - Vaclav Kotesovec, Jun 01 2018
E.g.f.: exp(x)*x*(12 + 21*x + 5*x^2)/6. - Stefano Spezia, Feb 21 2024
MAPLE
[n*(n+1)*(5*n+1)/6$n=1..40]; # Muniru A Asiru, Jan 01 2019
MATHEMATICA
Table[Range[x].Range[x+1, 2x], {x, 40}] (* or *) LinearRecurrence[{4, -6, 4, -1}, {2, 11, 32, 70}, 40] (* Harvey P. Dale, Jun 01 2018 *)
PROG
(PARI) a(n) = n*(n+1)*(5*n+1)/6;
(Magma) [n*(n+1)*(5*n+1)/6 : n in [1..40]]; // Vincenzo Librandi, Jan 01 2019
(GAP) a:=List([1..40], n->n*(n+1)*(5*n+1)/6);; Print(a); # Muniru A Asiru, Jan 01 2019
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Barry E. Williams, Dec 16 1999
EXTENSIONS
More terms from James A. Sellers, Jan 19 2000
STATUS
approved