login
A331617
E.g.f.: exp(1 / (1 - arctan(x)) - 1).
4
1, 1, 3, 11, 49, 265, 1683, 12035, 95169, 832337, 7998467, 83033403, 922112305, 10978263257, 139956480467, 1889161216179, 26798589518593, 401123509624737, 6346168059440515, 105040097140558699, 1805102151607613361, 32421358229074354601
OFFSET
0,3
COMMENTS
a(53) is negative. - Vaclav Kotesovec, Jan 26 2020
LINKS
FORMULA
a(0) = 1; a(n) = Sum_{k=1..n} binomial(n-1,k-1) * A191700(k) * a(n-k).
MATHEMATICA
nmax = 21; CoefficientList[Series[Exp[1/(1 - ArcTan[x]) - 1], {x, 0, nmax}], x] Range[0, nmax]!
A191700[0] = 1; A191700[n_] := A191700[n] = Sum[Binomial[n, k] If[OddQ[k], (-1)^Boole[IntegerQ[(k + 1)/4]] (k - 1)!, 0] A191700[n - k], {k, 1, n}]; a[0] = 1; a[n_] := a[n] = Sum[Binomial[n - 1, k - 1] A191700[k] a[n - k], {k, 1, n}]; Table[a[n], {n, 0, 21}]
PROG
(PARI) seq(n)={Vec(serlaplace(exp(1/(1 - atan(x + O(x*x^n))) - 1)))} \\ Andrew Howroyd, Jan 22 2020
KEYWORD
sign
AUTHOR
Ilya Gutkovskiy, Jan 22 2020
STATUS
approved