This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A356661 #13 Aug 21 2022 09:26:47 %S A356661 1,4,12,60,240,1860,10080,95760,766080,8210160,79833600,1100484000, %T A356661 12454041600,188172784800,2683799838720,44951306400000, %U A356661 711374856192000,13745322470880000,243290200817664000,5142812718440517120,103294640229580800000,2351280996859354560000 %N A356661 a(n) = n! * Sum_{d|n} 1/d^(n/d - 1). %F A356661 a(p) = 2 * p! for prime p. %F A356661 E.g.f.: Sum_{k>=1} x^k/(1 - x^k/k). %t A356661 a[n_] := n! * DivisorSum[n, 1/#^(n/# - 1) &]; Array[a, 22] (* _Amiram Eldar_, Aug 21 2022 *) %o A356661 (PARI) a(n) = n!*sumdiv(n, d, 1/d^(n/d-1)); %o A356661 (PARI) my(N=30, x='x+O('x^N)); Vec(serlaplace(sum(k=1, N, x^k/(1-x^k/k)))) %Y A356661 Cf. A087905, A087909, A098558, A356662. %K A356661 nonn %O A356661 1,2 %A A356661 _Seiichi Manyama_, Aug 21 2022