OFFSET
1,3
COMMENTS
See A368366 for further information.
LINKS
Paolo Xausa, Table of n, a(n) for n = 1..30
MATHEMATICA
A368367[n_] := With[{f = Range[n]!}, Total[f]^n - n^n*Apply[Times, f]];
Array[A368367, 10] (* Paolo Xausa, Jan 29 2024 *)
PROG
(Python)
from math import prod, factorial
def A368367(n): return sum(factorial(i) for i in range(1, n+1))**n-n**n*prod(i**(n-i+1) for i in range(2, n+1)) # Chai Wah Wu, Jan 25 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Jan 24 2024
STATUS
approved