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.

Previous Showing 61-70 of 182 results. Next

A368601 Number of ways to choose a set of n nonempty subsets of {1..n} such that it is possible to choose a different element from each.

Original entry on oeis.org

1, 1, 3, 32, 1201, 151286, 62453670, 84707326890, 384641855115279
Offset: 0

Views

Author

Gus Wiseman, Jan 01 2024

Keywords

Comments

The axiom of choice says that, given any set of nonempty sets Y, it is possible to choose a set containing an element from each. The strict version requires this set to have the same cardinality as Y, meaning no element is chosen more than once.

Examples

			The a(2) = 3 set-systems:
  {{1},{2}}
  {{1},{1,2}}
  {{2},{1,2}}
Non-isomorphic representatives of the a(3) = 32 set-systems:
  {{1},{2},{3}}
  {{1},{2},{1,3}}
  {{1},{2},{1,2,3}}
  {{1},{1,2},{1,3}}
  {{1},{1,2},{2,3}}
  {{1},{1,2},{1,2,3}}
  {{1},{2,3},{1,2,3}}
  {{1,2},{1,3},{2,3}}
  {{1,2},{1,3},{1,2,3}}
		

Crossrefs

For a unique choice we have A003024, any length A367904 (ranks A367908).
Sets of n nonempty subsets of {1..n} are counted by A136556.
For any length we have A367902, ranks A367906, no singletons A367770.
The complement is A368600, any length A367903 (see also A367907, A367769).
A000372 counts antichains, covering A006126, nonempty A014466.
A003465 counts covering set-systems, unlabeled A055621.
A058891 counts set-systems, unlabeled A000612.
A059201 counts covering T_0 set-systems.
A323818 counts covering connected set-systems, unlabeled A323819.

