cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-10 of 18 results. Next

A134954 Number of "hyperforests" on n labeled nodes, i.e., hypergraphs that have no cycles, assuming that each edge contains at least two vertices.

Original entry on oeis.org

1, 1, 2, 8, 55, 562, 7739, 134808, 2846764, 70720278, 2021462055, 65365925308, 2359387012261, 94042995460130, 4102781803365418, 194459091322828280, 9950303194613104995, 546698973373090998382, 32101070021048906407183, 2006125858248695722280564
Offset: 0

Views

Author

Don Knuth, Jan 26 2008

Keywords

Comments

The part of the name "assuming that each edge contains at least two vertices" is ambiguous. It may mean that not all n vertices have to be covered by some edge of the hypergraph, i.e., it is not necessarily a spanning hyperforest. However it is common to represent uncovered vertices as singleton edges, as in my example. - Gus Wiseman, May 20 2018

Examples

			From _Gus Wiseman_, May 20 2018: (Start)
The a(3) = 8 labeled spanning hyperforests are the following:
{{1,2,3}}
{{1,3},{2,3}}
{{1,2},{2,3}}
{{1,2},{1,3}}
{{3},{1,2}}
{{2},{1,3}}
{{1},{2,3}}
{{1},{2},{3}}
(End)

References

  • D. E. Knuth: The Art of Computer Programming, Volume 4, Generating All Combinations and Partitions Fascicle 3, Section 7.2.1.4. Generating all partitions. Page 38, Algorithm H. - Washington Bomfim, Sep 25 2008

Links

Crossrefs

Cf. A030019, A035053, A048143, A054921, A134955, A134956, A134957, A144959, A242817, A304716, A304717, A304867, A304911, A304912.

Programs

  • Maple
    b:= proc(n) option remember; add(Stirling2(n-1,i) *n^(i-1), i=0..n-1) end: B:= proc(n) x-> add(b(k) *x^k/k!, k=0..n) end: a:= n-> coeff(series(exp(B(n)(x)), x, n+1), x,n) *n!: seq(a(n), n=0..30);  # Alois P. Heinz, Sep 09 2008
  • Mathematica
    b[n_] := b[n] = Sum[StirlingS2[n-1, i]*n^(i-1), {i, 0, n-1}]; B[n_][x_] := Sum[b[k] *x^k/k!, {k, 0, n}]; a[0]=1; a[n_] := SeriesCoefficient[ Exp[B[n][x]], {x, 0, n}] *n!; Table[a[n], {n, 0, 30}] (* Jean-François Alcover, Feb 13 2015, after Alois P. Heinz *)

Formula

Exponential transform of A030019. - N. J. A. Sloane, Jan 30 2008
Binomial transform of A304911. - Gus Wiseman, May 20 2018
a(n) = Sum of n!*Product_{k=1..n} (A030019(k)/k!)^c_k / (c_k)! over all the partitions of n, c_1 + 2c_2 + ... + nc_n; c_1, c_2, ..., c_n >= 0. - Washington Bomfim, Sep 25 2008
a(n) ~ exp((n+1)/LambertW(1)) * n^(n-2) / (sqrt(1+LambertW(1)) * exp(2*n+2) * (LambertW(1))^n). - Vaclav Kotesovec, Jul 26 2014

A304867 Number of non-isomorphic hypertrees of weight n.

Original entry on oeis.org

1, 1, 1, 1, 2, 2, 5, 6, 13, 20, 41, 70, 144, 266, 545, 1072, 2210, 4491, 9388, 19529, 41286, 87361, 186657, 399927, 862584, 1866461, 4058367, 8852686, 19384258, 42570435, 93783472, 207157172, 458805044, 1018564642, 2266475432, 5053991582, 11292781891, 25280844844
Offset: 0

Views

Author

Gus Wiseman, May 20 2018

Keywords

Comments

A hypertree E is a connected antichain of finite sets satisfying Sum_{e in E} (|e| - 1) = |U(E)| - 1. The weight of a hypertree is the sum of cardinalities of its elements. Weight is generally not the same as number of vertices (see A035053).
From Kevin Ryde, Feb 25 2020: (Start)
a(n), except at n=1, is the number of free trees of n edges (so n+1 vertices) where any two leaves are an even distance apart. All trees are bipartite graphs and this condition is equivalent to all leaves being in the same bipartite half. The diameter of a tree is always between two leaves so these trees have even diameter (A000676).
The correspondence between hypertrees and these free trees is described for instance by Bacher (start of section 1.2). In such a free tree, call a vertex "even" if it is an even distance from a leaf. The hypertree vertices are these even vertices. Each hyperedge is the set of vertices surrounding an odd vertex, so hypertree weight is the total number of edges in the free tree.
(End)

