A303547 -id:A303547 - cp's OEIS frontend

A303546 Number of non-isomorphic aperiodic multiset partitions of weight n.

1, 3, 9, 29, 90, 285, 909, 2984, 9935, 34113, 119368, 428923, 1574223, 5915235, 22699730, 89000042, 356058539, 1453069854, 6044132793, 25612564200, 110503626702, 485161228675, 2166488899641, 9835209480533, 45370059225227

Offset: 1

Gus Wiseman, Apr 26 2018

nonn

A multiset is aperiodic if its multiplicities are relatively prime. For this sequence neither the parts nor their multiset union are required to be aperiodic, only the multiset of parts.

			Non-isomorphic representatives of the a(3) = 9 aperiodic multiset partitions are:
  {{1,1,1}}, {{1,2,2}}, {{1,2,3}},
  {{1},{1,1}}, {{1},{2,2}}, {{1},{2,3}}, {{2},{1,2}},
  {{1},{2},{2}}, {{1},{2},{3}}.

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

Cf. A000740, A000837, A001055, A007716, A007916, A049311, A100953, A255906, A269134, A275024, A293993, A301700, A303386, A303431, A303547.

a(n) = Sum_{d|n} mu(d) * A007716(n/d).

A321283 Number of non-isomorphic multiset partitions of weight n in which the part sizes are relatively prime.

Original entry on oeis.org

1, 1, 2, 7, 21, 84, 214, 895, 2607, 9591, 31134, 119313, 400950, 1574123, 5706112, 22572991, 86933012, 356058243, 1427784135, 6044132304, 25342935667, 110414556330, 481712291885, 2166488898387, 9784077216457, 45369658599779, 211869746691055, 1011161497851296, 4871413403219085

Offset: 0

Gus Wiseman, Nov 06 2018

nonn

Also the number of nonnegative integer matrices up to row and column permutations with sum of elements equal to n and no zero rows or columns, in which the row sums are relatively prime.
Also the number of non-isomorphic multiset partitions of weight n in which the multiset union of the parts is aperiodic, where a multiset is aperiodic if its multiplicities are relatively prime.
The weight of a multiset partition is the sum of sizes of its parts. Weight is generally not the same as number of vertices.

			Non-isomorphic representatives of the a(1) = 1 through a(4) = 21 multiset partitions with relatively prime part-sizes:
  {{1}}  {{1},{1}}  {{1},{1,1}}    {{1},{1,1,1}}
         {{1},{2}}  {{1},{2,2}}    {{1},{1,2,2}}
                    {{1},{2,3}}    {{1},{2,2,2}}
                    {{2},{1,2}}    {{1},{2,3,3}}
                    {{1},{1},{1}}  {{1},{2,3,4}}
                    {{1},{2},{2}}  {{2},{1,2,2}}
                    {{1},{2},{3}}  {{3},{1,2,3}}
                                   {{1},{1},{1,1}}
                                   {{1},{1},{2,2}}
                                   {{1},{1},{2,3}}
                                   {{1},{2},{1,2}}
                                   {{1},{2},{2,2}}
                                   {{1},{2},{3,3}}
                                   {{1},{2},{3,4}}
                                   {{1},{3},{2,3}}
                                   {{2},{2},{1,2}}
                                   {{1},{1},{1},{1}}
                                   {{1},{1},{2},{2}}
                                   {{1},{2},{2},{2}}
                                   {{1},{2},{3},{3}}
                                   {{1},{2},{3},{4}}
Non-isomorphic representatives of the a(1) = 1 through a(4) = 21 multiset partitions with aperiodic multiset union:
  {{1}}  {{1,2}}    {{1,2,2}}      {{1,2,2,2}}
         {{1},{2}}  {{1,2,3}}      {{1,2,3,3}}
                    {{1},{2,2}}    {{1,2,3,4}}
                    {{1},{2,3}}    {{1},{2,2,2}}
                    {{2},{1,2}}    {{1,2},{2,2}}
                    {{1},{2},{2}}  {{1},{2,3,3}}
                    {{1},{2},{3}}  {{1,2},{3,3}}
                                   {{1},{2,3,4}}
                                   {{1,2},{3,4}}
                                   {{1,3},{2,3}}
                                   {{2},{1,2,2}}
                                   {{3},{1,2,3}}
                                   {{1},{1},{2,3}}
                                   {{1},{2},{2,2}}
                                   {{1},{2},{3,3}}
                                   {{1},{2},{3,4}}
                                   {{1},{3},{2,3}}
                                   {{2},{2},{1,2}}
                                   {{1},{2},{2},{2}}
                                   {{1},{2},{3},{3}}
                                   {{1},{2},{3},{4}}

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

