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 11-20 of 102 results. Next

A279375 Number of set partitions of strict integer partitions of n that have distinct block-sums.

Original entry on oeis.org

1, 1, 1, 3, 3, 5, 9, 12, 16, 24, 39, 49, 70, 94, 127, 202, 247, 340, 450, 606, 772, 1169, 1407, 1920, 2454, 3267, 4089, 5469, 7293, 9222, 11884, 15291, 19417, 24890, 31469, 39662, 52619, 64764, 82502, 103576, 131169, 162726, 206015, 254233, 318464, 406262, 499210, 620593, 773673, 957073, 1181593
Offset: 0

Views

Author

Gus Wiseman, Dec 11 2016

Keywords

Comments

Also twice partitioned numbers where all partitions are strict. Also triangles of weight n in the multisystem of strict partitions. Strict partitions are an example of a multisystem that is neither transitive nor partitive nor contractible but is decomposable; see link for details.

Examples

			The a(6)=9 set partitions of strict integer partitions of 6 are: ((6)), ((51)), ((5)(1)), ((42)), ((4)(2)), ((321)), ((32)(1)), ((31)(2)), ((3)(2)(1)). The set partition ((3)(21)) is not counted because its blocks do not have distinct sums.
		

Crossrefs

Programs

  • Mathematica
    nn=20;sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    Table[Total[Length[Select[sps[Sort[#]],UnsameQ@@Total/@#&]]&/@Select[IntegerPartitions[n],UnsameQ@@#&]],{n,nn}]

A296119 Number of ways to choose a strict factorization of each factor in a factorization of n.

Original entry on oeis.org

1, 1, 1, 2, 1, 3, 1, 4, 2, 3, 1, 7, 1, 3, 3, 7, 1, 7, 1, 7, 3, 3, 1, 16, 2, 3, 4, 7, 1, 12, 1, 12, 3, 3, 3, 21, 1, 3, 3, 16, 1, 12, 1, 7, 7, 3, 1, 33, 2, 7, 3, 7, 1, 16, 3, 16, 3, 3, 1, 34, 1, 3, 7, 23, 3, 12, 1, 7, 3, 12, 1, 50, 1, 3, 7, 7, 3, 12, 1, 33, 7, 3
Offset: 1

Views

Author

Gus Wiseman, Dec 05 2017

Keywords

Examples

			The a(24) = 16 twice-factorizations:
(2)*(2)*(2)*(3),
(2)*(2)*(2*3), (2)*(2)*(6), (2)*(3)*(4),
(2)*(2*6), (2)*(3*4), (2)*(12), (3)*(2*4), (3)*(8), (4)*(2*3), (4)*(6),
(2*3*4), (2*12), (3*8), (4*6), (24).
		

Crossrefs

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Sum[Times@@(Length[Select[facs[#],UnsameQ@@#&]]&/@fac),{fac,facs[n]}],{n,100}]
  • PARI
    A045778(n, m=n) = ((n<=m) + sumdiv(n, d, if((d>1)&&(d<=m)&&(dA045778(n/d, d-1))));
    A296119(n, m=n) = if(1==n, 1, sumdiv(n, d, if((d>1)&&(d<=m), A045778(d)*A296119(n/d, d)))); \\ Antti Karttunen, Oct 08 2018

Formula

Dirichlet g.f.: 1/Product_{n > 1}(1 - A045778(n)/n^s).

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

A072706 Number of unimodal partitions/compositions of n into distinct terms.

Original entry on oeis.org

1, 1, 1, 3, 3, 5, 9, 11, 15, 21, 33, 39, 55, 69, 93, 127, 159, 201, 261, 327, 411, 537, 653, 819, 1011, 1257, 1529, 1899, 2331, 2829, 3441, 4179, 5031, 6093, 7305, 8767, 10575, 12573, 14997, 17847, 21223, 25089, 29757, 35055, 41379, 48801, 57285, 67131
Offset: 0

Views

Author

Henry Bottomley, Jul 04 2002

Keywords

Comments

Also the number of ways to partition a strict integer partition of n into two unordered blocks. - Gus Wiseman, Dec 31 2019

Examples

			a(6)=9 since 6 can be written as 1+2+3, 1+3+2, 1+5, 2+3+1, 2+4, 3+2+1, 4+2, 5+1, or 6, but not for example 1+4+1 (which does not have distinct terms) nor 2+1+3 (which is not unimodal).
From _Joerg Arndt_, Mar 25 2014: (Start)
The a(10) = 33 such compositions of 10 are:
01:  [ 1 2 3 4 ]
02:  [ 1 2 4 3 ]
03:  [ 1 2 7 ]
04:  [ 1 3 4 2 ]
05:  [ 1 3 6 ]
06:  [ 1 4 3 2 ]
07:  [ 1 4 5 ]
08:  [ 1 5 4 ]
09:  [ 1 6 3 ]
10:  [ 1 7 2 ]
11:  [ 1 9 ]
12:  [ 2 3 4 1 ]
13:  [ 2 3 5 ]
14:  [ 2 4 3 1 ]
15:  [ 2 5 3 ]
16:  [ 2 7 1 ]
17:  [ 2 8 ]
18:  [ 3 4 2 1 ]
19:  [ 3 5 2 ]
20:  [ 3 6 1 ]
21:  [ 3 7 ]
22:  [ 4 3 2 1 ]
23:  [ 4 5 1 ]
24:  [ 4 6 ]
25:  [ 5 3 2 ]
26:  [ 5 4 1 ]
27:  [ 6 3 1 ]
28:  [ 6 4 ]
29:  [ 7 2 1 ]
30:  [ 7 3 ]
31:  [ 8 2 ]
32:  [ 9 1 ]
33:  [ 10 ]
(End)
		

Crossrefs

The non-strict version is A001523.

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n>i*(i+1)/2, 0, `if`(n=0, 1,
          expand(b(n, i-1) +`if`(i>n, 0, x*b(n-i, i-1)))))
        end:
    a:= n->(p->add(coeff(p, x, i)*ceil(2^(i-1)), i=0..degree(p)))(b(n$2)):
    seq(a(n), n=0..100);  # Alois P. Heinz, Mar 25 2014
  • Mathematica
    b[n_, i_] := b[n, i] = If[n > i*(i + 1)/2, 0, If[n == 0, 1, Expand[b[n, i - 1] + If[i > n, 0, x*b[n - i, i - 1]]]]]; a[n_] := Function[{p}, Sum[Coefficient[p, x, i]*Ceiling[2^(i - 1)], {i, 0, Exponent[p, x]}]][b[n, n]]; Table[a[n], {n, 0, 100}] (* Jean-François Alcover, Jan 16 2015, after Alois P. Heinz *)
    Table[If[n==0,1,Sum[2^(Length[ptn]-1),{ptn,Select[IntegerPartitions[n],UnsameQ@@#&]}]],{n,0,15}] (* Gus Wiseman, Dec 31 2019 *)
  • PARI
    N=66; q='q+O('q^N); Vec( 1 + sum(n=1, N, 2^(n-1)*q^(n*(n+1)/2) / prod(k=1, n, 1-q^k ) ) ) \\ Joerg Arndt, Mar 25 2014

Formula

a(n) = sum_k A072705(n, k) = A032020(n)-A072707(k) = A032302(n)/2 (n>0).
G.f.: 1/2*(1+Product_{k>0} (1+2*x^k)). - Vladeta Jovovic, Jun 24 2003
G.f.: 1 + sum(n>=1, 2^(n-1)*q^(n*(n+1)/2) / prod(k=1..n, 1-q^k ) ). [Joerg Arndt, Jan 20 2014]
a(n) ~ c^(1/4) * exp(2*sqrt(c*n)) / (4*sqrt(3*Pi)*n^(3/4)), where c = -polylog(2, -2) = A266576 = 1.436746366883680946362902023893583354... - Vaclav Kotesovec, Sep 22 2019

A279784 Twice partitioned numbers where the latter partitions are constant.

Original entry on oeis.org

1, 1, 3, 5, 12, 18, 40, 60, 121, 186, 344, 524, 955, 1432, 2484, 3756, 6352, 9493, 15750, 23414, 38128, 56513, 90406, 133312, 211194, 309657, 484214, 708267, 1097159, 1597290, 2454245, 3560444, 5430091, 7854174, 11894335, 17151394, 25838413, 37145198, 55648059
Offset: 0

Views

Author

Gus Wiseman, Dec 18 2016

Keywords

Comments

Also number of ways to choose a divisor of each part of an integer partition of n.

Examples

			The a(4)=12 twice-partitions are:
((4)), ((3)(1)), ((2)(2)), ((22)),
((2)(1)(1)), ((2)(11)), ((11)(2)),
((1)(1)(1)(1)), ((11)(1)(1)), ((11)(11)), ((111)(1)), ((1111)).
		

Crossrefs

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n=0 or i=1, 1,
          b(n, i-1)+`if`(i>n, 0, numtheory[tau](i)*b(n-i, i)))
        end:
    a:= n-> b(n$2):
    seq(a(n), n=0..50);  # Alois P. Heinz, Dec 20 2016
  • Mathematica
    nn=20;CoefficientList[Series[Product[1/(1-DivisorSigma[0,n]x^n),{n,nn}],{x,0,nn}],x]

Formula

G.f.: exp(Sum_{k>=1} Sum_{j>=1} d(j)^k*x^(j*k)/k), where d(j) is the number of the divisors of j (A000005). - Ilya Gutkovskiy, Jul 17 2018
From Vaclav Kotesovec, Jul 28 2018: (Start)
a(n) ~ c * 2^(n/2), where
c = 203.986136154799274492709451797084688042886818134781591... if n is even and
c = 201.491703180375661735217350021245093454724452720559762... if n is odd.
In closed form, a(n) ~ ((2 + sqrt(2)) * Product_{k>=3} (1/(1 - tau(k) / 2^(k/2))) + (-1)^n * (2 - sqrt(2)) * Product_{k>=3} (1/(1 - (-1)^k * tau(k) / 2^(k/2)))) * 2^(n/2 - 1), where tau() is A000005. (End)

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.

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

A381431 Heinz number of the section-sum partition of the prime indices of n.

Original entry on oeis.org

1, 2, 3, 4, 5, 5, 7, 8, 9, 7, 11, 10, 13, 11, 11, 16, 17, 15, 19, 14, 13, 13, 23, 20, 25, 17, 27, 22, 29, 13, 31, 32, 17, 19, 17, 25, 37, 23, 19, 28, 41, 17, 43, 26, 33, 29, 47, 40, 49, 35, 23, 34, 53, 45, 19, 44, 29, 31, 59, 26, 61, 37, 39, 64, 23, 19, 67, 38
Offset: 1

Views

Author

Gus Wiseman, Feb 26 2025

Keywords

Comments

The image first differs from A320340, A364347, A350838 in containing a(150) = 65.
The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). This gives a bijective correspondence between positive integers and integer partitions.
The section-sum partition (A381436) of a multiset or partition y is defined as follows: (1) determine and remember the sum of all distinct parts, (2) remove one instance of each distinct part, (3) repeat until no parts are left. The remembered values comprise the section-sum partition. For example, starting with (3,2,2,1,1) we get (6,3).
Equivalently, the k-th part of the section-sum partition is the sum of all (distinct) parts that appear at least k times. Compare to the definition of the conjugate of a partition, where we count parts >= k.
The conjugate of a section-sum partition is a Look-and-Say partition; see A048767, union A351294, count A239455.

Examples

			Prime indices of 180 are (3,2,2,1,1), with section-sum partition (6,3), so a(180) = 65.
The terms together with their prime indices begin:
   1: {}
   2: {1}
   3: {2}
   4: {1,1}
   5: {3}
   5: {3}
   7: {4}
   8: {1,1,1}
   9: {2,2}
   7: {4}
  11: {5}
  10: {1,3}
  13: {6}
  11: {5}
  11: {5}
  16: {1,1,1,1}
		

Crossrefs

The conjugate is A048767, union A351294, complement A351295, fix A048768 (count A217605).
Taking length instead of sum in the definition gives A238745, conjugate A181819.
Partitions of this type are counted by A239455, complement A351293.
The union is A381432, complement A381433.
Values appearing only once are A381434, more than once A381435.
These are the Heinz numbers of rows of A381436, conjugate A381440.
Greatest prime index of each term is A381437, counted by A381438.
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, counted by A001222.
A122111 represents conjugation in terms of Heinz numbers.
Set multipartitions: A050320, A089259, A116540, A270995, A296119, A318360, A318361.
Partition ideals: A300383, A317141, A381078, A381441, A381452, A381454.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    egs[y_]:=If[y=={},{},Table[Total[Select[Union[y],Count[y,#]>=i&]],{i,Max@@Length/@Split[y]}]];
    Table[Times@@Prime/@egs[prix[n]],{n,100}]

Formula

A122111(a(n)) = A048767(n).

A293243 Numbers that cannot be written as a product of distinct squarefree numbers.

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, Oct 03 2017

Keywords

Comments

First differs from A212164 at a(441).
Numbers n such that A050326(n) = 0. - Felix Fröhlich, Oct 04 2017
Includes A246547, and all numbers of the form p^a*q^b where p and q are primes, a >= 1 and b >= 3. - Robert Israel, Oct 10 2017
Also numbers whose prime indices cannot be partitioned into a set of sets. For example, the prime indices of 90 are {1,2,2,3}, and we have sets of sets: {{2},{1,2,3}}, {{1,2},{2,3}}, {{1},{2},{2,3}}, {{2},{3},{1,2}}, so 90 is not in the sequence. - Gus Wiseman, Apr 28 2025

Examples

			120 is not in the sequence because 120 = 2*6*10. 3600 is not in the sequence because 3600 = 2*6*10*30.
		

Crossrefs

These are the zeros of A050326.
Multiset partitions of this type (set of sets) are counted by A050342.
Twice-partitions of this type (set of sets) are counted by A279785, see also A358914.
Normal multisets of this type are counted by A292432, A292444, A381996, A382214.
The case of a unique choice is A293511, counted by A382079.
For distinct block-sums instead of blocks see A381806, A381990, A381992, A382075.
Partitions of this type are counted by A382078.
The complement is A382200, counted by A382077.
A001055 counts factorizations, strict A045778.
A050320 counts factorizations into squarefree numbers.
A050345 counts factorizations partitioned into into distinct sets.
A317141 counts coarsenings of prime indices, refinements A300383.

Programs

  • Maple
    N:= 1000: # to get all terms <= N
    A:= Vector(N):
    A[1]:= 1:
    for n from 2 to N do
      if numtheory:-issqrfree(n) then
          S:= [$1..N/n]; T:= n*S; A[T]:= A[T]+A[S]
        fi;
    od:
    select(t -> A[t]=0, [$1..N]); # Robert Israel, Oct 10 2017
  • Mathematica
    nn=500;
    sqfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[sqfacs[n/d],Min@@#>d&]],{d,Select[Rest[Divisors[n]],SquareFreeQ]}]];
    Select[Range[nn],Length[sqfacs[#]]===0&]

A302494 Products of distinct primes of squarefree index.

Original entry on oeis.org

1, 2, 3, 5, 6, 10, 11, 13, 15, 17, 22, 26, 29, 30, 31, 33, 34, 39, 41, 43, 47, 51, 55, 58, 59, 62, 65, 66, 67, 73, 78, 79, 82, 83, 85, 86, 87, 93, 94, 101, 102, 109, 110, 113, 118, 123, 127, 129, 130, 134, 137, 139, 141, 143, 145, 146, 149, 155, 157, 158, 163
Offset: 1

Views

Author

Gus Wiseman, Apr 08 2018

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n.

Examples

			Entry A302242 describes a correspondence between positive integers and multiset multisystems. In this case it gives the following sequence of set systems.
01: {}
02: {{}}
03: {{1}}
05: {{2}}
06: {{},{1}}
10: {{},{2}}
11: {{3}}
13: {{1,2}}
15: {{1},{2}}
17: {{4}}
22: {{},{3}}
26: {{},{1,2}}
29: {{1,3}}
30: {{},{1},{2}}
31: {{5}}
33: {{1},{3}}
34: {{},{4}}
39: {{1},{1,2}}
		

Crossrefs

Programs

  • Mathematica
    Select[Range[100],Or[#===1,SquareFreeQ[#]&&And@@SquareFreeQ/@PrimePi/@FactorInteger[#][[All,1]]]&]
  • PARI
    is(n) = if(bigomega(n)!=omega(n), return(0), my(f=factor(n)[, 1]~); for(k=1, #f, if(!issquarefree(primepi(f[k])) && primepi(f[k])!=1, return(0)))); 1 \\ Felix Fröhlich, Apr 10 2018

A358914 Number of twice-partitions of n into distinct strict partitions.

Original entry on oeis.org

1, 1, 1, 3, 4, 7, 13, 20, 32, 51, 83, 130, 206, 320, 496, 759, 1171, 1786, 2714, 4104, 6193, 9286, 13920, 20737, 30865, 45721, 67632, 99683, 146604, 214865, 314782, 459136, 668867, 972425, 1410458, 2040894, 2950839, 4253713, 6123836, 8801349, 12627079
Offset: 0

Views

Author

Gus Wiseman, Dec 11 2022

Keywords

Comments

A twice-partition of n (A063834) is a sequence of integer partitions, one of each part of an integer partition of n.

Examples

			The a(1) = 1 through a(6) = 13 twice-partitions:
  ((1))  ((2))  ((3))     ((4))      ((5))      ((6))
                ((21))    ((31))     ((32))     ((42))
                ((2)(1))  ((3)(1))   ((41))     ((51))
                          ((21)(1))  ((3)(2))   ((321))
                                     ((4)(1))   ((4)(2))
                                     ((21)(2))  ((5)(1))
                                     ((31)(1))  ((21)(3))
                                                ((31)(2))
                                                ((3)(21))
                                                ((32)(1))
                                                ((41)(1))
                                                ((3)(2)(1))
                                                ((21)(2)(1))
		

Crossrefs

The unordered version is A050342, non-strict A261049.
This is the distinct case of A270995.
The case of strictly decreasing sums is A279785.
The case of constant sums is A279791.
For distinct instead of weakly decreasing sums we have A336343.
This is the twice-partition case of A358913.
A001970 counts multiset partitions of integer partitions.
A055887 counts sequences of partitions.
A063834 counts twice-partitions.
A330462 counts set systems by total sum and length.
A358830 counts twice-partitions with distinct lengths.

Programs

  • Mathematica
    twiptn[n_]:=Join@@Table[Tuples[IntegerPartitions/@ptn],{ptn,IntegerPartitions[n]}];
    Table[Length[Select[twiptn[n],UnsameQ@@#&&And@@UnsameQ@@@#&]],{n,0,10}]
  • PARI
    seq(n,k)={my(u=Vec(eta(x^2 + O(x*x^n))/eta(x + O(x*x^n))-1)); Vec(prod(k=1, n, my(c=u[k]); sum(j=0, min(c,n\k), x^(j*k)*c!/(c-j)!,  O(x*x^n))))} \\ Andrew Howroyd, Dec 31 2022

Extensions

Terms a(26) and beyond from Andrew Howroyd, Dec 31 2022
Previous Showing 11-20 of 102 results. Next