OFFSET
0,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..300
FORMULA
E.g.f.: exp(1 - cosh(x)) (even powers only). - Ilya Gutkovskiy, Jan 27 2020
EXAMPLE
a(5) = 254. Consider the shapes [10], [8, 2], [6, 4], [6, 2, 2], [4, 4, 2], [4, 2, 2, 2], [2, 2, 2, 2, 2]. Computing the number of associated set partitions gives -1 + 45 + 210 - 630 - 1575 + 3150 - 945 = 254.
MAPLE
b:= proc(n, t) option remember; `if`(n=0, 1-2*t, add(
b(n-2*j, 1-t)*binomial(n-1, 2*j-1), j=1..n/2))
end:
a:= n-> b(2*n, 0):
seq(a(n), n=0..20); # Alois P. Heinz, Jun 28 2016
MATHEMATICA
Table[Sum[BellY[2 n, k, Mod[Range[2 n], 2] - 1], {k, 0, 2 n}], {n, 0, 20}] (* Vladimir Reshetnikov, Nov 09 2016 *)
PROG
(Sage)
def ComplementaryBell(m, n):
shapes = ([x*m for x in p] for p in Partitions(n))
return sum((-1)^len(s)*SetPartitions(sum(s), s).cardinality() for s in shapes)
A260884 = lambda n: ComplementaryBell(2, n)
[A260884(n) for n in (0..19)]
(PARI) a(n) = {polcoef(serlaplace(exp(1 - cosh(x + O(x^(2*n + 1))))), 2*n)} \\ Andrew Howroyd, Jan 27 2020
CROSSREFS
KEYWORD
sign
AUTHOR
Peter Luschny, Sep 18 2015
STATUS
approved