Examples

			Non-isomorphic representatives of the a(6) = 5 hypertrees are the following:
  {{1,2,3,4,5,6}}
  {{1,2},{1,3,4,5}}
  {{1,2,3},{1,4,5}}
  {{1,2},{1,3},{1,4}}
  {{1,2},{1,3},{2,4}}
Non-isomorphic representatives of the a(7) = 6 hypertrees are the following:
  {{1,2,3,4,5,6,7}}
  {{1,2},{1,3,4,5,6}}
  {{1,2,3},{1,4,5,6}}
  {{1,2},{1,3},{1,4,5}}
  {{1,2},{1,3},{2,4,5}}
  {{1,3},{2,4},{1,2,5}}
From _Kevin Ryde_, Feb 25 2020: (Start)
a(6) = 5 hypertrees of weight 6 and their corresponding free trees of 6 edges (7 vertices).  Each * is an "odd" vertex (odd distance to a leaf).  Each hyperedge is the set of "even" vertices surrounding an odd.
  {1,2,3,4,5,6}       3   2
                       \ /
                      4-*-1      (star 7)
                       / \
                      5   6
  .
  {1,2},{1,3,4,5}               /-3
                      2--*--1--*--4
                                \-5
  .
  {1,2,3},{1,4,5}     2-\       /-4
                         *--1--*
                      3-/       \-5
  .
  {1,2},{1,3},{1,4}    /-*--2
                      1--*--3
                       \-*--4
  .
  {1,2},{2,4},{1,3}   3--*--1--*--2--*--4   (path 7)
(End)

Links

Crossrefs

Cf. A007716, A030019, A035053, A048143, A054921, A134955, A134957, A144959, A304911, A304912, A318601.

Programs

  • 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]&];
