%I #12 Jul 12 2017 03:18:22
%S 1,6,4,21,25,10,56,90,65,20,126,245,240,135,35,252,560,665,510,245,56,
%T 462,1134,1540,1435,945,406,84,792,2100,3150,3360,2695,1596,630,120,
%U 1287,3630,5880,6930,6370,4606,2520,930,165,2002,5940,10230
%N Third accumulation array of the polygonal number array (A086270), by antidiagonals.
%C The chain of accumulation arrays is A144257->A086270->A185732->A185733->A184734.
%H G. C. Greubel, <a href="/A185734/b185734.txt">Table of n, a(n) for the first 50 rows, flattened</a>
%F T(n,k) = C(k+3,4)*C(n+2,3)*(k*n-n+3*k+17)/20, k>=1, n>=1.
%F T(n,k) = Sum_{j=1..n} Sum_{l=1..k} A185733(j,l), by definition.
%e Northwest corner:
%e 1....6......21.....56....126
%e 4....25.....90....245....560
%e 10...65....240....665...1540
%e 20...135...510...1435...3360
%t f[n_,k_]:= k*(1+k)*(2+k)*n*(1+n)*(10+2*k-n+k*n)/144;
%t TableForm[Table[f[n,k],{n,1,10},{k,1,15}]] (* array A185733 *)
%t Table[f[n-k+1,k],{n,14},{k,n,1,-1}]//Flatten (* A185733 *)
%t s[n_,k_]:=Sum[f[i,j],{i,1,n},{j,1,k}]; (* accumulation array of {f(n,k)} *)
%t FullSimplify[s[n,k]] (* the formula for A185734 *)
%t TableForm[Table[s[n,k],{n,1,10},{k,1,15}]] (* array A185734 *)
%t Table[s[n-k+1,k],{n,14},{k,n,1,-1}]//Flatten (* A185734 *)
%Y Cf. A144112, A144257, A086270, A185732, A185733.
%K nonn,tabl,easy
%O 1,2
%A _Clark Kimberling_, Feb 02 2011