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.

Previous Showing 41-50 of 182 results. Next

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.

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

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.

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

A048109 Numbers having equally many squarefree and nonsquarefree divisors; number of unitary divisors of n (A034444) = number of non-unitary divisors of n (A048105).

Original entry on oeis.org

8, 24, 27, 40, 54, 56, 88, 104, 120, 125, 135, 136, 152, 168, 184, 189, 232, 248, 250, 264, 270, 280, 296, 297, 312, 328, 343, 344, 351, 375, 376, 378, 408, 424, 440, 456, 459, 472, 488, 513, 520, 536, 552, 568, 584, 594, 616, 621, 632, 664, 680, 686, 696
Offset: 1

Views

Author

Keywords

Comments

For these terms the number of divisors should be a special power of two because ud(n) = 2^r and nud(n) = ud(n). In particular the exponent of 2 is 1+A001221(n), the number of distinct prime factors + 1. Thus this is a subsequence of A036537 where A000005(A036537(n)) = 2^s; here s = 1+A001221(n).
Let us introduce a function D(n) = sigma_0(n)/2^(alpha(1)+...+alpha(r)), sigma_0(n) number of divisors of n (A000005), prime factorization of n = p(1)^alpha(1) * ... * p(r)^alpha(r), alpha(1)+...+alpha(r) is sequence (A001222). This function splits the set of positive integers into subsets, according to the value of D(n). Squarefree numbers (A005117) has D(n)=1, other numbers are "deviated" from the squarefree ideal and have 0 < D(n) < 1. So for D(n)=1/2 we have A048109, D(n)=3/4 we have A060687. - Ctibor O. Zizka, Sep 21 2008
Integers n such that there are exactly 3 Abelian groups of order n. That is, n such that A000688(n)=3. In other words, in the prime factorization of n there is exactly one prime with exponent of 3 and the others have exponent of 1. - Geoffrey Critzer, Jun 09 2015
The asymptotic density of this sequence is (6/Pi^2) * Sum_{k>=1} 1/(prime(k)^2*(prime(k)+1)) = (1/zeta(2)) * Sum_{k>=3} (-1)^(k+1) * P(k) = 0.0741777413672596019212880156082745910562809066233004356300970463709875..., where P is the prime zeta function. - Amiram Eldar, Jul 11 2020

Examples

			n = 88 = 2*2*2*11 has 8 divisors, of which 4 are unitary divisors (because of 2 distinct prime factors) and 4 are nonunitary divisors: U={1,88,11,8} and NU = {2,44,4,22}.
		

Crossrefs

