OFFSET
0,5
COMMENTS
Number of partitions of n-set into blocks congruent to 1 mod 3.
FORMULA
E.g.f.: exp(exp(x)/3 - 2*sin(Pi/6 - sqrt(3)*x/2) / (3*exp(x/2))). - Vaclav Kotesovec, Apr 15 2020
a(0) = 1; a(n) = Sum_{k=0..floor((n-1)/3)} binomial(n-1,3*k) * a(n-3*k-1). - Seiichi Manyama, Sep 22 2023
MATHEMATICA
nmax = 27; CoefficientList[Series[Exp[Sum[x^(3 k + 1)/(3 k + 1)!, {k, 0, nmax}]], {x, 0, nmax}], x] Range[0, nmax]!
a[n_] := a[n] = Sum[Boole[MemberQ[{1}, Mod[k, 3]]] Binomial[n - 1, k - 1] a[n - k], {k, 1, n}]; a[0] = 1; Table[a[n], {n, 0, 27}]
nmax = 30; CoefficientList[Series[Exp[Exp[x]/3 - 2*Sin[Pi/6 - Sqrt[3]*x/2] / (3*Exp[x/2])], {x, 0, nmax}], x] * Range[0, nmax]! (* Vaclav Kotesovec, Apr 15 2020 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Apr 08 2020
STATUS
approved