OFFSET
1,3
MAPLE
f:= proc(n) local d; mul(abs(combinat:-stirling1(n, d)), d=numtheory:-divisors(n)) end proc:
map(f, [$1..30]); # Robert Israel, Jul 15 2019
MATHEMATICA
Table[Product[Abs[StirlingS1[n, d]], {d, Divisors[n]}], {n, 1, 16}]
PROG
(PARI) a(n) = my(d=divisors(n)); prod(k=1, #d, abs(stirling(n, d[k], 1))); \\ Michel Marcus, Jul 16 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Jul 15 2019
STATUS
approved