A304912 Number of non-isomorphic spanning hyperforests of weight n.

1, 1, 2, 3, 6, 9, 18, 29, 56, 97, 186, 337, 657, 1238, 2442, 4768, 9569, 19174, 39151, 80154, 166211, 346239, 727853, 1537611, 3270710, 6989669, 15018389, 32405378, 70230238, 152772075, 333552711, 730632928, 1605459844, 3537861659, 7817447580, 17317397837

Offset

0,3

Comments

A spanning hyperforest is an antichain of finite nonempty sets, which cover a set of n vertices, whose connected components are hypertrees (see A304867). The weight of a hypertree is the sum of cardinalities of its elements. Weight is generally not the same as number of vertices (see A134957).

Links

Andrew Howroyd, Table of n, a(n) for n = 0..500

Formula

Euler transform of A304867.

Example

The a(6) = 18 spanning hyperforests are the following:
  {{1,2,3,4,5,6}}
  {{1},{2,3,4,5,6}}
  {{1,2},{3,4,5,6}}
  {{1,5},{2,3,4,5}}
  {{1,2,3},{4,5,6}}
  {{1,2,5},{3,4,5}}
  {{1},{2},{3,4,5,6}}
  {{1},{2,3},{4,5,6}}
  {{1},{2,5},{3,4,5}}
  {{1,2},{3,4},{5,6}}
  {{1,2},{3,5},{4,5}}
  {{1,3},{2,4},{3,4}}
  {{1,4},{2,4},{3,4}}
  {{1},{2},{3},{4,5,6}}
  {{1},{2},{3,4},{5,6}}
  {{1},{2},{3,5},{4,5}}
  {{1},{2},{3},{4},{5,6}}
  {{1},{2},{3},{4},{5},{6}}

Mathematica

etr[p_] := Module[{b}, b[n_] := b[n] = If[n == 0, 1, Sum[Sum[d*p[d], {d, Divisors[j]}]*b[n - j], {j, 1, n}]/n]; b];
EulerT[v_List] := With[{q = etr[v[[#]]&]}, q /@ Range[Length[v]]];
ser[v_] := Sum[v[[i]] x^(i - 1), {i, 1, Length[v]}] + O[x]^Length[v];
c[n_] := Module[{v = {1}}, For[i = 1, i <= Ceiling[n/2], i++, v = Join[{1}, EulerT[Join[{0}, EulerT[v]]]]]; v];
seq[n_] := Module[{u = c[n]}, x*ser[EulerT[u]]*(1 - x*ser[u]) + (1 - x)* ser[u] + x + O[x]^n // CoefficientList[#, x]& // Rest // EulerT // Prepend[#, 1]&];
seq[36] (* Jean-François Alcover, Feb 09 2020, after Andrew Howroyd *)

Prog

(PARI) EulerT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, n, 1/n))))-1, -#v)}
c(n)={my(v=[1]); for(i=2, ceil(n/2), v=concat([1], EulerT(concat([0], EulerT(v))))); v}
seq(n)={my(u=c(n)); concat([1], EulerT(Vec(x*Ser(EulerT(u))*(1-x*Ser(u)) + (1 - x)*(Ser(u) - 1)+ O(x*x^n))))} \\ Andrew Howroyd, Aug 29 2018

Crossrefs

Cf. A007716, A035053, A048143, A054921, A134954, A134955, A134957, A144959, A286520, A293993, A293994, A304911, A304867.

Sequence in context: A191398 A066313 A224958 * A018499 A107847 A059966

Adjacent sequences: A304909 A304910 A304911 * A304913 A304914 A304915

Keyword

nonn

Author

Gus Wiseman, May 20 2018

Extensions

Terms a(10) and beyond from Andrew Howroyd, Aug 29 2018

Status

approved