%I #10 Dec 05 2023 09:01:23
%S 1,4,35,462,8136,179112,4731786,145838844,5137045848,203566459392,
%T 8963064065088,434109674396736,22936702911358608,1312878755037640320,
%U 80928769156102447920,5344960170283958863008,376543135663291116638208,28184733661095459402610176
%N Expansion of e.g.f. 1/(1 - x + 3*log(1 - x)).
%F a(0) = 1; a(n) = n * a(n-1) + 3 * Sum_{k=1..n} (k-1)! * binomial(n,k) * a(n-k).
%o (PARI) a_vector(n) = my(v=vector(n+1)); v[1]=1; for(i=1, n, v[i+1]=i*v[i]+3*sum(j=1, i, (j-1)!*binomial(i, j)*v[i-j+1])); v;
%Y Cf. A052820, A367922.
%Y Cf. A367846.
%K nonn
%O 0,2
%A _Seiichi Manyama_, Dec 05 2023
|