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

A003182 Dedekind numbers: inequivalent monotone Boolean functions of n or fewer variables, or antichains of subsets of an n-set.

Original entry on oeis.org

2, 3, 5, 10, 30, 210, 16353, 490013148, 1392195548889993358, 789204635842035040527740846300252680
Offset: 0

Views

Author

Keywords

Comments

NP-equivalence classes of unate Boolean functions of n or fewer variables.
Also the number of simple games with n players in minimal winning form up to isomorphism. - Fabián Riquelme, Mar 13 2018
The labeled case is A000372. - Gus Wiseman, Feb 23 2019
First differs from A348260(n + 1) at a(5) = 210, A348260(6) = 233. - Gus Wiseman, Nov 28 2021
Pawelski & Szepietowski show that a(n) = A001206(n) (mod 2) and that a(9) = 6 (mod 210). - Charles R Greathouse IV, Feb 16 2023

Examples

			From _Gus Wiseman_, Feb 20 2019: (Start)
Non-isomorphic representatives of the a(0) = 2 through a(3) = 10 antichains:
  {}    {}     {}         {}
  {{}}  {{}}   {{}}       {{}}
        {{1}}  {{1}}      {{1}}
               {{1,2}}    {{1,2}}
               {{1},{2}}  {{1},{2}}
                          {{1,2,3}}
                          {{1},{2,3}}
                          {{1},{2},{3}}
                          {{1,3},{2,3}}
                          {{1,2},{1,3},{2,3}}
(End)
		

References

  • I. Anderson, Combinatorics of Finite Sets. Oxford Univ. Press, 1987, p. 38.
  • Arocha, Jorge Luis (1987) "Antichains in ordered sets" [ In Spanish ]. Anales del Instituto de Matematicas de la Universidad Nacional Autonoma de Mexico 27: 1-21.
  • J. Berman, Free spectra of 3-element algebras, in R. S. Freese and O. C. Garcia, editors, Universal Algebra and Lattice Theory (Puebla, 1982), Lect. Notes Math. Vol. 1004, 1983.
  • G. Birkhoff, Lattice Theory. American Mathematical Society, Colloquium Publications, Vol. 25, 3rd ed., Providence, RI, 1967, p. 63.
  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 273.
  • M. A. Harrison, Introduction to Switching and Automata Theory. McGraw Hill, NY, 1965, p. 188.
  • D. E. Knuth, The Art of Computer Programming, Vol. 4A, Section 7.1.1, p. 79.
  • W. F. Lunnon, The IU function: the size of a free distributive lattice, pp. 173-181 of D. J. A. Welsh, editor, Combinatorial Mathematics and Its Applications. Academic Press, NY, 1971.
  • Saburo Muroga, Threshold Logic and Its Applications. Wiley, NY, 1971, p. 38, Table 2.3.2. - Row 13.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • D. H. Wiedemann, personal communication.

Crossrefs

Formula

a(n) = A306505(n) + 1. - Gus Wiseman, Jul 02 2019

Extensions

a(7) added by Timothy Yusun, Sep 27 2012
a(8) from Pawelski added by Michel Marcus, Sep 01 2021
a(9) from Pawelski added by Michel Marcus, May 11 2023

A009998 Triangle in which j-th entry in i-th row is (j+1)^(i-j).

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 4, 3, 1, 1, 8, 9, 4, 1, 1, 16, 27, 16, 5, 1, 1, 32, 81, 64, 25, 6, 1, 1, 64, 243, 256, 125, 36, 7, 1, 1, 128, 729, 1024, 625, 216, 49, 8, 1, 1, 256, 2187, 4096, 3125, 1296, 343, 64, 9, 1, 1, 512, 6561, 16384, 15625, 7776, 2401, 512, 81, 10, 1
Offset: 0

Views

Author

Keywords

Comments

Read as a square array this is the Hilbert transform of triangle A123125 (see A145905 for the definition of this term). For example, the fourth row of A123125 is (0,1,4,1) and the expansion (x + 4*x^2 + x^3)/(1-x)^4 = x + 8*x^2 + 27*x^3 + 64*x^4 + ... generates the entries in the fourth row of this array read as a square. - Peter Bala, Oct 28 2008