Programs

  • Maple
    filter:= proc(n) local F;
      F:= ifactors(n)[2];
      mul(t[2]+1,t=F) = 2^(1+nops(F))
    end proc;
    select(filter, [$1..1000]); # Robert Israel, Jun 09 2015
  • Mathematica
    Position[Table[FiniteAbelianGroupCount[n], {n, 1, 1000}],3] // Flatten (* Geoffrey Critzer, Jun 09 2015 *)
  • PARI
    is(n)=select(e->e>1, factor(n)[,2])==[3]~ \\ Charles R Greathouse IV, Jun 10 2015
    
  • PARI
    isok(n) = sumdiv(n, d, issquarefree(d)) == sumdiv(n, d, !issquarefree(d)); \\ Michel Marcus, Jun 24 2015
    
  • Python
    from math import isqrt
    from sympy import mobius, primerange
    def A048109(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 g(x): return sum(mobius(k)*(x//k**2) for k in range(1, isqrt(x)+1))
        def f(x): return int(n+x+sum(sum(-g(x//p**j) if j&1 else g(x//p**j) for j in range(3,x.bit_length())) for p in primerange(isqrt(x)+1)))
        return bisection(f,n,n) # Chai Wah Wu, Feb 24 2025

Formula

Numbers k such that d(k) = 2^(omega(k)+1) or A000005(k) = 2^(A001221(k) + 1) = 2 * A034444(k).

Extensions

New name based on comment by Ivan Neretin, Jun 19 2015

A381716 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 constant blocks with distinct sums.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Mar 10 2025

Keywords

Comments

First differs from A381635 at a(1728) = 4, A381635(1728) = 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.

Examples

			The prime indices of 1728 are {1,1,1,1,1,1,2,2,2}, with multiset partitions into constant multisets with distinct sums:
  {{1,1,1,1,1,1},{2,2},{2}}
  {{1,1,1,1,1},{1},{2,2,2}}
  {{1,1,1,1,1},{1},{2,2},{2}}
  {{1,1,1,1},{1,1},{2,2,2}}
  {{1,1,1},{1,1},{1},{2,2,2}}
with block-sums: {1,5,6}, {2,4,6}, {1,2,3,6}, {1,2,4,5}, so a(1728) = 4.
		

Crossrefs

Without distinct sums we have A000688, after sums A381455 (upper), A381453 (lower).
More on multiset partitions into constant blocks: A006171, A279784, A295935.
For strict instead of constant we have A381633, before sums A381634.
Before taking sums we had A381635.
Positions of 0 are A381636.
For distinct blocks instead of sums we have A381715.
A001055 counts multiset partitions, see A317141 (upper), A300383 (lower).
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798.

Programs

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

A381717 Number of integer partitions of n that cannot be partitioned into constant multisets with distinct block-sums.

Original entry on oeis.org

0, 0, 0, 0, 1, 0, 0, 1, 3, 2, 3, 6, 7, 10, 15, 15, 28, 37, 47, 64, 71, 97, 139, 173, 215, 273, 361, 439, 551, 691, 853, 1078, 1325, 1623, 2046, 2458, 2998, 3697, 4527, 5472, 6590, 7988, 9590, 11598, 13933, 16560, 19976, 23822, 28420, 33797, 40088, 47476, 56369, 66678
Offset: 0

Views

Author

Gus Wiseman, Mar 16 2025

Keywords

Comments

Conjecture: Also the number of integer partitions of n having no permutation with all distinct run-sums, ranked by zeros of A382876. In other words, a partition has a permutation with all distinct run-sums iff it has a multiset partition into constant blocks with all distinct block-sums, where the run-sums of a sequence are obtained by splitting it into maximal runs and taking their sums.

Examples

			For y = (3,2,2,1) we have the multiset partition {{3},{2,2},{1}}, so y is not counted under a(8).
For y = (3,2,1,1,1) there are 3 multiset partitions into constant multisets:
  {{3},{2},{1,1,1}}
  {{3},{2},{1,1},{1}}
  {{3},{2},{1},{1},{1}}
but none of these has distinct block-sums, so y is counted under a(8).
For y = (3,3,1,1,1,1,1,1) we have multiset partitions:
  {{1},{3,3},{1,1,1,1,1}}
  {{1,1},{3,3},{1,1,1,1}}
  {{1},{1,1},{3,3},{1,1,1}}
so y is not counted under a(12).
The a(4) = 1 through a(13) = 10 partitions:
  211  .  .  3211  422    4221  6211   4322     633      5422
                   4211   5211  33211  7211     8211     6331
                   32111        42211  43211    43221    9211
                                       422111   44211    54211
                                       431111   53211    63211
                                       3221111  432111   333211
                                                4221111  432211
                                                         532111
                                                         4321111
                                                         42211111
		

Crossrefs

Twice-partitions of this type (constant with distinct) are counted by A279786.
Multiset partitions of this type are ranked by A326535 /\ A355743.
These partitions are ranked by A381636, zeros of A381635.
For strict instead of constant blocks we have A381990, see A381806, A381633, A382079.
For equal instead of distinct block-sums we have A381993.
A000041 counts integer partitions, strict A000009.
A000688 counts factorizations into prime powers, see A381455, A381453.
A001055 counts factorizations, strict A045778, see A317141, A300383.
A050361 counts factorizations into distinct prime powers.

Programs

  • Mathematica
    mce[y_]:=Table[ConstantArray[y[[1]],#]&/@ptn,{ptn,IntegerPartitions[Length[y]]}];
    Table[Length[Select[IntegerPartitions[n],Select[Join@@@Tuples[mce/@Split[#]],UnsameQ@@Total/@#&]=={}&]],{n,0,30}]

Extensions

a(37)-a(53) from Robert Price, Mar 31 2025

A212172 Row n of table represents second signature of n: list of exponents >= 2 in canonical prime factorization of n, in nonincreasing order, or 0 if no such exponent exists.

Original entry on oeis.org

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

Views

Author

Matthew Vandermast, Jun 03 2012

Keywords

Comments

Length of row n equals A056170(n) if A056170(n) is positive, or 1 if A056170(n) = 0.
The multiset of exponents >=2 in the prime factorization of n completely determines a(n) for over 20 sequences in the database (see crossreferences). It also determines the fractions A034444(n)/A000005(n) and A037445(n)/A000005(n).
For squarefree numbers, this multiset is { } (the empty multiset). The use of 0 in the table to represent each n with no exponents >=2 in its prime factorization accords with the usual OEIS practice of using 0 to represent nonexistent elements when possible. In comments, the second signature of squarefree numbers will be represented as { }.
For each second signature {S}, there exist values of j and k such that, if the second signature of n is {S}, then A085082(n) is congruent to j modulo k. These values are nontrivial unless {S} = { }. Analogous (but not necessarily identical) values of j and k also exist for each second signature with respect to A088873 and A181796.
Each sequence of integers with a given second signature {S} has a positive density, unlike the analogous sequences for prime signatures. The highest of these densities is 6/Pi^2 = 0.607927... for A005117 ({S} = { }).

Examples

			First rows of table read: 0; 0; 0; 2; 0; 0; 0; 3; 2; 0; 0; 2;...
12 = 2^2*3 has positive exponents 2 and 1 in its canonical prime factorization (1s are often left implicit as exponents). Since only exponents that are 2 or greater appear in a number's second signature, 12's second signature is {2}.
30 = 2*3*5 has no exponents greater than 1 in its prime factorization. The multiset of its exponents >= 2 is { } (the empty multiset), represented in the table with a 0.
72 = 2^3*3^2 has positive exponents 3 and 2 in its prime factorization, as does 108 = 2^2*3^3. Rows 72 and 108 both read {3,2}.
		

Crossrefs

A181800 gives first integer of each second signature. Also see A212171, A212173-A212181, A212642-A212644.
Functions determined by exponents >=2 in the prime factorization of n:
Additive: A046660, A056170.
Other: A007424, A051903 (for n > 1), A056626, A066301, A071325, A072411, A091050, A107078, A185102 (for n > 1), A212180.
Sequences that contain all integers of a specific second signature: A005117 (second signature { }), A060687 ({2}), A048109 ({3}).

Programs

  • Magma
    &cat[IsEmpty(e)select [0]else Reverse(Sort(e))where e is[pe[2]:pe in Factorisation(n)|pe[2]gt 1]:n in[1..102]]; // Jason Kimberley, Jun 13 2012
  • Mathematica
    row[n_] := Select[ FactorInteger[n][[All, 2]], # >= 2 &] /. {} -> 0 /. {k__} -> Sequence[k]; Table[row[n], {n, 1, 100}] (* Jean-François Alcover, Apr 16 2013 *)

Formula

For nonsquarefree n, row n is identical to row A057521(n) of table A212171.

A355743 Numbers whose prime indices are all prime-powers.

Original entry on oeis.org

1, 3, 5, 7, 9, 11, 15, 17, 19, 21, 23, 25, 27, 31, 33, 35, 41, 45, 49, 51, 53, 55, 57, 59, 63, 67, 69, 75, 77, 81, 83, 85, 93, 95, 97, 99, 103, 105, 109, 115, 119, 121, 123, 125, 127, 131, 133, 135, 147, 153, 155, 157, 159, 161, 165, 171, 175, 177, 179, 187
Offset: 1

Views

Author

Gus Wiseman, Jul 24 2022

Keywords

Comments

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 MM-numbers of multiset partitions into constant multisets, where the multiset of multisets with MM-number n is formed by taking the multiset of prime indices of each part of the multiset of prime indices of n. For example, the prime indices of 78 are {1,2,6}, so the multiset of multisets with MM-number 78 is {{},{1},{1,2}}.

Examples

			The terms together with their prime indices begin:
   1: {}
   3: {2}
   5: {3}
   7: {4}
   9: {2,2}
  11: {5}
  15: {2,3}
  17: {7}
  19: {8}
  21: {2,4}
  23: {9}
  25: {3,3}
  27: {2,2,2}
  31: {11}
  33: {2,5}
  35: {3,4}
  41: {13}
  45: {2,2,3}
		

Crossrefs

The multiplicative version is A000688, strict A050361, coprime A354911.
The case of only primes (not all prime-powers) is A076610, strict A302590.
Allowing prime index 1 gives A302492.
These are the products of elements of A302493.
Requiring n to be a prime-power gives A302601.
These are the positions of 1's in A355741.
The squarefree case is A356065.
The complement is A356066.
A001222 counts prime-power divisors.
A023894 counts ptns into prime-powers, strict A054685, with 1's A023893.
A034699 gives maximal prime-power divisor.
A246655 lists the prime-powers (A000961 includes 1), towers A164336.
A355742 chooses a prime-power divisor of each prime index.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],And@@PrimePowerQ/@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.
		

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.

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

A355742 Number of ways to choose a sequence of prime-power divisors, one of each prime index of n. Product of bigomega over the prime indices of n, with multiplicity.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Jul 20 2022

Keywords

Comments

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.

Examples

			The prime indices of 49 are {4,4}, and the a(49) = 4 choices are: (2,2), (2,4), (4,2), (4,4).
The prime indices of 777 are {2,4,12}, and the a(777) = 6 choices are: (2,2,2), (2,2,3), (2,2,4), (2,4,2), (2,4,3), (2,4,4).
		

Crossrefs

The unordered version is A001970, row-sums of A061260.
Positions of 1's are A076610, just primes A355743.
Positions of 0's are A299174.
Allowing all divisors (not just primes) gives A355731, firsts A355732.
Choosing only prime factors (not prime-powers) gives A355741.
Counting multisets of primes gives A355744.
The case of weakly increasing primes A355745, all divisors A355735.
A000688 counts factorizations into prime powers.
A001414 adds up distinct prime factors, counted by A001221.
A003963 multiplies together the prime indices of n.
A056239 adds up prime indices, row sums of A112798, counted by A001222.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Table[Times@@PrimeOmega/@primeMS[n],{n,100}]

Formula

Totally multiplicative with a(prime(k)) = A001222(k).

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

Views

Author

Gus Wiseman, Mar 08 2025

Keywords

Comments

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

Examples

			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)
		

Crossrefs

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.

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[#],SameQ@@#&]&/@prix[n]]]],{n,nn}]

Formula

a(A002110(n)) = A381807(n).
Previous Showing 41-50 of 182 results. Next