[go: up one dir, main page]

login
A100454
a(n) = sum of n-th column in array in A100452.
2
1, 7, 24, 58, 109, 188, 307, 444, 641, 885, 1149, 1493, 1936, 2358, 2975, 3645, 4267, 5102, 6057, 6941, 8124, 9395, 10458, 12140, 13561, 15336, 17110, 19204, 21124, 23596, 26219, 28587, 31254, 34593, 37252, 40545, 44524, 47451, 51724, 55853, 60068, 64152, 69801, 73657, 79372
OFFSET
1,2
LINKS
MATHEMATICA
t[1, n_]:= n^2; (* t = A100452 *)
t[m_, n_]/; 1<m<=n:= t[m, n]= (n-m+1)*Floor[(t[m-1, n] -1)/(n-m+1)];
t[_, _]=0;
A100454[n_]:= A100454[n]= Sum[t[n-k+1, n], {k, n}];
Table[A100454[n], {n, 60}] (* G. C. Greubel, Apr 07 2023 *)
PROG
(Magma)
function t(n, k) // t = A100452
if k eq 1 then return n^2;
else return (n-k+1)*Floor((t(n, k-1) -1)/(n-k+1));
end if;
end function;
A100454:= func< n | (&+[t(n, n-k+1): k in [1..n]]) >;
[A100454(n): n in [1..60]]; // G. C. Greubel, Apr 07 2023
(SageMath)
def t(n, k): # t = A100452
if (k==1): return n^2
else: return (n-k+1)*((t(n, k-1) -1)//(n-k+1))
def A100454(n): return sum(t(n, n-k+1) for k in range(1, n+1))
[A100454(n) for n in range(1, 61)] # G. C. Greubel, Apr 07 2023
CROSSREFS
Cf. A100452.
Sequence in context: A146298 A079671 A212511 * A081436 A024205 A008779
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Nov 22 2004
EXTENSIONS
Terms a(26) onward added by G. C. Greubel, Apr 07 2023
STATUS
approved