Programs

  • Mathematica
    Table[Length[Select[Subsets[Rest[Subsets[Range[n]]], {n}],Length[Select[Tuples[#], UnsameQ@@#&]]>0&]],{n,0,3}]
  • Python
    from itertools import combinations, product, chain
    def v(c):
        for elements in product(*c):
            if len(set(elements)) == len(elements):
                return True
        return False
    def a(n):
        if n == 0:
            return 1
        subsets = list(chain.from_iterable(combinations(range(1, n + 1), r) for r in
    range(1, n + 1)))
        cs = combinations(subsets, n)
        c = sum(1 for c in cs if v(c))
        return c
    [print(a(n)) for n in range(7)] # Robert P. P. McKone, Jan 02 2024

Formula

a(n) + A368600(n) = A136556(n).

Extensions

a(6) from Robert P. P. McKone, Jan 02 2024
a(7)-a(8) from Christian Sievers, Jul 25 2024

A253317 Indices in A261283 where records occur.

Original entry on oeis.org

0, 1, 2, 3, 8, 9, 10, 11, 128, 129, 130, 131, 136, 137, 138, 139, 32768, 32769, 32770, 32771, 32776, 32777, 32778, 32779, 32896, 32897, 32898, 32899, 32904, 32905, 32906, 32907, 2147483648, 2147483649, 2147483650, 2147483651, 2147483656, 2147483657
Offset: 1

Views

Author

Philippe Beaudoin, Dec 30 2014

Keywords

Comments

From Gus Wiseman, Dec 29 2023: (Start)
These are numbers whose binary indices are all powers of 2, where a binary index of n (row n of A048793) is any position of a 1 in its reversed binary expansion. For example, the terms together with their binary expansions and binary indices begin:
0: 0 ~ {}
1: 1 ~ {1}
2: 10 ~ {2}
3: 11 ~ {1,2}
8: 1000 ~ {4}
9: 1001 ~ {1,4}
10: 1010 ~ {2,4}
11: 1011 ~ {1,2,4}
128: 10000000 ~ {8}
129: 10000001 ~ {1,8}
130: 10000010 ~ {2,8}
131: 10000011 ~ {1,2,8}
136: 10001000 ~ {4,8}
137: 10001001 ~ {1,4,8}
138: 10001010 ~ {2,4,8}
139: 10001011 ~ {1,2,4,8}
For powers of 3 we have A368531.
(End)

Crossrefs

Cf. A053644 (most significant bit).
A048793 lists binary indices, length A000120, sum A029931.
A070939 gives length of binary expansion.
A096111 gives product of binary indices.

Programs

  • Maple
    a := proc(n) local k, A:
    A := [seq(0,i=1..n)]: A[1]:=0:
    for k from 1 to n-1 do
       A[k+1] := A[k-2^ilog2(k)+1]+2^(2^ilog2(k)-1): od:
    return A[n]: end proc: # Lorenzo Sauras Altuzarra, Dec 18 2019
    # second Maple program:
    a:= n-> (l-> add(l[i+1]*2^(2^i-1), i=0..nops(l)-1))(Bits[Split](n-1)):
    seq(a(n), n=1..38);  # Alois P. Heinz, Dec 13 2023
  • Mathematica
    Nest[Append[#1, #1[[-#2]] + 2^(#2 - 1)] & @@ {#, 2^(IntegerLength[Length[#], 2] - 1)} &, {0, 1}, 36] (* Michael De Vlieger, May 08 2020 *)
  • PARI
    a(n)={if(n<=1, 0, my(t=1<Andrew Howroyd, Dec 20 2019

Formula

a(1) = 0 and a(n) = a(n-A053644(n-1)) + 2^(A053644(n-1)-1). - Lorenzo Sauras Altuzarra, Dec 18 2019
a(n) = A358126(n-1) / 2. - Tilman Piesk, Dec 18 2022
a(2^n+1) = 2^(2^n-1) = A058891(n+1). - Gus Wiseman, Dec 29 2023
a(2^n) = A072639(n). - Gus Wiseman, Dec 29 2023
G.f.: 1/(1-x) * Sum_{k>=0} (2^(-1+2^k))*x^2^k/(1+x^2^k). - John Tyler Rascoe, May 22 2024

Extensions

Corrected reference in name from A253315 to A261283. - Tilman Piesk, Dec 18 2022

A279944 Number of positions in the free pure symmetric multifunction in one symbol with j-number n.

Original entry on oeis.org

1, 3, 5, 5, 7, 7, 9, 4, 7, 9, 11, 6, 9, 11, 13, 7, 8, 11, 13, 15, 9, 10, 13, 15, 9, 17, 6, 11, 12, 15, 17, 6, 11, 19, 8, 9, 13, 14, 17, 19, 8, 13, 21, 10, 11, 15, 16, 19, 11, 21, 10, 15, 23, 12, 13, 17, 18, 21, 13, 23, 12, 17, 25, 7, 14, 15, 19, 20, 23, 15, 25, 14, 19, 27, 9, 16, 17, 21, 22, 25, 9, 17, 27, 16, 21, 29, 11, 18, 19, 23, 24, 27, 11, 19, 29, 18, 23, 31, 13, 11
Offset: 1

Views

Author

Gus Wiseman, Dec 24 2016

Keywords

Comments

A free pure symmetric multifunction in one symbol f in PSM(x) is either (case 1) f = the symbol x, or (case 2) f = an expression of the form h[g_1,...,g_k] where h is in PSM(x), each of the g_i for i=1..(k>0) is in PSM(x), and for i < j we have g_i <= g_j under a canonical total ordering of PSM(x), such as the Mathematica ordering of expressions. For a positive integer n we define a free pure symmetric multifunction j(n) by: j(1)=x; j(n>1) = j(h)[j(g_1),...,j(g_k)] where n = r(h)^(p(g_1)*...*p(g_k)-1). Here r(n) is the n-th number that is not a perfect power (A007916) and p(n) is the n-th prime number (A000040). See example. Then a(n) is the number of brackets [...] plus the number of x's in j(n).

Examples

			The first 20 free pure symmetric multifunctions in x are:
j(1)  = j(1)            = x
j(2)  = j(1)[j(1)]      = x[x]
j(3)  = j(2)[j(1)]      = x[x][x]
j(4)  = j(1)[j(2)]      = x[x[x]]
j(5)  = j(3)[j(1)]      = x[x][x][x]
j(6)  = j(4)[j(1)]      = x[x[x]][x]
j(7)  = j(5)[j(1)]      = x[x][x][x][x]
j(8)  = j(1)[j(1),j(1)] = x[x,x]
j(9)  = j(2)[j(2)]      = x[x][x[x]]
j(10) = j(6)[j(1)]      = x[x[x]][x][x]
j(11) = j(7)[j(1)]      = x[x][x][x][x][x]
j(12) = j(8)[j(1)]      = x[x,x][x]
j(13) = j(9)[j(1)]      = x[x][x[x]][x]
j(14) = j(10)[j(1)]     = x[x[x]][x][x][x]
j(15) = j(11)[j(1)]     = x[x][x][x][x][x][x]
j(16) = j(1)[j(3)]      = x[x[x][x]]
j(17) = j(12)[j(1)]     = x[x,x][x][x]
j(18) = j(13)[j(1)]     = x[x][x[x]][x][x]
j(19) = j(14)[j(1)]     = x[x[x]][x][x][x][x]
j(20) = j(15)[j(1)]     = x[x][x][x][x][x][x][x].
		

Crossrefs

Cf. A279984 (numbers j(n)[x]=j(prime(n))), A277576 (numbers j(n)=x[x][x][x]...), A058891 (numbers j(n)=x[x,...,x]), A279969 (numbers j(n)=x[x[...[x]]]).

Programs

  • Mathematica
    nn=100;
    radQ[n_]:=If[n===1,False,SameQ[GCD@@FactorInteger[n][[All,2]],1]];
    rad[n_]:=rad[n]=If[n===0,1,NestWhile[#+1&,rad[n-1]+1,Not[radQ[#]]&]];
    Set@@@Array[radPi[rad[#]]==#&,nn];
    jfac[n_]:=With[{g=GCD@@FactorInteger[n+1][[All,2]]},JIX[radPi[Power[n+1,1/g]],Flatten[Cases[FactorInteger[g+1],{p_,k_}:>ConstantArray[PrimePi[p],k]]]]];
    diwt[n_]:=If[n===1,1,Apply[1+diwt[#1]+Total[diwt/@#2]&,jfac[n-1]]];
    Array[diwt,nn]

Formula

a(A007916(h)^(A000040(g_1)*...*A000040(g_k)-1)) = 1 + a(h) + a(g_1) + ... + a(g_k).

A326363 Number of maximal intersecting antichains of subsets of {1..n}.

Original entry on oeis.org

1, 2, 4, 6, 21, 169, 11749, 12160648
Offset: 0

Views

Author

Gus Wiseman, Jul 01 2019

Keywords

Comments

A set system (set of sets) is an antichain if no element is a subset of any other, and is intersecting if no two element are disjoint.

Examples

			The a(1) = 1 through a(4) = 21 maximal intersecting antichains:
  {}   {}    {}            {}
  {1}  {1}   {1}           {1}
       {2}   {2}           {2}
       {12}  {3}           {3}
             {123}         {4}
             {12}{13}{23}  {1234}
                           {12}{13}{23}
                           {12}{14}{24}
                           {13}{14}{34}
                           {23}{24}{34}
                           {12}{134}{234}
                           {13}{124}{234}
                           {14}{123}{234}
                           {23}{124}{134}
                           {24}{123}{134}
                           {34}{123}{124}
                           {12}{13}{14}{234}
                           {12}{23}{24}{134}
                           {13}{23}{34}{124}
                           {14}{24}{34}{123}
                           {123}{124}{134}{234}
		

Crossrefs

The case with nonempty, non-singleton edges is A326362.
Antichains of nonempty, non-singleton sets are A307249.
Minimal covering antichains are A046165.
Maximal intersecting antichains are A007363.
Maximal antichains of nonempty sets are A326359.

Programs

  • Mathematica
    stableSets[u_,Q_]:=If[Length[u]==0,{{}},With[{w=First[u]},Join[stableSets[DeleteCases[u,w],Q],Prepend[#,w]&/@stableSets[DeleteCases[u,r_/;r==w||Q[r,w]||Q[w,r]],Q]]]];
    fasmax[y_]:=Complement[y,Union@@(Most[Subsets[#]]&/@y)];
    Table[Length[fasmax[stableSets[Subsets[Range[n],{0,n}],Or[Intersection[#1,#2]=={},SubsetQ[#1,#2]]&]]],{n,0,5}]
    (* 2nd program *)
    n = 2^6; g = CompleteGraph[n]; i = 0;
    While[i < n, i++; j = i; While[j < n, j++; If[BitAnd[i, j] == 0 || BitAnd[i, j] == i || BitAnd[i, j] == j, g = EdgeDelete[g, i <-> j]]]];
    sets = FindClique[g, Infinity, All];
    Length[sets] (* Elijah Beregovsky, May 06 2020 *)

Formula

For n > 1, a(n) = A007363(n + 1) + 1 = A326362(n) + n + 1.

Extensions

a(7) from Elijah Beregovsky, May 06 2020

A367916 Number of sets of nonempty subsets of {1..n} with the same number of edges as covered vertices.

Original entry on oeis.org

1, 2, 6, 45, 1376, 161587, 64552473, 85987037645, 386933032425826, 6005080379837219319, 328011924848834642962619, 64153024576968812343635391868, 45547297603829979923254392040011994, 118654043008142499115765307533395739785599
Offset: 0

Views

Author

Gus Wiseman, Dec 08 2023

Keywords

Examples

			The a(0) = 1 through a(2) = 6 set-systems:
  {}  {}     {}
      {{1}}  {{1}}
             {{2}}
             {{1},{2}}
             {{1},{1,2}}
             {{2},{1,2}}
		

Crossrefs

The covering case is A054780.
For graphs we have A367862, covering A367863, unlabeled A006649.
These set-systems have ranks A367917.
A000372 counts antichains, covering A006126, nonempty A014466.
A003465 counts set-systems covering {1..n}, unlabeled A055621.
A058891 counts set-systems, unlabeled A000612.
A059201 counts covering T_0 set-systems.
A136556 counts set-systems on {1..n} with n edges.

Programs

  • Mathematica
    Table[Length[Select[Subsets[Rest[Subsets[Range[n]]]], Length[Union@@#]==Length[#]&]],{n,0,3}]
  • PARI
    \\ Here b(n) is A054780(n).
    b(n) = sum(k=0, n, (-1)^(n-k) * binomial(n,k) * binomial(2^k-1, n))
    a(n) = sum(k=0, n, binomial(n,k) * b(k)) \\ Andrew Howroyd, Dec 29 2023

Formula

Binomial transform of A054780.

A368598 Number of non-isomorphic n-element sets of singletons or pairs of elements of {1..n}, or unlabeled loop-graphs with n edges and up to n vertices.

Original entry on oeis.org

1, 1, 2, 6, 17, 52, 173, 585, 2064, 7520, 28265, 109501, 437394, 1799843, 7629463, 33302834, 149633151, 691702799, 3287804961, 16058229900, 80533510224, 414384339438, 2185878202630, 11811050484851, 65318772618624, 369428031895444, 2135166786135671, 12601624505404858
Offset: 0

Views

Author

Gus Wiseman, Jan 05 2024

Keywords

Comments

It doesn't matter for this sequence whether we use loops such as {x,x} or half-loops such as {x}.

Examples

			Non-isomorphic representatives of the a(0) = 1 through a(4) = 17 set-systems:
  {}  {{1}}  {{1},{2}}    {{1},{2},{3}}        {{1},{2},{3},{4}}
             {{1},{1,2}}  {{1},{2},{1,2}}      {{1},{2},{3},{1,2}}
                          {{1},{2},{1,3}}      {{1},{2},{3},{1,4}}
                          {{1},{1,2},{1,3}}    {{1},{2},{1,2},{1,3}}
                          {{1},{1,2},{2,3}}    {{1},{2},{1,2},{3,4}}
                          {{1,2},{1,3},{2,3}}  {{1},{2},{1,3},{1,4}}
                                               {{1},{2},{1,3},{2,3}}
                                               {{1},{2},{1,3},{2,4}}
                                               {{1},{3},{1,2},{2,4}}
                                               {{1},{1,2},{1,3},{1,4}}
                                               {{1},{1,2},{1,3},{2,3}}
                                               {{1},{1,2},{1,3},{2,4}}
                                               {{1},{1,2},{2,3},{3,4}}
                                               {{2},{1,2},{1,3},{1,4}}
                                               {{4},{1,2},{1,3},{2,3}}
                                               {{1,2},{1,3},{1,4},{2,3}}
                                               {{1,2},{1,3},{2,4},{3,4}}
		

Crossrefs

For any number of edges of any size we have A000612, covering A055621.
For any number of edges we have A000666, A054921, A322700.
The labeled version is A014068.
Counting by weight gives A320663, or A339888 with loops {x,x}.
The covering case is A368599.
For edges of any size we have A368731, covering A368186.
Row sums of A368836.
A000085 counts set partitions into singletons or pairs.
A001515 counts length-n set partitions into singletons or pairs.
A100861 counts set partitions into singletons or pairs by number of pairs.
A111924 counts set partitions into singletons or pairs by length.

Programs

  • Mathematica
    brute[m_]:=First[Sort[Table[Sort[Sort /@ (m/.Rule@@@Table[{(Union@@m)[[i]],p[[i]]}, {i,Length[p]}])], {p,Permutations[Range[Length[Union@@m]]]}]]];
    Table[Length[Union[brute /@ Subsets[Subsets[Range[n],{1,2}],{n}]]],{n,0,5}]
  • PARI
    a(n) = polcoef(G(n, O(x*x^n)), n) \\ G defined in A070166. - Andrew Howroyd, Jan 09 2024

Formula

a(n) = A070166(n, n). - Andrew Howroyd, Jan 09 2024

Extensions

Terms a(7) and beyond from Andrew Howroyd, Jan 09 2024

A368835 Number of unlabeled n-edge loop-graphs with at most n vertices such that it is not possible to choose a different vertex from each edge.

Original entry on oeis.org

0, 0, 0, 1, 5, 23, 98, 394, 1560, 6181, 24655, 99701, 410513, 1725725, 7423757, 32729320, 148027044, 687188969, 3275077017, 16022239940, 80431483586, 414094461610, 2185052929580, 11808696690600, 65312048149993, 369408792148714, 2135111662435080, 12601466371445619
Offset: 0

Views

Author

Gus Wiseman, Jan 13 2024

Keywords

Examples

			Non-isomorphic representatives of the a(4) = 5 loop-graphs:
  {{1,1},{2,2},{3,3},{1,2}}
  {{1,1},{2,2},{1,2},{1,3}}
  {{1,1},{2,2},{1,2},{3,4}}
  {{1,1},{2,2},{1,3},{2,3}}
  {{1,1},{1,2},{1,3},{2,3}}
		

Crossrefs

The case of a unique choice is A000081, row sums of A106234.
The labeled version is A368596, covering A368730.
Without the choice condition we have A368598.
The complement is A368984, row sums of A368926.
A000085, A100861, A111924 count set partitions into singletons or pairs.
A006125 counts graphs, unlabeled A000088.
A006129 counts covering graphs, unlabeled A002494.
A014068 counts loop-graphs, unlabeled A000666.
A058891 counts set-systems (without singletons A016031), unlabeled A000612.
A322661 counts labeled covering half-loop-graphs, connected A062740.

Programs

  • Mathematica
    Table[Length[Union[sysnorm /@ Select[Subsets[Subsets[Range[n],{1,2}],{n}],Select[Tuples[#], UnsameQ@@#&]=={}&]]],{n,0,5}]

Formula

a(n) = A368598(n) - A368984(n). - Andrew Howroyd, Jan 14 2024

Extensions

a(8) onwards from Andrew Howroyd, Jan 14 2024

A370638 Number of subsets of {1..n} such that a unique set can be obtained by choosing a different binary index of each element.

Original entry on oeis.org

1, 2, 4, 6, 12, 19, 30, 45, 90, 147, 230, 343, 504, 716, 994, 1352, 2704, 4349, 6469, 9162, 12585, 16862, 22122, 28617, 36653, 46431, 58075, 72097, 88456, 107966, 130742, 157647, 315294, 494967, 704753, 950080, 1234301, 1565165, 1945681, 2387060, 2890368, 3470798
Offset: 0

Views

Author

Gus Wiseman, Mar 09 2024

Keywords

Comments

A binary index of n is any position of a 1 in its reversed binary expansion. The binary indices of n are row n of A048793.

Examples

			The set {3,4} has binary indices {{1,2},{3}}, with two choices {1,3}, {2,3}, so is not counted under a(4).
The a(0) = 1 through a(5) = 19 subsets:
  {}  {}   {}     {}     {}       {}
      {1}  {1}    {1}    {1}      {1}
           {2}    {2}    {2}      {2}
           {1,2}  {1,2}  {4}      {4}
                  {1,3}  {1,2}    {1,2}
                  {2,3}  {1,3}    {1,3}
                         {1,4}    {1,4}
                         {2,3}    {1,5}
                         {2,4}    {2,3}
                         {1,2,4}  {2,4}
                         {1,3,4}  {4,5}
                         {2,3,4}  {1,2,4}
                                  {1,2,5}
                                  {1,3,4}
                                  {1,3,5}
                                  {2,3,4}
                                  {2,3,5}
                                  {2,4,5}
                                  {3,4,5}
		

Crossrefs

Set systems of this type are counted by A367904, ranks A367908.
A version for MM-numbers of multisets is A368101.
For prime indices we have A370584.
This is the unique version of A370636, complement A370637.
The maximal case is A370640, differences A370641.
Factorizations of this type are counted by A370645.
The case A370818 is the restriction to A000225.
A048793 lists binary indices, A000120 length, A272020 reverse, A029931 sum.
A058891 counts set-systems, A003465 covering, A323818 connected.
A070939 gives length of binary expansion.
A096111 gives product of binary indices.

Programs

  • Mathematica
    bpe[n_]:=Join@@Position[Reverse[IntegerDigits[n,2]],1];
    Table[Length[Select[Subsets[Range[n]],Length[Union[Sort /@ Select[Tuples[bpe/@#],UnsameQ@@#&]]]==1&]],{n,0,10}]

Formula

a(2^n - 1) = A370818(n).

Extensions

More terms from Jinyuan Wang, Mar 28 2025

A368599 Number of non-isomorphic n-element sets of singletons or pairs of elements of {1..n} with union {1..n}, or unlabeled loop-graphs with n edges covering n vertices.

Original entry on oeis.org

1, 1, 2, 5, 13, 34, 97, 277, 825, 2486, 7643, 23772, 74989, 238933, 769488, 2500758, 8199828, 27106647, 90316944, 303182461, 1025139840, 3490606305, 11967066094, 41302863014, 143493606215, 501772078429, 1765928732426, 6254738346969, 22294413256484, 79968425399831
Offset: 0

Views

Author

Gus Wiseman, Jan 06 2024

Keywords

Comments

It doesn't matter for this sequence whether we use loops such as {x,x} or half-loops such as {x}.

Examples

			The a(0) = 1 through a(4) = 13 set-systems:
  {}  {{1}}  {{1},{2}}    {{1},{2},{3}}        {{1},{2},{3},{4}}
             {{1},{1,2}}  {{1},{2},{1,3}}      {{1},{2},{3},{1,4}}
                          {{1},{1,2},{1,3}}    {{1},{2},{1,2},{3,4}}
                          {{1},{1,2},{2,3}}    {{1},{2},{1,3},{1,4}}
                          {{1,2},{1,3},{2,3}}  {{1},{2},{1,3},{2,4}}
                                               {{1},{2},{1,3},{3,4}}
                                               {{1},{1,2},{1,3},{1,4}}
                                               {{1},{1,2},{1,3},{2,4}}
                                               {{1},{1,2},{2,3},{2,4}}
                                               {{1},{1,2},{2,3},{3,4}}
                                               {{1},{2,3},{2,4},{3,4}}
                                               {{1,2},{1,3},{1,4},{2,3}}
                                               {{1,2},{1,3},{2,4},{3,4}}
		

Crossrefs

For any number of edges we have A000666, A054921, A322700.
For any number of edges of any size we have A055621, non-covering A000612.
For edges of any size we have A368186, covering case of A368731.
The labeled version is A368597, covering case of A014068.
This is the covering case of A368598.
A000085 counts set partitions into singletons or pairs.
A001515 counts length-n set partitions into singletons or pairs.
A100861 counts set partitions into singletons or pairs by number of pairs.
A111924 counts set partitions into singletons or pairs by length.

Programs

  • Mathematica
    brute[m_]:=First[Sort[Table[Sort[Sort /@ (m/.Rule@@@Table[{(Union@@m)[[i]],p[[i]]}, {i,Length[p]}])], {p,Permutations[Range[Length[Union@@m]]]}]]];
    Table[Length[Union[brute /@ Select[Subsets[Subsets[Range[n],{1,2}],{n}], Union@@#==Range[n]&]]],{n,0,5}]
  • PARI
    a(n) = polcoef(G(n, O(x*x^n)) - if(n, G(n-1, O(x*x^n))), n) \\ G defined in A070166. - Andrew Howroyd, Jan 09 2024

Formula

a(n) = A070166(n,n) - A070166(n-1,n) for n > 0. - Andrew Howroyd, Jan 09 2024

Extensions

Terms a(7) and beyond from Andrew Howroyd, Jan 09 2024

A054780 Number of n-covers of a labeled n-set.

Original entry on oeis.org

1, 1, 3, 32, 1225, 155106, 63602770, 85538516963, 386246934638991, 6001601072676524540, 327951891446717800997416, 64149416776011080449232990868, 45546527789182522411309599498741023, 118653450898277491435912500458608964207578
Offset: 0

Views

Author

Vladeta Jovovic, May 21 2000

Keywords

Comments

Also, number of n X n rational {0,1}-matrices with no zero rows or columns and with all rows distinct, up to permutation of rows.

Examples

			From _Gus Wiseman_, Dec 19 2023: (Start)
Number of ways to choose n nonempty sets with union {1..n}. For example, the a(3) = 32 covers are:
  {1}{2}{3}  {1}{2}{13}  {1}{2}{123}  {1}{12}{123}  {12}{13}{123}
             {1}{2}{23}  {1}{3}{123}  {1}{13}{123}  {12}{23}{123}
             {1}{3}{12}  {1}{12}{13}  {1}{23}{123}  {13}{23}{123}
             {1}{3}{23}  {1}{12}{23}  {2}{12}{123}
             {2}{3}{12}  {1}{13}{23}  {2}{13}{123}
             {2}{3}{13}  {2}{3}{123}  {2}{23}{123}
                         {2}{12}{13}  {3}{12}{123}
                         {2}{12}{23}  {3}{13}{123}
                         {2}{13}{23}  {3}{23}{123}
                         {3}{12}{13}  {12}{13}{23}
                         {3}{12}{23}
                         {3}{13}{23}
(End)
		

Crossrefs

Main diagonal of A055154.
Covers with any number of edges are counted by A003465, unlabeled A055621.
Connected graphs of this type are counted by A057500, unlabeled A001429.
This is the covering case of A136556.
The case of graphs is A367863, covering case of A116508, unlabeled A006649.
Binomial transform is A367916.
These set-systems have ranks A367917.
The unlabeled version is A368186.
A006129 counts covering graphs, connected A001187, unlabeled A002494.
A046165 counts minimal covers, ranks A309326.

Programs

  • Mathematica
    Join[{1}, Table[Sum[StirlingS1[n+1, k+1]*(2^k - 1)^n, {k, 0, n}]/n!, {n, 1, 15}]] (* Vaclav Kotesovec, Jun 04 2022 *)
    Table[Length[Select[Subsets[Rest[Subsets[Range[n]]],{n}],Union@@#==Range[n]&]],{n,0,4}] (* Gus Wiseman, Dec 19 2023 *)
  • PARI
    a(n) = sum(k=0, n, (-1)^k*binomial(n, k)*binomial(2^(n-k)-1, n)) \\ Andrew Howroyd, Jan 20 2024

Formula

a(n) = Sum_{k=0..n} (-1)^k*binomial(n, k)*binomial(2^(n-k)-1, n).
a(n) = (1/n!)*Sum_{k=0..n} Stirling1(n+1, k+1)*(2^k-1)^n.
G.f.: Sum_{n>=0} log(1+(2^n-1)*x)^n/((1+(2^n-1)*x)*n!). - Paul D. Hanna and Vladeta Jovovic, Jan 16 2008
a(n) ~ 2^(n^2) / n!. - Vaclav Kotesovec, Jun 04 2022
Inverse binomial transform of A367916. - Gus Wiseman, Dec 19 2023
Previous Showing 61-70 of 182 results. Next