(PARI) x='x+O('x^30); Vec(serlaplace(1/(4-sum(k=1, 3, exp(k*x))))) \\ G. C. Greubel, Oct 09 2018
(Magma) m:=30; R<x>:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!(1/(4-Exp(x)-Exp(2*x)-Exp(3*x)))); [Factorial(n-1)*b[n]: n in [1..m]]; // G. C. Greubel, Oct 09 2018