Cf. A000740, A000837, A007716, A007916, A100953, A301700, A303386, A303431, A303546, A303547, A320800-A320810.

\\ See links in A339645 for combinatorial species functions.
seq(n)={my(A=symGroupSeries(n)); NumUnlabeledObjsSeq(sCartProd(sExp(A), 1 + sum(d=1, n, moebius(d) * (-1 + sExp(O(x*x^n) + sum(i=1, n\d, polcoef(A,i*d)*x^(i*d)))) )))} \\ Andrew Howroyd, Jan 17 2023

\\ faster self contained program.
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, my(g=gcd(t, q[j])); g*x^(q[j]/g)) + O(x*x^k), -k))}
a(n)={if(n==0, 1, my(s=0); forpart(q=n, my(u=vector(n, t, K(q, t, n\t))); s+=permcount(q)*polcoef(sum(d=1, n, moebius(d)*exp(sum(t=1, n\d, sum(i=1, n\(t*d), u[t][i*d]*x^(i*d*t))/t, O(x*x^n)) )), n)); s/n!)} \\ Andrew Howroyd, Jan 17 2023

a(n) = A007716(n) - A320810(n). - Andrew Howroyd, Jan 17 2023

Terms a(11) and beyond from Andrew Howroyd, Jan 17 2023

A320810 Number of non-isomorphic multiset partitions of weight n whose part-sizes have a common divisor > 1.

Original entry on oeis.org

0, 2, 3, 12, 7, 84, 15, 410, 354, 3073, 56, 28300, 101, 210036, 126839, 2070047, 297, 25295952, 490, 269662769, 89071291, 3449056162, 1255, 51132696310, 400625539, 713071048480, 145126661415, 11351097702297, 4565, 199926713003444, 6842, 3460838122540969

Offset: 1

Gus Wiseman, Nov 15 2018

nonn

Also the number of nonnegative integer matrices up to row and column permutations with sum of elements equal to n and no zero rows or columns, in which the column sums are not relatively prime.
Also the number of non-isomorphic multiset partitions of weight n in which the multiset union of the parts is periodic, where a multiset is periodic if its multiplicities have a common divisor > 1.
The weight of a multiset partition is the sum of sizes of its parts. Weight is generally not the same as number of vertices.

			Non-isomorphic representatives of the a(2) = 1 through a(5) = 7 multiset partitions whose part-sizes have a common divisor:
  {{1,1}}  {{1,1,1}}  {{1,1,1,1}}    {{1,1,1,1,1}}
  {{1,2}}  {{1,2,2}}  {{1,1,2,2}}    {{1,1,2,2,2}}
           {{1,2,3}}  {{1,2,2,2}}    {{1,2,2,2,2}}
                      {{1,2,3,3}}    {{1,2,2,3,3}}
                      {{1,2,3,4}}    {{1,2,3,3,3}}
                      {{1,1},{1,1}}  {{1,2,3,4,4}}
                      {{1,1},{2,2}}  {{1,2,3,4,5}}
                      {{1,2},{1,2}}
                      {{1,2},{2,2}}
                      {{1,2},{3,3}}
                      {{1,2},{3,4}}
                      {{1,3},{2,3}}
Non-isomorphic representatives of the a(2) = 1 through a(5) = 7 multiset partitions with periodic multiset union:
  {{1,1}}    {{1,1,1}}      {{1,1,1,1}}        {{1,1,1,1,1}}
  {{1},{1}}  {{1},{1,1}}    {{1,1,2,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,2,2}}      {{1},{1,1},{1,1}}
                            {{1,1},{2,2}}      {{1},{1},{1},{1,1}}
                            {{1,2},{1,2}}      {{1},{1},{1},{1},{1}}
                            {{1},{1},{1,1}}
                            {{1},{1},{2,2}}
                            {{1},{2},{1,2}}
                            {{1},{1},{1},{1}}
                            {{1},{1},{2},{2}}

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

Cf. A007716, A018783, A047966, A120733, A303546, A303547, A305563, A319149, A319162, A319164, A319810.

