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

A000688 Number of Abelian groups of order n; number of factorizations of n into prime powers.

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, 3, 1, 2, 1, 2, 1, 1, 1, 7, 1, 2, 2, 4, 1, 1, 1, 3, 1, 1, 1
Offset: 1

Views

Author

Keywords

Comments

Equivalently, number of Abelian groups with n conjugacy classes. - Michael Somos, Aug 10 2010
a(n) depends only on prime signature of n (cf. A025487). So a(24) = a(375) since 24 = 2^3*3 and 375 = 3*5^3 both have prime signature (3, 1).
Also number of rings with n elements that are the direct product of fields; these are the commutative rings with n elements having no nilpotents; likewise the commutative rings where for every element x there is a k > 0 such that x^(k+1) = x. - Franklin T. Adams-Watters, Oct 20 2006
Range is A033637.
a(n) = 1 if and only if n is from A005117 (squarefree numbers). See the Ahmed Fares comment there, and the formula for n>=2 below. - Wolfdieter Lang, Sep 09 2012
Also, from a theorem of Molnár (see [Molnár]), the number of (non-isomorphic) abelian groups of order 2*n + 1 is equal to the number of non-congruent lattice Z-tilings of R^n by crosses, where a "cross" is a unit cube in R^n for which at each facet is attached another unit cube (Z, R are the integers and reals, respectively). (Cf. [Horak].) - L. Edson Jeffery, Nov 29 2012
Zeta(k*s) is the Dirichlet generating function of the characteristic function of numbers which are k-th powers (k=1 in A000012, k=2 in A010052, k=3 in A010057, see arXiv:1106.4038 Section 3.1). The infinite product over k (here) is the number of representations n=product_i (b_i)^(e_i) where all exponents e_i are distinct and >=1. Examples: a(n=4)=2: 4^1 = 2^2. a(n=8)=3: 8^1 = 2^1*2^2 = 2^3. a(n=9)=2: 9^1 = 3^2. a(n=12)=2: 12^1 = 3*2^2. a(n=16)=5: 16^1 = 2*2^3 = 4^2 = 2^2*4^1 = 2^4. If the e_i are the set {1,2} we get A046951, the number of representations as a product of a number and a square. - R. J. Mathar, Nov 05 2016
See A060689 for the number of non-abelian groups of order n. - M. F. Hasler, Oct 24 2017
Kendall & Rankin prove that the density of {n: a(n) = m} exists for each m. - Charles R Greathouse IV, Jul 14 2024

Examples

			a(1) = 1 since the trivial group {e} is the only group of order 1, and it is Abelian; alternatively, since the only factorization of 1 into prime powers is the empty product.
