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

A296150 Triangle whose n-th row is the integer partition with Heinz number n.

Original entry on oeis.org

1, 2, 1, 1, 3, 2, 1, 4, 1, 1, 1, 2, 2, 3, 1, 5, 2, 1, 1, 6, 4, 1, 3, 2, 1, 1, 1, 1, 7, 2, 2, 1, 8, 3, 1, 1, 4, 2, 5, 1, 9, 2, 1, 1, 1, 3, 3, 6, 1, 2, 2, 2, 4, 1, 1, 10, 3, 2, 1, 11, 1, 1, 1, 1, 1, 5, 2, 7, 1, 4, 3, 2, 2, 1, 1, 12, 8, 1, 6, 2, 3, 1, 1, 1, 13, 4
Offset: 1

Views

Author

Gus Wiseman, Feb 05 2018

Keywords

Comments

Same as A112798 with rows reversed. Row lengths are A001222. Rows sums are A056239.
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Examples

			Sequence of partitions begins: (), (1), (2), (11), (3), (21), (4), (111), (22), (31), (5), (211), (6), (41), (32), (1111), (7), (221).
		

Crossrefs

Programs

  • Maple
    f := n -> op(map(numtheory:-pi, sort(map(`$`@op, ifactors(n)[2]), `>`))):
    map(f, [$1..100]); # Robert Israel, Feb 09 2018
  • Mathematica
    Table[If[n===1,{},Join@@Cases[FactorInteger[n]//Reverse,{p_,k_}:>Table[PrimePi[p],{k}]]],{n,50}]

A317141 In the ranked poset of integer partitions ordered by refinement, number of integer partitions coarser (greater) than or equal to the integer partition with Heinz number n.

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 3, 2, 2, 1, 4, 1, 2, 2, 5, 1, 4, 1, 4, 2, 2, 1, 6, 2, 2, 3, 4, 1, 5, 1, 7, 2, 2, 2, 8, 1, 2, 2, 7, 1, 5, 1, 4, 4, 2, 1, 10, 2, 4, 2, 4, 1, 7, 2, 7, 2, 2, 1, 9, 1, 2, 4, 11, 2, 5, 1, 4, 2, 5, 1, 12, 1, 2, 4, 4, 2, 5, 1, 11, 5, 2, 1, 10, 2
Offset: 1

Views

Author

Gus Wiseman, Jul 22 2018

Keywords

Comments

The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).

Examples

			The a(24) = 6 partitions coarser than or equal to (2111) are (2111), (311), (221), (32), (41), (5), with Heinz numbers 24, 20, 18, 15, 14, 11.
		

Crossrefs

Programs

  • Maple
    g:= l-> `if`(l=[], {[]}, (t-> map(sort, map(x->
            [seq(subsop(i=x[i]+t, x), i=1..nops(x)),
            [x[], t]][], g(subsop(-1=[][], l)))))(l[-1])):
    a:= n-> nops(g(map(i-> numtheory[pi](i[1])$i[2], ifactors(n)[2]))):
    seq(a(n), n=1..100);  # Alois P. Heinz, Jul 22 2018
  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    mps[set_]:=Union[Sort[Sort/@(#/.x_Integer:>set[[x]])]&/@sps[Range[Length[set]]]];
    ptncaps[ptn_]:=Union[Sort/@Apply[Plus,mps[ptn],{2}]];
    Table[Length[ptncaps[primeMS[n]]],{n,100}]

A299202 Moebius function of the multiorder of integer partitions indexed by their Heinz numbers.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Feb 05 2018

Keywords

Comments

By convention, mu() = 0.
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Examples

			Heinz number of (2,1,1) is 12, so mu(2,1,1) = a(12) = 2.
		

Crossrefs

Programs

  • Mathematica
    nn=120;
    ptns=Table[If[n===1,{},Join@@Cases[FactorInteger[n]//Reverse,{p_,k_}:>Table[PrimePi[p],{k}]]],{n,nn}];
    tris=Join@@Map[Tuples[IntegerPartitions/@#]&,ptns];
    mu[y_]:=mu[y]=If[Length[y]===1,1,-Sum[Times@@mu/@t,{t,Select[tris,And[Length[#]>1,Sort[Join@@#,Greater]===y]&]}]];
    mu/@ptns

Formula

mu(y) = Sum_{g(t)=y} (-1)^d(t), where the sum is over all enriched p-trees (A289501, A299203) whose multiset of leaves is the integer partition y, and d(t) is the number of non-leaf nodes in t.

A299200 Number of twice-partitions whose domain is the integer partition with Heinz number n.

Original entry on oeis.org

1, 1, 2, 1, 3, 2, 5, 1, 4, 3, 7, 2, 11, 5, 6, 1, 15, 4, 22, 3, 10, 7, 30, 2, 9, 11, 8, 5, 42, 6, 56, 1, 14, 15, 15, 4, 77, 22, 22, 3, 101, 10, 135, 7, 12, 30, 176, 2, 25, 9, 30, 11, 231, 8, 21, 5, 44, 42, 297, 6, 385, 56, 20, 1, 33, 14, 490, 15, 60, 15, 627, 4
Offset: 1

Views

Author

Gus Wiseman, Feb 05 2018

Keywords

Comments

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Examples

			The a(15) = 6 twice-partitions: (3)(2), (3)(11), (21)(2), (21)(11), (111)(2), (111)(11).
		

Crossrefs

Programs

  • Maple
    with(numtheory): with(combinat):
    a:= n-> mul(numbpart(pi(i[1]))^i[2], i=ifactors(n)[2]):
    seq(a(n), n=1..82);  # Alois P. Heinz, Jan 14 2021
  • Mathematica
    Table[Times@@Cases[FactorInteger[n],{p_,k_}:>PartitionsP[PrimePi[p]]^k],{n,100}]
  • PARI
    a(n) = {my(f = factor(n)); for (k=1, #f~, f[k, 1] = numbpart(primepi(f[k, 1]));); factorback(f);} \\ Michel Marcus, Feb 26 2018

Formula

Multiplicative with a(prime(n)) = A000041(n).

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

A357982 Replace prime(k) with A000009(k) in the prime factorization of n.

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, 4, 4, 1, 2, 10, 2, 12, 1, 3, 5, 4, 1, 15, 6, 4, 2, 18, 2, 22, 3, 2, 8, 27, 1, 4, 4, 5, 4, 32, 1, 6, 2, 6, 10, 38, 2, 46, 12, 2, 1, 8, 3, 54, 5, 8, 4, 64, 1, 76, 15, 4, 6, 6, 4, 89, 2, 1
Offset: 1

Views

Author

Gus Wiseman, Oct 25 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. This sequence gives the number of ways to choose a strict partition of each prime index of n.
The indices i, where a(i) = 1, form A003586, and the indices j, where a(j) > 1, form A059485. - Ivan N. Ianakiev, Oct 27 2022

Examples

			The a(121) = 9 twice-partitions are: (5)(5), (5)(41), (5)(32), (41)(5), (41)(41), (41)(32), (32)(5), (32)(41), (32)(32).
		

Crossrefs

Other multiplicative sequences: A003961, A357852, A064988, A064989, A357980.
The non-strict version is A299200.
A horizontal version is A357978, non-strict A357977.
A000040 lists the primes.
A056239 adds up prime indices, row-sums of A112798.

Programs

  • Mathematica
    Table[Times@@Cases[FactorInteger[n],{p_,k_}:>PartitionsQ[PrimePi[p]]^k],{n,100}]
  • PARI
    f9(n) = polcoeff( prod( k=1, n, 1 + x^k, 1 + x * O(x^n)), n); \\ A000009
    a(n) = my(f=factor(n)); for (k=1, #f~, f[k,1] = f9(primepi(f[k,1]))); factorback(f); \\ Michel Marcus, Oct 26 2022

A317142 Number of refinement-ordered pairs of strict integer partitions of n.

Original entry on oeis.org

1, 1, 1, 3, 3, 5, 9, 12, 16, 24, 37, 47, 68, 90, 123, 180, 228, 307, 408, 540, 694, 970, 1207, 1598, 2048, 2669, 3357, 4382, 5599, 7109, 8990, 11428, 14330, 18144, 22652, 28343, 35746, 44269, 55094, 68384, 84780, 104477, 129360, 158682, 195323, 240177, 293704
Offset: 0

Views

Author

Gus Wiseman, Jul 22 2018

Keywords

Comments

If x and y are strict partitions of the same integer and it is possible to produce x by further partitioning the parts of y, flattening, and sorting, then x <= y.
This sequence is dominated by A294617 (set partitions of strict partitions).

Examples

			The a(9) = 24 refinement-ordered pairs:
    (9)<=(9)
  (5,4)<=(9)   (5,4)<=(5,4)
  (6,3)<=(9)   (6,3)<=(6,3)
  (7,2)<=(9)   (7,2)<=(7,2)
  (8,1)<=(9)   (8,1)<=(8,1)
(4,3,2)<=(9) (4,3,2)<=(5,4) (4,3,2)<=(6,3) (4,3,2)<=(7,2) (4,3,2)<=(4,3,2)
(5,3,1)<=(9) (5,3,1)<=(5,4) (5,3,1)<=(6,3) (5,3,1)<=(8,1) (5,3,1)<=(5,3,1)
(6,2,1)<=(9) (6,2,1)<=(6,3) (6,2,1)<=(7,2) (6,2,1)<=(8,1) (6,2,1)<=(6,2,1)
		

Crossrefs

Programs

  • Mathematica
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    mps[set_]:=Union[Sort[Sort/@(#/.x_Integer:>set[[x]])]&/@sps[Range[Length[set]]]];
    Table[Sum[Length[Union[Select[Sort/@Map[Total,mps[ptn],{2}],UnsameQ@@#&]]],{ptn,Select[IntegerPartitions[n],UnsameQ@@#&]}],{n,30}]

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

A299203 Number of enriched p-trees whose multiset of leaves is the integer partition with Heinz number n.

Original entry on oeis.org

0, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 4, 1, 1, 1, 5, 1, 3, 1, 3, 1, 1, 1, 11, 1, 1, 2, 3, 1, 5, 1, 12, 1, 1, 1, 15, 1, 1, 1, 11, 1, 4, 1, 3, 3, 1, 1, 38, 1, 3, 1, 3, 1, 9, 1, 9, 1, 1, 1, 21, 1, 1, 4, 34, 1, 4, 1, 3, 1, 5, 1, 54, 1, 1, 3, 3, 1, 4, 1, 33, 5, 1, 1, 23, 1, 1, 1, 9, 1, 20, 1, 3, 1, 1, 1, 117, 1, 3, 3, 12, 1, 4, 1, 9, 4, 1, 1, 57, 1, 4, 1, 34
Offset: 1

Views

Author

Gus Wiseman, Feb 05 2018

Keywords

Comments

By convention, a(1) = 0.
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Examples

			a(54) = 9: (((22)2)1), ((222)1), (((22)1)2), (((21)2)2), ((221)2), ((22)(21)), ((22)21), ((21)22), (2221).
a(40) = 11: ((31)(11)), (((31)1)1), ((3(11))1), ((311)1), (3((11)1)), (3(111)), (((11)1)3), ((111)3), ((31)11), (3(11)1), (3111).
a(36) = 15: ((22)(11)), ((2(11))2), (((11)2)2), (((21)1)2), ((211)2), (((22)1)1), (((21)2)1), ((221)1), ((21)(21)), (22(11)), (2(11)2), ((11)22), ((22)11), ((21)21), (2211).
		

Crossrefs

Programs

  • Mathematica
    nn=120;
    ptns=Table[If[n===1,{},Join@@Cases[FactorInteger[n]//Reverse,{p_,k_}:>Table[PrimePi[p],{k}]]],{n,nn}];
    tris=Join@@Map[Tuples[IntegerPartitions/@#]&,ptns];
    qci[y_]:=qci[y]=If[Length[y]===1,1,Sum[Times@@qci/@t,{t,Select[tris,And[Length[#]>1,Sort[Join@@#,Greater]===y]&]}]];
    qci/@ptns

A061260 G.f.: Product_{k>=1} (1-y*x^k)^(-numbpart(k)), where numbpart(k) = number of partitions of k, cf. A000041.

Original entry on oeis.org

1, 2, 1, 3, 2, 1, 5, 6, 2, 1, 7, 11, 6, 2, 1, 11, 23, 15, 6, 2, 1, 15, 40, 32, 15, 6, 2, 1, 22, 73, 67, 37, 15, 6, 2, 1, 30, 120, 134, 79, 37, 15, 6, 2, 1, 42, 202, 255, 172, 85, 37, 15, 6, 2, 1, 56, 320, 470, 348, 187, 85, 37, 15, 6, 2, 1, 77, 511, 848, 697, 397, 194, 85, 37, 15, 6, 2, 1
Offset: 1

Views

Author

Vladeta Jovovic, Apr 23 2001

Keywords

Comments

Multiset transformation of A000041. - R. J. Mathar, Apr 30 2017
Number of orderless twice-partitions of n of length k. A twice-partition of n is a choice of a partition of each part in a partition of n. The T(5,3) = 6 orderless twice-partitions: (3)(1)(1), (21)(1)(1), (111)(1)(1), (2)(2)(1), (2)(11)(1), (11)(11)(1). - Gus Wiseman, Mar 23 2018

Examples

			:  1;
:  2,   1;
:  3,   2,   1;
:  5,   6,   2,   1;
:  7,  11,   6,   2,  1;
: 11,  23,  15,   6,  2,  1;
: 15,  40,  32,  15,  6,  2,  1;
: 22,  73,  67,  37, 15,  6,  2, 1;
: 30, 120, 134,  79, 37, 15,  6, 2, 1;
: 42, 202, 255, 172, 85, 37, 15, 6, 2, 1;
		

Crossrefs

Row sums: A001970, first column: A000041.
T(2,n) gives A061261,

Programs

  • Maple
    b:= proc(n, i, p) option remember; `if`(p>n, 0, `if`(n=0, 1,
          `if`(min(i, p)<1, 0, add(b(n-i*j, i-1, p-j)*binomial(
           combinat[numbpart](i)+j-1, j), j=0..min(n/i, p)))))
        end:
    T:= (n, k)-> b(n$2, k):
    seq(seq(T(n, k), k=1..n), n=1..14);  # Alois P. Heinz, Apr 13 2017
  • Mathematica
    b[n_, i_, p_] := b[n, i, p] = If[p > n, 0, If[n == 0, 1, If[Min[i, p] < 1, 0, Sum[b[n - i*j, i - 1, p - j]*Binomial[PartitionsP[i] + j - 1, j], {j, 0, Min[n/i, p]}]]]];
    T[n_, k_] := b[n, n, k];
    Table[T[n, k], {n, 1, 14}, {k, 1, n}] // Flatten (* Jean-François Alcover, May 17 2018, after Alois P. Heinz *)
Showing 1-10 of 32 results. Next