Examples

			Triangle begins:
  1;
  1,  1;
  1,  2,  1;
  1,  4,  3,  1;
  1,  8,  9,  4,  1;
  1, 16, 27, 16,  5,  1;
  1, 32, 81, 64, 25,  6,  1;
  ...
From _Gus Wiseman_, May 01 2021: (Start)
The rows of the triangle are obtained by reading antidiagonals upward in the following table of A(k,n) = n^k, with offset k = 0, n = 1:
         n=1:     n=2:     n=3:     n=4:     n=5:     n=6:
   k=0:   1        1        1        1        1        1
   k=1:   1        2        3        4        5        6
   k=2:   1        4        9       16       25       36
   k=3:   1        8       27       64      125      216
   k=4:   1       16       81      256      625     1296
   k=5:   1       32      243     1024     3125     7776
   k=6:   1       64      729     4096    15625    46656
   k=7:   1      128     2187    16384    78125   279936
   k=8:   1      256     6561    65536   390625  1679616
   k=9:   1      512    19683   262144  1953125 10077696
  k=10:   1     1024    59049  1048576  9765625 60466176
(End)
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 24.

Crossrefs

Row sums give A026898.
Column n = 2 of the array is A000079.
Column n = 3 of the array is A000244.
Row k = 2 of the array is A000290.
Row k = 3 of the array is A000578.
Diagonal n = k of the array is A000312.
Diagonal n = k + 1 of the array is A000169.
Diagonal n = k + 2 of the array is A000272.
The transpose of the array is A009999.
The numbers of divisors of the entries are A343656 (row sums: A343657).
A007318 counts k-sets of elements of {1..n}.
A059481 counts k-multisets of elements of {1..n}.

Programs

  • Haskell
    a009998 n k = (k + 1) ^ (n - k)
    a009998_row n = a009998_tabl !! n
    a009998_tabl = map reverse a009999_tabl
    -- Reinhard Zumkeller, Feb 02 2014
    
  • Maple
    E := (n,x) -> `if`(n=0,1,x*(1-x)*diff(E(n-1,x),x)+E(n-1,x)*(1+(n-1)*x));
    G := (n,x) -> E(n,x)/(1-x)^(n+1);
    A009998 := (n,k) -> coeff(series(G(n-k,x),x,18),x,k);
    seq(print(seq(A009998(n,k),k=0..n)),n=0..6);
    # Peter Luschny, Aug 02 2010
  • Mathematica
    Flatten[Table[(j+1)^(i-j),{i,0,20},{j,0,i}]] (* Harvey P. Dale, Dec 25 2012 *)
  • PARI
    T(i,j)=(j+1)^(i-j) \\ Charles R Greathouse IV, Feb 06 2017

Formula

T(n,n) = 1; T(n,k) = (k+1)*T(n-1,k) for k=0..n-1. - Reinhard Zumkeller, Feb 02 2014
T(n,m) = (m+1)*Sum_{k=0..n-m}((n+1)^(k-1)*(n-m)^(n-m-k)*(-1)^(n-m-k)*binomial(n-m-1,k-1)). - Vladimir Kruchinin, Sep 12 2015

Extensions

a(62) corrected to 512 by T. D. Noe, Dec 20 2007

A305001 Number of labeled antichains of finite sets spanning n vertices without singletons.

Original entry on oeis.org

1, 0, 1, 5, 87, 6398, 7745253, 2414573042063, 56130437190053518791691, 286386577668298410118121281898931424413687
Offset: 0

Views

Author

Gus Wiseman, May 23 2018

Keywords

Comments

From Gus Wiseman, Jul 03 2019: (Start)
Also the number of antichains covering n vertices and having empty intersection (meaning there is no vertex in common to all the edges). For example, the a(3) = 5 antichains are:
{{3},{1,2}}
{{2},{1,3}}
{{1},{2,3}}
{{1},{2},{3}}
{{1,2},{1,3},{2,3}}
(End)