a(p) = 1 for any prime p, since the only factorization into prime powers is p = p^1, and (in view of Lagrange's theorem) there is only one group of prime order p; it is isomorphic to (Z/pZ,+) and thus Abelian.
From _Wolfdieter Lang_, Jul 22 2011: (Start)
a(8) = 3 because 8 = 2^3, hence a(8) = pa(3) = A000041(3) = 3 from the partitions (3), (2, 1) and (1, 1, 1), leading to the 3 factorizations of 8: 8, 4*2 and 2*2*2.
a(36) = 4 because 36 = 2^2*3^2, hence a(36) = pa(2)*pa(2) = 4 from the partitions (2) and (1, 1), leading to the 4 factorizations of 36: 2^2*3^2, 2^2*3^1*3^1, 2^1*2^1*3^2 and 2^1*2^1*3^1*3^1.
(End)
		

References

  • Steven R. Finch, Mathematical Constants, Cambridge, 2003, pp. 274-278.
  • D. S. Mitrinovic et al., Handbook of Number Theory, Kluwer, Section XIII.12, p. 468.
  • J. S. Rose, A Course on Group Theory, Camb. Univ. Press, 1978, see p. 7.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • A. Speiser, Die Theorie der Gruppen von endlicher Ordnung, 4. Auflage, Birkhäuser, 1956.

Crossrefs

Cf. A080729 (Dgf at s=2), A369634 (Dgf at s=3).

Programs

  • Haskell
    a000688 = product . map a000041 . a124010_row
    -- Reinhard Zumkeller, Aug 28 2014
    
  • Maple
    with(combinat): readlib(ifactors): for n from 1 to 120 do ans := 1: for i from 1 to nops(ifactors(n)[2]) do ans := ans*numbpart(ifactors(n)[2][i][2]) od: printf(`%d,`,ans): od: # James Sellers, Dec 07 2000
  • Mathematica
    f[n_] := Times @@ PartitionsP /@ Last /@ FactorInteger@n; Array[f, 107] (* Robert G. Wilson v, Sep 22 2006 *)
    Table[FiniteAbelianGroupCount[n], {n, 200}] (* Requires version 7.0 or later. - Vladimir Joseph Stephan Orlovsky, Jul 01 2011 *)
  • PARI
    A000688(n)=local(f);f=factor(n);prod(i=1,matsize(f)[1],numbpart(f[i,2])) \\ Michael B. Porter, Feb 08 2010
    
  • PARI
    a(n)=my(f=factor(n)[,2]); prod(i=1,#f,numbpart(f[i])) \\ Charles R Greathouse IV, Apr 16 2015
    
  • Python
    from sympy import factorint, npartitions
    from math import prod
    def A000688(n): return prod(map(npartitions,factorint(n).values())) # Chai Wah Wu, Jan 14 2022
  • Sage
    def a(n):
        F=factor(n)
        return prod([number_of_partitions(F[i][1]) for i in range(len(F))])
    # Ralf Stephan, Jun 21 2014
    

Formula

Multiplicative with a(p^k) = number of partitions of k = A000041(k); a(mn) = a(m)a(n) if (m, n) = 1.
a(2n) = A101872(n).
a(n) = Product_{j = 1..N(n)} A000041(e(j)), n >= 2, if
n = Product_{j = 1..N(n)} prime(j)^e(j), N(n) = A001221(n). See the Richert reference, quoting A. Speiser's book on finite groups (in German, p. 51 in words). - Wolfdieter Lang, Jul 23 2011
In terms of the cycle index of the symmetric group: Product_{q=1..m} [z^{v_q}] Z(S_v) 1/(1-z) where v is the maximum exponent of any prime in the prime factorization of n, v_q are the exponents of the prime factors, and Z(S_v) is the cycle index of the symmetric group on v elements. - Marko Riedel, Oct 03 2014
Dirichlet g.f.: Sum_{n >= 1} a(n)/n^s = Product_{k >= 1} zeta(ks) [Kendall]. - Álvar Ibeas, Nov 05 2014
a(n)=2 for all n in A054753 and for all n in A085987. a(n)=3 for all n in A030078 and for all n in A065036. a(n)=4 for all n in A085986. a(n)=5 for all n in A030514 and for all n in A178739. a(n)=6 for all n in A143610. - R. J. Mathar, Nov 05 2016
A050360(n) = a(A025487(n)). a(n) = A050360(A101296(n)). - R. J. Mathar, May 26 2017
a(n) = A000001(n) - A060689(n). - M. F. Hasler, Oct 24 2017
From Amiram Eldar, Nov 01 2020: (Start)
a(n) = a(A057521(n)).
Asymptotic mean: lim_{n->oo} (1/n) * Sum_{k=1..n} a(k) = A021002. (End)
a(n) = A005361(n) except when n is a term of A046101, since A000041(x) = x for x <= 3. - Miles Englezou, Feb 17 2024
Inverse Moebius transform of A188585: a(n) = Sum_{d|n} A188585(d). - Amiram Eldar, Jun 10 2025

A055887 Number of ordered partitions of partitions.

Original entry on oeis.org

1, 1, 3, 8, 22, 59, 160, 431, 1164, 3140, 8474, 22864, 61697, 166476, 449210, 1212113, 3270684, 8825376, 23813776, 64257396, 173387612, 467856828, 1262431711, 3406456212, 9191739970, 24802339472, 66924874539, 180585336876, 487278670744, 1314838220172
Offset: 0

Views

Author

Christian G. Bower, Jun 09 2000

Keywords

Comments

Jordan matrices are upper bidiagonal matrices such that (A) the diagonal entries are in sorted order, (B) there are only 1's and 0's on the superdiagonal, (C) for each superdiagonal 1, the two diagonal entries to the left and below it must be equal. Let J(N) be the number of N X N Jordan matrices where the diagonal values are, without loss of generality, taken to be a prefix of some fixed strictly increasing sequence x_1, x_2, x_3, ... If Jordan blocks sorted by eigenvalue with ties broken by block size during the sorting, then J(1, 2, 3, ...) is this sequence. - Warren D. Smith, Jan 28 2002
Number of compositions of n into parts k >= 1 where there are A000041(k) sorts of part k. - Joerg Arndt, Sep 30 2012
Also number of chains of multisets that partition a normal multiset of weight n, where a multiset is normal if it spans an initial interval of positive integers. - Gus Wiseman, Oct 28 2015
From Gus Wiseman, Jul 31 2022: (Start)
Also the number of ways to choose a multiset partition into constant multisets of a multiset of length n covering an initial interval of positive integers. This interpretation involves only multisets, not sequences. For example, the a(1) = 1 through a(3) = 8 multiset partitions are:
{{1}} {{1,1}} {{1,1,1}}
{{1},{1}} {{1},{1,1}}
{{1},{2}} {{1},{2,2}}
{{2},{1,1}}
{{1},{1},{1}}
{{1},{1},{2}}
{{1},{2},{2}}
{{1},{2},{3}}
Factorizations into prime powers, are counted by A000688.
The strongly normal case is A063834.
The strongly normal strict case is A270995.
Twice-partitions of type PPR are counted by A279784, factorizations A295935.
The strict case is A304969.
(End)

Examples

			The a(4) = 22 chains of multisets, where notation x-y means "y is a submultiset of x", are: (o-o-o-o) (oo-o-o) (oo-oo) (ooo-o) (oooo) (oe-o-o) (ooe-o) (oooe) (oe-oe) (ooe-e) (oee-o) (ooee) (oei-o) (ooei) (oe-e-e) (oee-e) (oeee) (oei-e) (oeei) (oei-i) (oeii) (oeis).
From _Gus Wiseman_, Jul 31 2022: (Start)
a(n) is the number of ways to choose an integer partition of each part of an integer composition of n. The a(0) = 1 through a(3) = 8 choices are:
  ()  ((1))  ((2))     ((3))
             ((11))    ((21))
             ((1)(1))  ((111))
                       ((1)(2))
                       ((2)(1))
                       ((1)(11))
                       ((11)(1))
                       ((1)(1)(1))
(End)
		

Crossrefs

Row sums of A060642.
Cf. A326346.
The unordered version is A001970, row-sums of A061260.
A000041 counts integer partitions, strict A000009.
A011782 counts integer compositions.
A072233 counts partitions by sum and length.

Programs

  • Maple
    with(combstruct); SeqSetSetU := [T, {T=Sequence(S), S=Set(U,card >= 1), U=Set(Z,card >=1)},unlabeled];
    P := (x) -> product( 1/(1-x^k), k=1..20 ) - 1; F := (x) -> series( 1/(1-P(x)) - 1, x, 21 ); # F(x) is g.f. for this sequence # Warren D. Smith, Jan 28 2002
    A055887rec:= proc(n::integer) local k; option remember; with(combinat): if n = 0 then 1 else add(numbpart(k) *procname(n - k), k=1..n); end if; end proc: seq (A055887rec(n), n=0..10); # Thomas Wieder, Nov 26 2007
  • Mathematica
    a = 1/Product[(1 - x^k), {k, 1, \[Infinity]}] - 1; CoefficientList[Series[1/(1 - a), {x, 0, 20}], x] (* Geoffrey Critzer, Dec 23 2010 *)
    (1/(2 - 1/QPochhammer[x]) + O[x]^30)[[3]] (* Vladimir Reshetnikov, Sep 22 2016 *)
    Table[Sum[Times@@PartitionsP/@c,{c,Join@@Permutations/@IntegerPartitions[n]}],{n,0,10}] (* Gus Wiseman, Jul 31 2022 *)
  • PARI
    Vec(1/(2-1/eta(x+O(x^66)))) \\ Joerg Arndt, Sep 30 2012

Formula

Invert transform of partitions numbers A000041.
Let p(k) be the number of integer partitions of k. Furthermore, set a(0)=1. Then a(n) = Sum_{k=1..n} p(k)*a(n-k). - Thomas Wieder, Nov 26 2007
G.f.: 1/( 1 - Sum_{k>=1} p(k)*x^k ) where p(k) = A000041(k) is the number of integer partitions of k. - Joerg Arndt, Sep 30 2012
a(n) ~ c * d^n, where d = 2.698329106474211231263998666188376330713465125913986356769... (see A246828) and c = 0.414113793172792357745578049739573823627306487211379286647... - Vaclav Kotesovec, Mar 29 2014

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.

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

A381636 Numbers whose prime indices cannot be partitioned into constant blocks with distinct sums.

Original entry on oeis.org

12, 60, 63, 84, 120, 126, 132, 156, 204, 228, 252, 276, 300, 315, 325, 348, 372, 420, 444, 492, 504, 516, 560, 564, 588, 630, 636, 650, 660, 693, 708, 720, 732, 780, 804, 819, 840, 852, 876, 924, 931, 948, 975, 996, 1008, 1020, 1068, 1071, 1092, 1140, 1164
Offset: 1

Views

Author

Gus Wiseman, Mar 10 2025

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 numbers that cannot be written as a product of prime powers > 1 with distinct sums of prime indices (A056239).
Contains no squarefree numbers.
Conjecture: These are the zeros of A382876.

Examples

			The prime indices of 300 are {1,1,2,3,3}, with partitions into constant blocks:
  {{2},{1,1},{3,3}}
  {{1},{1},{2},{3,3}}
  {{2},{3},{3},{1,1}}
  {{1},{1},{2},{3},{3}}
but none of these has distinct block-sums, so 300 is in the sequence.
The terms together with their prime indices begin:
   12: {1,1,2}
   60: {1,1,2,3}
   63: {2,2,4}
   84: {1,1,2,4}
  120: {1,1,1,2,3}
  126: {1,2,2,4}
  132: {1,1,2,5}
  156: {1,1,2,6}
  204: {1,1,2,7}
  228: {1,1,2,8}
  252: {1,1,2,2,4}
  276: {1,1,2,9}
  300: {1,1,2,3,3}
		

Crossrefs

More on multiset partitions into constant blocks: A006171, A279784, A295935.
These are the positions of 0 in A381635, after taking block-sums A381716.
Partitions of this type are counted by A381717.
For strict instead of constant blocks we have A381806, zeros of A381633.
For equal instead of distinct block-sums we have A381871.
A000688 counts multiset partitions into constant, see A381455 (upper), A381453 (lower).
A001055 counts multiset partitions, see A317141 (upper), A300383 (lower).
A050361 counts multiset partitions into distinct constant blocks, after sums A381715.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798.

Programs

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

A381635 Number of ways to partition 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 09 2025

Keywords

Comments

First differs from A381716 at a(1728) = 5, A381716(1728) = 4.
Also the number of factorizations on n into prime powers > 1 with distinct sums of prime indices (A056239).

Examples

			The a(432) = 3 multiset partitions:
  {{2,2,2},{1,1,1,1}}
  {{1},{1,1,1},{2,2,2}}
  {{1},{2},{2,2},{1,1,1}}
Note {{2},{2,2},{1,1,1,1}} is not included, as it does not have distinct block-sums.
		

Crossrefs

Without distinct block-sums we have A000688, after sums A381455 (upper), A381453 (lower).
For distinct blocks instead of sums we have A050361, after sums A381715.
For strict instead of constant we have A381633 (zeros A381806), after sums A381634.
Positions of 0 are A381636.
Taking block-sums (and sorting) gives A381716.
Other multiset partitions of prime indices:
More on multiset partitions into constant blocks: A006171, A279784, A295935.
A001055 counts multiset partitions, 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.

Programs

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

A382525 Number of times n appears in A048767 (rank of Look-and-Say partition of prime indices). Number of ordered set partitions whose block-sums are the prime signature of n.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Apr 05 2025

Keywords

Comments

The Look-and-Say partition of a multiset or partition y is obtained by interchanging parts with multiplicities. Hence, the multiplicity of k in the Look-and-Say partition of y is the sum of all parts that appear exactly k times. For example, starting with (3,2,2,1,1) we get (2,2,2,1,1,1), the multiset union of ((1,1,1),(2,2),(2)).
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.
Also the number of ways to choose a set of disjoint strict integer partitions, one of each nonzero multiplicity in the prime factorization of n.

Examples

			The a(27) = 2 partitions with Look-and-Say partition (2,2,2) are: (3,3), (2,2,1,1).
The prime indices of 3456 are {1,1,1,1,1,1,1,2,2,2}, and the partitions with Look-and-Say partition (2,2,2,1,1,1,1,1,1,1) are:
  (7,3,3)
  (7,2,2,1,1)
  (6,3,3,1)
  (5,3,3,2)
  (4,3,3,2,1)
  (4,3,2,2,1,1)
so a(3456) = 6.
		

Crossrefs

Positions of positive terms are A351294, conjugate A381432.
Positions of 0 are A351295, conjugate A381433.
Positions of 1 are A381540, conjugate A381434.
Positions of terms > 1 are A381541, conjugate A381435.
Positions of first appearances are A382775.
A000670 counts ordered set partitions.
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.
A239455 counts Look-and-Say partitions, complement A351293.
A381436 lists the section-sum partition of prime indices, ranks A381431.
A381440 lists the Look-and-Say partition of prime indices, ranks A048767.

Programs

  • Mathematica
    stp[y_]:=Select[Tuples[Select[IntegerPartitions[#],UnsameQ@@#&]&/@y],UnsameQ@@Join@@#&];
    Table[Length[stp[Last/@FactorInteger[n]]],{n,100}]

Formula

a(2^n) = A000009(n).
a(prime(n)) = 1.

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

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.
Showing 1-10 of 58 results. Next