A063834 -id:A063834 - cp's OEIS frontend

A289501 Number of enriched p-trees of weight n.

1, 1, 2, 4, 12, 32, 112, 352, 1296, 4448, 16640, 59968, 231168, 856960, 3334400, 12679424, 49991424, 192890880, 767229952, 2998427648, 12015527936, 47438950400, 191117033472, 760625733632, 3082675150848, 12346305839104, 50223511928832, 202359539335168

Offset: 0

Gus Wiseman, Jul 07 2017

nonn

An enriched p-tree of weight n is either (case 1) the number n itself, or (case 2) a sequence of two or more enriched p-trees, having a weakly decreasing sequence of weights summing to n.

			The a(4) = 12 enriched p-trees are:
  4,
  (31), ((21)1), (((11)1)1), ((111)1),
  (22), (2(11)), ((11)2), ((11)(11)),
  (211), ((11)11),
  (1111).

Alois P. Heinz, Table of n, a(n) for n = 0..1588

Cf. A052337, A063834, A093637, A196545, A273873, A281145, A300660.

b:= proc(n, i) option remember; `if`(n=0, 1,
      `if`(i<1, 0, b(n, i-1)+a(i)*b(n-i, min(n-i, i))))
    end:
a:= n-> `if`(n=0, 1, 1+b(n, n-1)):
seq(a(n), n=0..30);  # Alois P. Heinz, Jul 07 2017

a[n_]:=a[n]=1+Sum[Times@@a/@y,{y,Rest[IntegerPartitions[n]]}];
Array[a,20]
(* Second program: *)
b[n_, i_] := b[n, i] = If[n == 0, 1,
     If[i<1, 0, b[n, i-1] + a[i] b[n-i, Min[n-i, i]]]];
a[n_] := If[n == 0, 1, 1 + b[n, n-1]];
a /@ Range[0, 30] (* Jean-François Alcover, May 09 2021, after Alois P. Heinz *)

seq(n)={my(v=vector(n)); for(n=1, n, v[n] = 1 + polcoef(1/prod(k=1, n-1, 1 - v[k]*x^k + O(x*x^n)), n)); concat([1], v)} \\ Andrew Howroyd, Aug 26 2018

O.g.f.: (1/(1-x) + Product_{i>0} 1/(1-a(i)*x^i))/2.

A279785 Number of ways to choose a strict partition of each part of a strict partition of n.

Original entry on oeis.org

1, 1, 1, 3, 4, 7, 11, 18, 28, 47, 71, 108, 166, 252, 382, 587, 869, 1282, 1938, 2832, 4153, 6148, 8962, 12965, 18913, 27301, 39380, 56747, 81226, 115907, 166358, 236000, 334647, 475517, 671806, 947552, 1335679, 1875175, 2630584, 3687589, 5150585, 7183548

Offset: 0

Gus Wiseman, Dec 18 2016

nonn

This sequence is obtained from the generalized Euler transform in A266964 by taking f(n) = -1, g(n) = -A000009(n). - Seiichi Manyama, Nov 14 2018

			The a(6)=11 twice-partitions are:
((6)), ((5)(1)), ((51)), ((4)(2)), ((42)), ((41)(1)),
((3)(2)(1)), ((31)(2)), ((32)(1)), ((321)), ((21)(2)(1)).

Seiichi Manyama, Table of n, a(n) for n = 0..10000 (terms 0..5000 from Alois P. Heinz)
Gus Wiseman, Sequences enumerating triangles of integer partitions

Cf. A000009, A063834, A270995, A279375, A327553, A327605.

with(numtheory):
g:= proc(n) option remember; `if`(n=0, 1, add(add(
      `if`(d::odd, d, 0), d=divisors(j))*g(n-j), j=1..n)/n)
    end:
b:= proc(n, i) option remember; `if`(n>i*(i+1)/2, 0,
      `if`(n=0, 1, b(n, i-1)+`if`(i>n, 0, g(i)*b(n-i, i-1))))
    end:
