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 85 results. Next

A327519 Number of factorizations of A305078(n - 1), the n-th connected number, into connected numbers > 1.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 1, 1, 2, 1, 2, 3, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 1, 4, 2, 1, 1, 1, 1, 5, 1, 2, 1, 2, 1, 1, 1, 1, 1, 2, 1, 2, 4, 2, 3, 1, 2, 1, 2, 1, 1, 4, 1, 1, 1, 2, 1, 1, 2, 4, 1, 1, 1, 2, 2, 7, 1, 1, 4, 1, 1, 2, 1, 2, 1, 1, 1, 1, 2, 2, 1, 1, 7, 2, 1
Offset: 1

Views

Author

Gus Wiseman, Sep 21 2019

Keywords

Comments

A number n with prime factorization n = prime(m_1)^s_1 * ... * prime(m_k)^s_k is connected if the simple labeled graph with vertex set {m_1,...,m_k} and edges between any two vertices with a common divisor greater than 1 is connected. Connected numbers are listed in A305078.

Examples

			The a(190) = 8 factorizations of 585 together with the corresponding multiset partitions of {2,2,3,6}:
  (3*3*5*13)  {{2},{2},{3},{6}}
  (3*3*65)    {{2},{2},{3,6}}
  (3*5*39)    {{2},{3},{2,6}}
  (3*195)     {{2},{2,3,6}}
  (5*9*13)    {{3},{2,2},{6}}
  (5*117)     {{3},{2,2,6}}
  (9*65)      {{2,2},{3,6}}
  (585)       {{2,2,3,6}}
		

Crossrefs

See link for additional cross-references.