\\ See links in A339645 for combinatorial species functions.
seq(n)={my(A=symGroupSeries(n));Vec(OgfSeries(sCartProd(sExp(A), -sum(d=2, n, moebius(d) * (-1 + sExp(O(x*x^n) + sum(i=1, n\d, polcoef(A,i*d)*x^(i*d)))) ))), -n)} \\ Andrew Howroyd, Jan 17 2023

a(n) = A007716(n) - A321283(n). - Andrew Howroyd, Jan 17 2023

Terms a(11) and beyond from Andrew Howroyd, Jan 17 2023

A320800 Number of non-isomorphic multiset partitions of weight n in which both the multiset union of the parts and the multiset union of the dual parts are aperiodic.

Original entry on oeis.org

1, 1, 1, 5, 14, 78, 157, 881, 2267, 9257, 28397

Offset: 0

Gus Wiseman, Nov 02 2018

The latter condition is equivalent to the parts having relatively prime sizes.
A multiset is aperiodic if its multiplicities are relatively prime.
The dual of a multiset partition has, for each vertex, one part consisting of the indices (or positions) of the parts containing that vertex, counted with multiplicity. For example, the dual of {{1,2},{2,2}} is {{1},{1,2,2}}.
The weight of a multiset partition is the sum of sizes of its parts. Weight is generally not the same as number of vertices.

			Non-isomorphic representatives of the a(1) = 1 through a(4) = 14 multiset partitions:
  {{1}}  {{1},{2}}  {{1},{2,2}}    {{1},{2,2,2}}
                    {{1},{2,3}}    {{1},{2,3,3}}
                    {{2},{1,2}}    {{1},{2,3,4}}
                    {{1},{2},{2}}  {{2},{1,2,2}}
                    {{1},{2},{3}}  {{3},{1,2,3}}
                                   {{1},{1},{2,3}}
                                   {{1},{2},{2,2}}
                                   {{1},{2},{3,3}}
                                   {{1},{2},{3,4}}
                                   {{1},{3},{2,3}}
                                   {{2},{2},{1,2}}
                                   {{1},{2},{2},{2}}
                                   {{1},{2},{3},{3}}
                                   {{1},{2},{3},{4}}

Cf. A000740, A000837, A007716, A007916, A100953, A301700, A303386, A303431, A303546, A303547, A316983, A320801-A320810.

A303709 Number of periodic factorizations of n using elements of A007916 (numbers that are not perfect powers).

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Apr 29 2018

nonn

A periodic factorization of n is a finite multiset of positive integers greater than 1 whose product is n and whose multiplicities have a common divisor greater than 1. Note that a factorization of a number that is not a perfect power (A007916) is always aperiodic (A303386), so the indices of nonzero entries of this sequence all lie at perfect powers (A001597).

			The a(900) = 5 periodic factorizations are (2*2*3*3*5*5), (2*2*15*15), (3*3*10*10), (5*5*6*6), (30*30).

Antti Karttunen, Table of n, a(n) for n = 1..100000

Cf. A000740, A000837, A001055, A001597, A007716, A007916, A052409, A052410, A281116, A303547, A303552, A303553, A303386, A303707, A303708, A303710.

