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