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.

A055621 Number of covers of an unlabeled n-set.

Original entry on oeis.org

1, 1, 4, 34, 1952, 18664632, 12813206150470528, 33758171486592987151274638874693632, 1435913805026242504952006868879460423801146743462225386100617731367239680
Offset: 0

Views

Author

Vladeta Jovovic, Jun 04 2000

Keywords

Examples

			There are 4 nonisomorphic covers of {1,2}, namely {{1},{2}}, {{1,2}}, {{1},{1,2}} and {{1},{2},{1,2}}.
From _Gus Wiseman_, Aug 14 2019: (Start)
Non-isomorphic representatives of the a(3) = 34 covers:
  {123}  {1}{23}    {1}{2}{3}      {1}{2}{3}{23}
         {13}{23}   {1}{3}{23}     {1}{2}{13}{23}
         {3}{123}   {2}{13}{23}    {1}{2}{3}{123}
         {23}{123}  {2}{3}{123}    {2}{3}{13}{23}
                    {3}{13}{23}    {1}{3}{23}{123}
                    {12}{13}{23}   {2}{3}{23}{123}
                    {1}{23}{123}   {3}{12}{13}{23}
                    {3}{23}{123}   {2}{13}{23}{123}
                    {13}{23}{123}  {3}{13}{23}{123}
                                   {12}{13}{23}{123}
.
  {1}{2}{3}{13}{23}     {1}{2}{3}{12}{13}{23}    {1}{2}{3}{12}{13}{23}{123}
  {1}{2}{3}{23}{123}    {1}{2}{3}{13}{23}{123}
  {2}{3}{12}{13}{23}    {2}{3}{12}{13}{23}{123}
  {1}{2}{13}{23}{123}
  {2}{3}{13}{23}{123}
  {3}{12}{13}{23}{123}
(End)
		

References

  • F. Bergeron, G. Labelle and P. Leroux, Combinatorial Species and Tree-Like Structures, Cambridge, 1998, p. 78 (2.3.39)

Crossrefs

Unlabeled set-systems are A000612 (partial sums).
The version with empty edges allowed is A003181.
The labeled version is A003465.
The T_0 case is A319637.
The connected case is A323819.
The T_1 case is A326974.

Programs

  • Maple
    b:= proc(n, i, l) `if`(n=0, 2^(w-> add(mul(2^igcd(t, l[h]),
          h=1..nops(l)), t=1..w)/w)(ilcm(l[])), `if`(i<1, 0,
          add(b(n-i*j, i-1, [l[], i$j])/j!/i^j, j=0..n/i)))
        end:
    a:= n-> `if`(n=0, 2, b(n$2, [])-b(n-1$2, []))/2:
    seq(a(n), n=0..8);  # Alois P. Heinz, Aug 14 2019
  • Mathematica
    b[n_, i_, l_] := b[n, i, l] = If[n==0, 2^Function[w, Sum[Product[2^GCD[t, l[[h]]], {h, 1, Length[l]}], {t, 1, w}]/w][If[l=={}, 1, LCM@@l]], If[i<1, 0, Sum[b[n-i*j, i-1, Join[l, Table[i, {j}]]]/j!/i^j, {j, 0, n/i}]]];
    a[n_] := If[n==0, 2, b[n, n, {}] - b[n-1, n-1, {}]]/2;
    a /@ Range[0, 8] (* Jean-François Alcover, Jan 31 2020, after Alois P. Heinz *)

Formula

a(n) = (A003180(n) - A003180(n-1))/2 = A000612(n) - A000612(n-1) for n>0.
Euler transform of A323819. - Gus Wiseman, Aug 14 2019

Extensions

More terms from David Moews (dmoews(AT)xraysgi.ims.uconn.edu) Jul 04 2002
a(0) = 1 prepended by Gus Wiseman, Aug 14 2019