OFFSET
0,4
FORMULA
a(n) = 2^binomial(n, 2) - A001187(n).
a(n) = n!*[x^n](g - log(g) - 1) where g = Sum_{n>=0} 2^binomial(n, 2) * x^n / n!. - Geoffrey Critzer, Nov 11 2011
MAPLE
upto := 18: g := add(2^binomial(n, 2) * x^n / n!, n = 0..upto+1):
ser := series(g - log(g) - 1, x, upto+1):
seq(n!*coeff(ser, x, n), n = 0..upto); # Peter Luschny, Feb 24 2023
MATHEMATICA
g=Sum[2^Binomial[n, 2]x^n/n!, {n, 0, 20}]; Range[0, 20]! CoefficientList[Series[g-Log[g]-1, {x, 0, 20}], x] (* Geoffrey Critzer, Nov 11 2011 *)
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Vladeta Jovovic, Apr 15 2000
EXTENSIONS
Edited and extended with a(0) = 0 by Peter Luschny, Feb 24 2023
STATUS
approved