Examples

			The a(3) = 5 antichains:
  {{1,2,3}}
  {{1,2},{1,3}}
  {{1,2},{2,3}}
  {{1,3},{2,3}}
  {{1,2},{1,3},{2,3}}
		

Crossrefs

The binomial transform is the non-covering case A307249.
The second binomial transform is A014466.

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]]]];
    Table[Length[Select[stableSets[Subsets[Range[n],{1,n}],SubsetQ],And[Union@@#==Range[n],#=={}||Intersection@@#=={}]&]],{n,0,5}] (* Gus Wiseman, Jul 03 2019 *)

Extensions

a(9) from A307249 - Dmitry I. Ignatov, Nov 27 2023

A046165 Number of minimal covers of n objects.

Original entry on oeis.org

1, 1, 2, 8, 49, 462, 6424, 129425, 3731508, 152424420, 8780782707, 710389021036, 80610570275140, 12815915627480695, 2855758994821922882, 892194474524889501292, 391202163933291014701953, 240943718535427829240708786, 208683398342300491409959279244
Offset: 0

Views

Author

Keywords

Comments

No edge of a minimal cover can be a subset of any other, so minimal covers are antichains, but the converse is not true. - Gus Wiseman, Jul 03 2019
a(n) is the number of undirected graphs on n nodes for which the intersection number and independence number are equal. See Proposition 2.3.7 and Theorem 2.3.3 of the Deligeorgaki et al. paper below. - Alex Markham, Oct 13 2022

Examples

			From _Gus Wiseman_, Jul 02 2019: (Start)
The a(1) = 1 through a(3) = 8 minimal covers:
  {{1}}  {{1,2}}    {{1,2,3}}
         {{1},{2}}  {{1},{2,3}}
                    {{2},{1,3}}
                    {{3},{1,2}}
                    {{1,2},{1,3}}
                    {{1,2},{2,3}}
                    {{1},{2},{3}}
                    {{1,3},{2,3}}
(End)
		

Crossrefs

Antichain covers are A006126.
Minimal covering simple graphs are A053530.
Maximal antichains are A326358.
Row sums of A035347 or of A282575.

Programs

  • Maple
    a:= n-> add(add((-1)^i* binomial(k,i) *(2^k-1-i)^n, i=0..k)/k!, k=0..n):
    seq(a(n), n=0..20);  # Alois P. Heinz, Aug 19 2008
  • Mathematica
    Table[Sum[Sum[Binomial[n,i]StirlingS2[i,k](2^k-k-1)^(n-i),{i,k,n}],{k,2,n}]+1,{n,1,20}] (* Geoffrey Critzer, Jun 27 2013 *)

Formula

E.g.f.: Sum_{n>=0} (exp(x)-1)^n*exp(x*(2^n-n-1))/n!. - Vladeta Jovovic, May 08 2004
a(n) = Sum_{k=1..n} Sum_{i=k..n} C(n,i)*Stirling2(i,k)*(2^k - k - 1)^(n - i). - Geoffrey Critzer, Jun 27 2013
a(n) ~ c * 2^(n^2/4 + n + 1/2) / sqrt(Pi*n), where c = JacobiTheta3(0,1/2) = EllipticTheta[3, 0, 1/2] = 2.1289368272118771586694585485449... if n is even, and c = JacobiTheta2(0,1/2) = EllipticTheta[2, 0, 1/2] = 2.1289312505130275585916134025753... if n is odd. - Vaclav Kotesovec, Mar 10 2014

Extensions

a(0)=1 prepended by Alois P. Heinz, Feb 18 2017

A343656 Array read by antidiagonals where A(n,k) is the number of divisors of n^k.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 3, 2, 1, 1, 4, 3, 3, 1, 1, 5, 4, 5, 2, 1, 1, 6, 5, 7, 3, 4, 1, 1, 7, 6, 9, 4, 9, 2, 1, 1, 8, 7, 11, 5, 16, 3, 4, 1, 1, 9, 8, 13, 6, 25, 4, 7, 3, 1, 1, 10, 9, 15, 7, 36, 5, 10, 5, 4, 1, 1, 11, 10, 17, 8, 49, 6, 13, 7, 9, 2, 1, 1, 12, 11, 19, 9, 64, 7, 16, 9, 16, 3, 6, 1
Offset: 1

Views

Author

Gus Wiseman, Apr 28 2021

Keywords

Comments

First differs from A343658 at A(4,2) = 5, A343658(4,2) = 6.
As a triangle, T(n,k) = number of divisors of k^(n-k).

Examples

			Array begins:
       k=0 k=1 k=2 k=3 k=4 k=5 k=6 k=7
  n=1:  1   1   1   1   1   1   1   1
  n=2:  1   2   3   4   5   6   7   8
  n=3:  1   2   3   4   5   6   7   8
  n=4:  1   3   5   7   9  11  13  15
  n=5:  1   2   3   4   5   6   7   8
  n=6:  1   4   9  16  25  36  49  64
  n=7:  1   2   3   4   5   6   7   8
  n=8:  1   4   7  10  13  16  19  22
  n=9:  1   3   5   7   9  11  13  15
Triangle begins:
  1
  1  1
  1  2  1
  1  3  2  1
  1  4  3  3  1
  1  5  4  5  2  1
  1  6  5  7  3  4  1
  1  7  6  9  4  9  2  1
  1  8  7 11  5 16  3  4  1
  1  9  8 13  6 25  4  7  3  1
  1 10  9 15  7 36  5 10  5  4  1
  1 11 10 17  8 49  6 13  7  9  2  1
  1 12 11 19  9 64  7 16  9 16  3  6  1
  1 13 12 21 10 81  8 19 11 25  4 15  2  1
For example, row n = 8 counts the following divisors:
  1  64  243  256  125  36  7  1
     32  81   128  25   18  1
     16  27   64   5    12
     8   9    32   1    9
     4   3    16        6
     2   1    8         4
     1        4         3
              2         2
              1         1
		

Crossrefs

Columns k=1..9 of the array give A000005, A048691, A048785, A344327, A344328, A344329, A343526, A344335, A344336.
Row n = 6 of the array is A000290.
Diagonal n = k of the array is A062319.
Array antidiagonal sums (row sums of the triangle) are A343657.
Dominated by A343658.
A000312 = n^n.
A007318 counts k-sets of elements of {1..n}.
A009998(n,k) = n^k (as an array, offset 1).
A059481 counts k-multisets of elements of {1..n}.

Programs

  • Mathematica
    Table[DivisorSigma[0,k^(n-k)],{n,10},{k,n}]
  • PARI
    A(n, k) = numdiv(n^k); \\ Seiichi Manyama, May 15 2021

Formula

A(n,k) = A000005(A009998(n,k)), where A009998(n,k) = n^k is the interpretation as an array.
A(n,k) = Sum_{d|n} k^omega(d). - Seiichi Manyama, May 15 2021

A007363 Maximal self-dual antichains on n points.

Original entry on oeis.org

0, 1, 3, 5, 20, 168, 11748, 12160647
Offset: 1

Views

Author

Keywords

Comments

From Gus Wiseman, Jul 02 2019: (Start)
If self-dual means (pairwise) intersecting, then a(n) is the number of maximal intersecting antichains of nonempty subsets of {1..(n - 1)}. A set of sets is an antichain if no part is a subset of any other, and is intersecting if no two parts are disjoint. For example, the a(2) = 1 through a(5) = 20 maximal intersecting antichains are:
{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}
(End)

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Intersecting antichains are A326372.
Intersecting antichains of nonempty sets are A001206.
Unlabeled intersecting antichains are A305857.
Maximal antichains of nonempty sets are A326359.
The case with empty edges allowed is A326363.

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],{1,n}],Or[Intersection[#1,#2]=={},SubsetQ[#1,#2]]&]]],{n,0,5}] (* Gus Wiseman, Jul 02 2019 *)
    (* 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]-1 (* Elijah Beregovsky, May 06 2020 *)

Formula

For n > 0, a(n) = A326363(n - 1) - 1 = A326362(n - 1) + n - 1. - Gus Wiseman, Jul 03 2019

Extensions

a(8) from Elijah Beregovsky, May 06 2020

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

A326359 Number of maximal antichains of nonempty subsets of {1..n}.

Original entry on oeis.org

1, 1, 2, 6, 28, 375, 31745, 123805913
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.

Examples

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

Crossrefs

Antichains of nonempty sets are A014466.
Minimal covering antichains are A046165.
Maximal intersecting antichains are A007363.
Maximal antichains of sets are A326358.

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],{1,n}],SubsetQ]]],{n,0,5}]

