cp's OEIS Frontend

A061718 a(n) = (n*(n+1)/2)^n.

%I A061718 #37 Jun 13 2025 08:20:05
%S A061718 1,9,216,10000,759375,85766121,13492928512,2821109907456,
%T A061718 756680642578125,253295162119140625,103510234140112521216,
%U A061718 50714860157241037295616,29345269354638035222576971
%N A061718 a(n) = (n*(n+1)/2)^n.
%C A061718 a(n) is the number of n X n matrices with nonnegative integer entries such that every row sum equals 2. - Sharon Sela (sharonsela(AT)hotmail.com), May 08 2002
%C A061718 Resultant of the polynomials P(n,x) and Q(n,x) where P(n,x)=sum(k=1,n,k*(-x)^k) and Q(n,x)=x^n-1. - _Benoit Cloitre_, Jan 26 2003
%C A061718 a(n) is also the number of positive-volume, axis-aligned, n-dimensional rectangular solids that have vertices in the set {0,1,...,n}^n. Proof: If (M_1,...,M_n) is the corner with the maximum coordinate values for such a solid, then there are (M_1)*...*(M_n) possibilities for the corner with the minimum coordinate values. The sum over all possibilities for M_1, ..., M_n can be factored into the product of n sums; each of the n sums simplifies to n(n+1)/2. - _Lee A. Newberg_, Aug 31 2009
%H A061718 Harry J. Smith, <a href="/A061718/b061718.txt">Table of n, a(n) for n = 1..100</a>
%F A061718 a(n) = Sum_{i=1..n,j=1..n,k=1..n,...} (i*j*k*...). E.g., a(2) = 9 because 1*1 + 1*2 + 2*1 + 2*2 = 9. - _Ben Paul Thurston_, Aug 15 2006
%F A061718 a(n) = [x^n] 1/(1 - (n*(n + 1)/2)*x). - _Ilya Gutkovskiy_, Oct 10 2017
%F A061718 a(n) ~ 2^(-n)*exp(1-1/(2*n))*n^(2*n). - _Stefano Spezia_, Jun 13 2025
%p A061718 a:=n->mul(sum(j, j=0..n),k=1..n): seq(a(n), n=1..13); # _Zerinvary Lajos_, Jun 02 2007
%p A061718 a:=n->mul(binomial(n+2,2), k=0..n): seq(a(n), n=0..12); # _Zerinvary Lajos_, Oct 02 2007
%t A061718 Table[((n(n+1))/2)^n,{n,20}] (* _Harvey P. Dale_, Dec 09 2022 *)
%t A061718 With[{nn=15},#[[1]]^#[[2]]&/@Thread[{Accumulate[Range[nn]],Range[nn]}]] (* _Harvey P. Dale_, Dec 25 2023 *)
%o A061718 (C++) /* e.g. n = 6: */
%o A061718 int main()
%o A061718 {
%o A061718     int sum = 0;
%o A061718     for(int i = 1; i < 7; i++)
%o A061718         for(int j = 1; j < 7; j++)
%o A061718             for(int k=1;k<7; k++)
%o A061718                 for(int l = 1; l < 7; l++)
%o A061718                     for(int m = 1; m < 7; m++)
%o A061718                         for(int n = 1; n < 7; n++)
%o A061718                             sum += i*j*k*l*m*n;
%o A061718     cout << sum << endl;
%o A061718     return 0;
%o A061718 } // _Ben Paul Thurston_, Aug 15 2006
%o A061718 (PARI) { for (n=1, 100, write("b061718.txt", n, " ", (n*(n + 1)/2)^n) ) } \\ _Harry J. Smith_, Jul 26 2009
%Y A061718 Cf. A066300.
%K A061718 easy,nonn
%O A061718 1,2
%A A061718 _Jason Earls_, Jun 20 2001