OFFSET
0,3
COMMENTS
Also the number of strict multiset partitions of {1, 1, 2, 2, 3, 3, ..., n, n}. For example, the a(2) = 5 strict multiset partitions of {1, 1, 2, 2} are (1122), (1)(122), (2)(112), (11)(22), (1)(2)(12). - Gus Wiseman, Jul 18 2018
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..300
FORMULA
EXAMPLE
From Gus Wiseman, Sep 02 2019: (Start)
These are set-systems covering {1..n} with vertex-degrees <= 2. For example, the a(3) = 40 covers are:
{123} {1}{23} {1}{2}{3} {1}{2}{3}{12}
{2}{13} {1}{2}{13} {1}{2}{3}{13}
{3}{12} {1}{2}{23} {1}{2}{3}{23}
{1}{123} {1}{3}{12} {1}{2}{13}{23}
{12}{13} {1}{3}{23} {1}{2}{3}{123}
{12}{23} {2}{3}{12} {1}{3}{12}{23}
{13}{23} {2}{3}{13} {2}{3}{12}{13}
{2}{123} {1}{12}{23}
{3}{123} {1}{13}{23}
{12}{123} {1}{2}{123}
{13}{123} {1}{3}{123}
{23}{123} {2}{12}{13}
{2}{13}{23}
{2}{3}{123}
{3}{12}{13}
{3}{12}{23}
{12}{13}{23}
{1}{23}{123}
{2}{13}{123}
{3}{12}{123}
(End)
MATHEMATICA
facs[n_]:=facs[n]=If[n<=1, {{}}, Join@@Table[Map[Prepend[#, d]&, Select[facs[n/d], Min@@#>=d&]], {d, Rest[Divisors[n]]}]];
Table[Length[Select[facs[Array[Prime, n, 1, Times]^2], UnsameQ@@#&]], {n, 0, 6}] (* Gus Wiseman, Jul 18 2018 *)
m = 20;
a094577[n_] := Sum[Binomial[n, k]*BellB[2 n - k], {k, 0, n}];
egf = Exp[(1 - Exp[x])/2]*Sum[a094577[n]*(x/2)^n/n!, {n, 0, m}] + O[x]^m;
CoefficientList[egf + O[x]^m, x]*Range[0, m-1]! (* Jean-François Alcover, May 13 2019 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Goran Kilibarda, Vladeta Jovovic, May 12 2004
STATUS
approved