A308037 a(n) = Sum_{d|n} Stirling2(n,d).

1, 2, 2, 9, 2, 123, 2, 1830, 3027, 43038, 2, 2023728, 2, 49337473, 213142023, 2313595723, 2, 216927216877, 2, 6712023695345, 82312699558575, 366282502967439, 2, 113350450913387211, 2436684974110753, 1850568574287104493, 106563274551407600878, 231678790379913209098, 2

Offset

1,2

Links

Table of n, a(n) for n=1..29.

Formula

a(n) = 2 <=> n is prime <=> n in { A000040 }. - Alois P. Heinz, May 10 2019

Maple

a:= n-> add(Stirling2(n, d), d=numtheory[divisors](n)):
seq(a(n), n=1..30);  # Alois P. Heinz, May 10 2019

Mathematica

a[n_] := a[n] = Sum[StirlingS2[n, d], {d, Divisors[n]}]; Table[a[n], {n, 1, 29}]

Prog

(PARI) a(n) = sumdiv(n, d, stirling(n, d, 2)); \\ Michel Marcus, May 10 2019

Crossrefs

Cf. A000040, A008277, A056045, A096308.

Sequence in context: A327859 A343824 A182106 * A011403 A113554 A085846

Adjacent sequences: A308034 A308035 A308036 * A308038 A308039 A308040

Keyword

nonn

Author

Ilya Gutkovskiy, May 10 2019

Status

approved