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

A321455 Number of ways to factor n into factors > 1 all having the same sum of prime indices.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Nov 10 2018

Keywords

Comments

Also the number of multiset partitions of the multiset of prime indices of n with equal block-sums.
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. The sum of prime indices of n is A056239(n).

Examples

			The a(1440) = 6 factorizations into factors all having the same sum of prime indices:
  (10*12*12)
  (5*6*6*8)
  (9*10*16)
  (30*48)
  (36*40)
  (1440)
The a(900) = 5 multiset partitions with equal block-sums:
  {{1,1,2,2,3,3}}
  {{3,3},{1,1,2,2}}
  {{1,2,3},{1,2,3}}
  {{1,3},{1,3},{2,2}}
  {{3},{3},{1,2},{1,2}}

Links

Crossrefs

Positions of 1's are A321453. Positions of terms > 1 are A321454.
Cf. A001055, A035470, A056239, A279787, A305551, A321469, A322794, A326515, A326516, A326518, A326534.

Programs

  • Mathematica
    hwt[n_]:=Total[Cases[FactorInteger[n],{p_,k_}:>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[Select[facs[n],SameQ@@hwt/@#&]],{n,100}]
  • PARI
    A056239(n) = if(1==n, 0, my(f=factor(n)); sum(i=1, #f~, f[i, 2] * primepi(f[i, 1])));
all_have_same_sum_of_pis(facs) = if(!#facs, 1, (#Set(apply(A056239,facs)) == 1));
A321455(n, m=n, facs=List([])) = if(1==n, all_have_same_sum_of_pis(facs), my(s=0, newfacs); fordiv(n, d, if((d>1)&&(d<=m), newfacs = List(facs); listput(newfacs,d); s += A321455(n/d, d, newfacs))); (s)); \\ Antti Karttunen, Jan 20 2025

Extensions

Data section extended to a(108) by Antti Karttunen, Jan 20 2025

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

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 2, 1, 1, 2, 3, 1, 4, 2, 2, 1, 5, 1, 6, 2, 2, 3, 8, 1, 3, 4, 1, 2, 10, 2, 12, 1, 3, 5, 4, 1, 15, 6, 4, 2, 18, 2, 22, 3, 2, 8, 27, 1, 3, 3, 5, 4, 32, 1, 6, 2, 6, 10, 38, 2, 46, 12, 2, 1, 8, 3, 54, 5, 8, 4, 64, 1, 76, 15, 3, 6, 6, 4, 89, 2, 1
Offset: 1

Views

Author

Gus Wiseman, Mar 08 2025

Keywords

Comments

First differs from A357982 at a(25) = 3, A357982(25) = 4.
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}.

Examples

			The a(25) = 3 multisets are: {3,3}, {1,2,3}, {1,1,2,2}.

Crossrefs

For constant instead of strict partitions see A381453, A355733, A381455, A000688.
Positions of 1 are A003586.
The upper version is A381078, before sums A050320.
For distinct block-sums see A381634, A381633, A381806.
Multiset partitions of prime indices:
- For multiset partitions (A001055) see A317141 (upper), A300383 (lower).
- For strict multiset partitions (A045778) see A381452.
- 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 sets of constant multisets with distinct sums (A381635, zeros A381636) see A381716.
More on set systems: A050342, A116539, A296120, A318361.
More on set multipartitions: A089259, A116540, A270995, A296119, A318360.
More on set multipartitions with distinct sums: A279785, A381717, A381718.
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.
A358914 counts twice-partitions into distinct strict partitions.
Cf. A000720, A001222, A002846, A005117, A066328, A213242, A213385, A213427, A293511, A299200, A299201, A299202, A300385, A317142.

Programs

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

Formula

a(A002110(n)) = A381808(n).

A381633 Number of ways to partition the prime indices of n into sets with distinct 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, 4, 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, 3, 1, 2, 1, 0, 2, 5, 1, 1, 2, 4, 1, 0, 1, 2, 1, 1, 2, 5, 1, 0, 0, 2, 1, 4, 2, 2, 2
Offset: 1

Views

Author

Gus Wiseman, Mar 09 2025

Keywords

Comments

First differs from A050326 at 30, 60, 70, 90, ...
First differs from A339742 at 42, 66, 78, 84, ...
First differs from A381634 at a(210) = 12, A381634(210) = 10.
Also the number of factorizations on n into squarefree numbers > 1 with distinct sums of prime indices.
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, sum A056239.

Examples

			The A050320(60) = 6 ways to partition {1,1,2,3} into sets are:
  {{1},{1,2,3}}
  {{1,2},{1,3}}
  {{1},{1},{2,3}}
  {{1},{2},{1,3}}
  {{1},{3},{1,2}}
  {{1},{1},{2},{3}}
Of these, only the following have distinct block-sums:
  {{1},{1,2,3}}
  {{1,2},{1,3}}
  {{1},{2},{1,3}}
So a(60) = 3.

Crossrefs

Without distinct block-sums we have A050320, after sums A381078 (lower A381454).
For distinct blocks instead of sums we have A050326, after sums A381441, see A358914.
Taking block-sums (and sorting) gives A381634.
For constant instead of strict blocks we have A381635, see A381716, A381636.
Positions of 0 are A381806, superset of A293243.
Positions of 1 are A381870, superset of A293511.
More on set multipartitions with distinct sums: A279785, A381717, A381718.
More on set multipartitions: A089259, A116540, A270995, A296119, A318360.
A000041 counts integer partitions, strict A000009.
A001055 count multiset partitions of prime indices, see A317141 (upper), A300383 (lower).
A003963 gives product of prime indices.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798.
A265947 counts refinement-ordered pairs of integer partitions.
Cf. A000720, A001222, A002846, A005117, A213242, A299202, A300385, A317142, A317143, A321469.

Programs

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

A275780 Number of set partitions of [n] into blocks with distinct element sums.

Original entry on oeis.org

1, 1, 2, 4, 12, 43, 160, 668, 3098, 15465, 83100, 477651, 2914505, 18795814, 127790544, 911448954, 6808162094, 53067398065, 430956571977, 3636314065247, 31841519540324, 288664242344692, 2706949104147162, 26205222185730884, 261681461422075548, 2691088457402830312
Offset: 0

Views

Author

Alois P. Heinz, Aug 08 2016

Keywords

Examples

			a(3) = 4: 123, 13|2, 1|23, 1|2|3.
a(4) = 12: 1234, 123|4, 124|3, 12|34, 134|2, 13|24, 1|234, 1|23|4, 14|2|3, 1|24|3, 1|2|34, 1|2|3|4.

Links

Crossrefs

Cf. A000110, A007837, A035470, A038041, A275781, A278339, A300335, A321469, A326512, A326513, A326519, A326535.

Programs

  • Mathematica
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
Table[Length[Select[sps[Range[n]],UnsameQ@@Total/@#&]],{n,0,10}] (* Gus Wiseman, Jul 13 2019 *)

Formula

a(n) = A000110(n) - A275781(n).

Extensions

a(17)-a(25) from Christian Sievers, Aug 20 2024

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

Views

Author

Gus Wiseman, Mar 06 2025

Keywords

Comments

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

Examples

			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.

Crossrefs

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.

Programs

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

Formula

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

Views

Author

Gus Wiseman, Mar 26 2025

Keywords

Comments

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.

Examples

			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

Crossrefs

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.

Programs

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

Extensions

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

Views

Author

Gus Wiseman, Mar 12 2025

Keywords

Comments

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.

Examples

			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.

Crossrefs

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.

Programs

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

A326519 Number of normal multiset partitions of weight n where each part has a different sum.

Original entry on oeis.org

1, 1, 3, 11, 51, 259, 1461, 9133, 62348, 459547, 3632419
Offset: 0

Views

Author

Gus Wiseman, Jul 12 2019

Keywords

Comments

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

Examples

			The a(0) = 1 through a(3) = 11 normal multiset partitions:
  {}  {{1}}  {{1,1}}    {{1,1,1}}
             {{1,2}}    {{1,1,2}}
             {{1},{2}}  {{1,2,2}}
                        {{1,2,3}}
                        {{1},{1,1}}
                        {{1},{1,2}}
                        {{1},{2,2}}
                        {{1},{2,3}}
                        {{2},{1,2}}
                        {{2},{1,3}}
                        {{1},{2},{3}}

Links

Crossrefs

Cf. A038041, A255906, A275780, A317583, A321469, A326517, A326518, A326520, A326521, A326535.

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

Extensions

a(8)-a(10) from Robert Price, Apr 03 2025

A326535 MM-numbers of multiset partitions where each part has a different sum.

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 10, 11, 13, 14, 15, 17, 19, 21, 22, 23, 26, 29, 30, 31, 33, 34, 37, 38, 39, 41, 42, 43, 46, 47, 51, 53, 55, 57, 58, 59, 61, 62, 65, 66, 67, 69, 71, 73, 74, 77, 78, 79, 82, 83, 85, 86, 87, 89, 91, 93, 94, 95, 97, 101, 102, 103, 106, 107, 109
Offset: 1

Views

Author

Gus Wiseman, Jul 12 2019

Keywords

Comments

First differs from A298540 in having 187.
These are numbers where each prime index has a different sum of prime indices. 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. The multiset multisystem with MM-number n is obtained by taking the multiset of prime indices of each prime index of n. For example, the prime indices of 78 are {1,2,6}, so the multiset multisystem with MM-number 78 is {{},{1},{1,2}}.

Examples

			The sequence of multiset partitions where each part has a different sum, preceded by their MM-numbers, begins:
   1: {}
   2: {{}}
   3: {{1}}
   5: {{2}}
   6: {{},{1}}
   7: {{1,1}}
  10: {{},{2}}
  11: {{3}}
  13: {{1,2}}
  14: {{},{1,1}}
  15: {{1},{2}}
  17: {{4}}
  19: {{1,1,1}}
  21: {{1},{1,1}}
  22: {{},{3}}
  23: {{2,2}}
  26: {{},{1,2}}
  29: {{1,3}}
  30: {{},{1},{2}}
  31: {{5}}

Links

Crossrefs

Cf. A038041, A112798, A275780, A302242, A320324, A321469, A326519, A326533, A326534, A326536, A326537.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Select[Range[100],UnsameQ@@Total/@primeMS/@primeMS[#]&]

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

Views

Author

Gus Wiseman, Mar 05 2025

Keywords

Comments

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

Examples

			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.

Links

Crossrefs

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.

Programs

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

Formula

a(A002110(n)) = A066723(n).
Showing 1-10 of 42 results. Next

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.