A304716 -id:A304716 - cp's OEIS frontend

A305078 Heinz numbers of connected integer partitions.

2, 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

Offset: 1

Gus Wiseman, May 24 2018

nonn

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. This sequence lists all Heinz numbers of multisets S such that G(S) is a connected graph.

			The sequence of all connected multiset multisystems (see A302242, A112798) begins:
   2: {{}}
   3: {{1}}
   5: {{2}}
   7: {{1,1}}
   9: {{1},{1}}
  11: {{3}}
  13: {{1,2}}
  17: {{4}}
  19: {{1,1,1}}
  21: {{1},{1,1}}
  23: {{2,2}}
  25: {{2},{2}}
  27: {{1},{1},{1}}
  29: {{1,3}}
  31: {{5}}
  37: {{1,1,2}}
  39: {{1},{1,2}}
  41: {{6}}
  43: {{1,4}}
  47: {{2,3}}
  49: {{1,1},{1,1}}
  53: {{1,1,1,1}}
  57: {{1},{1,1,1}}
  59: {{7}}
  61: {{1,2,2}}
  63: {{1},{1},{1,1}}
  65: {{2},{1,2}}
  67: {{8}}
  71: {{1,1,3}}
  73: {{2,4}}
  79: {{1,5}}
  81: {{1},{1},{1},{1}}
  83: {{9}}
  87: {{1},{1,3}}
  89: {{1,1,1,2}}
  91: {{1,1},{1,2}}
  97: {{3,3}}

Madeline Locus Dawsey, Tyler Russell and Dannie Urban, Polynomials Associated to Integer Partitions, arXiv:2108.00943 [math.NT], 2021.