seq[40] (* Jean-François Alcover, Feb 08 2020, after Andrew Howroyd *)
  • PARI
    EulerT(v)={Vec(exp(x*Ser(dirmul(v,vector(#v,n,1/n))))-1, -#v)}
c(n)={my(v=[1]); for(i=1, ceil(n/2), v=concat([1], EulerT(concat([0], EulerT(v))))); v}
seq(n)={my(u=c(n)); Vec(x*Ser(EulerT(u))*(1-x*Ser(u)) + (1 - x)*Ser(u) + x + O(x*x^n))} \\ Andrew Howroyd, Aug 29 2018

Formula

a(n) = Sum_{k=1..floor(n/2)} A318601(n+1-k, k). - Andrew Howroyd, Aug 29 2018

Extensions

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

A134955 Number of "hyperforests" on n unlabeled nodes, i.e., hypergraphs that have no cycles, assuming that each edge contains at least two vertices.

Original entry on oeis.org

1, 1, 2, 4, 9, 20, 50, 128, 351, 1009, 3035, 9464, 30479, 100712, 340072, 1169296, 4082243, 14438577, 51643698, 186530851, 679530937, 2494433346, 9219028889, 34280914106, 128179985474, 481694091291, 1818516190252, 6894350122452
Offset: 0

Views

Author

Don Knuth, Jan 26 2008

Keywords

Comments

A hyperforest is an antichain of finite nonempty sets (edges) whose connected components are hypertrees. It is spanning if all vertices are covered by some edge. However, it is common to represent uncovered vertices as singleton edges. For example, {{1,2},{1,4}} and {{3},{1,2},{1,4}} may represent the same hyperforest, the former being free of singletons (see example 2) and the latter being spanning (see example 1). This is different from a hyperforest with singleton edges allowed, which is one whose non-singleton edges only are required to form an antichain. For example, {{1},{2},{1,3},{2,3}} is a hyperforest with singleton edges allowed. - Gus Wiseman, May 22 2018
Equivalently, number of block graphs on n nodes, that is, graphs where every block is a complete graph. These graphs can be characterized as induced-diamond-free chordal graphs. - Falk Hüffner, Jul 25 2019

Examples

			From _Gus Wiseman_, May 20 2018: (Start)
Non-isomorphic representatives of the a(4) = 9 spanning hyperforests are the following:
  {{1,2,3,4}}
  {{1},{2,3,4}}
  {{1,2},{3,4}}
  {{1,4},{2,3,4}}
  {{1},{2},{3,4}}
  {{1},{2,4},{3,4}}
  {{1,3},{2,4},{3,4}}
  {{1,4},{2,4},{3,4}}
  {{1},{2},{3},{4}}
Non-isomorphic representatives of the a(4) = 9 hyperforests spanning up to 4 vertices without singleton edges are the following:
  {}
  {{1,2}}
  {{1,2,3}}
  {{1,2,3,4}}
  {{1,2},{3,4}}
  {{1,3},{2,3}}
  {{1,4},{2,3,4}}
  {{1,3},{2,4},{3,4}}
  {{1,4},{2,4},{3,4}}
(End)

References

  • D. E. Knuth: The Art of Computer Programming, Volume 4, Generating All Combinations and Partitions Fascicle 3, Section 7.2.1.4. Generating all partitions. Page 38, Algorithm H. - Washington Bomfim, Sep 25 2008

Links

Crossrefs

Cf. A030019, A035053 (hypertrees), A048143, A054921, A134954 (labeled case), A134955, A134957, A144959, A245566, A304716, A304717, A304867, A304911.

Programs

  • Maple
    with(numtheory): etr:= proc(p) local b; b:=proc(n) option remember; `if`(n=0,1, add(add(d*p(d), d=divisors(j)) *b(n-j), j=1..n)/n) end end: b:= etr(B): c:= etr(b): B:= n-> if n=0 then 0 else c(n-1) fi: C:= etr(B): aa:= proc(n) option remember; B(n)+C(n) -add(B(k)*C(n-k), k=0..n) end: a:= etr(aa): seq(a(n), n=0..27); # Alois P. Heinz, Sep 09 2008
  • 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]; b = etr[B]; c = etr[b]; B[n_] := If[n == 0, 0, c[n-1]]; CC = etr[B]; aa[n_] := aa[n] = B[n]+CC[n]-Sum[B[k]*CC[n-k], {k, 0, n}]; a = etr[aa]; Table[a[n], {n, 0, 27}] (* Jean-François Alcover, Feb 13 2015, after Alois P. Heinz*)
  • PARI
    \\ here b is A007563 as vector
EulerT(v)={Vec(exp(x*Ser(dirmul(v,vector(#v,n,1/n))))-1, -#v)}
b(n)={my(v=[1]);for(i=2, n, v=concat([1], EulerT(EulerT(v)))); v}
seq(n)={my(u=b(n)); concat([1], EulerT(Vec(x*Ser(EulerT(u))*(1-x*Ser(u)))))} \\ Andrew Howroyd, May 22 2018

Formula

Euler transform of A035053. - N. J. A. Sloane, Jan 30 2008
a(n) = Sum of prod_{k=1}^n\,{A035053(k) + c_k -1 /choose c_k} over all the partitions of n, c_1 + 2c_2 + ... + nc_n; c_1, c_2, ..., c_n >= 0. - Washington Bomfim, Sep 25 2008
a(n) ~ c * d^n / n^(5/2), where d = 4.189610958393826965527036454524... (see A245566), c = 0.36483930544... . - Vaclav Kotesovec, Jul 26 2014

A134957 Number of hyperforests with n unlabeled vertices: analog of A134955 when edges of size 1 are allowed (with no two equal edges).

Original entry on oeis.org

1, 2, 6, 20, 75, 310, 1422, 7094, 37877, 213610, 1256422, 7641700, 47735075, 304766742, 1981348605, 13079643892, 87480944764, 591771554768, 4042991170169, 27864757592632, 193549452132550, 1353816898675732, 9529263306483357, 67457934248821368, 480019516988969011
Offset: 0

Views

Author

Don Knuth, Jan 26 2008

Keywords

Examples

			From _Gus Wiseman_, May 20 2018: (Start)
Non-isomorphic representatives of the a(3) = 20 hyperforests are the following:
  {}
  {{1}}
  {{1,2}}
  {{1,2,3}}
  {{1},{2}}
  {{1},{2,3}}
  {{2},{1,2}}
  {{3},{1,2,3}}
  {{1,3},{2,3}}
  {{1},{2},{3}}
  {{1},{2},{1,2}}
  {{1},{3},{2,3}}
  {{2},{3},{1,2,3}}
  {{2},{1,3},{2,3}}
  {{3},{1,3},{2,3}}
  {{1,2},{1,3},{2,3}}
  {{1},{2},{3},{2,3}}
  {{1},{2},{3},{1,2,3}}
  {{1},{2},{1,3},{2,3}}
  {{2},{3},{1,3},{2,3}}
  {{3},{1,2},{1,3},{2,3}}
  {{1},{2},{3},{1,3},{2,3}}
  {{2},{3},{1,2},{1,3},{2,3}}
  {{1},{2},{3},{1,2},{1,3},{2,3}}
(End)

Links

Crossrefs

Cf. A030019, A035053, A048143, A054921, A134955, A134957, A134959, A144959, A304867, A304911.

Programs

  • 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];
b[n_] := Module[{v = {1}}, For[i = 2, i <= n, i++, v = Join[{1}, EulerT[EulerT[2 v]]]]; v];
seq[n_] := Module[{u = 2 b[n]}, Join[{1}, EulerT[ser[EulerT[u]]*(1 - x*ser[u]) + O[x]^n // CoefficientList[#, x]&]]];
seq[24] (* Jean-François Alcover, Feb 10 2020, after Andrew Howroyd *)
  • PARI
    \\ here b(n) is A318494 as vector
EulerT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, n, 1/n))))-1, -#v)}
b(n)={my(v=[1]); for(i=2, n, v=concat([1], EulerT(EulerT(2*v)))); v}
seq(n)={my(u=2*b(n)); concat([1], EulerT(Vec(Ser(EulerT(u))*(1-x*Ser(u)))))} \\ Andrew Howroyd, Aug 27 2018

Formula

Euler transform of A134959. - Gus Wiseman, May 20 2018

Extensions

Terms a(7) and beyond from Andrew Howroyd, Aug 27 2018

A144958 Number of unlabeled acyclic graphs covering n vertices.

Original entry on oeis.org

1, 0, 1, 1, 3, 4, 10, 17, 39, 77, 176, 381, 891, 2057, 4941, 11915, 29391, 73058, 184236, 468330, 1202349, 3108760, 8097518, 21218776, 55925742, 148146312, 394300662, 1053929982, 2828250002, 7617271738, 20584886435, 55802753243
Offset: 0

Views

Author

Washington Bomfim, Sep 27 2008

Keywords

Comments

a(n) is the number of forests with n unlabeled nodes without isolated vertices. This follows from the fact that for n>0 A005195(n-1) counts the forests with one or more isolated nodes.
The labeled version is A105784. The connected case is A000055. This is the covering case of A005195. - Gus Wiseman, Apr 29 2024

Examples

			From _Gus Wiseman_, Apr 29 2024: (Start)
Edge-sets of non-isomorphic representatives of the a(0) = 1 through a(5) = 4 forests:
  {}    .    {12}    {13,23}    {12,34}       {12,35,45}
                                {13,24,34}    {13,24,35,45}
                                {14,24,34}    {14,25,35,45}
                                              {15,25,35,45}
(End)

Links

Crossrefs

The connected case is A000055.
This is the covering case of A005195, labeled A001858.
The labeled version is A105784.
For triangles instead of cycles we have A372169, non-covering A006785.
Unique cycle: A372191 (lab A372195), non-covering A236570 (lab A372193).
A006125 counts simple graphs, unlabeled A000088.
A006129 counts covering graphs, unlabeled A002494.
Cf. A000272, A053530, A054548, A137917, A144959, A372174.

Programs

  • Mathematica
    brute[m_]:=First[Sort[Table[Sort[Sort/@(m/.Rule@@@Table[{i,p[[i]]},{i,Length[p]}])],{p,Permutations[Union@@m]}]]];
cyc[y_]:=Select[Join@@Table[Select[Join@@Permutations/@Subsets[Union@@y,{k}],And@@Table[MemberQ[Sort/@y,Sort[{#[[i]],#[[If[i==k,1,i+1]]]}]],{i,k}]&],{k,3,Length[y]}],Min@@#==First[#]&];
Table[Length[Union[Union[brute/@Select[Subsets[Subsets[Range[n],{2}]],Union@@#==Range[n]&&Length[cyc[#]]==0&]]]],{n,0,5}] (* Gus Wiseman, Apr 29 2024 *)
  • PARI
    EulerT(v)={Vec(exp(x*Ser(dirmul(v,vector(#v,n,1/n))))-1, -#v)}
TreeGf(N)={my(A=vector(N, j, 1)); for (n=1, N-1, A[n+1] = 1/n * sum(k=1, n, sumdiv(k, d, d*A[d]) * A[n-k+1] ) ); x*Ser(A)}
seq(n)={my(t=TreeGf(n), v=EulerT(Vec(t - t^2/2 + subst(t,x,x^2)/2))); concat([1,0], vector(#v-1, i, v[i+1]-v[i]))} \\ Andrew Howroyd, Aug 01 2024

Formula

a(n) = A005195(n) - A005195(n-1).

Extensions

Name changed and 1 prepended by Gus Wiseman, Apr 29 2024.

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

Original entry on oeis.org

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

Views

Author

Gus Wiseman, May 20 2018

Keywords

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).

Examples

			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}}

Links

Crossrefs

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

Programs

  • 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 *)
  • 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

Formula

Euler transform of A304867.

Extensions

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

A304911 Number of labeled hyperforests spanning n vertices without singleton edges.

Original entry on oeis.org

1, 0, 1, 4, 32, 351, 5057, 90756, 1956971, 49366904, 1427680932, 46590895869, 1694163054597, 67938488277050, 2978980898086377, 141801848209013050, 7282651452378019772, 401410357608479625207, 23635996827115264290005
Offset: 0

Views

Author

Gus Wiseman, May 20 2018

Keywords

Examples

			The a(3) = 4 hyperforests are {{1,2,3}}, {{1,3},{2,3}}, {{1,2},{2,3}}, {{1,2},{1,3}}.

Crossrefs

Cf. A006126, A030019, A035053, A048143, A054921, A134954, A134955, A134957, A144959, A293993, A293994, A304716, A304717, A304867.

Formula

E.g.f.: exp(A030019(x) - x - 1) where A030019(x) is the e.g.f. of A030019.

A134959 Number of spanning hypertrees with n unlabeled vertices: analog of A035053 when edges of size 1 are allowed (with no two equal edges).

Original entry on oeis.org

1, 2, 3, 10, 35, 150, 707, 3700, 20470, 119260, 719341, 4466316, 28367118, 183620874, 1207563011, 8049914664, 54295152117, 369981325578, 2544017965638, 17633790542978, 123108792874528, 865045359778662, 6114040341515978, 43443726772579152, 310195170229429300
Offset: 0

Views

Author

Don Knuth, Jan 26 2008

Keywords

Examples

			From _Gus Wiseman_, May 20 2018: (Start)
Non-isomorphic representatives of the a(3) = 10 hypertrees are the following:
  {{1,2,3}}
  {{3},{1,2,3}}
  {{1,3},{2,3}}
  {{2},{3},{1,2,3}}
  {{2},{1,3},{2,3}}
  {{3},{1,3},{2,3}}
  {{1},{2},{3},{1,2,3}}
  {{1},{2},{1,3},{2,3}}
  {{2},{3},{1,3},{2,3}}
  {{1},{2},{3},{1,3},{2,3}}
(End)

Links

Crossrefs

Cf. A030019, A035053, A048143, A054921, A134955, A134957, A144959, A304867, A304911, A318494.

Programs

  • 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];
b[n_] := Module[{v = {1}}, For[i = 2, i <= n, i++, v = Join[{1}, EulerT[EulerT[2 v]]]]; v];
seq[n_] := Module[{u = 2 b[n]}, 1 + x*ser[EulerT[u]]*(1 - x*ser[u]) + O[x]^n // CoefficientList[#, x]&];
seq[25] (* Jean-François Alcover, Feb 10 2020, after Andrew Howroyd *)
  • PARI
    \\ here b(n) is A318494 as vector
EulerT(v)={Vec(exp(x*Ser(dirmul(v,vector(#v,n,1/n))))-1, -#v)}
b(n)={my(v=[1]); for(i=2, n, v=concat([1], EulerT(EulerT(2*v)))); v}
seq(n)={my(u=2*b(n)); Vec(1 + x*Ser(EulerT(u))*(1-x*Ser(u)))} \\ Andrew Howroyd, Aug 27 2018

Formula

Inverse Euler transform of A134957. - Gus Wiseman, May 20 2018

Extensions

Terms a(7) and beyond from Andrew Howroyd, Aug 27 2018

A134956 Number of hyperforests with n labeled vertices: analog of A134954 when edges of size 1 are allowed (with no two equal edges).

Original entry on oeis.org

1, 2, 8, 64, 880, 17984, 495296, 17255424, 728771584, 36208782336, 2069977144320, 133869415030784, 9664049202221056, 770400218809384960, 67219977066339008512, 6372035504466437079040, 652103070162164448952320, 71656927837957783339925504
Offset: 0

Views

Author

Don Knuth, Jan 26 2008

Keywords

Examples

			From _Gus Wiseman_, May 21 2018: (Start)
The a(2) = 8 hyperforests are the following:
  {{1},{2},{1,2}}
  {{1},{1,2}}
  {{2},{1,2}}
  {{1,2}}
  {{1},{2}}
  {{1}}
  {{2}}
  {}
(End)

References

  • D. E. Knuth: The Art of Computer Programming, Volume 4, Generating All Combinations and Partitions Fascicle 3, Section 7.2.1.4. Generating all partitions. Page 38, Algorithm H. - Washington Bomfim, Sep 25 2008

Links

Crossrefs

Cf. A134958. - Washington Bomfim, Sep 25 2008
Cf. A030019, A035053, A048143, A054921, A134954, A134955, A134957, A144959, A304716, A304717, A304867, A304911.

Programs

  • Maple
    with(combinat): p:= proc(n) option remember; add(stirling2(n-1, i) *n^(i-1), i=0..n-1) end: g:= proc(n) option remember; p(n) +add(binomial(n-1, k-1) *p(k) *g(n-k), k=1..n-1) end: a:= n-> `if`(n=0, 1, 2^n * g(n)): seq(a(n), n=0..30); # Alois P. Heinz, Oct 07 2008
  • Mathematica
    p[n_] := p[n] = Sum[ StirlingS2[n-1, i]*n^(i-1), {i, 0, n-1}]; g[n_] := g[n] = p[n] + Sum[Binomial[n-1, k-1]*p[k]*g[n-k], {k, 1, n-1}]; a[n_] := If[n == 0, 1, 2^n* g[n]]; Table[a[n], {n, 0, 30}] (* Jean-François Alcover, Feb 13 2015, after Alois P. Heinz *)

Formula

Equals 2^n*A134954(n).
a(n) = Sum of n!prod_{k=1}^n\{ frac{ A134958(k)^{c_k} }{ k!^{c_k} c_k! } } over all the partitions of n, c_1 + 2c_2 + ... + nc_n; c_1, c_2, ..., c_n >= 0. - Washington Bomfim, Sep 25 2008

A304386 Number of unlabeled hypertrees (connected antichains with no cycles) spanning up to n vertices and allowing singleton edges.

Original entry on oeis.org

1, 2, 5, 15, 50, 200, 907, 4607, 25077, 144337, 863678, 5329994, 33697112, 217317986, 1424880997, 9474795661, 63769947778, 433751273356, 2977769238994, 20611559781972, 143720352656500, 1008765712435162, 7122806053951140, 50566532826530292, 360761703055959592
Offset: 0

Views

Author

Gus Wiseman, May 21 2018

Keywords

Examples

			Non-isomorphic representatives of the a(3) = 15 hypertrees are the following:
  {}
  {{1}}
  {{1,2}}
  {{1,2,3}}
  {{2},{1,2}}
  {{1,3},{2,3}}
  {{3},{1,2,3}}
  {{1},{2},{1,2}}
  {{3},{1,2},{2,3}}
  {{3},{1,3},{2,3}}
  {{2},{3},{1,2,3}}
  {{1},{2},{3},{1,2,3}}
  {{2},{3},{1,2},{1,3}}
  {{2},{3},{1,3},{2,3}}
  {{1},{2},{3},{1,3},{2,3}}

Links

Crossrefs

Cf. A030019, A035053, A048143, A054921, A134954, A134955, A134957, A134959, A144959, A304386, A304717, A304867, A304911, A304912, A304968, A304970.

Programs

  • PARI
    \\ here b(n) is A318494 as vector
EulerT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, n, 1/n))))-1, -#v)}
b(n)={my(v=[1]); for(i=2, n, v=concat([1], EulerT(EulerT(2*v)))); v}
seq(n)={my(u=2*b(n)); Vec(1 + x*Ser(EulerT(u))*(1-x*Ser(u))/(1-x))} \\ Andrew Howroyd, Aug 27 2018

Formula

Partial sums of b(1) = 1, b(n) = A134959(n) otherwise.

Extensions

Terms a(7) and beyond from Andrew Howroyd, Aug 27 2018
Showing 1-10 of 18 results. Next

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.