Formula

For n > 0, a(n) = A326358(n) - 1.

Extensions

a(6) from Andrew Howroyd, Aug 14 2019
a(7) from Dmitry I. Ignatov, Oct 12 2021

A326361 Number of maximal intersecting antichains of sets covering n vertices with no singletons.

Original entry on oeis.org

1, 1, 1, 2, 12, 133, 11386, 12143511
Offset: 0

Views

Author

Gus Wiseman, Jul 01 2019

Keywords

Comments

Covering means there are no isolated vertices. A set system (set of sets) is an antichain if no part is a subset of any other, and is intersecting if no two parts are disjoint.

Examples

			The a(4) = 12 antichains:
  {{1,2,3,4}}
  {{1,2},{1,3,4},{2,3,4}}
  {{1,3},{1,2,4},{2,3,4}}
  {{1,4},{1,2,3},{2,3,4}}
  {{2,3},{1,2,4},{1,3,4}}
  {{2,4},{1,2,3},{1,3,4}}
  {{3,4},{1,2,3},{1,2,4}}
  {{1,2},{1,3},{1,4},{2,3,4}}
  {{1,2},{2,3},{2,4},{1,3,4}}
  {{1,3},{2,3},{3,4},{1,2,4}}
  {{1,4},{2,4},{3,4},{1,2,3}}
  {{1,2,3},{1,2,4},{1,3,4},{2,3,4}}
		