Cf. A001221, A048143, A056239, A112798, A286518, A286520, A290103, A302242, A303837, A304118, A304714, A304716, A305052, A305055, A305079.

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[Union[Append[Delete[s,List/@c[[1]]],LCM@@s[[c[[1]]]]]]]]];
Select[Range[300],Length[zsm[primeMS[#]]]==1&]

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

Gus Wiseman, May 24 2018

nonn

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

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

Antti Karttunen, Table of n, a(n) for n = 1..16384
Antti Karttunen, Data supplement: n, a(n) computed for n = 1..100000
Index entries for sequences computed from indices in prime factorization
Index entries for sequences related to Heinz numbers

Cf. A001221, A048143, A056239, A112798, A286518, A286520, A290103, A302242, A303837, A304118, A304714, A304716, A305052, A305055, A305078, A305501.

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

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

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

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

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

Don Knuth, Jan 26 2008

nonn

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

			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)

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

Alois P. Heinz, Table of n, a(n) for n = 0..370
N. J. A. Sloane, Transforms

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

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

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

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

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

Gus Wiseman, Jul 23 2019

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 set-system {{1,2},{1,4},{3}} with BII-number 268 has two connected components, so a(268) = 2.

John Tyler Rascoe, Table of n, a(n) for n = 0..10000
John Tyler Rascoe, Log scatterplot of a(n), n=0..32906.

Positions of 0's and 1's are A326749.
Cf. A000120, A001187, A029931, A048143, A048793, A070939, A072639, A304716, A305078, A305079 (same for MM-numbers), A323818, A326031, A326702.
Ranking sequences using BII-numbers: A309314 (hyperforests), A326701 (set partitions), A326703 (chains), A326704 (antichains), A326750 (clutters), A326751 (blobs), A326752 (hypertrees), A326754 (covers).

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

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

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

A259936 Number of ways to express the integer n as a product of its unitary divisors (A034444).

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 2, 1, 1, 2, 1, 2, 2, 2, 1, 2, 1, 2, 1, 2, 1, 5, 1, 1, 2, 2, 2, 2, 1, 2, 2, 2, 1, 5, 1, 2, 2, 2, 1, 2, 1, 2, 2, 2, 1, 2, 2, 2, 2, 2, 1, 5, 1, 2, 2, 1, 2, 5, 1, 2, 2, 5, 1, 2, 1, 2, 2, 2, 2, 5, 1, 2, 1, 2, 1, 5, 2, 2, 2, 2, 1, 5, 2, 2, 2, 2, 2, 2, 1, 2, 2, 2, 1, 5, 1, 2, 5

Offset: 1

Geoffrey Critzer, Jul 09 2015

nonn

Equivalently, a(n) is the number of ways to express the cyclic group Z_n as a direct sum of its Hall subgroups.  A Hall subgroup of a finite group G is a subgroup whose order is coprime to its index.
a(n) is the number of ways to partition the set of distinct prime factors of n.
Also the number of singleton or pairwise coprime factorizations of n. - Gus Wiseman, Sep 24 2019

			a(60) = 5 because we have: 60 = 4*3*5 = 4*15 = 3*20 = 5*12.
For n = 36, its unitary divisors are 1, 4, 9, 36. From these we obtain 36 either as 1*36 or 4*9, thus a(36) = 2. - _Antti Karttunen_, Oct 21 2017

Alois P. Heinz, Table of n, a(n) for n = 1..20000
Wikipedia, Hall subgroup
Index entries for sequences computed from exponents in factorization of n

Cf. A000110, A001055, A001221, A034444, A089233, A258466, A281116, A285572.
Differs from A050320 for the first time at n=36.
Differs from A354870 for the first time at n=210, where a(210) = 15, while A354870(210) = 12.
Cf. A304716, A302569, A304711, A305079.
Related classes of factorizations:
- No conditions: A001055
- Strict: A045778
- Constant: A089723
- Distinct multiplicities: A255231
- Singleton or coprime: A259936
- Relatively prime: A281116
- Aperiodic: A303386
- Stable (indivisible): A305149
- Connected: A305193
- Strict relatively prime: A318721
- Uniform: A319269
- Intersecting: A319786
- Constant or distinct factors coprime: A327399
- Constant or relatively prime: A327400
- Coprime: A327517
- Not relatively prime: A327658
- Distinct factors coprime: A327695

map(combinat:-bell @ nops @ numtheory:-factorset, [$1..100]); # Robert Israel, Jul 09 2015

Table[BellB[PrimeNu[n]], {n, 1, 75}]
(* second program *)
facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
Table[Length[Select[facs[n],Length[#]==1||CoprimeQ@@#&]],{n,100}] (* Gus Wiseman, Sep 24 2019 *)

a(n) = my(t=omega(n), x='x, m=contfracpnqn(matrix(2, t\2, y, z, if( y==1, -z*x^2, 1 - (z+1)*x)))); polcoeff(1/(1 - x + m[2, 1]/m[1, 1]) + O(x^(t+1)), t) \\ Charles R Greathouse IV, Jun 30 2017

a(n) = A000110(A001221(n)).

a(n > 1) = A327517(n) + 1. - Gus Wiseman, Sep 24 2019

Incorrect comment removed by Antti Karttunen, Jun 11 2022

A305052 z-density of the integer partition with Heinz number n. Clutter density of the n-th multiset multisystem (A302242).

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

Offset: 1

Gus Wiseman, May 24 2018

sign

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
The z-density of a multiset S of positive integers is Sum_{s in S} (omega(s) - 1) - omega(lcm(S)) where omega = A001221 is number of distinct prime factors.
First nonnegative entry after a(1) = 0 is a(169) = 0.

			The 1105th multiset multisystem is {{2},{1,2},{4}} with clutter density -2, so a(1105) = -2.
The 5429th multiset multisystem is {{1,2,2},{1,1,1,2}} with clutter density 0, so a(5429) = 0.
The 11837th multiset multisystem is {{1,1},{1,1,1},{1,1,1,2}} with clutter density -1, so a(11837) = -1.
The 42601th multiset multisystem is {{1,2},{1,3},{1,2,3}} with clutter density 1, so a(42601) = 1.

Cf. A001221, A030019, A048143, A056239, A112798, A285572, A286518, A286520, A290103, A302242, A303837, A304118, A304714, A304716, A304911, A305001.

zens[n_]:=If[n==1,0,Total@Cases[FactorInteger[n],{p_,k_}:>k*(PrimeNu[PrimePi[p]]-1)]-PrimeNu[LCM@@Cases[FactorInteger[n],{p_,k_}:>PrimePi[p]]]];
Array[zens,100]

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

Gus Wiseman, Dec 04 2018

nonn

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.

			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.

Wikipedia, k-edge-connected graph

Cf. A001222, A003963, A013922, A054921, A056239, A095983, A112798, A302242, A304714, A304716, A305078, A305079, A322335, A322336, A322337.

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

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

Don Knuth, Jan 26 2008

nonn

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

			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)

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

Alois P. Heinz, Table of n, a(n) for n = 0..1000 (first 201 terms from T. D. Noe)
N. J. A. Sloane, Transforms
Wikipedia, Block graph

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

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

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

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

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

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

Gus Wiseman, Nov 18 2019

nonn

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

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

Connected numbers are A305078.
Stable numbers are A316476.
Clutters (of sets) are A048143.
Cf. A056239, A112798, A289509, A302242, A302494, A304716, A318991, A319837, A320275, A320456, A328514, A329553, A329555.

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]&]

Equals {1} followed by the intersection of A305078 and A316476.

A144959 A134955(n) - A134955(n-1). Number of hyperforests spanning n unlabeled nodes without isolated vertices.

Original entry on oeis.org

1, 0, 1, 2, 5, 11, 30, 78, 223, 658, 2026, 6429, 21015, 70233, 239360, 829224, 2912947, 10356334, 37205121, 134887153, 493000086, 1814902409, 6724595543, 25061885217, 93899071368, 353514105817, 1336822098961, 5075833932200

Offset: 0

Washington Bomfim, Sep 27 2008

nonn

a(n) is the number of hyperforests with n unlabeled nodes without isolated vertices. This follows from the fact that for n>0 A134955(n-1) counts the hyperforests of order n with one or more isolated nodes.

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

Cf. A030019, A035053, A048143, A054921, A134954, A134955, A134957, A144958 (unlabeled forests without isolated vertices), A144959, A304716, A304717, A304867, A304911.

etr[p_] := 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[0] = 0; b[n_] := b[n] = etr[etr[b]][n-1];
c[1] = 0; c[n_] := b[n] + etr[b][n] - Sum[b[k]*etr[b][n-k], {k, 0, n}];
a = etr[c];
Table[a[n], {n, 0, 27}] (* Jean-François Alcover, Jul 12 2018, after Alois P. Heinz's code for A035053 *)

\\ 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(concat([0], Vec(Ser(EulerT(u))*(1-x*Ser(u))-1))))} \\ Andrew Howroyd, May 22 2018

Euler transform of b(1) = 0, b(n > 1) = A035053(n). - Gus Wiseman, May 21 2018

cp's OEIS Frontend

A305078 Heinz numbers of connected integer partitions.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

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

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

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

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

Formula

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

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Python

Formula

A259936 Number of ways to express the integer n as a product of its unitary divisors (A034444).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

Extensions

A305052 z-density of the integer partition with Heinz number n. Clutter density of the n-th multiset multisystem (A302242).

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

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

Views

Author

Keywords

Comments

Examples

Links

Crossrefs