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

A286518 Number of finite connected sets of positive integers greater than one with least common multiple n.

Original entry on oeis.org

1, 1, 1, 2, 1, 4, 1, 4, 2, 4, 1, 20, 1, 4, 4, 8, 1, 20, 1, 20, 4, 4, 1, 88, 2, 4, 4, 20, 1, 96, 1, 16, 4, 4, 4, 196, 1, 4, 4, 88, 1, 96, 1, 20, 20, 4, 1, 368, 2, 20, 4, 20, 1, 88, 4, 88, 4, 4, 1, 1824, 1, 4, 20, 32, 4, 96, 1, 20, 4, 96, 1, 1688, 1, 4, 20, 20, 4, 96, 1, 368, 8, 4, 1, 1824, 4, 4, 4, 88, 1, 1824, 4, 20
Offset: 1

Views

Author

Gus Wiseman, Jul 24 2017

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 that are not relatively prime. 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.
a(n) depends only on prime signature of n (cf. A025487). - Antti Karttunen, Feb 17 2024

Examples

			The a(6)=4 sets are: {6}, {2,6}, {3,6}, {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]]]]]]]]];
    Table[Length[Select[Subsets[Rest[Divisors[n]]],zsm[#]==={n}&]],{n,2,20}]
  • PARI
    isconnected(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)); };
    A286518aux(n, parts, from=1, ss=List([])) = { my(k = #parts, s=0, newss); if(lcm(Vec(ss))==n && isconnected(ss), s++); for(i=from, k, newss = List(ss); listput(newss, parts[i]); s += A286518aux(n, parts, i+1, newss)); (s) };
    A286518(n) = if(1==n, n, A286518aux(n, divisors(n))); \\ Antti Karttunen, Feb 17 2024

Formula

From Antti Karttunen, Feb 17 2024: (Start)
a(n) <= A069626(n).
It seems that a(n) >= A318670(n), for all n > 1.
(End)

Extensions

Term a(1)=1 prepended and more terms added by Antti Karttunen, Feb 17 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

Views

Author

Geoffrey Critzer, Jul 09 2015

Keywords

Comments

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

Examples

			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
		

Crossrefs

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

Programs

  • Maple
    map(combinat:-bell @ nops @ numtheory:-factorset, [$1..100]); # Robert Israel, Jul 09 2015
  • Mathematica
    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 *)
  • PARI
    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

Formula

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

Extensions

Incorrect comment removed by Antti Karttunen, Jun 11 2022

A320632 Numbers k such that there exists a pair of factorizations of k into factors > 1 where no factor of one divides any factor of the other.

Original entry on oeis.org

36, 60, 72, 84, 90, 100, 108, 120, 126, 132, 140, 144, 150, 156, 168, 180, 196, 198, 200, 204, 210, 216, 220, 225, 228, 234, 240, 252, 260, 264, 270, 276, 280, 288, 294, 300, 306, 308, 312, 315, 324, 330, 336, 340, 342, 348, 350, 360, 364, 372, 378, 380, 390
Offset: 1

Views

Author

Gus Wiseman, Dec 09 2018

Keywords

Comments

Positions of nonzero terms in A322437 or A322438.
Mats Granvik has conjectured that these are all the positive integers k such that sigma_0(k) - 2 > (bigomega(k) - 1) * omega(k), where sigma_0 = A000005, omega = A001221, and bigomega = A001222. - Gus Wiseman, Nov 12 2019
Numbers with more semiprime divisors than prime divisors. - Wesley Ivan Hurt, Jun 10 2021

Examples

			An example of such a pair for 36 is (4*9)|(6*6).
		

Crossrefs

The following are additional cross-references relating to Granvik's conjecture.
bigomega(n) * omega(n) is A113901(n).
(bigomega(n) - 1) * omega(n) is A307409(n).
sigma_0(n) - bigomega(n) * omega(n) is A328958(n).
sigma_0(n) - 2 - (omega(n) - 1) * nu(n) is A328959(n).

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Select[Range[100],Select[Subsets[facs[#],{2}],And[!Or@@Divisible@@@Tuples[#],!Or@@Divisible@@@Reverse/@Tuples[#]]&]!={}&]
  • PARI
    factorizations(n, m=n, f=List([]), z=List([])) = if(1==n, listput(z,Vec(f)); z, my(newf); fordiv(n, d, if((d>1)&&(d<=m), newf = List(f); listput(newf,d); z = factorizations(n/d, d, newf, z))); (z));
    is_ndf_pair(fac1,fac2) = { for(i=1,#fac1,for(j=1,#fac2,if(!(fac1[i]%fac2[j])||!(fac2[j]%fac1[i]),return(0)))); (1); };
    has_at_least_one_ndfpair(z) = { for(i=1,#z,for(j=i+1,#z,if(is_ndf_pair(z[i],z[j]),return(1)))); (0); };
    isA320632(n) = has_at_least_one_ndfpair(Vec(factorizations(n))); \\ Antti Karttunen, Dec 10 2020

A319786 Number of factorizations of n where no two factors are relatively prime.

Original entry on oeis.org

1, 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, 4, 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, 7, 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, 4, 1, 1, 1, 4, 1
Offset: 1

Views

Author

Gus Wiseman, Sep 27 2018

Keywords

Comments

First differs from A305193 at a(36) = 4, A305193(36) = 5.
a(n) depends only on prime signature of n (cf. A025487). - Antti Karttunen, Nov 07 2018

Examples

			The a(48) = 7 factorizations are (2*2*2*6), (2*2*12), (2*4*6), (2*24), (4*12), (6*8), (48).
		

Crossrefs

Programs

  • Mathematica
    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],!Or@@CoprimeQ@@@Subsets[#,{2}]&]],{n,100}]
  • PARI
    A319786(n, m=n, facs=List([])) = if(1==n, (1!=gcd(Vec(facs))), my(s=0, newfacs); fordiv(n, d, if((d>1)&&(d<=m), newfacs = List(facs); listput(newfacs,d); s += A319786(n/d, d, newfacs))); (s)); \\ Antti Karttunen, Nov 07 2018

Extensions

More terms from Antti Karttunen, Nov 07 2018

A322435 Number of pairs of factorizations of n into factors > 1 where no factor of the second divides any factor of the first.

Original entry on oeis.org

1, 0, 0, 1, 0, 1, 0, 2, 1, 1, 0, 3, 0, 1, 1, 5, 0, 3, 0, 3, 1, 1, 0, 7, 1, 1, 2, 3, 0, 4, 0, 7, 1, 1, 1, 15, 0, 1, 1, 7, 0, 4, 0, 3, 3, 1, 0, 16, 1, 3, 1, 3, 0, 7, 1, 7, 1, 1, 0, 18, 0, 1, 3, 16, 1, 4, 0, 3, 1, 4, 0, 32, 0, 1, 3, 3, 1, 4, 0, 16, 5, 1, 0, 18, 1
Offset: 1

Views

Author

Gus Wiseman, Dec 08 2018

Keywords

Examples

			The a(36) = 15 pairs of factorizations:
  (2*2*3*3)|(4*9)
  (2*2*3*3)|(6*6)
  (2*2*3*3)|(36)
    (2*2*9)|(6*6)
    (2*2*9)|(36)
    (2*3*6)|(4*9)
    (2*3*6)|(36)
     (2*18)|(36)
    (3*3*4)|(6*6)
    (3*3*4)|(36)
     (3*12)|(36)
      (4*9)|(6*6)
      (4*9)|(36)
      (6*6)|(4*9)
      (6*6)|(36)
		

Crossrefs

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Length[Select[Tuples[facs[n],2],!Or@@Divisible@@@Tuples[#]&]],{n,100}]

A322437 Number of unordered pairs of factorizations of n into factors > 1 where no factor of one divides any factor of the other.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Dec 08 2018

Keywords

Comments

First differs from A322438 at a(144) = 3, A322438(144) = 4.
From Antti Karttunen, Dec 11 2020: (Start)
Zeros occur on numbers that are either of the form p^k, or q * p^k, or p*q*r, for some primes p, q, r, and exponent k >= 0. [Note also that in all these cases, when x > 1, A307408(x) = 2+A307409(x) = 2 + (A001222(x) - 1)*A001221(x) = A000005(x)].
Proof:
It is easy to see that for such numbers it is not possible to obtain two such distinct factorizations, that no factor of the other would not divide some factor of the other.
Conversely, the complement set of above is formed of such composites n that have at least one unitary divisor that is either of the form
(1) p^x * q^y, with x, y >= 2,
or
(2) p^x * q^y * r^z, with x >= 2, and y, z >= 1,
or
(3) p^x * q^y * r^z * s^w, with x, y, z, w >= 1,
where p, q, r, s are distinct primes. Let's indicate with C the remaining portion of k coprime to p, q, r and s (which could be also 1). Then in case (1) we can construct two factorizations, the first having factors (p*q*C) and (p^(x-1) * q^(y-1)), and the second having factors (p^x * C) and (q^y) that are guaranteed to satisfy the condition that no factor in the other factorization divides any of the factors of the other factorization. For case (2) pairs like {(p * q^y * C), (p^(x-1) * r^z)} and {(p^x * C), (q^y * r^z)}, and for case (3) pairs like {(p^x * q^y * C), (r^z * s^w)} and {(p^x * r^z * C), (q^y * s^w)} offer similar examples, therefore a(n) > 0 for all such cases.
(End)

Examples

			The a(120) = 2 pairs of such factorizations:
   (6*20)|(8*15)
   (8*15)|(10*12)
The a(144) = 3 pairs of factorizations:
   (6*24)|(9,16)
   (8*18)|(12*12)
   (9*16)|(12*12)
The a(210) = 3 pairs of factorizations:
   (6*35)|(10*21)
   (6*35)|(14*15)
  (10*21)|(14*15)
[Note that 210 is the first squarefree number obtaining nonzero value]
The a(240) = 4 pairs of factorizations:
   (6*40)|(15*16)
   (8*30)|(12*20)
  (10*24)|(15*16)
  (12*20)|(15*16)
The a(1728) = 14 pairs of factorizations:
    (6*6*48)|(27*64)
   (6*12*24)|(27*64)
     (6*288)|(27*64)
    (8*8*27)|(12*12*12)
  (12*12*12)|(27*64)
  (12*12*12)|(32*54)
    (12*144)|(27*64)
    (12*144)|(32*54)
    (16*108)|(24*72)
     (18*96)|(27*64)
     (24*72)|(27*64)
     (24*72)|(32*54)
     (27*64)|(36*48)
     (32*54)|(36*48)
		

Crossrefs

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Length[Select[Subsets[facs[n],{2}],And[!Or@@Divisible@@@Tuples[#],!Or@@Divisible@@@Reverse/@Tuples[#]]&]],{n,100}]
  • PARI
    factorizations(n, m=n, f=List([]), z=List([])) = if(1==n, listput(z,Vec(f)); z, my(newf); fordiv(n, d, if((d>1)&&(d<=m), newf = List(f); listput(newf,d); z = factorizations(n/d, d, newf, z))); (z));
    is_ndf_pair(fac1,fac2) = { for(i=1,#fac1,for(j=1,#fac2,if(!(fac1[i]%fac2[j])||!(fac2[j]%fac1[i]),return(0)))); (1); };
    number_of_ndfpairs(z) = sum(i=1,#z,sum(j=i+1,#z,is_ndf_pair(z[i],z[j])));
    A322437(n) = number_of_ndfpairs(Vec(factorizations(n))); \\ Antti Karttunen, Dec 10 2020

Formula

For n > 0, a(A002110(n)) = A322441(n)/2 = A339626(n). - Antti Karttunen, Dec 10 2020

Extensions

Data section extended up to a(120) and more examples added by Antti Karttunen, Dec 10 2020

A322436 Number of pairs of factorizations of n into factors > 1 where no factor of the second properly divides any factor of the first.

Original entry on oeis.org

1, 1, 1, 3, 1, 3, 1, 5, 3, 3, 1, 8, 1, 3, 3, 11, 1, 8, 1, 8, 3, 3, 1, 18, 3, 3, 5, 8, 1, 12, 1, 15, 3, 3, 3, 31, 1, 3, 3, 18, 1, 12, 1, 8, 8, 3, 1, 39, 3, 8, 3, 8, 1, 18, 3, 18, 3, 3, 1, 42, 1, 3, 8, 33, 3, 12, 1, 8, 3, 12, 1, 67, 1, 3, 8, 8, 3, 12, 1, 39, 11
Offset: 1

Views

Author

Gus Wiseman, Dec 08 2018

Keywords

Examples

			The a(12) = 8 pairs of factorizations:
  (2*2*3)|(2*2*3)
  (2*2*3)|(2*6)
  (2*2*3)|(3*4)
  (2*2*3)|(12)
    (2*6)|(12)
    (3*4)|(3*4)
    (3*4)|(12)
     (12)|(12)
		

Crossrefs

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    divpropQ[x_,y_]:=And[x!=y,Divisible[x,y]];
    Table[Length[Select[Tuples[facs[n],2],!Or@@divpropQ@@@Tuples[#]&]],{n,100}]

A322438 Number of unordered pairs of factorizations of n into factors > 1 where no factor of one properly divides any factor of the other.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Dec 08 2018

Keywords

Comments

First differs from A322437 at a(144) = 4, A322437(144) = 3.
First differs from A379958 at a(120) = 2, A379958(120) = 1.

Examples

			The a(240) = 5 pairs of factorizations::
  (4*4*15)|(4*6*10)
    (6*40)|(15*16)
    (8*30)|(12*20)
   (10*24)|(15*16)
   (12*20)|(15*16)
		

Crossrefs

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    divpropQ[x_,y_]:=And[x!=y,Divisible[x,y]];
    Table[Length[Select[Subsets[facs[n],{2}],And[!Or@@divpropQ@@@Tuples[#],!Or@@divpropQ@@@Reverse/@Tuples[#]]&]],{n,100}]
  • PARI
    factorizations(n, m=n, f=List([]), z=List([])) = if(1==n, listput(z,Vec(f)); z, my(newf); fordiv(n, d, if((d>1)&&(d<=m), newf = List(f); listput(newf,d); z = factorizations(n/d, d, newf, z))); (z));
    is_proper_ndf_pair(fac1,fac2) = { for(i=1,#fac1,for(j=1,#fac2,if((fac1[i]!=fac2[j]) && (!(fac1[i]%fac2[j]) || !(fac2[j]%fac1[i])),return(0)))); (1); };
    number_of_proper_ndfpairs(z) = sum(i=1,#z,sum(j=i+1,#z,is_proper_ndf_pair(z[i],z[j])));
    A322438(n) = number_of_proper_ndfpairs(Vec(factorizations(n))); \\ Antti Karttunen, Jan 24 2025

Extensions

Data section extended up to a(144) by Antti Karttunen, Jan 24 2025

A305253 Number of connected factorizations of n into factors greater than 1 whose distinct factors are pairwise indivisible.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, May 28 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 whose distinct factors are pairwise indivisible and such that G(S) is a connected graph.

Examples

			The a(360) = 8 factorizations: (360), (4*90), (10*36), (12*30), (15*24), (18*20), (4*6*15), (6*6*10).
		

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]]}]];
    sacs[n_]:=Select[facs[n],Function[f,Length[zsm[f]]==1&&Select[Tuples[Union[f],2],UnsameQ@@#&&Divisible@@#&]=={}]]
    Table[Length[sacs[n]],{n,500}]
  • 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)); };
    A305253aux(n, m, facs) = if(1==n, is_connected(Vec(facs)), my(s=0, newfacs); fordiv(n, d, if((d>1)&&(d<=m)&&factorback(apply(x -> (x==d)||(x%d),Vec(facs))), newfacs = List(facs); listput(newfacs,d); s += A305253aux(n/d, d, newfacs))); (s));
    A305253(n) = if(1==n,0,A305253aux(n, n, List([]))); \\ Antti Karttunen, Dec 06 2018

Formula

a(n) <= A305193(n) <= A001055(n). - Antti Karttunen, Dec 06 2018

Extensions

Definition clarified by Gus Wiseman, more terms from Antti Karttunen, Dec 06 2018

A321270 Number of connected multiset partitions of a multiset whose multiplicities are the prime indices of n.

Original entry on oeis.org

1, 1, 2, 1, 3, 2, 5, 1, 5, 4, 7, 3, 11, 7, 10, 1, 15, 9, 22, 7, 19, 12, 30, 5, 22, 19, 28, 14, 42, 22, 56, 1, 33, 30, 42, 20, 77, 45
Offset: 1

Views

Author

Gus Wiseman, Nov 01 2018

Keywords

Comments

This multiset (row n of A305936) is generally not the same as the multiset of prime indices of n. For example, the prime indices of 12 are {1,1,2}, while a multiset whose multiplicities are {1,1,2} is {1,1,2,3}.

Examples

			The a(2) = 1 through a(12) = 3 connected multiset partitions:
  {{1}}  {{11}}    {{12}}  {{111}}      {{112}}    {{1111}}
         {{1}{1}}          {{1}{11}}    {{1}{12}}  {{1}{111}}
                           {{1}{1}{1}}             {{11}{11}}
                                                   {{1}{1}{11}}
                                                   {{1}{1}{1}{1}}
.
  {{123}}  {{1122}}      {{1112}}      {{11111}}          {{1123}}
           {{1}{122}}    {{1}{112}}    {{1}{1111}}        {{1}{123}}
           {{12}{12}}    {{11}{12}}    {{11}{111}}        {{12}{13}}
           {{2}{112}}    {{1}{1}{12}}  {{1}{1}{111}}
           {{1}{2}{12}}                {{1}{11}{11}}
                                       {{1}{1}{1}{11}}
                                       {{1}{1}{1}{1}{1}}
The a(18) = 9, a(27) = 28, and a(36) = 20 connected multiset partitions of {1,1,2,2,3}, {1,1,2,2,3,3}, and {1,1,2,2,3,4} respectively:
  {{1,1,2,2,3}}      {{1,1,2,2,3,3}}        {{1,1,2,2,3,4}}
  {{1},{1,2,2,3}}    {{1},{1,2,2,3,3}}      {{1},{1,2,2,3,4}}
  {{1,2},{1,2,3}}    {{1,1,2},{2,3,3}}      {{1,1,2},{2,3,4}}
  {{1,3},{1,2,2}}    {{1,1,3},{2,2,3}}      {{1,2},{1,2,3,4}}
  {{2},{1,1,2,3}}    {{1,2},{1,2,3,3}}      {{1,2,2},{1,3,4}}
  {{2,3},{1,1,2}}    {{1,2,2},{1,3,3}}      {{1,2,3},{1,2,4}}
  {{1},{1,2},{2,3}}  {{1,2,3},{1,2,3}}      {{1,3},{1,2,2,4}}
  {{1},{2},{1,2,3}}  {{1,3},{1,2,2,3}}      {{1,4},{1,2,2,3}}
  {{2},{1,2},{1,3}}  {{2},{1,1,2,3,3}}      {{2},{1,1,2,3,4}}
                     {{2,3},{1,1,2,3}}      {{2,3},{1,1,2,4}}
                     {{3},{1,1,2,2,3}}      {{2,4},{1,1,2,3}}
                     {{1},{1,2},{2,3,3}}    {{1},{1,2},{2,3,4}}
                     {{1},{1,3},{2,2,3}}    {{1},{2},{1,2,3,4}}
                     {{1},{2},{1,2,3,3}}    {{1,2},{1,3},{2,4}}
                     {{1,2},{1,3},{2,3}}    {{1,2},{1,4},{2,3}}
                     {{1},{2,3},{1,2,3}}    {{1},{2,3},{1,2,4}}
                     {{1},{3},{1,2,2,3}}    {{1},{2,4},{1,2,3}}
                     {{2},{1,2},{1,3,3}}    {{2},{1,2},{1,3,4}}
                     {{2},{1,3},{1,2,3}}    {{2},{1,3},{1,2,4}}
                     {{2},{2,3},{1,1,3}}    {{2},{1,4},{1,2,3}}
                     {{2},{3},{1,1,2,3}}
                     {{3},{1,2},{1,2,3}}
                     {{3},{1,3},{1,2,2}}
                     {{3},{2,3},{1,1,2}}
                     {{1},{2},{1,3},{2,3}}
                     {{1},{2},{3},{1,2,3}}
                     {{1},{3},{1,2},{2,3}}
                     {{2},{3},{1,2},{1,3}}
		

Crossrefs

Showing 1-10 of 12 results. Next