a:= n-> b(n$2):
seq(a(n), n=0..70);  # Alois P. Heinz, Dec 20 2016

nn=20;CoefficientList[Series[Product[(1+PartitionsQ[k]x^k),{k,nn}],{x,0,nn}],x]
(* Second program: *)
g[n_] := g[n] = If[n==0, 1, Sum[Sum[If[OddQ[d], d, 0], {d, Divisors[j]}]* g[n - j], {j, 1, n}]/n]; b[n_, i_] := b[n, i] = If[n > i*(i + 1)/2, 0, If[n==0, 1, b[n, i-1] + If[i>n, 0, g[i]*b[n-i, i-1]]]]; a[n_] := b[n, n]; Table[a[n], {n, 0, 70}] (* Jean-François Alcover, Feb 07 2017, after Alois P. Heinz *)

G.f.: Product_{k>0} (1 + A000009(k)*x^k). - Seiichi Manyama, Nov 14 2018

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

Christian G. Bower, Jun 09 2000

nonn

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)

			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)

Alois P. Heinz, Table of n, a(n) for n = 0..2320
N. J. A. Sloane, Transforms
N. J. A. Sloane and Thomas Wieder, The Number of Hierarchical Orderings, arXiv:math/0307064 [math.CO], 2003; Order 21 (2004), 83-89, DOI:10.1007/s11083-004-9460-9.
Eric Weisstein's World of Mathematics, q-Pochhammer Symbol.

Cf. A055888, A083355, A095975, A246828.
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.
Cf. A001055, A050361, A063834, A279784, A355743.

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

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

Vec(1/(2-1/eta(x+O(x^66)))) \\ Joerg Arndt, Sep 30 2012

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
G.f.: 1/(2 - 1/(x)inf), where (x)_inf is the q-Pochhammer symbol. - _Vladimir Reshetnikov, Sep 22 2016

A133494 Diagonal of the array of iterated differences of A047848.

Original entry on oeis.org

1, 1, 3, 9, 27, 81, 243, 729, 2187, 6561, 19683, 59049, 177147, 531441, 1594323, 4782969, 14348907, 43046721, 129140163, 387420489, 1162261467, 3486784401, 10460353203, 31381059609, 94143178827, 282429536481, 847288609443, 2541865828329, 7625597484987, 22876792454961, 68630377364883

Offset: 0

Paul Barry, Paul Curtz, Dec 23 2007

a(n) is the number of ways to choose a composition C, and then choose a composition of each part of C. - Geoffrey Critzer, Mar 19 2012
a(n) is the top left entry of the n-th power of the 3 X 3 matrix [1, 1, 1; 1, 1, 1; 1, 1, 1]. - R. J. Mathar, Feb 03 2014
a(n) is the reptend length of 1/3^(n+1) in decimal. - Jianing Song, Nov 14 2018
Also the number of pairs of integer compositions, the first summing to n and the second with sum equal to the length of the first. If an integer composition is regarded as an arrow from sum to length, these are composable pairs, and the obvious composition operation founds a category of integer compositions. For example, we have (2,1,1,4) . (1,2,1) . (1,2) = (2,6), where dots represent the composition operation. The version without empty compositions is A000244. Composable triples are counted by 1 followed by A000302. The unordered version is A022811. - Gus Wiseman, Jul 14 2022

			From _Gus Wiseman_, Jul 15 2020: (Start)
The a(0) = 1 through a(3) = 9 ways to choose a composition of each part of a composition:
  ()  (1)  (2)      (3)
           (1,1)    (1,2)
           (1),(1)  (2,1)
                    (1,1,1)
                    (1),(2)
                    (2),(1)
                    (1),(1,1)
                    (1,1),(1)
                    (1),(1),(1)
(End)

G. C. Greubel, Table of n, a(n) for n = 0..1000
Bishal Deb, Cyclic sieving phenomena via combinatorics of continued fractions, arXiv:2508.13709 [math.CO], 2025. See p. 42.
Index entries for linear recurrences with constant coefficients, signature (3).

The strict version is A336139.
Splittings of partitions are A323583.
Multiset partitions of partitions are A001970.
Partitions of each part of a partition are A063834.
Compositions of each part of a partition are A075900.
Strict partitions of each part of a strict partition are A279785.
Compositions of each part of a strict partition are A304961.
Strict compositions of each part of a composition are A307068.
Compositions of each part of a strict composition are A336127.
Cf. A000012, A000244, A000302, A001906, A006951, A011782, A022811, A071919.
Cf. A072706, A078008, A112467, A182105, A316245, A336128, A336135.
Cf. A339006, A355382, A355384, A355387.

[n eq 0 select 1 else 3^(n-1): n in [0..30]]; // G. C. Greubel, Nov 20 2023

a:= n-> ceil(3^(n-1)):
seq(a(n), n=0..30);  # Alois P. Heinz, Jul 26 2020

CoefficientList[Series[(1 - 2 x)/(1 - 3 x), {x, 0, 50}], x] (* Vladimir Joseph Stephan Orlovsky, Jun 21 2011 *)
Join[{1}, 3^(Range[0, 30])] (* G. C. Greubel, Nov 20 2023 *)

a(n)=max(1,3^(n-1)) \\ Charles R Greathouse IV, Jul 07 2011

Vec((1-2*x)/(1-3*x) + O(x^100)) \\  Altug Alkan, Oct 30 2015

