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

A305078 Heinz numbers of connected integer partitions.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, May 24 2018

Keywords

Comments

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.

Examples

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

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

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

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

Views

Author

Gus Wiseman, May 24 2018

Keywords

Comments

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.

Examples

			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.
		

Crossrefs

Programs

  • Mathematica
    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]

A304118 Number of z-blobs with least common multiple n > 1.

Original entry on oeis.org

0, 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, 2, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 7, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 3, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 7, 1, 1, 1, 1
Offset: 1

Views

Author

Gus Wiseman, May 19 2018

Keywords

Comments

Given a finite set S of positive integers greater than 1, let G(S) be the simple labeled graph with vertex set S and edges between any two vertices that have a common divisor greater than 1. For example, G({6,14,15,35}) is a 4-cycle. A set S is said to be connected if G(S) is a connected graph. The clutter density of S is defined to be Sum_{s in S} (omega(s) - 1) - omega(LCM(S)), where omega = A001221 and LCM is least common multiple. A z-blob is a finite connected set S of pairwise indivisible positive integers greater than 1 such that no cap of S with at least two edges has clutter density -1.
If n is squarefree with k prime factors, then a(n) = A275307(k).

Examples

			The a(60) = 7 z-blobs together with the corresponding multiset systems (see A112798, A302242) are the following.
        (60): {{1,1,2,3}}
     (12,30): {{1,1,2},{1,2,3}}
     (20,30): {{1,1,3},{1,2,3}}
   (6,15,20): {{1,2},{2,3},{1,1,3}}
  (10,12,15): {{1,3},{1,1,2},{2,3}}
  (12,15,20): {{1,1,2},{2,3},{1,1,3}}
  (12,20,30): {{1,1,2},{1,1,3},{1,2,3}}
