A050326 -id:A050326 - cp's OEIS frontend

A085987 Product of exactly four primes, three of which are distinct (p^2qr).

60, 84, 90, 126, 132, 140, 150, 156, 198, 204, 220, 228, 234, 260, 276, 294, 306, 308, 315, 340, 342, 348, 350, 364, 372, 380, 414, 444, 460, 476, 490, 492, 495, 516, 522, 525, 532, 550, 558, 564, 572, 580, 585, 620, 636, 644, 650, 666, 693, 708, 726

Offset: 1

Alford Arnold, Jul 08 2003

nonn

A014613 is completely determined by A030514, A065036, A085986, A085987 and A046386 since p(4) = 5. (cf. A000041). More generally, the first term of sequences which completely determine the k-almost primes can be found in A036035 (a resorted version of A025487).

A050326(a(n)) = 4. - Reinhard Zumkeller, May 03 2013

			a(1) = 60 since 60 = 2*2*3*5 and has three distinct prime factors.

Cf. A001248, A006881, A030078, A054753, A007304, A050997, A046387, A036035, A086974.

Subsequence of A014613, A307341, A178212.

f[n_]:=Sort[Last/@FactorInteger[n]]=={1,1,2}; Select[Range[2000], f] (* Vladimir Joseph Stephan Orlovsky, May 03 2011 *)
pefp[{a_,b_,c_}]:={a^2 b c,a b^2 c,a b c^2}; Module[{upto=800},Select[ Flatten[ pefp/@Subsets[Prime[Range[PrimePi[upto/6]]],{3}]]//Union,#<= upto&]] (* Harvey P. Dale, Oct 02 2018 *)

list(lim)=my(v=List(),t,x,y,z);forprime(p=2,lim^(1/4),t=lim\p^2;forprime(q=p+1,sqrtint(t),forprime(r=q+1,t\q,x=p^2*q*r;y=p*q^2*r;listput(v,x);if(y<=lim,listput(v,y);z=p*q*r^2;if(z<=lim,listput(v,z))))));vecsort(Vec(v)) \\ Charles R Greathouse IV, Jul 15 2011

is(n)=vecsort(factor(n)[,2]~)==[1,1,2] \\ Charles R Greathouse IV, Oct 19 2015

from math import isqrt
from sympy import primepi, primerange, integer_nthroot
def A085987(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        kmin = kmax >> 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x): return n+x+sum((t:=primepi(s:=isqrt(y:=x//r**2)))+(t*(t-1)>>1)-sum(primepi(y//k) for k in primerange(1, s+1)) for r in primerange(isqrt(x)+1))+sum(primepi(x//p**3) for p in primerange(integer_nthroot(x,3)[0]+1))-primepi(integer_nthroot(x,4)[0])
    return bisection(f,n,n) # Chai Wah Wu, Mar 27 2025

More terms from Reinhard Zumkeller, Jul 25 2003

A381441 Number of multisets that can be obtained by partitioning the prime indices of n into a set of sets (set system) and taking their sums.

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Mar 06 2025

nonn

First differs from A050326 at a(210) = 13, A050326(210) = 15. This comes from the set systems {{3},{1,2,4}} and {{1,2},{3,4}}, and from {{4},{1,2,3}} and {{1,3},{2,4}}.
Also the number of multisets that can be obtained by taking the sums of prime indices of each factor in a strict factorization of n into squarefree numbers > 1.
A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.
A multiset partition can be regarded as an arrow in the ranked poset of integer partitions. For example, we have {{1},{1,2},{1,3},{1,2,3}}: {1,1,1,1,2,2,3,3} -> {1,3,4,6}, or (33221111) -> (6431) (depending on notation).
Sets of sets are generally not transitive. For example, we have arrows: {{1},{1,2}}: {1,1,2} -> {1,3} and {{1,3}}: {1,3} -> {4}, but there is no set of sets {1,1,2} -> {4}.

			The prime indices of 60 are {1,1,2,3}, with partitions into sets of sets:
  {{1},{1,2,3}}
  {{1,2},{1,3}}
  {{1},{2},{1,3}}
  {{1},{3},{1,2}}
with block-sums: {1,6}, {3,4}, {1,2,4}, {1,3,3}, which are all different, so a(60) = 4.

Before taking sums we had A050326, non-strict A050320.
Positions of 0 are A293243.
Positions of 1 are A293511.
This is the strict version of A381078 (lower A381454).
For distinct block-sums (instead of blocks) we have A381634, before sums A381633.
Other multiset partitions of prime indices:
- For multisets of constant multisets (A000688) see A381455 (upper), A381453 (lower).
- For multiset partitions (A001055) see A317141 (upper), A300383 (lower).
- For strict multiset partitions (A045778) see A381452.
- For sets of constant multisets (A050361) see A381715.
- For strict multiset partitions with distinct sums (A321469) see A381637.
- For sets of constant multisets with distinct sums (A381635) see A381716, A381636.
More on set systems: A050342, A116539, A279785, A296120, A318361.
A000041 counts integer partitions, strict A000009.
A000040 lists the primes, differences A001223.
A003963 gives product of prime indices.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798.
A122111 represents conjugation in terms of Heinz numbers.
A265947 counts refinement-ordered pairs of integer partitions.
Cf. A000720, A001222, A002846, A005117, A066328, A213242, A213385, A213427, A299202, A300385, A317142.

prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
Table[Length[Union[Sort[Total/@prix/@#]&/@Select[facs[n],UnsameQ@@#&&And@@SquareFreeQ/@#&]]],{n,100}]

a(A002110(n)) = A066723(n).

A381718 Number of normal multiset partitions of weight n into sets with distinct sums.

Original entry on oeis.org

1, 1, 2, 6, 23, 106, 549, 3184, 20353, 141615, 1063399, 8554800, 73281988, 665141182, 6369920854, 64133095134, 676690490875, 7462023572238, 85786458777923, 1025956348473929, 12739037494941490

Offset: 0

Gus Wiseman, Mar 26 2025

nonn

We call a multiset or multiset partition normal iff it covers an initial interval of positive integers. The weight of a multiset partition is the sum of sizes of its blocks.

			The a(1) = 1 through a(3) = 6 multiset partitions:
  {{1}}  {{1,2}}    {{1,2,3}}
         {{1},{2}}  {{1},{1,2}}
                    {{1},{2,3}}
                    {{2},{1,2}}
                    {{2},{1,3}}
                    {{1},{2},{3}}
The a(4) = 23 factorizations:
  2*3*6  5*30    3*30    2*30    210
         10*15   6*15    6*10    2*105
         2*5*15  2*3*15  2*3*10  3*70
         3*5*10                  5*42
                                 7*30
                                 6*35
                                 10*21
                                 2*3*35
                                 2*5*21
                                 2*7*15
                                 3*5*14
                                 2*3*5*7

For distinct blocks instead of sums we have A116539, see A050326.
Without distinct sums we have A116540 (normal set multipartitions).
Twice-partitions of this type are counted by A279785.
Without strict blocks we have A326519.
Factorizations of this type are counted by A381633.
For constant instead of strict blocks we have A382203.
For distinct sizes instead of sums we have A382428, non-strict blocks A326517.
For equal instead of distinct block-sums we have A382429, non-strict blocks A326518.
A000670 counts patterns, ranked by A055932 and A333217, necklace A019536.
A001055 count factorizations, strict A045778.
Normal multiset partitions: A034691, A035310, A255906.
Set multipartitions: A089259, A270995, A296119, A318360.
Cf. A000110, A007716, A050320, A255903, A275780, A317532, A317583, A321469, A381635, A382204, A382214, A382216.

allnorm[n_Integer]:=Function[s,Array[Count[s,y_/;y<=#]+1&,n]]/@Subsets[Range[n-1]+1];
sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
mps[mset_]:=Union[Sort[Sort/@(#/.x_Integer:>mset[[x]])]&/@sps[Range[Length[mset]]]];
Table[Length[Join@@(Select[mps[#],UnsameQ@@Total/@#&&And@@UnsameQ@@@#&]&/@allnorm[n])],{n,0,5}]

a(10)-a(11) from Robert Price, Mar 31 2025

a(12)-a(20) from Christian Sievers, Apr 05 2025

A381806 Numbers that cannot be written as a product of squarefree numbers with distinct sums of prime indices.

Original entry on oeis.org

4, 8, 9, 16, 24, 25, 27, 32, 40, 48, 49, 54, 56, 64, 72, 80, 81, 88, 96, 104, 108, 112, 121, 125, 128, 135, 136, 144, 152, 160, 162, 169, 176, 184, 189, 192, 200, 208, 216, 224, 232, 240, 243, 248, 250, 256, 272, 288, 289, 296, 297, 304, 320, 324, 328, 336

Offset: 1

Gus Wiseman, Mar 12 2025

nonn

First differs from A212164 in having 3600.
First differs from A293243 in having 18000.
A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.
Also numbers whose prime indices cannot be partitioned into a multiset of sets with distinct sums.

			There are 4 factorizations of 18000 into squarefree numbers:
  (2*2*3*5*10*30)
  (2*2*5*6*10*15)
  (2*2*10*15*30)
  (2*5*6*10*30)
but none of these has all distinct sums of prime indices, so 18000 is in the sequence.

Strongly normal multisets of this type are counted by A292444.
These are the zeros in A381633, see A050320, A321469, A381078, A381634.
For distinct blocks see A050326, A293243, A293511, A358914, A381441.
For more on set multipartitions see A089259, A116540, A270995, A296119, A318360.
For more on set multipartitions with distinct sums see A279785, A381718.
For constant instead of strict blocks we have A381636, see A381635, A381716.
Partitions of this type are counted by A381990, complement A381992.
The complement is A382075.
A001055 counts multiset partitions, strict A045778.
A003963 gives product of prime indices.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798.
A317141 counts coarsenings of prime indices, refinements A300383.
Cf. A000688, A000720, A001222, A005117, A299202, A300385, A381454, A381870.

hwt[n_]:=Total[Cases[FactorInteger[n],{p_,k_}:>PrimePi[p]*k]];
sqfics[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[sqfics[n/d],Min@@#>=d&]],{d,Select[Rest[Divisors[n]],SquareFreeQ]}]]
Select[Range[nn],Length[Select[sqfics[#],UnsameQ@@hwt/@#&]]==0&]

A381992 Number of integer partitions of n that can be partitioned into sets with distinct sums.

Original entry on oeis.org

1, 1, 1, 2, 3, 5, 6, 9, 13, 17, 25, 33, 44, 59, 77, 100, 134, 170, 217, 282, 360, 449, 571, 719, 899, 1122, 1391, 1727, 2136, 2616, 3209, 3947, 4800, 5845, 7094, 8602, 10408, 12533, 15062, 18107, 21686, 25956, 30967, 36936, 43897, 52132, 61850, 73157, 86466, 101992, 120195

Offset: 0

Gus Wiseman, Mar 16 2025

nonn

Also the number of integer partitions of n whose Heinz number belongs to A382075 (can be written as a product of squarefree numbers with distinct sums of prime indices).

			There are 6 ways to partition (3,2,2,1) into sets:
  {{2},{1,2,3}}
  {{1,2},{2,3}}
  {{1},{2},{2,3}}
  {{2},{2},{1,3}}
  {{2},{3},{1,2}}
  {{1},{2},{2},{3}}
Of these, 3 have distinct block sums:
  {{2},{1,2,3}}
  {{1,2},{2,3}}
  {{1},{2},{2,3}}
so (3,2,2,1) is counted under a(8).
The a(1) = 1 through a(8) = 13 partitions:
  (1)  (2)  (3)    (4)      (5)      (6)        (7)        (8)
            (2,1)  (3,1)    (3,2)    (4,2)      (4,3)      (5,3)
                   (2,1,1)  (4,1)    (5,1)      (5,2)      (6,2)
                            (2,2,1)  (3,2,1)    (6,1)      (7,1)
                            (3,1,1)  (4,1,1)    (3,2,2)    (3,3,2)
                                     (2,2,1,1)  (3,3,1)    (4,2,2)
                                                (4,2,1)    (4,3,1)
                                                (5,1,1)    (5,2,1)
                                                (3,2,1,1)  (6,1,1)
                                                           (3,2,2,1)
                                                           (3,3,1,1)
                                                           (4,2,1,1)
                                                           (3,2,1,1,1)

More on set multipartitions: A089259, A116540, A270995, A296119, A318360.
Twice-partitions of this type are counted by A279785.
Multiset partitions of this type are counted by A381633, zeros of A381634.
For constant instead of strict blocks see A381717, A381636, A381635, A381716, A381991.
Normal multiset partitions of this type are counted by A381718, see A116539.
The complement is counted by A381990, ranked by A381806.
These partitions are ranked by A382075.
For distinct blocks instead of sums we have A382077, complement A382078.
For a unique choice we have A382079.
A000041 counts integer partitions, strict A000009.
A050320 counts multiset partitions of prime indices into sets.
A050326 counts multiset partitions of prime indices into distinct sets.
A265947 counts refinement-ordered pairs of integer partitions.
A382201 lists MM-numbers of sets with distinct sums.
Cf. A002846, A047966, A213427, A279786, A299202, A317142, A381870.
Cf. A293243, A293511, A358914, A381078, A381441, A381454.

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]]]];
Table[Length[Select[IntegerPartitions[n],Length[Select[mps[#], And@@UnsameQ@@@#&&UnsameQ@@Total/@#&]]>0&]],{n,0,10}]

a(21)-a(50) from Bert Dobbelaere, Mar 29 2025

A381078 Number of multisets that can be obtained by partitioning the prime indices of n into a multiset of sets (set multipartition) and taking their sums.

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, 3, 1, 2, 2, 2, 1, 5, 1, 2, 2, 2, 1, 2, 1, 2, 2, 2, 1, 2, 2, 2, 2, 2, 1, 6, 1, 2, 2, 1, 2, 5, 1, 2, 2, 5, 1, 3, 1, 2, 2, 2, 2, 5, 1, 2, 1, 2, 1, 6, 2, 2, 2

Offset: 1

Gus Wiseman, Mar 05 2025

nonn

First differs from A050320 at a(210) = 13, A050320(210) = 15. This comes from the set multipartitions {{3},{1,2,4}} and {{1,2},{3,4}}, and from {{4},{1,2,3}} and {{1,3},{2,4}}.
Also the number of multisets that can be obtained by taking the sums of prime indices of each factor in a factorization of n into squarefree numbers > 1.
A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.
A multiset partition can be regarded as an arrow in the ranked poset of integer partitions. For example, we have {{1},{1,2},{1,3},{1,2,3}}: {1,1,1,1,2,2,3,3} -> {1,3,4,6}, or (33221111) -> (6431) (depending on notation).
Set multipartitions are generally not transitive. For example, we have arrows: {{1},{1,2}}: {1,1,2} -> {1,3} and {{1,3}}: {1,3} -> {4}, but there is no set multipartition {1,1,2} -> {4}.

			The prime indices of 60 are {1,1,2,3}, with set multipartitions:
  {{1},{1,2,3}}
  {{1,2},{1,3}}
  {{1},{1},{2,3}}
  {{1},{2},{1,3}}
  {{1},{3},{1,2}}
  {{1},{1},{2},{3}}
with block-sums: {1,6}, {3,4}, {1,1,5}, {1,2,4}, {1,3,3}, {1,1,2,3}, which are all different multisets, so a(60) = 6.

Robert Price, Table of n, a(n) for n = 1..1000

Before taking sums we had A050320, strict A050326 (zeros A293243), distinct sums A381633.
For distinct blocks we have A381441.
The lower version is A381454.
For distinct block-sums we have A381634.
Other multiset partitions of prime indices:
- For multisets of constant multisets (A000688) see A381455 (upper), A381453 (lower).
- For multiset partitions (A001055) see A317141 (upper), A300383 (lower).
- For strict multiset partitions (A045778) see A381452.
- For sets of constant multisets (A050361) see A381717.
- For strict multiset partitions with distinct sums (A321469) see A381637.
- For sets of constant multisets with distinct sums (A381635) see A381716, A381636.
More on set multipartitions: A089259, A116540, A270995, A296119, A318360.
A000041 counts integer partitions, strict A000009.
A000040 lists the primes.
A003963 gives product of prime indices.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798.
A122111 represents conjugation in terms of Heinz numbers.
A265947 counts refinement-ordered pairs of integer partitions.
Cf. A000720, A001222, A002846, A005117, A025487, A066328, A213242, A213385, A213427, A299201, A299202, A300385.

hwt[n_]:=Total[Cases[FactorInteger[n],{p_,k_}:>PrimePi[p]*k]];
sqfacs[n_]:=If[n<=1,{{}},Join@@Table[(Prepend[#,d]&)/@Select[sqfacs[n/d],Min@@#>=d&],{d,Select[Rest[Divisors[n]],SquareFreeQ]}]];
Table[Length[Union[Sort[hwt/@#]&/@sqfacs[n]]],{n,100}]

a(A002110(n)) = A066723(n).

A381453 Number of multisets that can be obtained by choosing a constant integer partition of each prime index of n and taking the multiset union.

Original entry on oeis.org

1, 1, 2, 1, 2, 2, 3, 1, 3, 2, 2, 2, 4, 3, 4, 1, 2, 3, 4, 2, 6, 2, 3, 2, 3, 4, 4, 3, 4, 4, 2, 1, 4, 2, 6, 3, 6, 4, 8, 2, 2, 6, 4, 2, 6, 3, 4, 2, 6, 3, 4, 4, 5, 4, 4, 3, 8, 4, 2, 4, 6, 2, 8, 1, 8, 4, 2, 2, 6, 6, 6, 3, 4, 6, 6, 4, 6, 8, 4, 2, 5, 2, 2, 6, 4, 4, 8

Offset: 1

Gus Wiseman, Mar 08 2025

nonn

First differs from A355733 and A355735 at a(21) = 6, A355733(21) = A355735(21) = 5.
A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.
A multiset partition can be regarded as an arrow in the ranked poset of integer partitions. For example, we have {{1},{1,2},{1,3},{1,2,3}}: {1,1,1,1,2,2,3,3} -> {1,3,4,6}, or (33221111) -> (6431) (depending on notation).
Multisets of constant multisets are generally not transitive. For example, we have arrows: {{1,1},{2}}: {1,1,2} -> {2,2} and {{2,2}}: {2,2} -> {4}, but there is no multiset of constant multisets {1,1,2} -> {4}.

			The a(21) = 6 multisets are: {2,4}, {1,1,4}, {2,2,2}, {1,1,2,2}, {2,1,1,1,1}, {1,1,1,1,1,1}.
The a(n) partitions for n = 1, 3, 7, 13, 53, 21 (G = 16):
  ()  (2)   (4)     (6)       (G)                 (42)
      (11)  (22)    (33)      (88)                (411)
            (1111)  (222)     (4444)              (222)
                    (111111)  (22222222)          (2211)
                              (1111111111111111)  (21111)
                                                  (111111)

Robert Price, Table of n, a(n) for n = 1..300

Positions of 1 are A000079.
The strict case is A008966.
Before sorting we had A355731.
Choosing divisors instead of constant multisets gives A355733.
The upper version is A381455, before taking sums A000688.
Multiset partitions of prime indices:
- For multiset partitions (A001055) see A317141 (upper), A300383 (lower).
- For strict multiset partitions (A045778) see A381452.
- For set multipartitions (A050320) see A381078 (upper), A381454 (lower).
- For set systems (A050326, zeros A293243) see A381441 (upper).
- For sets of constant multisets (A050361) see A381715.
- For strict multiset partitions with distinct sums (A321469) see A381637.
- For set systems with distinct sums (A381633, zeros A381806) see A381634.
- For sets of constant multisets with distinct sums (A381635, zeros A381636) see A381716.
More on multiset partitions into constant blocks: A006171, A279784, A295935.
A000041 counts integer partitions, strict A000009.
A000040 lists the primes.
A003963 gives product of prime indices.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798.
A122111 represents conjugation in terms of Heinz numbers.
A265947 counts refinement-ordered pairs of integer partitions.
Cf. A000720, A000961, A001222, A002577, A018818, A213242, A213385, A213427, A275870, A299200, A300273, A300385.

prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Table[Length[Union[Sort/@Join@@@Tuples[Select[IntegerPartitions[#],SameQ@@#&]&/@prix[n]]]],{n,nn}]

a(A002110(n)) = A381807(n).

A381455 Number of multisets that can be obtained by taking the sum of each block of a multiset partition of the prime indices of n into a multiset of constant multisets.

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, 3, 2, 1, 3, 2, 1, 1, 1, 7, 1, 1, 1, 4, 1, 1, 1, 3, 1, 1, 1, 2, 2, 1, 1, 5, 2, 2, 1, 2, 1, 3, 1, 3, 1, 1, 1, 2, 1, 1, 2, 11, 1, 1, 1, 2, 1, 1, 1, 6, 1, 1, 2, 2, 1, 1, 1, 5, 5, 1, 1, 2, 1, 1, 1

Offset: 1

Gus Wiseman, Mar 06 2025

nonn

First differs from A000688 at a(144) = 9, A000688(144) = 10.
First differs from A295879 at a(128) = 15, A295879(128) = 13.
Also the number of multisets that can be obtained by taking the sums of prime indices of each factor in a factorization of n into prime powers > 1.
A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.
A multiset partition can be regarded as an arrow in the ranked poset of integer partitions. For example, we have {{1},{1,2},{1,3},{1,2,3}}: {1,1,1,1,2,2,3,3} -> {1,3,4,6}, or (33221111) -> (6431) (depending on notation).
Multisets of constant multisets are generally not transitive. For example, we have arrows: {{1,1},{2}}: {1,1,2} -> {2,2} and {{2,2}}: {2,2} -> {4}, but there is no multiset of constant multisets {1,1,2} -> {4}.

			The prime indices of 36 are {1,1,2,2}, with the following 4 partitions into a multiset of constant multisets:
  {{1,1},{2,2}}
  {{1},{1},{2,2}}
  {{2},{2},{1,1}}
  {{1},{1},{2},{2}}
with block-sums: {2,4}, {1,1,4}, {2,2,2}, {1,1,2,2}, which are all different, so a(36) = 4.
The prime indices of 144 are {1,1,1,1,2,2}, with the following 10 partitions into a multiset of constant multisets:
  {{2,2},{1,1,1,1}}
  {{1},{2,2},{1,1,1}}
  {{2},{2},{1,1,1,1}}
  {{1,1},{1,1},{2,2}}
  {{1},{1},{1,1},{2,2}}
  {{1},{2},{2},{1,1,1}}
  {{2},{2},{1,1},{1,1}}
  {{1},{1},{1},{1},{2,2}}
  {{1},{1},{2},{2},{1,1}}
  {{1},{1},{1},{1},{2},{2}}
with block-sums: {4,4}, {1,3,4}, {2,2,4}, {2,2,4}, {1,1,2,4}, {1,2,2,3}, {2,2,2,2}, {1,1,1,1,4}, {1,1,2,2,2}, {1,1,1,1,2,2}, of which 9 are distinct, so a(144) = 9.
The a(n) partitions for n = 4, 8, 16, 32, 36, 64, 72, 128:
  (2)   (3)    (4)     (5)      (42)    (6)       (43)     (7)
  (11)  (21)   (22)    (32)     (222)   (33)      (322)    (43)
        (111)  (31)    (41)     (411)   (42)      (421)    (52)
               (211)   (221)    (2211)  (51)      (2221)   (61)
               (1111)  (311)            (222)     (4111)   (322)
                       (2111)           (321)     (22111)  (331)
                       (11111)          (411)              (421)
                                        (2211)             (511)
                                        (3111)             (2221)
                                        (21111)            (3211)
                                        (111111)           (4111)
                                                           (22111)
                                                           (31111)
                                                           (211111)
                                                           (1111111)

Robert Price, Table of n, a(n) for n = 1..1000

Before taking sums we had A000688.
Positions of 1 are A005117.
There is a chain from the prime indices of n to a singleton iff n belongs to A300273.
The lower version is A381453.
For distinct blocks we have A381715, before sum A050361.
For distinct block-sums we have A381716, before sums A381635 (zeros A381636).
Other multiset partitions of prime indices:
- For multiset partitions (A001055) see A317141 (upper), A300383 (lower).
- For strict multiset partitions (A045778) see A381452.
- For set multipartitions (A050320) see A381078 (upper), A381454 (lower).
- For set systems (A050326) see A381441 (upper).
- For strict multiset partitions with distinct sums (A321469) see A381637.
- For set systems with distinct sums (A381633) see A381634, A293243.
More on multiset partitions into constant blocks: A006171, A279784, A295935.
A000041 counts integer partitions, strict A000009.
A000040 lists the primes.
A003963 gives product of prime indices.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798.
A122111 represents conjugation in terms of Heinz numbers.
A265947 counts refinement-ordered pairs of integer partitions.
Cf. A000720, A001222, A002577, A002846, A018818, A213242, A213427, A275870, A289078, A299200, A299202, A300385.

hwt[n_]:=Total[Cases[FactorInteger[n],{p_,k_}:>PrimePi[p]*k]];
sqfacs[n_]:=If[n<=1,{{}},Join@@Table[(Prepend[#,d]&)/@Select[sqfacs[n/d],Min@@#>=d&],{d,Select[Rest[Divisors[n]],PrimePowerQ]}]];
Table[Length[Union[Sort[hwt/@#]&/@sqfacs[n]]],{n,100}]

a(s) = 1 for any squarefree number s.

a(p^k) = A000041(k) for any prime p.

A339618 Heinz numbers of non-graphical integer partitions of even numbers.

Original entry on oeis.org

3, 7, 9, 10, 13, 19, 21, 22, 25, 28, 29, 30, 34, 37, 39, 43, 46, 49, 52, 53, 55, 57, 61, 62, 63, 66, 70, 71, 75, 76, 79, 82, 84, 85, 87, 88, 89, 91, 94, 100, 101, 102, 107, 111, 113, 115, 116, 117, 118, 121, 129, 130, 131, 133, 134, 136, 138, 139, 146, 147

Offset: 1

Gus Wiseman, Dec 18 2020

nonn

An integer partition is graphical if it comprises the multiset of vertex-degrees of some graph. Graphical partitions are counted by A000569.
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), giving a bijective correspondence between positive integers and integer partitions.
The following are equivalent characteristics for any positive integer n:
(1) the multiset of prime indices of n can be partitioned into distinct strict pairs (a set of edges);
(2) n can be factored into distinct squarefree semiprimes;
(3) the unordered prime signature of n is graphical.

			The sequence of terms together with their prime indices begins:
      3: {2}         43: {14}        79: {22}
      7: {4}         46: {1,9}       82: {1,13}
      9: {2,2}       49: {4,4}       84: {1,1,2,4}
     10: {1,3}       52: {1,1,6}     85: {3,7}
     13: {6}         53: {16}        87: {2,10}
     19: {8}         55: {3,5}       88: {1,1,1,5}
     21: {2,4}       57: {2,8}       89: {24}
     22: {1,5}       61: {18}        91: {4,6}
     25: {3,3}       62: {1,11}      94: {1,15}
     28: {1,1,4}     63: {2,2,4}    100: {1,1,3,3}
     29: {10}        66: {1,2,5}    101: {26}
     30: {1,2,3}     70: {1,3,4}    102: {1,2,7}
     34: {1,7}       71: {20}       107: {28}
     37: {12}        75: {2,3,3}    111: {2,12}
     39: {2,6}       76: {1,1,8}    113: {30}
For example, there are three possible multigraphs with degrees (1,1,3,3):
  {{1,2},{1,2},{1,2},{3,4}}
  {{1,2},{1,2},{1,3},{2,4}}
  {{1,2},{1,2},{1,4},{2,3}}.
Since none of these is a graph, the Heinz number 100 belongs to the sequence.

Eric Weisstein's World of Mathematics, Graphical partition.

A181819 applied to A320894 gives this sequence.
A300061 is a superset.
A339617 counts these partitions.
A320922 ranks the complement, counted by A000569.
A006881 lists squarefree semiprimes.
A320656 counts factorizations into squarefree semiprimes.
A339659 counts graphical partitions of 2n into k parts.
The following count vertex-degree partitions and give their Heinz numbers:
- A058696 counts partitions of 2n (A300061).
- A000070 counts non-multigraphical partitions of 2n (A339620).
- A209816 counts multigraphical partitions (A320924).
- A339655 counts non-loop-graphical partitions of 2n (A339657).
- A339656 counts loop-graphical partitions (A339658).
- A339617 counts non-graphical partitions of 2n (A339618 [this sequence]).
- A000569 counts graphical partitions (A320922).
The following count partitions of even length and give their Heinz numbers:
- A027187 has no additional conditions (A028260).
- A096373 cannot be partitioned into strict pairs (A320891).
- A338914 can be partitioned into strict pairs (A320911).
- A338915 cannot be partitioned into distinct pairs (A320892).
- A338916 can be partitioned into distinct pairs (A320912).
- A339559 cannot be partitioned into distinct strict pairs (A320894).
- A339560 can be partitioned into distinct strict pairs (A339561).
Cf. A001358, A007717, A050326, A320923, A338899.

strs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[strs[n/d],Min@@#>d&]],{d,Select[Divisors[n],And[SquareFreeQ[#],PrimeOmega[#]==2]&]}]];
nrmptn[n_]:=Join@@MapIndexed[Table[#2[[1]],{#1}]&,If[n==1,{},Flatten[Cases[FactorInteger[n]//Reverse,{p_,k_}:>Table[PrimePi[p],{k}]]]]];
Select[Range[100],EvenQ[Length[nrmptn[#]]]&&strs[Times@@Prime/@nrmptn[#]]=={}&]

Equals A300061 \ A320922.

For all n, A181821(a(n)) and A304660(a(n)) belong to A320894.

A179643 Products of exactly 2 distinct squares of primes and a different prime (p^2 * q^2 * r).

Original entry on oeis.org

180, 252, 300, 396, 450, 468, 588, 612, 684, 700, 828, 882, 980, 1044, 1100, 1116, 1300, 1332, 1452, 1476, 1548, 1575, 1692, 1700, 1900, 1908, 2028, 2124, 2156, 2178, 2196, 2205, 2300, 2412, 2420, 2450, 2475, 2548, 2556, 2628, 2844, 2900, 2925, 2988

Offset: 1

Vladimir Joseph Stephan Orlovsky, Jul 21 2010

nonn

A050326(a(n)) = 5, subsequence of A225228. - Reinhard Zumkeller, May 03 2013

			180 = 2^2 * 3^2 * 5, 252 = 2^2 * 3^2 * 7, 300 = 2^2 * 3 * 5^2, ...

T. D. Noe, Table of n, a(n) for n = 1..1000
Will Nicholes, List of Prime Signatures
Index to sequences related to prime signature

Cf. A006881, A007304, A065036, A085986, A085987, A178739, A179642.

Cf. A179695.

f[n_]:=Sort[Last/@FactorInteger[n]]=={1,2,2}; Select[Range[3000], f]

list(lim)=my(v=List(),t);forprime(p=2,sqrt(lim\12),forprime(q=p+1,sqrt(lim\p^2\2),t=(p*q)^2;forprime(r=2,lim\t,if(p==r||q==r,next);listput(v,t*r))));vecsort(Vec(v)) \\ Charles R Greathouse IV, Jul 19 2011

from math import isqrt
from sympy import primepi, primerange, integer_nthroot
def A179643(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        kmin = kmax >> 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x): return n+x+sum((t:=primepi(s:=isqrt(y:=isqrt(x//r))))+(t*(t-1)>>1)-sum(primepi(y//k) for k in primerange(1, s+1)) for r in primerange(x+1))+sum(primepi(isqrt(x//p**3)) for p in primerange(integer_nthroot(x,3)[0]+1))-primepi(integer_nthroot(x,5)[0])
    return bisection(f,n,n) # Chai Wah Wu, Mar 27 2025

cp's OEIS Frontend

A085987 Product of exactly four primes, three of which are distinct (p^2*q*r).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Python

Extensions

A381441 Number of multisets that can be obtained by partitioning the prime indices of n into a set of sets (set system) and taking their sums.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A381718 Number of normal multiset partitions of weight n into sets with distinct sums.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Extensions

A381806 Numbers that cannot be written as a product of squarefree numbers with distinct sums of prime indices.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A381992 Number of integer partitions of n that can be partitioned into sets with distinct sums.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Extensions

A381078 Number of multisets that can be obtained by partitioning the prime indices of n into a multiset of sets (set multipartition) and taking their sums.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A381453 Number of multisets that can be obtained by choosing a constant integer partition of each prime index of n and taking the multiset union.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A381455 Number of multisets that can be obtained by taking the sum of each block of a multiset partition of the prime indices of n into a multiset of constant multisets.

Views

Author

Keywords

A085987 Product of exactly four primes, three of which are distinct (p^2qr).