OFFSET
0,3
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..200
FORMULA
a(n) = n! * Sum_{k=1..n} (-1)^(k-1) * 3^(n-k)/k.
a(n) = 3 * n * a(n-1) - (-1)^n * (n-1)!.
a(n) = (2*n+1) * a(n-1) + 3 * (n-1)^2 * a(n-2).
a(n) ~ log(4/3) * 3^n * n!. - Vaclav Kotesovec, May 23 2025
MATHEMATICA
a[n_]:= n! * Sum[(-1)^(k-1)*3^(n-k)/k, {k, 1, n}]; Table[a[n], {n, 0, 18}] (* Vincenzo Librandi, May 22 2025 *)
PROG
(PARI) a(n) = n!*sum(k=1, n, (-1)^(k-1)*3^(n-k)/k);
(Magma) [0] cat [n le 1 select 1 else 3 * n * Self(n-1) - (-1)^n * Factorial(n-1): n in [1..20]]; // Vincenzo Librandi, May 22 2025
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Seiichi Manyama, May 22 2025
STATUS
approved
