OFFSET
1,2
LINKS
Robert Israel, Table of n, a(n) for n = 1..447
MAPLE
N:= 30:
V:= Vector(N):
for d from 1 to N do
for k from 1 to floor(N/d) do
n:= k*d; V[n]:= V[n] + n!/d!*binomial(n-1, d-1);
od od:
convert(V, list); # Robert Israel, Jan 01 2020
MATHEMATICA
a[n_] := n! Sum[Binomial[n - 1, d - 1]/d!, {d, Divisors[n]}]; Table[a[n], {n, 1, 22}]
PROG
(PARI) a(n) = n! * sumdiv(n, d, binomial(n-1, d-1) / d!); \\ Michel Marcus, Nov 26 2019
(Magma) [Factorial(n)*( &+[ Binomial(n-1, d-1)/Factorial(d):d in Divisors(n)]): n in [1..22]]; // Marius A. Burtea, Jan 02 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Nov 26 2019
STATUS
approved