Crossrefs

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[Select[stableSets[Subsets[Range[n]],Or[Intersection[#1,#2]=={},SubsetQ[#1,#2]]&],Union@@#==Range[n]&]]],{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 = Select[FindClique[g, Infinity, All], BitOr @@ # == n - 1 &];
    Length[sets] (* Elijah Beregovsky, May 05 2020 *)

Extensions

a(6)-a(7) from Elijah Beregovsky, May 05 2020

A326362 Number of maximal intersecting antichains of nonempty, non-singleton subsets of {1..n}.

Original entry on oeis.org

1, 1, 1, 2, 16, 163, 11742, 12160640
Offset: 0

Views

Author

Gus Wiseman, Jul 01 2019

Keywords

Comments

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

Examples

			The a(4) = 16 maximal intersecting antichains:
  {{1,2,3,4}}
  {{1,2},{1,3},{2,3}}
  {{1,2},{1,4},{2,4}}
  {{1,3},{1,4},{3,4}}
  {{2,3},{2,4},{3,4}}
  {{1,2},{1,3,4},{2,3,4}}
  {{1,3},{1,2,4},{2,3,4}}
  {{1,4},{1,2,3},{2,3,4}}
  {{2,3},{1,2,4},{1,3,4}}
  {{2,4},{1,2,3},{1,3,4}}
  {{3,4},{1,2,3},{1,2,4}}
  {{1,2},{1,3},{1,4},{2,3,4}}
  {{1,2},{2,3},{2,4},{1,3,4}}
  {{1,3},{2,3},{3,4},{1,2,4}}
  {{1,4},{2,4},{3,4},{1,2,3}}
  {{1,2,3},{1,2,4},{1,3,4},{2,3,4}}
		

Crossrefs

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],{2,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]-Log[2,n]-1 (* Elijah Beregovsky, May 06 2020 *)

Formula

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

Extensions

a(7) from Elijah Beregovsky, May 06 2020
Showing 1-10 of 18 results. Next