OFFSET
1,4
LINKS
Andrew Howroyd, Table of n, a(n) for n = 1..200
FORMULA
Inverse Moebius transform of A326888.
PROG
(PARI) a(n)={sum(k=1, n, -k!*(stirling((n+1)\2, k, 2) + stirling(n\2+1, k, 2))/4 + k!*sumdiv(n, d, eulerphi(d)*stirling(n/d, k, 2))/(2*n))} \\ Andrew Howroyd, Sep 13 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Andrew Howroyd, Sep 13 2019
STATUS
approved