The a(120) = 14 z-blobs together with the corresponding multiset systems are the following.
       (120): {{1,1,1,2,3}}
     (24,30): {{1,1,1,2},{1,2,3}}
     (24,60): {{1,1,1,2},{1,1,2,3}}
     (30,40): {{1,2,3},{1,1,1,3}}
     (40,60): {{1,1,1,3},{1,1,2,3}}
   (6,15,40): {{1,2},{2,3},{1,1,1,3}}
  (10,15,24): {{1,3},{2,3},{1,1,1,2}}
  (12,15,40): {{1,1,2},{2,3},{1,1,1,3}}
  (12,30,40): {{1,1,2},{1,2,3},{1,1,1,3}}
  (15,20,24): {{2,3},{1,1,3},{1,1,1,2}}
  (15,24,40): {{2,3},{1,1,1,2},{1,1,1,3}}
  (20,24,30): {{1,1,3},{1,1,1,2},{1,2,3}}
  (24,30,40): {{1,1,1,2},{1,2,3},{1,1,1,3}}
  (24,40,60): {{1,1,1,2},{1,1,1,3},{1,1,2,3}}
		

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]]]]]]]]];
    zensity[s_]:=Total[(PrimeNu[#]-1&)/@s]-PrimeNu[LCM@@s];
    zreeQ[s_]:=And[Length[s]>=2,zensity[s]==-1];
    zlobQ[s_]:=Apply[And,Composition[Not,zreeQ]/@Apply[LCM,zptns[s],{2}]];
    zswell[s_]:=Union[LCM@@@Select[Subsets[s],Length[zsm[#]]==1&]];
    zkernels[s_]:=Table[Select[s,Divisible[w,#]&],{w,zswell[s]}];
    zptns[s_]:=Select[stableSets[zkernels[s],Length[Intersection[#1,#2]]>0&],Union@@#==s&];
    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[If[n==1,0,Length[Select[Rest[Subsets[Rest[Divisors[n]]]],And[zsm[#]=={n},Select[Tuples[#,2],UnsameQ@@#&&Divisible@@#&]=={},zlobQ[#]]&]]],{n,100}]

A303838 Number of z-forests with least common multiple n > 1.

Original entry on oeis.org

0, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 3, 1, 2, 2, 1, 1, 3, 1, 3, 2, 2, 1, 4, 1, 2, 1, 3, 1, 8, 1, 1, 2, 2, 2, 5, 1, 2, 2, 4, 1, 8, 1, 3, 3, 2, 1, 5, 1, 3, 2, 3, 1, 4, 2, 4, 2, 2, 1, 16, 1, 2, 3, 1, 2, 8, 1, 3, 2, 8, 1, 7, 1, 2, 3, 3, 2, 8, 1, 5, 1, 2, 1, 16, 2, 2
Offset: 1

Views

Author

Gus Wiseman, May 19 2018

Keywords

Comments

Given a finite set S of positive integers greater than 1, let G(S) be the simple labeled graph with vertex set S and edges between any two vertices that have a common divisor greater than 1. For example, G({6,14,15,35}) is a 4-cycle. A set S is said to be connected if G(S) is a connected graph. The clutter density of S is defined to be Sum_{s in S} (omega(s) - 1) - omega(LCM(S)), where omega = A001221 and LCM is least common multiple. A z-forest is a finite set of pairwise indivisible positive integers greater than 1 such that all connected components are z-trees, meaning they have clutter density -1.
This is a generalization to multiset systems of the usual definition of hyperforest (viz. hypergraph F such that two distinct hyperedges of F intersect in at most a common vertex and such that every cycle of F is contained in a hyperedge).
If n is squarefree with k prime factors, then a(n) = A134954(k).
Differs from A324837 at positions {1, 180, 210, ...}. For example, a(210) = 55, A324837(210) = 49.

Examples

			The a(60) = 16 z-forests together with the corresponding multiset systems (see A112798, A302242) are the following.
       (60): {{1,1,2,3}}
     (3,20): {{2},{1,1,3}}
     (4,15): {{1,1},{2,3}}
     (4,30): {{1,1},{1,2,3}}
     (5,12): {{3},{1,1,2}}
     (6,20): {{1,2},{1,1,3}}
    (10,12): {{1,3},{1,1,2}}
    (12,15): {{1,1,2},{2,3}}
    (12,20): {{1,1,2},{1,1,3}}
    (15,20): {{2,3},{1,1,3}}
    (3,4,5): {{2},{1,1},{3}}
   (3,4,10): {{2},{1,1},{1,3}}
    (4,5,6): {{1,1},{3},{1,2}}
   (4,6,10): {{1,1},{1,2},{1,3}}
   (4,6,15): {{1,1},{1,2},{2,3}}
  (4,10,15): {{1,1},{1,3},{2,3}}
		

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]]]]]]]]];
    zensity[s_]:=Total[(PrimeNu[#]-1&)/@s]-PrimeNu[LCM@@s];
    Table[Length[Select[Rest[Subsets[Rest[Divisors[n]]]],Function[s,LCM@@s==n&&And@@Table[zensity[Select[s,Divisible[m,#]&]]==-1,{m,zsm[s]}]&&Select[Tuples[s,2],UnsameQ@@#&&Divisible@@#&]=={}]]],{n,100}]

A305193 Number of connected factorizations of n.

Original entry on oeis.org

0, 1, 1, 2, 1, 1, 1, 3, 2, 1, 1, 2, 1, 1, 1, 5, 1, 2, 1, 2, 1, 1, 1, 4, 2, 1, 3, 2, 1, 1, 1, 7, 1, 1, 1, 5, 1, 1, 1, 4, 1, 1, 1, 2, 2, 1, 1, 7, 2, 2, 1, 2, 1, 4, 1, 4, 1, 1, 1, 3, 1, 1, 2, 11, 1, 1, 1, 2, 1, 1, 1, 10, 1, 1, 2, 2, 1, 1, 1, 7, 5, 1, 1, 3, 1, 1, 1, 4, 1, 3, 1, 2, 1, 1, 1, 12, 1, 2, 2, 5, 1, 1, 1, 4, 1
Offset: 1

Views

Author

Gus Wiseman, May 27 2018

Keywords

Comments

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 counts factorizations S such that G(S) is a connected graph.
a(n) depends only on prime signature of n (cf. A025487). - Antti Karttunen, Nov 07 2018

Examples

			The a(72) = 10 factorizations:
(72),
(2*2*18), (2*3*12), (2*6*6), (3*4*6),
(2*36), (3*24), (4*18), (6*12),
(2*2*3*6).
		

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]]]]]]]]];
    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[zsm[#]]==1&]],{n,100}]
  • PARI
    is_connected(facs) = { my(siz=length(facs)); if(1==siz,1,my(m=matrix(siz,siz,i,j,(gcd(facs[i],facs[j])!=1))^siz); for(n=1,siz,if(0==vecmin(m[n,]),return(0))); (1)); };
    A305193aux(n, m, facs) = if(1==n, is_connected(Set(facs)), my(s=0, newfacs); fordiv(n, d, if((d>1)&&(d<=m), newfacs = List(facs); listput(newfacs,d); s += A305193aux(n/d, d, newfacs))); (s)); \\ Antti Karttunen, Nov 07 2018
    A305193(n) = if(1==n,0,A305193aux(n, n, List([]))); \\ Antti Karttunen, Nov 07 2018

Extensions

More terms from Antti Karttunen, Nov 07 2018

A304382 Number of z-trees summing to n. Number of connected strict integer partitions of n with pairwise indivisible parts and clutter density -1.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 3, 2, 2, 1, 3, 2, 4, 3, 5, 2, 5, 4, 6, 3, 7, 6, 8, 4, 9, 8, 13, 9, 15, 8, 14, 12, 16, 12, 20, 20, 24, 15, 27, 20, 33, 27, 35
Offset: 1

Views

Author

Gus Wiseman, May 21 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.
The clutter density of a multiset S of positive integers is Sum_{s in S} (omega(s) - 1) - omega(LCM(S)).

Examples

			The a(30) = 8 z-trees together with the corresponding multiset systems are the following.
       (30): {{1,2,3}}
     (26,4): {{1,6},{1,1}}
     (22,8): {{1,5},{1,1,1}}
     (21,9): {{2,4},{2,2}}
    (16,14): {{1,1,1,1},{1,4}}
   (15,9,6): {{2,3},{2,2},{1,2}}
  (14,10,6): {{1,4},{1,3},{1,2}}
  (12,10,8): {{1,1,2},{1,3},{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]]]]]]]]];
    zensity[s_]:=Total[(PrimeNu[#]-1&)/@s]-PrimeNu[LCM@@s];
    zreeQ[s_]:=And[Length[s]>=2,zensity[s]==-1];
    strConnAnti[n_]:=Select[IntegerPartitions[n],UnsameQ@@#&&Length[zsm[#]]==1&&Select[Tuples[#,2],UnsameQ@@#&&Divisible@@#&]=={}&];
    Table[Length[Select[strConnAnti[n],Length[#]==1||zreeQ[#]&]],{n,20}]

A317073 Number of antichains of multisets with multiset-join a normal multiset of size n.

Original entry on oeis.org

1, 1, 3, 16, 198, 9890, 8592538
Offset: 0

Views

Author

Gus Wiseman, Jul 20 2018

Keywords

Comments

An antichain of multisets is a finite set of finite nonempty multisets, none of which is a submultiset of any other. A multiset is normal if it spans an initial interval of positive integers. The multiset-join of a set of multisets has the same vertices with multiplicities equal to the maxima of the multiplicities in the edges.

Examples

			The a(3) = 16 antichains of multisets:
  (111),
  (122), (12)(22), (1)(22),
  (112), (11)(12), (2)(11),
  (123), (13)(23), (12)(23), (12)(13), (12)(13)(23), (3)(12), (2)(13), (1)(23), (1)(2)(3).
		

Crossrefs

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]]]];
    multijoin[mss__]:=Join@@Table[Table[x,{Max[Count[#,x]&/@{mss}]}],{x,Union[mss]}]
    submultisetQ[M_,N_]:=Or[Length[M]==0,MatchQ[{Sort[List@@M],Sort[List@@N]},{{x_,Z___},{_,x_,W___}}/;submultisetQ[{Z},{W}]]];
    allnorm[n_]:=Function[s,Array[Count[s,y_/;y<=#]+1&,n]]/@Subsets[Range[n-1]+1];
    auu[m_]:=Select[stableSets[Union[Rest[Subsets[m]]],submultisetQ],multijoin@@#==m&];
    Table[Length[Join@@Table[auu[m],{m,allnorm[n]}]],{n,5}]

Extensions

a(6) from Robert Price, Jun 21 2021

A317077 Number of connected multiset partitions of normal multisets of size n.

Original entry on oeis.org

1, 1, 3, 8, 28, 110, 519, 2749, 16317, 106425, 755425, 5781956, 47384170, 413331955, 3818838624, 37213866876, 381108145231, 4088785729738, 45829237977692, 535340785268513, 6502943193997922, 81984445333355812, 1070848034863526547, 14467833457108560375, 201894571410270034773
Offset: 0

Views

Author

Gus Wiseman, Jul 20 2018

Keywords

Comments

A multiset is normal if it spans an initial interval of positive integers.

Examples

			The a(3) = 8 connected multiset partitions are (111), (1)(11), (1)(1)(1), (122), (2)(12), (112), (1)(12), (123).
		

Crossrefs

Programs

  • Mathematica
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    mps[set_]:=Union[Sort[Sort/@(#/.x_Integer:>set[[x]])]&/@sps[Range[Length[set]]]];
    csm[s_]:=With[{c=Select[Tuples[Range[Length[s]],2],And[OrderedQ[#],UnsameQ@@#,Length[Intersection@@s[[#]]]>0]&]},If[c=={},s,csm[Union[Append[Delete[s,List/@c[[1]]],multijoin@@s[[c[[1]]]]]]]]];
    allnorm[n_]:=Function[s,Array[Count[s,y_/;y<=#]+1&,n]]/@Subsets[Range[n-1]+1];
    Length/@Table[Join@@Table[Select[mps[m],Length[csm[#]]==1&],{m,allnorm[n]}],{n,8}]
  • PARI
    EulerT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, n, 1/n))))-1, -#v)}
    Connected(v)={my(u=vector(#v));for(n=1, #u, u[n]=v[n] - sum(k=1, n-1, binomial(n-1,k)*v[k]*u[n-k]));u}
    seq(n)={my(u=vector(n, k, x*Ser(EulerT(vector(n,i,binomial(i+k-1,i)))))); Vec(1+vecsum(Connected(vector(n, k, sum(i=1, k, (-1)^(k-i)*binomial(k,i)*u[i])))))} \\ Andrew Howroyd, Jan 16 2023

Extensions

Terms a(9) and beyond from Andrew Howroyd, Jan 16 2023

A324837 Number of minimal subsets of {1...n} with least common multiple n.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Mar 19 2019

Keywords

Comments

Note that the elements must be pairwise indivisible divisors of n.
Differs from A303838 at positions {1, 180, 210, ...}. For example, a(210) = 49, A303838(210) = 55. - Gus Wiseman, Apr 01 2019

Examples

			The a(30) = 8 subsets are: {30}, {2,15}, {3,10}, {5,6}, {6,10}, {6,15}, {10,15}, {2,3,5}.
		

Crossrefs

Programs

  • Mathematica
    minim[s_]:=Complement[s,First/@Select[Tuples[s,2],UnsameQ@@#&&SubsetQ@@#&]];
    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[minim[Select[Rest[stableSets[Divisors[n],Divisible]],LCM@@#==n&]]],{n,100}]
Showing 1-10 of 26 results. Next