%I #18 Jan 07 2025 01:59:51
%S 61,1385,12284,68060,281210,948002,2749340,7097948,16700255,36419955,
%T 74551048,144631240,267951892,476948260,819683560,1365672424,
%U 2213323585,3499318141,5410278500,8197124100
%N Fourth column of triangle A060058.
%H Indranil Ghosh, <a href="/A060061/b060061.txt">Table of n, a(n) for n = 0..5000</a>
%F a(n) = Sum_{j3=1..n+1} j3^2*Sum_{j2=1..j3+1} j2^2*Sum_{j1=1..j2+1} j1^2.
%F a(n) = A060058(n+3, 3) = binomial(n+6, 6)*(280*n^3+2436*n^2+5906*n+3843)/(7*9).
%F G.f.: (61+775*x+1179*x^2+225*x^3)/(1-x)^10 = p(3, x)/(1-x)^(3*3+1) with p(3, x)=sum(A060063(3, m)*x^m, m=0..3).
%t Table[Binomial[n+6,6]*(280*n^3+2436*n^2+5906n+3843)/63,{n,0,19}] (* _Indranil Ghosh_, Feb 21 2017 *)
%o (Python)
%o import math
%o def C(n, r):
%o f=math.factorial
%o return f(n)//f(r)//f(n-r)
%o def A060061(n):
%o return (C(n+6, 6)*(280*n**3+2436*n**2+5906*n+3843))//63 # _Indranil Ghosh_, Feb 21 2017
%Y Cf. A000330, A060060.
%K nonn,easy
%O 0,1
%A _Wolfdieter Lang_, Mar 16 2001