radQ[n_]:=Or[n===1,GCD@@FactorInteger[n][[All,2]]===1];
facsr[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facsr[n/d],Min@@#>=d&]],{d,Select[Rest[Divisors[n]],radQ]}]];
Table[Length[Select[facsr[n],GCD@@Length/@Split[#]!=1&]],{n,200}]

gcd_of_multiplicities(lista) = { my(u=length(lista)); if(u<2, u, my(g=0, pe = lista[1], j=1); for(i=2,u,if(lista[i]==pe, j++, g = gcd(j,g); j=1; pe = lista[i])); gcd(g,j)); }; \\ the supplied lista (newfacs) should be monotonic
A303709(n, m=n, facs=List([])) = if(1==n, (1!=gcd_of_multiplicities(facs)), my(s=0, newfacs); fordiv(n, d, if((d>1)&&(d<=m)&&!ispower(d), newfacs = List(facs); listput(newfacs,d); s += A303709(n/d, d, newfacs))); (s)); \\ Antti Karttunen, Dec 06 2018

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

Changed a(1) to 1 by Gus Wiseman, Dec 06 2018

A320802 Number of non-isomorphic aperiodic multiset partitions of weight n whose dual is also an aperiodic multiset partition.

Original entry on oeis.org

1, 1, 2, 8, 26, 89, 274, 908, 2955, 9926, 34021, 119367, 428612, 1574222, 5914324, 22699632, 88997058, 356058538, 1453059643, 6044132792, 25612530061, 110503625785, 485161109305, 2166488899640, 9835209048655, 45370059225137, 212582814591083, 1011306624492831

Offset: 0

Gus Wiseman, Nov 06 2018

nonn

Also the number of nonnegative integer matrices with sum of entries equal to n and no zero rows or columns where the multiset of rows and the multiset of columns are both aperiodic, up to row and column permutations.
A multiset is aperiodic if its multiplicities are relatively prime.
The dual of a multiset partition has, for each vertex, one part consisting of the indices (or positions) of the parts containing that vertex, counted with multiplicity. For example, the dual of {{1,2},{2,2}} is {{1},{1,2,2}}.
The weight of a multiset partition is the sum of sizes of its parts. Weight is generally not the same as number of vertices.
Also the number of non-isomorphic aperiodic multiset partitions of weight n whose parts have relatively prime periods, where the period of a multiset is the GCD of its multiplicities.

			Non-isomorphic representatives of the a(1) = 1 through a(4) = 26 multiset partitions:
  {{1}}  {{1,1}}    {{1,1,1}}      {{1,1,1,1}}
         {{1},{2}}  {{1,2,2}}      {{1,2,2,2}}
                    {{1},{1,1}}    {{1,2,3,3}}
                    {{1},{2,2}}    {{1},{1,1,1}}
                    {{1},{2,3}}    {{1},{1,2,2}}
                    {{2},{1,2}}    {{1,1},{2,2}}
                    {{1},{2},{2}}  {{1},{2,2,2}}
                    {{1},{2},{3}}  {{1,2},{2,2}}
                                   {{1},{2,3,3}}
                                   {{1,2},{3,3}}
                                   {{1},{2,3,4}}
                                   {{1,3},{2,3}}
                                   {{2},{1,2,2}}
                                   {{3},{1,2,3}}
                                   {{1},{1},{1,1}}
                                   {{1},{1},{2,2}}
                                   {{1},{1},{2,3}}
                                   {{1},{2},{1,2}}
                                   {{1},{2},{2,2}}
                                   {{1},{2},{3,3}}
                                   {{1},{2},{3,4}}
                                   {{1},{3},{2,3}}
                                   {{2},{2},{1,2}}
                                   {{1},{2},{2},{2}}
                                   {{1},{2},{3},{3}}
                                   {{1},{2},{3},{4}}

Jinyuan Wang, Table of n, a(n) for n = 0..50

Cf. A000740, A000837, A007716, A007916, A100953, A301700, A303386, A303431, A303546, A303547, A316983, A320800-A320810.

Second Moebius transform of A007716, or Moebius transform of A303546, where the Moebius transform of a sequence b is a(n) = Sum_{d|n} mu(d) * b(n/d).

a(26)-a(27) from Jinyuan Wang, Jun 27 2020

A321390 Third Moebius transform of A007716. Number of non-isomorphic aperiodic multiset partitions of weight n whose parts have relatively prime periods and whose dual is also an aperiodic multiset partition.

Original entry on oeis.org

1, 1, 1, 7, 24, 88, 265, 907, 2929, 9918, 33931, 119366, 428314, 1574221, 5913415, 22699536, 88994103, 356058537, 1453049451, 6044132791, 25612496016, 110503624870, 485160989937, 2166488899639, 9835208617114, 45370059225048

Offset: 0

Gus Wiseman, Nov 08 2018

nonn

The Moebius transform c of a sequence b is c(n) = Sum_{d|n} mu(d) * b(n/d).
Also the number of nonnegative integer matrices with sum of entries equal to n and no zero rows or columns where the multiset of rows and the multiset of columns are both aperiodic and the nonzero entries are relatively prime, up to row and column permutations.
A multiset is aperiodic if its multiplicities are relatively prime. The period of a multiset is the GCD of its multiplicities.
The dual of a multiset partition has, for each vertex, one part consisting of the indices (or positions) of the parts containing that vertex, counted with multiplicity. For example, the dual of {{1,2},{2,2}} is {{1},{1,2,2}}.
The weight of a multiset partition is the sum of sizes of its parts. Weight is generally not the same as number of vertices.

			Non-isomorphic representatives of the a(1) = 1 through a(4) = 24 multiset partitions:
  {{1}}  {{1},{2}}  {{1,2,2}}      {{1,2,2,2}}
                    {{1},{1,1}}    {{1,2,3,3}}
                    {{1},{2,2}}    {{1},{1,1,1}}
                    {{1},{2,3}}    {{1},{1,2,2}}
                    {{2},{1,2}}    {{1},{2,2,2}}
                    {{1},{2},{2}}  {{1,2},{2,2}}
                    {{1},{2},{3}}  {{1},{2,3,3}}
                                   {{1,2},{3,3}}
                                   {{1},{2,3,4}}
                                   {{1,3},{2,3}}
                                   {{2},{1,2,2}}
                                   {{3},{1,2,3}}
                                   {{1},{1},{1,1}}
                                   {{1},{1},{2,2}}
                                   {{1},{1},{2,3}}
                                   {{1},{2},{1,2}}
                                   {{1},{2},{2,2}}
                                   {{1},{2},{3,3}}
                                   {{1},{2},{3,4}}
                                   {{1},{3},{2,3}}
                                   {{2},{2},{1,2}}
                                   {{1},{2},{2},{2}}
                                   {{1},{2},{3},{3}}
                                   {{1},{2},{3},{4}}

Cf. A000740, A000837, A007716, A007916, A100953, A301700, A303386, A303431, A303546, A303547, A316983, A320800-A320810.

A108572 Number of partitions of n which, as multisets, are nontrivial repetitions of a multiset.

Original entry on oeis.org

0, 0, 0, 1, 0, 3, 0, 4, 2, 7, 0, 13, 0, 15, 8, 21, 0, 37, 0, 44, 16, 56, 0, 93, 6, 101, 29, 137, 0, 217, 0, 230, 57, 297, 20, 450, 0, 490, 102, 643, 0, 918, 0, 1004, 202, 1255, 0, 1783, 14, 1992, 298, 2438, 0, 3364, 61, 3734, 491, 4565, 0, 6251, 0, 6842, 818

Offset: 1

Len Smiley, Jul 25 2005

nonn

The singleton and the all-ones partitions are ignored, so that a(n)=0 if n is prime. If a partition is listed as m_1^am_2^bm_3^c..., then it is counted exactly when gcd(a,b,c,...)>1. These are equinumerous (conjugate) with those partitions for which gcd(m_1,m_2,...)>1 (less 1, the singleton), hence the formula.

			a(25) = 6: 1^(15)2^5 = 5{1, 1, 1, 2}, 1^52^(10) = 5{1, 2, 2}, 1^(10)3^5 = 5{3, 1, 1}, 2^53^5 = 5{3, 2}, 1^44^4 = 5{4, 1}, 5^5 = 5{5}.
Note that A000041(25)=P(25)=1958, only 6 of which satisfy the criterion.

Cf. A000837, A018783, A047966, A100953, A303386, A303547, A303553, A319149, A319162, A319164, A319810.

with(combinat):PartMulti:=proc(n::nonnegint) local count,a,i,j,b,m,k,part_vec;
bigcount:=0; if isprime(n) then return(bigcount) else ps:=partition(n); b:=nops(ps);
for m from 2 to b-1 do p:=ps[m]; a:=nops(p); part_vec:=array(1..n);
for k from 1 to n do part_vec[k]:=0 od;
for i from 1 to a do j:=p[i]; part_vec[j]:=part_vec[j]+1 od;
g:=0; for j from 1 to n do g:=igcd(g,part_vec[j]) od;
if g>1 then bigcount:=bigcount+1 fi od; return(bigcount) end if end proc;
seq(PartMulti(q),q=1..49);

Table[Length[Select[IntegerPartitions[n],And[Length[#]1]&]],{n,20}] (* Gus Wiseman, Dec 06 2018 *)

a(n) = A018783(n)-1, n>1. - Vladeta Jovovic, Jul 28 2005

More terms from Gus Wiseman, Dec 06 2018

A303552 Number of periodic multisets of compositions of total weight n.

Original entry on oeis.org

0, 1, 1, 3, 1, 9, 1, 18, 7, 44, 1, 119, 1, 246, 48, 585, 1, 1470, 1, 3248, 250, 7535, 1, 18114, 42, 40593, 1373, 93726, 1, 218665, 1, 493735, 7539, 1127981, 285, 2587962, 1, 5841445, 40597, 13244166, 1, 30047413, 1, 67604050, 216745, 152258273, 1, 342747130

Offset: 1

Gus Wiseman, Apr 26 2018

nonn

A multiset is periodic if its multiplicities have a common divisor greater than 1.

			The a(6) = 9 periodic multisets of compositions are:
{1,1,1,1,1,1},
{1,1,2,2}, {1,1,11,11},
{2,2,2}, {11,11,11},
{3,3}, {21,21}, {12,12}, {111,111}.

Cf. A000740, A000837, A007716, A007916, A034691, A100953, A255906, A269134, A301700, A303386, A303431, A303547, A303551.

nn=60;
ser=Product[1/(1-x^n)^2^(n-1),{n,nn}]
Table[SeriesCoefficient[ser,{x,0,n}]-Sum[MoebiusMu[d]*SeriesCoefficient[ser,{x,0,n/d}],{d,Divisors[n]}],{n,1,nn}]

A303553 Number of periodic factorizations of n > 1 into positive factors greater than 1; a(1) = 1 by convention.

Original entry on oeis.org

1, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 2, 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, 4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1

Offset: 1

Gus Wiseman, Apr 26 2018

nonn

A multiset is periodic if its multiplicities have a common divisor greater than 1.

			The a(64)  = 4 periodic factorizations are (2*2*2*2*2*2), (2*2*4*4), (4*4*4), (8*8).
The a(144) = 4 periodic factorizations are (2*2*2*2*3*3), (2*2*6*6), (3*3*4*4), (12*12).
The a(256) = 5 periodic factorizations are (2*2*2*2*2*2*2*2), (2*2*2*2*4*4), (2*2*8*8), (4*4*4*4), (16*16).
The a(576) = 7 periodic factorizations are (2*2*2*2*2*2*3*3), (2*2*2*2*6*6), (2*2*3*3*4*4), (2*2*12*12), (3*3*8*8), (4*4*6*6), (24*24).

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

Cf. A000740, A000837, A001055, A007916, A018783(n) = a(2^n), A052409, A052410, A275870, A303386, A303431, A303547, A303552, A303709.

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],GCD@@Length/@Split[#]>1&]],{n,2,100}]

gcd_of_multiplicities(lista) = { my(u=length(lista)); if(u<2, u, my(g=0, pe = lista[1], j=1); for(i=2,u,if(lista[i]==pe, j++, g = gcd(j,g); j=1; pe = lista[i])); gcd(g,j)); }; \\ the supplied lista (newfacs) should be monotonic
A303553(n, m=n, facs=List([])) = if(1==n, (gcd_of_multiplicities(facs)!=1), my(s=0, newfacs); fordiv(n, d, if((d>1)&&(d<=m), newfacs = List(facs); listput(newfacs,d); s += A303553(n/d, d, newfacs))); (s)); \\ Antti Karttunen, Dec 06 2018

a(n) >= A303709(n). - Antti Karttunen, Dec 06 2018

a(1) = 1 prepended by Antti Karttunen, Dec 06 2018

cp's OEIS Frontend

A303546 Number of non-isomorphic aperiodic multiset partitions of weight n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Formula

A321283 Number of non-isomorphic multiset partitions of weight n in which the part sizes are relatively prime.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

PARI

PARI

Formula

Extensions

A320810 Number of non-isomorphic multiset partitions of weight n whose part-sizes have a common divisor > 1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

PARI

Formula

Extensions

A320800 Number of non-isomorphic multiset partitions of weight n in which both the multiset union of the parts and the multiset union of the dual parts are aperiodic.

Views

Author

Keywords

Comments

Examples

Crossrefs

A303709 Number of periodic factorizations of n using elements of A007916 (numbers that are not perfect powers).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A320802 Number of non-isomorphic aperiodic multiset partitions of weight n whose dual is also an aperiodic multiset partition.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Formula

Extensions

A321390 Third Moebius transform of A007716. Number of non-isomorphic aperiodic multiset partitions of weight n whose parts have relatively prime periods and whose dual is also an aperiodic multiset partition.

Views

Author

Keywords

Comments

Examples

Crossrefs

A108572 Number of partitions of n which, as multisets, are nontrivial repetitions of a multiset.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs