%I #7 Apr 08 2023 02:28:47
%S 1,7,24,58,109,188,307,444,641,885,1149,1493,1936,2358,2975,3645,4267,
%T 5102,6057,6941,8124,9395,10458,12140,13561,15336,17110,19204,21124,
%U 23596,26219,28587,31254,34593,37252,40545,44524,47451,51724,55853,60068,64152,69801,73657,79372
%N a(n) = sum of n-th column in array in A100452.
%H G. C. Greubel, <a href="/A100454/b100454.txt">Table of n, a(n) for n = 1..1000</a>
%t t[1, n_]:= n^2; (* t = A100452 *)
%t t[m_, n_]/; 1<m<=n:= t[m,n]= (n-m+1)*Floor[(t[m-1,n] -1)/(n-m+1)];
%t t[_, _]=0;
%t A100454[n_]:= A100454[n]= Sum[t[n-k+1,n], {k,n}];
%t Table[A100454[n], {n, 60}] (* _G. C. Greubel_, Apr 07 2023 *)
%o (Magma)
%o function t(n, k) // t = A100452
%o if k eq 1 then return n^2;
%o else return (n-k+1)*Floor((t(n, k-1) -1)/(n-k+1));
%o end if;
%o end function;
%o A100454:= func< n | (&+[t(n,n-k+1): k in [1..n]]) >;
%o [A100454(n): n in [1..60]]; // _G. C. Greubel_, Apr 07 2023
%o (SageMath)
%o def t(n, k): # t = A100452
%o if (k==1): return n^2
%o else: return (n-k+1)*((t(n, k-1) -1)//(n-k+1))
%o def A100454(n): return sum(t(n,n-k+1) for k in range(1,n+1))
%o [A100454(n) for n in range(1,61)] # _G. C. Greubel_, Apr 07 2023
%Y Cf. A100452.
%K nonn
%O 1,2
%A _N. J. A. Sloane_, Nov 22 2004
%E Terms a(26) onward added by _G. C. Greubel_, Apr 07 2023
|