Programs

  • Mathematica
    nn=100;
    zsm[s_]:=With[{c=Select[Subsets[Range[Length[s]],{2}],GCD@@s[[#]]>1&]},If[c=={},s,zsm[Sort[Append[Delete[s,List/@c[[1]]],LCM@@s[[c[[1]]]]]]]]];
    facsusing[s_,n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facsusing[Select[s,Divisible[n/d,#]&],n/d],Min@@#>=d&]],{d,Select[s,Divisible[n,#]&]}]];
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    y=Select[Range[nn],Length[zsm[primeMS[#]]]==1&];
    Table[Length[facsusing[y,n]],{n,y}]

A304716 Number of integer partitions of n whose distinct parts are connected.

Original entry on oeis.org

1, 2, 2, 3, 2, 5, 2, 6, 4, 9, 3, 15, 4, 18, 12, 25, 11, 41, 17, 54, 36, 72, 44, 113, 69, 145, 113, 204, 153, 302, 220, 394, 343, 541, 475, 771, 662, 1023, 968, 1398, 1314, 1929, 1822, 2566, 2565, 3440, 3446, 4677, 4688, 6187, 6407, 8216, 8544, 10975, 11436
Offset: 1

Views

Author

Gus Wiseman, May 17 2018

Keywords

Comments

Given a finite set S of positive integers greater than one, let G(S) be the simple labeled graph with vertex set S and edges between any two vertices with a common divisor. For example, G({6,14,15,35}) is a 4-cycle. A multiset S is said to be connected if G(S) is a connected graph.

Examples

			The a(12) = 15 connected integer partitions and their corresponding connected multiset multisystems (see A112798, A302242) are the following.
                     (12): {{1,1,2}}
                    (6 6): {{1,2},{1,2}}
                    (8 4): {{1,1,1},{1,1}}
                    (9 3): {{2,2},{2}}
                   (10 2): {{1,3},{1}}
                  (4 4 4): {{1,1},{1,1},{1,1}}
                  (6 3 3): {{1,2},{2},{2}}
                  (6 4 2): {{1,2},{1,1},{1}}
                  (8 2 2): {{1,1,1},{1},{1}}
                (3 3 3 3): {{2},{2},{2},{2}}
                (4 4 2 2): {{1,1},{1,1},{1},{1}}
                (6 2 2 2): {{1,2},{1},{1},{1}}
              (4 2 2 2 2): {{1,1},{1},{1},{1},{1}}
            (2 2 2 2 2 2): {{1},{1},{1},{1},{1},{1}}
(1 1 1 1 1 1 1 1 1 1 1 1): {{},{},{},{},{},{},{},{},{},{},{},{}}
		

Crossrefs

Programs

  • Mathematica
    zsm[s_]:=With[{c=Select[Tuples[Range[Length[s]],2],And[Less@@#,GCD@@s[[#]]]>1&]},If[c==={},s,zsm[Union[Append[Delete[s,List/@c[[1]]],LCM@@s[[c[[1]]]]]]]]];
    Table[Length[Select[IntegerPartitions[n],Length[zsm[Union[#]]]===1&]],{n,30}]

Formula

For n > 1, a(n) = A218970(n) + 1. - Gus Wiseman, Dec 04 2018

Extensions

Name changed to distinguish from A218970 by Gus Wiseman, Dec 04 2018

A305079 Number of connected components of the integer partition with Heinz number n.

Original entry on oeis.org

0, 1, 1, 2, 1, 2, 1, 3, 1, 2, 1, 3, 1, 2, 2, 4, 1, 2, 1, 3, 1, 2, 1, 4, 1, 2, 1, 3, 1, 3, 1, 5, 2, 2, 2, 3, 1, 2, 1, 4, 1, 2, 1, 3, 2, 2, 1, 5, 1, 2, 2, 3, 1, 2, 2, 4, 1, 2, 1, 4, 1, 2, 1, 6, 1, 3, 1, 3, 2, 3, 1, 4, 1, 2, 2, 3, 2, 2, 1, 5, 1, 2, 1, 3, 2, 2, 1
Offset: 1

Views

Author

Gus Wiseman, May 24 2018

Keywords

Comments

First differs from |A305052(n)| at a(169) = 1, A305052(169) = 0.
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
Given a finite multiset S of positive integers greater than one, let G(S) be the simple labeled graph with vertex set S and edges between any two vertices with a common divisor greater than 1. For example, G({6,14,15,35}) is a 4-cycle. If S is the integer partition with Heinz number n, a(n) is the number of connected components of G(S).

Examples

			The a(315) = 2 connected components of {2,2,3,4} are {{3},{2,2,4}}.
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n===1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    zsm[s_]:=With[{c=Select[Tuples[Range[Length[s]],2],And[Less@@#,GCD@@s[[#]]]>1&]},If[c=={},s,zsm[Sort[Append[Delete[s,List/@c[[1]]],LCM@@s[[c[[1]]]]]]]]];
    Table[Length[zsm[primeMS[n]]],{n,100}]
  • PARI
    zero_first_elem_and_connected_elems(ys) = { my(cs = List([ys[1]]), i=1); ys[1] = 0; while(i<=#cs, for(j=2,#ys,if(ys[j]&&(1!=gcd(cs[i],ys[j])), listput(cs,ys[j]); ys[j] = 0)); i++); (ys); };
    A007814(n) = valuation(n,2);
    A000265(n) = (n/2^A007814(n));
    A305079(n) = if(!(n%2),A007814(n)+A305079(A000265(n)), my(cs = apply(p -> primepi(p),factor(n)[,1]~), s=0); while(#cs, cs = select(c -> c, zero_first_elem_and_connected_elems(cs)); s++); (s)); \\ Antti Karttunen, Nov 10 2018

Formula

For all n, k > 0, we have a(2^n * k) = n + a(k).
For all x, y > 0, we have a(x * y) <= a(x) + a(y).
For x, y > 0 strongly coprime, we have a(x * y) = a(x) + a(y). Strongly coprime means every prime index of x is coprime to every prime index of y, where a prime index of n is a number m such that prime(m) divides n.
a(n) = A305501(A064989(n)) + A007814(n). - Antti Karttunen, Nov 10 2018

Extensions

Terms and Mathematica program corrected by Gus Wiseman, Nov 10 2018

A099947 Number of topologically connected set partitions of {1,...,n}.

Original entry on oeis.org

1, 1, 1, 1, 2, 6, 21, 85, 385, 1907, 10205, 58455, 355884, 2290536, 15518391, 110283179, 819675482, 6355429550, 51293023347, 430062712439, 3739408304962, 33665192703946, 313354708842791, 3011545611755271, 29847401178719637, 304713973031878687, 3201007359886598431
Offset: 0

Views

Author

N. J. A. Sloane, Nov 12 2004

Keywords

Comments

A set partition of {1,...,n} is topologically connected if the graph whose vertices are the blocks and whose edges are crossing pairs of blocks is connected, where two blocks cross each other if they are of the form {{...x...y...}, {...z...t...}} for some x < z < y < t or z < x < t < y. - Gus Wiseman, Feb 19 2019

Examples

			O.g.f.: A(x) = 1 + x + x^2 + x^3 + 2*x^4 + 6*x^5 + 21*x^6 + 85*x^7 +...
From _Paul D. Hanna_, Apr 16 2013: (Start)
The o.g.f. satisfies
(1) A(x) = 1 + x/A(x) + 2*x^2/A(x)^2 + 5*x^3/A(x)^3 + 15*x^4/A(x)^4 + 52*x^5/A(x)^5 + 203*x^6/A(x)^6 + ... + A000110(n)*x^n/A(x)^n + ...
(2) A(x) = 1 + x/(A(x)-x) + x^2/((A(x)-x)*(A(x)-2*x)) + x^3/((A(x)-x)*(A(x)-2*x)*(A(x)-3*x)) + x^4/((A(x)-x)*(A(x)-2*x)*(A(x)-3*x)*(A(x)-4*x)) + ... (End)
From _Gus Wiseman_, Feb 19 2019: (Start)
The a(1) = 1 through a(6) = 21 topologically connected set partitions:
  {{1}}  {{12}}  {{123}}  {{1234}}    {{12345}}    {{123456}}
                          {{13}{24}}  {{124}{35}}  {{1235}{46}}
                                      {{13}{245}}  {{124}{356}}
                                      {{134}{25}}  {{1245}{36}}
                                      {{135}{24}}  {{1246}{35}}
                                      {{14}{235}}  {{125}{346}}
                                                   {{13}{2456}}
                                                   {{134}{256}}
                                                   {{1345}{26}}
                                                   {{1346}{25}}
                                                   {{135}{246}}
                                                   {{1356}{24}}
                                                   {{136}{245}}
                                                   {{14}{2356}}
                                                   {{145}{236}}
                                                   {{146}{235}}
                                                   {{15}{2346}}
                                                   {{13}{25}{46}}
                                                   {{14}{25}{36}}
                                                   {{14}{26}{35}}
                                                   {{15}{24}{36}}
(End)
		

Crossrefs

Programs

  • Mathematica
    a[0] = 1; a[n_] := Module[{A = 1 + x}, For[i = 1, i <= n, i++, A = Sum[x^m/Product[A - k*x + x*O[x]^n, {k, 1, m}], {m, 0, n}]]; Coefficient[A, x^n]]; Table[a[n], {n, 0, 24}] (* Jean-François Alcover, Sep 13 2013, after Paul D. Hanna *)
    nn=8;
    nonXQ[stn_]:=!MatchQ[stn,{_,{_,x_,_,y_,_},_,{_,z_,_,t_,_},_}/;x_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    Solve[Table[BellB[n]==Sum[Product[a[Length[s]],{s,stn}],{stn,Select[sps[Range[n]],nonXQ]}],{n,nn}],Array[a,nn]] (* Gus Wiseman, Feb 19 2019 *)
  • PARI
    {a(n)=if(n<0, 0, polcoeff( x/serreverse(x*serlaplace(exp(exp(x+x*O(x^n))-1))), n))} /* Michael Somos, Sep 22 2005 */
    
  • PARI
    {a(n)=local(A=1+x); for(i=1, n, A=sum(m=0, n, x^m/prod(k=1, m, A - k*x +x*O(x^n)) )); polcoeff(A, n)} \\ Paul D. Hanna, Apr 16 2013

Formula

From Paul D. Hanna, Apr 16 2013: (Start)
O.g.f. A(x) satisfies
(1) A(x) = Sum_{n>=0} A000110(n)*x^n/A(x)^n, where A000110 are the Bell numbers.
(2) A(x) = Sum_{n>=0} x^n / Product_{k=1..n} (A(x) - k*x).
(3) A(x) = 1/(1 - x/(A(x) - 1*x/(1 - x/(A(x) - 2*x/(1 - x/(A(x) - 3*x/(1 - x/(A(x) - 4*x/(1 - x/(A(x) - ... )))))))))), a continued fraction. (End)
B(n) = Sum_p Product_{s in p} a(|s|) where p is a non-crossing set partition of {1,...,n} and B = A000110. In words, every set partition of {1,...,n} can be uniquely decomposed as a non-crossing set partition together with a topologically connected set partition of each block. - Gus Wiseman, Feb 19 2019

Extensions

Name edited by Gus Wiseman, Feb 19 2019

A326753 Number of connected components of the set-system with BII-number n.

Original entry on oeis.org

0, 1, 1, 2, 1, 1, 1, 1, 1, 2, 2, 3, 2, 2, 2, 2, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1
Offset: 0

Views

Author

Gus Wiseman, Jul 23 2019

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. We define the set-system with BII-number n to be obtained by taking the binary indices of each binary index of n. Every finite set of finite nonempty sets has a different BII-number. For example, 18 has reversed binary expansion (0,1,0,0,1), and since the binary indices of 2 and 5 are {2} and {1,3} respectively, the BII-number of {{2},{1,3}} is 18. Elements of a set-system are sometimes called edges.

Examples

			The set-system {{1,2},{1,4},{3}} with BII-number 268 has two connected components, so a(268) = 2.
		

Crossrefs

Positions of 0's and 1's are A326749.
Ranking sequences using BII-numbers: A309314 (hyperforests), A326701 (set partitions), A326703 (chains), A326704 (antichains), A326750 (clutters), A326751 (blobs), A326752 (hypertrees), A326754 (covers).

Programs

  • Mathematica
    bpe[n_]:=Join@@Position[Reverse[IntegerDigits[n,2]],1];
    csm[s_]:=With[{c=Select[Tuples[Range[Length[s]],2],And[OrderedQ[#],UnsameQ@@#,Length[Intersection@@s[[#]]]>0]&]},If[c=={},s,csm[Sort[Append[Delete[s,List/@c[[1]]],Union@@s[[c[[1]]]]]]]]];
    Table[Length[csm[bpe/@bpe[n]]],{n,0,100}]
  • Python
    from sympy.utilities.iterables import connected_components
    def bin_i(n): #binary indices
        return([(i+1) for i, x in enumerate(bin(n)[2:][::-1]) if x =='1'])
    def A326753(n):
        E,a = [],[bin_i(k) for k in bin_i(n)]
        m = len(a)
        for i in range(m):
            for j in a[i]:
                for k in range(m):
                    if j in a[k]:
                        E.append((i,k))
        return(len(connected_components((list(range(m)),E)))) # John Tyler Rascoe, Jul 16 2024

Formula

a(A072639(n)) = n. - John Tyler Rascoe, Jul 15 2024

A326786 Cut-connectivity of the set-system with BII-number n.

Original entry on oeis.org

0, 1, 1, 0, 2, 2, 2, 2, 1, 0, 0, 0, 0, 0, 0, 0, 2, 2, 0, 0, 1, 1, 1, 1, 2, 2, 0, 0, 1, 1, 1, 1, 2, 0, 2, 0, 1, 1, 1, 1, 2, 0, 2, 0, 1, 1, 1, 1, 1, 1, 1, 1, 3, 3, 3, 3, 1, 1, 1, 1, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3
Offset: 0

Views

Author

Gus Wiseman, Jul 25 2019

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. We define the set-system with BII-number n to be obtained by taking the binary indices of each binary index of n. Every finite set of finite nonempty sets has a different BII-number. For example, 18 has reversed binary expansion (0,1,0,0,1), and since the binary indices of 2 and 5 are {2} and {1,3} respectively, the BII-number of {{2},{1,3}} is 18.
Elements of a set-system are sometimes called edges. The cut-connectivity of a set-system is the minimum number of vertices that must be removed (together with any resulting empty or duplicate edges) to obtain a disconnected or empty set-system. Except for cointersecting set-systems (A326853), this is the same as vertex-connectivity (A327051).

Examples

			Positions of first appearances of each integer, together with the corresponding set-systems, are:
     0: {}
     1: {{1}}
     4: {{1,2}}
    52: {{1,2},{1,3},{2,3}}
  2868: {{1,2},{1,3},{2,3},{1,4},{2,4},{3,4}}
		

Crossrefs

Cf. A000120, A013922, A029931, A048793, A070939, A305078, A322388, A322389 (same for MM-numbers), A322390, A326031, A326701, A326749, A326753, A326787 (edge-connectivity), A327051 (vertex-connectivity).

Programs

  • Mathematica
    bpe[n_]:=Join@@Position[Reverse[IntegerDigits[n,2]],1];
    csm[s_]:=With[{c=Select[Tuples[Range[Length[s]],2],And[OrderedQ[#],UnsameQ@@#,Length[Intersection@@s[[#]]]>0]&]},If[c=={},s,csm[Sort[Append[Delete[s,List/@c[[1]]],Union@@s[[c[[1]]]]]]]]];
    vertConn[y_]:=If[Length[csm[bpe/@y]]!=1,0,Min@@Length/@Select[Subsets[Union@@bpe/@y],Function[del,Length[csm[DeleteCases[DeleteCases[bpe/@y,Alternatives@@del,{2}],{}]]]!=1]]];
    Table[vertConn[bpe[n]],{n,0,100}]

A318978 Heinz numbers of integer partitions with a common divisor > 1.

Original entry on oeis.org

3, 5, 7, 9, 11, 13, 17, 19, 21, 23, 25, 27, 29, 31, 37, 39, 41, 43, 47, 49, 53, 57, 59, 61, 63, 65, 67, 71, 73, 79, 81, 83, 87, 89, 91, 97, 101, 103, 107, 109, 111, 113, 115, 117, 121, 125, 127, 129, 131, 133, 137, 139, 147, 149, 151, 157, 159, 163, 167, 169
Offset: 1

Views

Author

Gus Wiseman, Sep 06 2018

Keywords

Comments

The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).
Is this the same as A305078 without the leading 2? - R. J. Mathar, Sep 08 2018

Examples

			The sequence of all integer partitions with a common divisor begins: (2), (3), (4), (2,2), (5), (6), (7), (8), (4,2), (9), (3,3), (2,2,2), (10), (11), (12), (6,2), (13), (14), (15), (4,4), (16), (8,2), (17), (18), (4,2,2), (6,3), (19), (20), (21), (22), (2,2,2,2), (23), (10,2), (24), (6,4), (25).
		

Crossrefs

Programs

  • Mathematica
    Select[Range[100],GCD@@PrimePi/@If[#==1,{},FactorInteger[#]][[All,1]]>1&]

A322338 Edge-connectivity of the integer partition with Heinz number n.

Original entry on oeis.org

0, 1, 1, 0, 1, 0, 1, 0, 2, 0, 1, 0, 1, 0, 0, 0, 1, 0, 1, 0, 2, 0, 1, 0, 2, 0, 3, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 2, 0, 1, 0, 1, 0, 0, 0, 1, 0, 2, 0, 0, 0, 1, 0, 0, 0, 2, 0, 1, 0, 1, 0, 3, 0, 2, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 4, 0, 1, 0, 0, 0, 2
Offset: 1

Views

Author

Gus Wiseman, Dec 04 2018

Keywords

Comments

The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).
The edge-connectivity of an integer partition is the minimum number of parts that must be removed so that the prime factorizations of the remaining parts form a disconnected (or empty) hypergraph.

Examples

			2093 is the Heinz number of (9,6,4), corresponding to the multiset partition {{1,1},{1,2},{2,2}}, which can be made disconnected by removing only the part {1,2}, so a(2093) = 1.
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    csm[s_]:=With[{c=Select[Tuples[Range[Length[s]],2],And[OrderedQ[#],UnsameQ@@#,Length[Intersection@@s[[#]]]>0]&]},If[c=={},s,csm[Sort[Append[Delete[s,List/@c[[1]]],Union@@s[[c[[1]]]]]]]]];
    Table[PrimeOmega[n]-Max@@PrimeOmega/@Select[Divisors[n],Length[csm[primeMS/@primeMS[#]]]!=1&],{n,100}]

A324173 Regular triangle read by rows where T(n,k) is the number of set partitions of {1,...,n} with k topologically connected components.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 1, 3, 1, 0, 2, 6, 6, 1, 0, 6, 15, 20, 10, 1, 0, 21, 51, 65, 50, 15, 1, 0, 85, 203, 252, 210, 105, 21, 1, 0, 385, 912, 1120, 938, 560, 196, 28, 1, 0, 1907, 4527, 5520, 4620, 2898, 1302, 336, 36, 1, 0, 10205, 24370, 29700, 24780, 15792, 7812, 2730, 540, 45, 1
Offset: 0

Views

Author

Gus Wiseman, Feb 17 2019

Keywords

Comments

A set partition is crossing if it contains a pair of blocks of the form {{...x...y...}, {...z...t...}} where x < z < y < t or z < x < t < y.
The topologically connected components of a set partition correspond to the blocks of its minimal non-crossing coarsening.

Examples

			Triangle begins:
     1
     0     1
     0     1     1
     0     1     3     1
     0     2     6     6     1
     0     6    15    20    10     1
     0    21    51    65    50    15     1
     0    85   203   252   210   105    21     1
     0   385   912  1120   938   560   196    28     1
     0  1907  4527  5520  4620  2898  1302   336    36     1
     0 10205 24370 29700 24780 15792  7812  2730   540    45     1
Row n = 4 counts the following set partitions:
  {{1234}}    {{1}{234}}  {{1}{2}{34}}  {{1}{2}{3}{4}}
  {{13}{24}}  {{12}{34}}  {{1}{23}{4}}
              {{123}{4}}  {{12}{3}{4}}
              {{124}{3}}  {{1}{24}{3}}
              {{134}{2}}  {{13}{2}{4}}
              {{14}{23}}  {{14}{2}{3}}
		

Crossrefs

Programs

  • Mathematica
    croXQ[stn_]:=MatchQ[stn,{_,{_,x_,_,y_,_},_,{_,z_,_,t_,_},_}/;x0]&]},If[c=={},s,csm[Sort[Append[Delete[s,List/@c[[1]]],Union@@s[[c[[1]]]]]]]]];
    crosscmpts[stn_]:=csm[Union[Subsets[stn,{1}],Select[Subsets[stn,{2}],croXQ]]];
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    Table[Length[Select[sps[Range[n]],Length[crosscmpts[#]]==k&]],{n,0,8},{k,0,n}]

A329559 MM-numbers of multiset clutters (connected weak antichains of multisets).

Original entry on oeis.org

1, 2, 3, 5, 7, 9, 11, 13, 17, 19, 23, 25, 27, 29, 31, 37, 41, 43, 47, 49, 53, 59, 61, 67, 71, 73, 79, 81, 83, 89, 91, 97, 101, 103, 107, 109, 113, 121, 125, 127, 131, 137, 139, 149, 151, 157, 163, 167, 169, 173, 179, 181, 191, 193, 197, 199, 203, 211, 223, 227
Offset: 1

Views

Author

Gus Wiseman, Nov 18 2019

Keywords

Comments

A weak antichain of multisets is a multiset of multisets, none of which is a proper subset of any other.

Examples

			The sequence of terms tother with their corresponding clutters begins:
   1: {}              37: {{1,1,2}}            91: {{1,1},{1,2}}
   2: {{}}            41: {{6}}                97: {{3,3}}
   3: {{1}}           43: {{1,4}}             101: {{1,6}}
   5: {{2}}           47: {{2,3}}             103: {{2,2,2}}
   7: {{1,1}}         49: {{1,1},{1,1}}       107: {{1,1,4}}
   9: {{1},{1}}       53: {{1,1,1,1}}         109: {{10}}
  11: {{3}}           59: {{7}}               113: {{1,2,3}}
  13: {{1,2}}         61: {{1,2,2}}           121: {{3},{3}}
  17: {{4}}           67: {{8}}               125: {{2},{2},{2}}
  19: {{1,1,1}}       71: {{1,1,3}}           127: {{11}}
  23: {{2,2}}         73: {{2,4}}             131: {{1,1,1,1,1}}
  25: {{2},{2}}       79: {{1,5}}             137: {{2,5}}
  27: {{1},{1},{1}}   81: {{1},{1},{1},{1}}   139: {{1,7}}
  29: {{1,3}}         83: {{9}}               149: {{3,4}}
  31: {{5}}           89: {{1,1,1,2}}         151: {{1,1,2,2}}
		

Crossrefs

Connected numbers are A305078.
Stable numbers are A316476.
Clutters (of sets) are A048143.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    zsm[s_]:=With[{c=Select[Tuples[Range[Length[s]],2],And[Less@@#,GCD@@s[[#]]]>1&]},If[c=={},s,zsm[Sort[Append[Delete[s,List/@c[[1]]],LCM@@s[[c[[1]]]]]]]]];
    stableQ[u_,Q_]:=!Apply[Or,Outer[#1=!=#2&&Q[#1,#2]&,u,u,1],{0,1}];
    Select[Range[100],And[stableQ[primeMS[#],Divisible],Length[zsm[primeMS[#]]]<=1]&]

Formula

Equals {1} followed by the intersection of A305078 and A316476.
Showing 1-10 of 85 results. Next