A317748 -id:A317748 - cp's OEIS frontend

A281116 Number of factorizations of n>=2 into factors greater than 1 with no common divisor other than 1 (a(1)=0 by convention).

0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 2, 0, 1, 1, 0, 0, 2, 0, 2, 1, 1, 0, 3, 0, 1, 0, 2, 0, 4, 0, 0, 1, 1, 1, 5, 0, 1, 1, 3, 0, 4, 0, 2, 2, 1, 0, 5, 0, 2, 1, 2, 0, 3, 1, 3, 1, 1, 0, 8, 0, 1, 2, 0, 1, 4, 0, 2, 1, 4, 0, 9, 0, 1, 2, 2, 1, 4, 0, 5, 0, 1, 0, 8, 1, 1, 1, 3, 0, 8, 1, 2, 1, 1, 1, 7, 0, 2, 2, 5

Offset: 1

Gus Wiseman, Jan 15 2017

nonn

Let (e1, e2, ..., ek) be a prime-signature of n (that is, n = p^e1 * q^e2 * ... * r^ek for some primes, p, q, ..., r). Then a(n) is the number of ways of partitioning multiset {e1 x 1, e2 x 2, ..., ek x k} into multisets such that none of the numbers 1 .. k is present in all member multisets of that set partition. - Antti Karttunen, Sep 08 2018

			a(6)=1:  (2*3)
a(12)=2; (2*2*3)       (3*4)
a(24)=3: (2*2*2*3)     (2*3*4)     (3*8)
a(30)=4: (2*3*5)       (2*15)      (3*10)    (5*6)
a(36)=5: (2*2*3*3)     (2*2*9)     (2*3*6)   (3*3*4)   (4*9)
a(96)=7: (2*2*2*2*2*3) (2*2*2*3*4) (2*2*3*8) (2*3*4*4) (2*3*16) (3*4*8) (3*32).

Antti Karttunen, Table of n, a(n) for n = 1..65537
Index entries for sequences computed from exponents in factorization of n

Cf. A001055, A007916, A089233, A162247, A259936, A281113, A317751.

First column of A317748.

postfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[postfacs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
Table[Length[Select[postfacs[n],GCD@@#===1&]],{n,2,100}]

A281116(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 += A281116(n/d, d, newfacs))); (s)); \\ Antti Karttunen, Sep 08 2018

Term a(1) = 0 prepended by Antti Karttunen, Sep 08 2018

A317757 Number of non-isomorphic multiset partitions of size n such that the blocks have empty intersection.

Original entry on oeis.org

1, 0, 1, 4, 17, 56, 205, 690, 2446, 8506, 30429, 109449, 402486, 1501424, 5714194, 22132604, 87383864, 351373406, 1439320606, 6003166059, 25488902820, 110125079184, 483987225922, 2162799298162, 9823464989574, 45332196378784, 212459227340403, 1010898241558627, 4881398739414159

Offset: 0

Gus Wiseman, Aug 06 2018

nonn

			Non-isomorphic representatives of the a(4) = 17 multiset partitions:
  {1}{234},{2}{111},{2}{113},{11}{22},{11}{23},{12}{34},
  {1}{1}{22},{1}{1}{23},{1}{2}{11},{1}{2}{12},{1}{2}{13},{1}{2}{34},{2}{3}{11},
  {1}{1}{1}{2},{1}{1}{2}{2},{1}{1}{2}{3},{1}{2}{3}{4}.

Andrew Howroyd, Table of n, a(n) for n = 0..50

Cf. A007716, A035310, A255906, A281116, A317073, A317533.
Cf. A317748, A317751, A317752, A317755.
Cf. A319077, A319748, A319755, A319778, A319781, A319790.

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]]]];
strnorm[n_]:=Flatten[MapIndexed[Table[#2,{#1}]&,#]]&/@IntegerPartitions[n];
sysnorm[m_]:=If[Union@@m!=Range[Max@@Flatten[m]],sysnorm[m/.Rule@@@Table[{(Union@@m)[[i]],i},{i,Length[Union@@m]}]],First[Sort[sysnorm[m,1]]]];sysnorm[m_,aft_]:=If[Length[Union@@m]<=aft,{m},With[{mx=Table[Count[m,i,{2}],{i,Select[Union@@m,#>=aft&]}]},Union@@(sysnorm[#,aft+1]&/@Union[Table[Map[Sort,m/.{par+aft-1->aft,aft->par+aft-1},{0,1}],{par,First/@Position[mx,Max[mx]]}]])]];
Table[Length[Union[sysnorm/@Join@@Table[Select[mps[m],Intersection@@#=={}&],{m,strnorm[n]}]]],{n,6}]

EulerT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, n, 1/n))))-1, -#v)}
permcount(v) = {my(m=1, s=0, k=0, t); for(i=1, #v, t=v[i]; k=if(i>1&&t==v[i-1], k+1, 1); m*=t*k; s+=t); s!/m}
K(q, t, k)={EulerT(Vec(sum(j=1, #q, gcd(t, q[j])*x^lcm(t, q[j])) + O(x*x^k), -k))}
R(q, n)={vector(n, t, x*Ser(K(q, t, n)/t))}
a(n)={my(s=0); forpart(q=n, my(f=prod(i=1, #q, 1 - x^q[i]), u=R(q,n)); s+=permcount(q)*sum(k=0, n, my(c=polcoef(f,k)); if(c, c*polcoef(exp(sum(t=1, n\(k+1), x^(t*k)*u[t], O(x*x^n) ))/if(k,1-x^k,1), n))) ); s/n!} \\ Andrew Howroyd, May 30 2023

a(8)-a(10) from Gus Wiseman, Sep 27 2018

a(0)=1 prepended and terms a(11) and beyond from Andrew Howroyd, May 30 2023

A317752 Number of multiset partitions of normal multisets of size n such that the blocks have empty intersection.

Original entry on oeis.org

0, 1, 8, 49, 305, 1984, 13686, 100124, 776885, 6386677, 55532358, 509549386, 4921352952, 49899820572, 529807799836, 5876162077537, 67928460444139, 816764249684450, 10195486840926032, 131896905499007474, 1765587483656124106, 24419774819813602870

Offset: 1

Gus Wiseman, Aug 06 2018

nonn

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

			The a(3) = 8 multiset partitions with empty intersection:
  {{2},{1,1}}
  {{1},{2,2}}
  {{1},{2,3}}
  {{2},{1,3}}
  {{3},{1,2}}
  {{1},{1},{2}}
  {{1},{2},{2}}
  {{1},{2},{3}}

Andrew Howroyd, Table of n, a(n) for n = 1..100

Cf. A007716, A255903, A255906, A317073, A281116, A317077, A317532, A317533.

Cf. A317748, A317751, A317755, A317757, A317776.

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]]]];
allnorm[n_]:=Function[s,Array[Count[s,y_/;y<=#]+1&,n]]/@Subsets[Range[n-1]+1];
Table[Length[Join@@Table[Select[mps[m],Intersection@@#=={}&],{m,allnorm[n]}]],{n,6}]

P(n,k)={1/prod(i=1, n, (1 - x^i*y + O(x*x^n))^binomial(k+i-1, k-1))}
R(n,k)={my(p=P(n,k), q=p/(1-y+O(y*y^n))); Vec(sum(i=2, n, polcoef(p,i,y) + polcoef(q,i,y)*sum(j=1, n\i, (-1)^j*binomial(k,j)*x^(i*j))), -n)}
seq(n)={sum(k=2, n, R(n, k)*sum(r=k, n, binomial(r, k)*(-1)^(r-k)) )} \\ Andrew Howroyd, Feb 05 2021

Terms a(9) and beyond from Andrew Howroyd, Feb 05 2021

A317755 Number of multiset partitions of strongly normal multisets of size n such that the blocks have empty intersection.

Original entry on oeis.org

0, 1, 6, 30, 130, 629, 2930, 15019, 78224, 438626, 2548481

Offset: 1

Gus Wiseman, Aug 06 2018

A multiset is strongly normal if it spans an initial interval of positive integers with weakly decreasing multiplicities.

			The a(3) = 6 strongly normal multiset partitions with empty intersection:
  {{2},{1,1}}
  {{1},{2,3}}
  {{2},{1,3}}
  {{3},{1,2}}
  {{1},{1},{2}}
  {{1},{2},{3}}

Cf. A007716, A035310, A255906, A317073, A281116, A317077, A317078, A317532, A317533.

Cf. A317748, A317751, A317752, A317757, A317775.

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]]]];
strnorm[n_]:=Flatten[MapIndexed[Table[#2,{#1}]&,#]]&/@IntegerPartitions[n];
Table[Length[Join@@Table[Select[mps[m],Intersection@@#=={}&],{m,strnorm[n]}]],{n,6}]

a(10)-a(11) from Robert Price, May 08 2021

A318721 Number of strict relatively prime factorizations of n.

Original entry on oeis.org

1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 1, 1, 0, 0, 1, 0, 1, 1, 1, 0, 2, 0, 1, 0, 1, 0, 4, 0, 0, 1, 1, 1, 2, 0, 1, 1, 2, 0, 4, 0, 1, 1, 1, 0, 2, 0, 1, 1, 1, 0, 2, 1, 2, 1, 1, 0, 6, 0, 1, 1, 0, 1, 4, 0, 1, 1, 4, 0, 4, 0, 1, 1, 1, 1, 4, 0, 2, 0, 1, 0, 6, 1, 1, 1

Offset: 1

Gus Wiseman, Sep 02 2018

nonn

a(1) could be either 0 or 1.

			The a(60) = 6 factorizations are (3*20), (4*15), (5*12), (2*3*10), (2*5*6), (3*4*5).

Cf. A001055, A007716, A045778, A078374, A281116, A317748, A318720.

strfacs[n_]:=If[n<=1,{{}},Join@@Table[(Prepend[#1,d]&)/@Select[strfacs[n/d],Min@@#1>d&],{d,Rest[Divisors[n]]}]];
Table[Length[Select[strfacs[n],GCD@@#==1&]],{n,50}]

A318720 Numbers k such that there exists a strict relatively prime factorization of k in which no pair of factors is relatively prime.

Original entry on oeis.org

900, 1764, 1800, 2700, 3528, 3600, 4356, 4500, 4900, 5292, 5400, 6084, 6300, 7056, 7200, 8100, 8712, 8820, 9000, 9800, 9900, 10404, 10584, 10800, 11025, 11700, 12100, 12168, 12348, 12600, 12996, 13068, 13500, 14112, 14400, 14700, 15300, 15876, 16200, 16900

Offset: 1

Gus Wiseman, Sep 02 2018

nonn

From Amiram Eldar, Nov 01 2020: (Start)
Also, numbers with more than two non-unitary prime divisors, i.e., numbers k such that A056170(k) > 2, or equivalently, numbers divisible by the squares of three distinct primes.
The complement of the union of A005117, A190641 and A338539.
The asymptotic density of this sequence is 1 - 6/Pi^2 - (6/Pi^2)*A154945 - (3/Pi^2)*(A154945^2 - A324833) = 0.0033907041... (End)

			900 is in the sequence because the factorization 900 = (6*10*15) is relatively prime (since the GCD of (6,10,15) is 1) but each of the pairs (6,10), (6,15), (10,15) has a common divisor > 1. Larger examples are:
1800 = (6*15*20) = (10*12*15).
9900 = (6*10*165) = (6*15*110) = (10*15*66).
5400 = (6*20*45) = (10*12*45) = (10*15*36) = (15*18*20).
60 is not in the sequence because all its possible factorizations (4 * 15, 3 * 4 * 5, etc.) contain at least one pair that is coprime, if not more than one prime.

Amiram Eldar, Table of n, a(n) for n = 1..10000

Cf. A001055, A001221, A001222, A007716, A045778, A051185, A078374, A281116, A303140, A303283, A305843, A305854, A317748, A318715, A318717, A318721.

Cf. A005117, A036785, A056170, A154945, A190641, A324833, A338539.

strfacs[n_] := If[n <= 1, {{}}, Join@@Table[(Prepend[#1, d] &)/@Select[strfacs[n/d], Min@@#1 > d &], {d, Rest[Divisors[n]]}]]; Select[Range[10000], Function[n, Select[strfacs[n], And[GCD@@# == 1, And@@(GCD[##] > 1 &)@@@Select[Tuples[#, 2], Less@@# &]] &] != {}]]
Select[Range[20000], Count[FactorInteger[#][[;;,2]], ?(#1 > 1 &)] > 2 &] (* _Amiram Eldar, Nov 01 2020 *)

A317751 Number of divisors d of n such that there exists a factorization of n into factors > 1 with GCD d.

Original entry on oeis.org

0, 1, 1, 2, 1, 2, 1, 2, 2, 2, 1, 3, 1, 2, 2, 3, 1, 3, 1, 3, 2, 2, 1, 3, 2, 2, 2, 3, 1, 2, 1, 3, 2, 2, 2, 5, 1, 2, 2, 3, 1, 2, 1, 3, 3, 2, 1, 4, 2, 3, 2, 3, 1, 3, 2, 3, 2, 2, 1, 3, 1, 2, 3, 4, 2, 2, 1, 3, 2, 2, 1, 5, 1, 2, 3, 3, 2, 2, 1, 4, 3, 2, 1, 3, 2, 2, 2, 3, 1, 3, 2, 3, 2, 2, 2, 4, 1, 3, 3, 5, 1, 2, 1, 3, 2

Offset: 1

Gus Wiseman, Aug 06 2018

nonn

Also the number of distinct possible GCDs of factorizations of n into factors > 1.
Also the number of nonzero terms in row n of A317748.
a(prime^n) = A008619(n).
If n is squarefree and composite, a(n) = 2.

			The divisors of 36 that are possible GCDs of factorizations of 36 are {1, 2, 3, 6, 36}, so a(36) = 5.

Antti Karttunen, Table of n, a(n) for n = 1..65537
Index entries for sequences computed from exponents in factorization of n

Cf. A000005, A000837, A001055, A014963, A045778, A050370, A162247, A281116, A289509.

Cf. A317748, A317752, A317755, A317757.

facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
goc[n_,m_]:=Length[Select[facs[n],And[And@@(Divisible[#,m]&/@#),GCD@@(#/m)==1]&]];
Table[Length[Select[Divisors[n],goc[n,#]!=0&]],{n,100}]

A317751aux(n, m, facs, gcds) = if(1==n, setunion([gcd(Vec(facs))],gcds), my(newfacs); fordiv(n, d, if((d>1)&&(d<=m), newfacs = List(facs); listput(newfacs,d); gcds = setunion(gcds,A317751aux(n/d, d, newfacs, gcds)))); (gcds));
A317751(n) = if(1==n,0,length(A317751aux(n, n, List([]), Set([])))); \\ Antti Karttunen, Sep 08 2018

More terms from Antti Karttunen, Sep 08 2018

A318749 Number of pairwise relatively nonprime strict factorizations of n (no two factors are coprime).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 1, 1, 1, 3, 1, 1, 2, 2, 1, 1, 1, 3, 1, 1, 1, 3, 1, 1, 1, 3, 1, 1, 1, 2, 2, 1, 1, 5, 1, 2, 1, 2, 1, 3, 1, 3, 1, 1, 1, 3, 1, 1, 2, 4, 1, 1, 1, 2, 1, 1, 1, 5, 1, 1, 2, 2, 1, 1, 1, 5, 2, 1, 1, 3, 1, 1, 1, 3, 1, 3, 1, 2, 1, 1, 1, 7, 1, 2, 2, 3, 1, 1, 1, 3, 1

Offset: 1

Gus Wiseman, Sep 02 2018

nonn

a(n) depends only on prime signature of n (cf. A025487). - Antti Karttunen, Oct 08 2018

			The a(96) = 7 factorizations are (96), (2*48), (4*24), (6*16), (8*12), (2*4*12), (2*6*8).
The a(480) = 18 factorizations:
  (480)
  (2*240) (4*120) (6*80) (8*60) (10*48) (12*40) (16*30) (20*24)
  (2*4*60) (2*6*40) (2*8*30) (2*10*24) (2*12*20) (4*6*20) (4*10*12) (6*8*10)
  (2*4*6*10)

Antti Karttunen, Table of n, a(n) for n = 1..10000
Antti Karttunen, Data supplement: n, a(n) computed for n = 1..100000
Index entries for sequences computed from exponents in factorization of n

Cf. A001055, A045778, A051185, A281116, A303282, A305843, A305854, A317748, A318715, A318717, A318721.

strfacs[n_]:=If[n<=1,{{}},Join@@Table[(Prepend[#1,d]&)/@Select[strfacs[n/d],Min@@#1>d&],{d,Rest[Divisors[n]]}]];
Table[Length[Select[strfacs[n],And@@(GCD[##]>1&)@@@Select[Tuples[#,2],Less@@#&]&]],{n,50}]

A318749(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 += A318749(n/d, d-1, newfacs))); (s)); \\ Antti Karttunen, Oct 08 2018

More terms from Antti Karttunen, Oct 08 2018

cp's OEIS Frontend

A281116 Number of factorizations of n>=2 into factors greater than 1 with no common divisor other than 1 (a(1)=0 by convention).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Extensions

A317757 Number of non-isomorphic multiset partitions of size n such that the blocks have empty intersection.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Extensions

A317752 Number of multiset partitions of normal multisets of size n such that the blocks have empty intersection.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Extensions

A317755 Number of multiset partitions of strongly normal multisets of size n such that the blocks have empty intersection.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Extensions

A318721 Number of strict relatively prime factorizations of n.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A318720 Numbers k such that there exists a strict relatively prime factorization of k in which no pair of factors is relatively prime.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A317751 Number of divisors d of n such that there exists a factorization of n into factors > 1 with GCD d.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Extensions

A318749 Number of pairwise relatively nonprime strict factorizations of n (no two factors are coprime).

Views

Author