[(3^n + 2*int(n==0))//3 for n in range(31)] # G. C. Greubel, Nov 20 2023

Binomial transform of A078008. - Paul Curtz, Aug 04 2008
From R. J. Mathar, Nov 11 2008: (Start)
G.f.: (1 - 2*x)/(1 - 3*x).
a(n) = A000244(n-1), n > 0. (End)
From Philippe Deléham, Nov 13 2008: (Start)
a(n) = Sum_{k=0..n} A112467(n,k)*2^k.
a(n) = Sum_{k=0..n} A071919(n,k)*2^k. (End)
Let A(x) be the g.f. Then B(x) = x*A(x) satisfies B(x/(1-x)) = x/(1 - 2*B(x)). - Vladimir Kruchinin, Dec 05 2011
G.f.: 1/(1 - (Sum_{k>=1} (x/(1 - x))^k)). - Joerg Arndt, Sep 30 2012
For n > 0, a(n) = 2*(Sum_{k=0..n-1} a(k)) - 1 = 3^(n-1). - J. Conrad, Oct 29 2015
G.f.: 1 + x/(1 + x)*(1 + 4*x/(1 + 4*x)*(1 + 7*x/(1 + 7*x)*(1 + 10*x/(1 + 10*x)*(1 + .... - Peter Bala, May 27 2017
Invert transform of A011782(n) = 2^(n-1). Second invert transform of A000012. - Gus Wiseman, Jul 19 2020
a(n) = ceiling(3^(n-1)). - Alois P. Heinz, Jul 26 2020
From Elmo R. Oliveira, Mar 31 2025: (Start)
E.g.f.: (2 + exp(3*x))/3.
a(n) = 3*a(n-1) for n > 1. (End)

Definition clarified by R. J. Mathar, Nov 11 2008

A273873 Number of strict trees of weight n.

Original entry on oeis.org

1, 1, 2, 3, 6, 12, 28, 65, 166, 412, 1076, 2806, 7524, 20020, 54744, 148417, 410078, 1126732, 3144500, 8728570, 24555900, 68713420, 194469616, 548088278, 1559301428, 4418131108, 12628267512, 35957541462, 103150588492, 294924202032, 848878072440, 2435729999665

Offset: 1

Gus Wiseman, Jun 01 2016

nonn

A strict tree t is either (case 1) a positive integer t = n, or (case 2) a set t = {t1, t2, ..., tk} of two or more strict trees (i.e. branches) with distinct weights, where the weight of a strict tree in the second case is the sum of the weights of its branches; hence the multiset of weights is a strict integer partition of n. For example, {{{{{2,1},1},2},3},{4,{2,1}},{2,1},1} is a strict tree of weight 20.

			a(6) = 12: {6, (51), (42), ((41)1), (321), ((31)2), ((32)1), (((31)1)1), ((21)21), (((21)1)2), (((21)2)1), ((((21)1)1)1)}.

Alois P. Heinz, Table of n, a(n) for n = 1..2000
H. Gingold and A. Knopfmacher, Analytic properties of power product expansions, Canad. J. Math. 47(1995), 1219-1239
Gus Wiseman, Comcategories and Multiorders, (pdf version)
Gus Wiseman, All strict trees n=1..6.

Cf. A196545 (weakly ordered plane trees); A220418, A220420 (power product expansions); A271619, A063834 (twice partitioned numbers), A289501.

b:= proc(n, i) option remember; `if`(i*(i+1)/2 1+b(n, n-1):
seq(a(n), n=1..35);  # Alois P. Heinz, Jun 02 2016

STE[n_Integer?Positive]:=STE[n]=1+Plus@@Map[Function[ptn,Times@@STE/@ptn],Select[IntegerPartitions[n],And[Length[#]>1,UnsameQ@@#]&]];
Array[STE,30]
(* Second program: *)
b[n_, i_] := b[n, i] = If[i(i + 1)/2 < n, 0,
     If[n == 0, 1, b[n, i - 1] + b[n - i, Min[n - i, i - 1]] a[i]]];
a[n_] := If[n == 0, 1, 1 + b[n, n - 1]];
a /@ Range[35] (* Jean-François Alcover, May 09 2021, after Alois P. Heinz *)

seq(n)={my(v=vector(n)); for(n=1, n, v[n] = 1 + polcoef(prod(k=1, n-1, 1 + v[k]*x^k + O(x*x^n)), n)); v} \\ Andrew Howroyd, Aug 26 2018

Sum_{g(t)=y} (-1)^{d(t)} = mu(|y|<={y_1,...,y_k}), where mu is the Mobius function of the multiorder of integer partitions, g(t) is the multiset of leaves of a strict tree t, and d(t) is the number of branchings.

Strict trees are closely related to the coefficients appearing in a(i) = Sum_y c(y_1)*...*c(y_k) where Sum_i c(i)*x^i = Prod_i (1 + a(i)*x^i). The latter identity is the formal power product expansion (PPE) of an (ordinary) generating function.

A265947 Total size of all principal order ideals in the poset of integer partitions of n with the refinement order.

Original entry on oeis.org

1, 1, 3, 6, 14, 26, 55, 99, 192, 340, 619, 1063, 1873, 3129, 5308, 8718, 14385, 23116, 37346, 58949, 93294, 145131, 225623, 345833, 529976, 801675, 1211225, 1811558, 2703327, 3998289, 5901849, 8641160, 12623450, 18315370, 26503133, 38119289, 54691750, 78028166, 111041918, 157250528, 222105633

Offset: 0

Max Alekseyev, Dec 23 2015

nonn

a(n) is the number of refinement-ordered pairs of integer partitions of n. Every such pair (x,y) is a multiset union x and a multiset of sums y of some weakly ordered sequence of integer partitions, so this sequence is dominated by A063834 (twice partitioned numbers). - Gus Wiseman, May 01 2016

			a(4) = 14 ordered pairs of partitions: {(4,4), (4,22), (4,31), (4,211), (4,1111), (22,22), (22,211), (22,1111), (31,31), (31,211), (31,1111), (211,211), (211,1111), (1111,1111)}.

Jon Mark Perry et al., Counting refinements of partitions, Mathoverflow, 2015.

Cf. A001764, A002846, A213242, A213385, A213427, A063834.

def A265947(n):
    P = Posets.IntegerPartitions(n)
    return sum( len(P.order_ideal([p])) for p in P )

# Alternative:
def A265947(n):
    return Posets.IntegerPartitions(n).relations_number() # F. Chapoton, Feb 26 2020

A050361 Number of factorizations into distinct prime powers greater than 1.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1

Offset: 1

Christian G. Bower, Oct 15 1999

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

The number of unordered factorizations of n into 1 and exponentially odd prime powers, i.e., p^e where p is a prime and e is odd (A246551). - Amiram Eldar, Jun 12 2025

			From _Gus Wiseman_, Jul 30 2022: (Start)
The A000688(216) = 9 factorizations of 216 into prime powers are:
  (2*2*2*3*3*3)
  (2*2*2*3*9)
  (2*2*2*27)
  (2*3*3*3*4)
  (2*3*4*9)
  (2*4*27)
  (3*3*3*8)
  (3*8*9)
  (8*27)
Of these, the a(216) = 4 strict cases are:
  (2*3*4*9)
  (2*4*27)
  (3*8*9)
  (8*27)
(End)

Antti Karttunen, Table of n, a(n) for n = 1..100000 (first 10000 terms from Reinhard Zumkeller)
Index entries for sequences computed from exponents in factorization of n.

Cf. A000009, A050360, A050362, A050363, A050364, A101296, A246551.
Cf. A124010.
This is the strict case of A000688.
Positions of 1's are A004709, complement A046099.
The case of primes (instead of prime-powers) is A008966, non-strict A000012.
The non-strict additive version allowing 1's A023893, ranked by A302492.
The non-strict additive version is A023894, ranked by A355743.
The additive version (partitions) is A054685, ranked by A356065.
The additive version allowing 1's is A106244, ranked by A302496.
A001222 counts prime-power divisors.
A005117 lists all squarefree numbers.
A034699 gives maximal prime-power divisor.
A246655 lists all prime-powers (A000961 includes 1), towers A164336.
A296131 counts twice-factorizations of type PQR, non-strict A295935.
Cf. A001970, A002110, A025487, A055887, A063834, A076610, A085970, A279786, A302590, A302601, A354911, A355742.

a050361 = product . map a000009 . a124010_row
-- Reinhard Zumkeller, Aug 28 2014

A050361 := proc(n)
    local a,f;
    if n = 1 then
        1;
    else
        a := 1 ;
        for f in ifactors(n)[2] do
            a := a*A000009(op(2,f)) ;
        end do:
    end if;
end proc: # R. J. Mathar, May 25 2017

Table[Times @@ PartitionsQ[Last /@ FactorInteger[n]], {n, 99}] (* Arkadiusz Wesolowski, Feb 27 2017 *)

A000009(n,k=(n-!(n%2))) = if(!n,1,my(s=0); while(k >= 1, if(k<=n, s += A000009(n-k,k)); k -= 2); (s));
A050361(n) = factorback(apply(A000009,factor(n)[,2])); \\ Antti Karttunen, Nov 17 2019

Dirichlet g.f.: Product_{n is a prime power >1}(1 + 1/n^s).
Multiplicative with a(p^e) = A000009(e).
a(A002110(k))=1.
a(n) = A050362(A101296(n)). - R. J. Mathar, May 26 2017
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Product_{p prime} f(1/p) = 1.26020571070524171076..., where f(x) = (1-x) * Product_{k>=1} (1 + x^k). - Amiram Eldar, Oct 03 2023

A023893 Number of partitions of n into prime power parts (1 included); number of nonisomorphic Abelian subgroups of symmetric group S_n.

Original entry on oeis.org

1, 1, 2, 3, 5, 7, 10, 14, 20, 27, 36, 48, 63, 82, 105, 134, 171, 215, 269, 335, 415, 511, 626, 764, 929, 1125, 1356, 1631, 1953, 2333, 2776, 3296, 3903, 4608, 5427, 6377, 7476, 8744, 10205, 11886, 13818, 16032, 18565, 21463, 24768, 28536

Offset: 0

Olivier Gérard

nonn

			From _Gus Wiseman_, Jul 28 2022: (Start)
The a(0) = 1 through a(6) = 10 partitions:
  ()  (1)  (2)   (3)    (4)     (5)      (33)
           (11)  (21)   (22)    (32)     (42)
                 (111)  (31)    (41)     (51)
                        (211)   (221)    (222)
                        (1111)  (311)    (321)
                                (2111)   (411)
                                (11111)  (2211)
                                         (3111)
                                         (21111)
                                         (111111)
(End)

Seiichi Manyama, Table of n, a(n) for n = 0..10000
Index entries for sequences related to groups

Cf. A009490, A023894 (first differences), A062297 (number of Abelian subgroups).
Cf. A018819, A062051, A131995.
The multiplicative version (factorizations) is A000688.
Not allowing 1's gives A023894, strict A054685, ranked by A355743.
The version for just primes (not prime-powers) is A034891, strict A036497.
The strict version is A106244.
These partitions are ranked by A302492.
A000041 counts partitions, strict A000009.
A001222 counts prime-power divisors.
A072233 counts partitions by sum and length.
A246655 lists the prime-powers (A000961 includes 1), towers A164336.
Cf. A001970, A055887, A063834, A085970, A279784, A295935, A356065.

Table[Length[Select[IntegerPartitions[n],Count[Map[Length,FactorInteger[#]], 1] == Length[#] &]], {n, 0, 35}] (* Geoffrey Critzer, Oct 25 2015 *)
nmax = 50; Clear[P]; P[m_] := P[m] = Product[Product[1/(1-x^(p^k)), {k, 1, m}], {p, Prime[Range[PrimePi[nmax]]]}]/(1-x)+O[x]^nmax // CoefficientList[ #, x]&; P[1]; P[m=2]; While[P[m] != P[m-1], m++]; P[m] (* Jean-François Alcover, Aug 31 2016 *)

lista(m) = {x = t + t*O(t^m); gf = prod(k=1, m, if (isprimepower(k), 1/(1-x^k), 1))/(1-x); for (n=0, m, print1(polcoeff(gf, n, t), ", "));} \\ Michel Marcus, Mar 09 2013

from functools import lru_cache
from sympy import factorint
@lru_cache(maxsize=None)
def A023893(n):
    @lru_cache(maxsize=None)
    def c(n): return sum((p**(e+1)-p)//(p-1) for p,e in factorint(n).items())+1
    return (c(n)+sum(c(k)*A023893(n-k) for k in range(1,n)))//n if n else 1 # Chai Wah Wu, Jul 15 2024

G.f.: (Product_{p prime} Product_{k>=1} 1/(1-x^(p^k))) / (1-x).

A023894 Number of partitions of n into prime power parts (1 excluded).

Original entry on oeis.org

1, 0, 1, 1, 2, 2, 3, 4, 6, 7, 9, 12, 15, 19, 23, 29, 37, 44, 54, 66, 80, 96, 115, 138, 165, 196, 231, 275, 322, 380, 443, 520, 607, 705, 819, 950, 1099, 1268, 1461, 1681, 1932, 2214, 2533, 2898, 3305, 3768, 4285, 4872, 5530, 6267, 7094, 8022, 9060

Offset: 0

Olivier Gérard

nonn

			From _Gus Wiseman_, Jul 28 2022: (Start)
The a(0) = 1 through a(9) = 7 partitions:
  ()  .  (2)  (3)  (4)   (5)   (33)   (7)    (8)     (9)
                   (22)  (32)  (42)   (43)   (44)    (54)
                               (222)  (52)   (53)    (72)
                                      (322)  (332)   (333)
                                             (422)   (432)
                                             (2222)  (522)
                                                     (3222)
(End)

Seiichi Manyama, Table of n, a(n) for n = 0..10000
E. Grosswald, Partitions into prime powers

The multiplicative version (factorizations) is A000688, coprime A354911.
Allowing 1's gives A023893, strict A106244, ranked by A302492.
The strict version is A054685.
The version for just primes is ranked by A076610, squarefree A356065.
Twice-partitions of this type are counted by A279784, factorizations A295935.
These partitions are ranked by A355743.
A000041 counts partitions, strict A000009.
A001222 counts prime-power divisors.
A072233 counts partitions by sum and length.
A246655 lists the prime-powers (A000961 includes 1), towers A164336.
Cf. A001970, A055887, A063834, A085970.

Table[Length[Select[IntegerPartitions[n],And@@PrimePowerQ/@#&]],{n,0,30}] (* Gus Wiseman, Jul 28 2022 *)

is_primepower(n)= {ispower(n, , &n); isprime(n)}
lista(m) = {x = t + t*O(t^m); gf = prod(k=1, m, if (is_primepower(k), 1/(1-x^k), 1)); for (n=0, m, print1(polcoeff(gf, n, t), ", "));}
\\ Michel Marcus, Mar 09 2013

from functools import lru_cache
from sympy import factorint
@lru_cache(maxsize=None)
def A023894(n):
    @lru_cache(maxsize=None)
    def c(n): return sum((p**(e+1)-p)//(p-1) for p,e in factorint(n).items())
    return (c(n)+sum(c(k)*A023894(n-k) for k in range(1,n)))//n if n else 1 # Chai Wah Wu, Jul 15 2024

G.f.: Prod(p prime, Prod(k >= 1, 1/(1-x^(p^k))))

A271619 Twice partitioned numbers where the first partition is strict.

Original entry on oeis.org

1, 1, 2, 5, 8, 18, 34, 65, 109, 223, 386, 698, 1241, 2180, 3804, 6788, 11390, 19572, 34063, 56826, 96748, 163511, 272898, 452155, 755928, 1244732, 2054710, 3382147, 5534696, 8992209, 14733292, 23763685, 38430071, 62139578, 99735806, 160183001, 256682598

Offset: 0

Gus Wiseman, Apr 10 2016

nonn

Number of sequences of integer partitions of the parts of some strict partition of n.

This sequence is obtained from the generalized Euler transform in A266964 by taking f(n) = -1, g(n) = -A000041(n). - Seiichi Manyama, Nov 15 2018

			a(6)=34: {(6);(5)(1),(51);(4)(2),(42);(4)(11),(41)(1),(411);(33);(3)(2)(1),(31)(2),(32)(1),(321);(3)(11)(1),(31)(11),(311)(1),(3111);(22)(2),(222);(21)(2)(1),(22)(11),(211)(2),(221)(1),(2211);(21)(11)(1),(111)(2)(1),(211)(11),(1111)(2),(2111)(1),(21111);(111)(11)(1),(1111)(11),(11111)(1),(111111)}

Seiichi Manyama, Table of n, a(n) for n = 0..10000 (terms 0..5000 from Alois P. Heinz)

Cf. A000009, A000041, A063834 (twice partitioned numbers), A270995, A279785, A327552, A327607.

b:= proc(n, i) option remember; `if`(n>i*(i+1)/2, 0,
      `if`(n=0, 1, b(n, i-1) +`if`(i>n, 0,
       b(n-i, i-1)*combinat[numbpart](i))))
    end:
a:= n-> b(n$2):
seq(a(n), n=0..50);  # Alois P. Heinz, Apr 11 2016

With[{n = 50}, CoefficientList[Series[Product[(1 + PartitionsP[i] x^i), {i, 1, n}], {x, 0, n}], x]]

G.f.: Product_{i>=1} (1 + A000041(i) * x^i).

cp's OEIS Frontend

A289501 Number of enriched p-trees of weight n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A279785 Number of ways to choose a strict partition of each part of a strict partition of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A055887 Number of ordered partitions of partitions.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A133494 Diagonal of the array of iterated differences of A047848.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

PARI

SageMath

Formula

Extensions

A273873 Number of strict trees of weight n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A265947 Total size of all principal order ideals in the poset of integer partitions of n with the refinement order.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Sage

Sage