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

A353832 Heinz number of the multiset of run-sums of the prime indices of n.

Original entry on oeis.org

1, 2, 3, 3, 5, 6, 7, 5, 7, 10, 11, 9, 13, 14, 15, 7, 17, 14, 19, 15, 21, 22, 23, 15, 13, 26, 13, 21, 29, 30, 31, 11, 33, 34, 35, 21, 37, 38, 39, 25, 41, 42, 43, 33, 35, 46, 47, 21, 19, 26, 51, 39, 53, 26, 55, 35, 57, 58, 59, 45, 61, 62, 49, 13, 65, 66, 67, 51, 69, 70, 71, 35, 73, 74, 39, 57, 77, 78, 79, 35, 19
Offset: 1

Views

Author

Gus Wiseman, May 23 2022

Keywords

Comments

The sequence of runs of a sequence consists of its maximal consecutive constant subsequences when read left-to-right. For example, the runs of (2,2,1,1,1,3,2,2) are (2,2), (1,1,1), (3), (2,2), with sums (4,3,3,4).
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.
This sequence represents the transformation f(P) described by Kimberling at A237685.

Examples

			The prime indices of 1260 are {1,1,2,2,3,4}, with run-sums (2,4,3,4), and the multiset {2,3,4,4} has Heinz number 735, so a(1260) = 735.
		

Crossrefs

The number of distinct prime factors of a(n) is A353835, weak A353861.
The version for compositions is A353847, listed A353932.
The greatest prime factor of a(n) has index A353862, least A353931.
A001222 counts prime factors, distinct A001221.
A056239 adds up prime indices, row sums of A112798 and A296150.
A300273 ranks collapsible partitions, counted by A275870.
A353833 ranks partitions with all equal run-sums, counted by A304442.
A353838 ranks partitions with all distinct run-sums, counted by A353837.
A353840-A353846 pertain to partition run-sum trajectory.
A353851 counts compositions w/ all equal run-sums, ranked by A353848.
A353864 counts rucksack partitions, ranked by A353866.
A353865 counts perfect rucksack partitions, ranked by A353867.
Cf. A005811, A047966, A071625, A073093, A181819, A182850, A182857, A304660, A323014, A353834, A353839, A353841 (1 + iterations needed to reach a squarefree number).

Programs

  • Mathematica
    Table[Times@@Prime/@Cases[If[n==1,{},FactorInteger[n]],{p_,k_}:>PrimePi[p]*k],{n,100}]
  • PARI
    pis_to_runs(n) = { my(runs=List([]), f=factor(n)); for(i=1,#f~,while(f[i,2], listput(runs,primepi(f[i,1])); f[i,2]--)); (runs); };
    A353832(n) = if(1==n,n,my(pruns = pis_to_runs(n), m=1, runsum=pruns[1]); for(i=2,#pruns,if(pruns[i] == pruns[i-1], runsum += pruns[i], m *= prime(runsum); runsum = pruns[i])); (m*prime(runsum))); \\ Antti Karttunen, Jan 20 2025

Formula

A001222(a(n)) = A001221(n).
A001221(a(n)) = A353835(n).
A061395(a(n)) = A353862(n).

Extensions

More terms from Antti Karttunen, Jan 20 2025

A239312 Number of condensed integer partitions of n.

Original entry on oeis.org

1, 1, 1, 2, 3, 3, 5, 6, 9, 10, 14, 16, 23, 27, 33, 41, 51, 62, 75, 93, 111, 134, 159, 189, 226, 271, 317, 376, 445, 520, 609, 714, 832, 972, 1129, 1304, 1520, 1753, 2023, 2326, 2692, 3077, 3540, 4050, 4642, 5298, 6054, 6887, 7854, 8926, 10133, 11501, 13044
Offset: 0

Views

Author

Clark Kimberling, Mar 15 2014

Keywords

Comments

Suppose that p is a partition of n. Let x(1), x(2), ..., x(k) be the distinct parts of p, and let m(i) be the multiplicity of x(i) in p. Let c(p) be the partition {m(1)*x(1), m(2)*x(2), ..., x(k)*m(k)} of n. Call a partition q of n a condensed partition of n if q = c(p) for some partition p of n. Then a(n) is the number of distinct condensed partitions of n. Note that c(p) = p if and only if p has distinct parts and that condensed partitions can have repeated parts.
Also the number of integer partitions of n such that it is possible to choose a different divisor of each part. For example, the partition (6,4,4,1) has choices (3,2,4,1), (3,4,2,1), (6,2,4,1), (6,4,2,1) so is counted under a(15). - Gus Wiseman, Mar 12 2024

Examples

			a(5) = 3 gives the number of partitions of 5 that result from condensations as shown here: 5 -> 5, 41 -> 41, 32 -> 32, 311 -> 32, 221 -> 41, 2111 -> 32, 11111 -> 5.
From _Gus Wiseman_, Mar 12 2024: (Start)
The a(1) = 1 through a(9) = 10 condensed partitions:
  (1)  (2)  (3)    (4)    (5)    (6)      (7)      (8)      (9)
            (2,1)  (2,2)  (3,2)  (3,3)    (4,3)    (4,4)    (5,4)
                   (3,1)  (4,1)  (4,2)    (5,2)    (5,3)    (6,3)
                                 (5,1)    (6,1)    (6,2)    (7,2)
                                 (3,2,1)  (3,2,2)  (7,1)    (8,1)
                                          (4,2,1)  (3,3,2)  (4,3,2)
                                                   (4,2,2)  (4,4,1)
                                                   (4,3,1)  (5,2,2)
                                                   (5,2,1)  (5,3,1)
                                                            (6,2,1)
(End)
		

Crossrefs

The strict case is A000009.
These partitions have ranks A368110, complement A355740.
The complement is counted by A370320.
The version for prime factors (not all divisors) is A370592, ranks A368100.
The complement for prime factors is A370593, ranks A355529.
For a unique choice we have A370595, ranks A370810.
For multiple choices we have A370803, ranks A370811.
The case without ones is A370805, complement A370804.
The version for factorizations is A370814, complement A370813.
A000005 counts divisors.
A000041 counts integer partitions.
A237685 counts partitions of depth 1, or A353837 if we include depth 0.
A355731 counts choices of a divisor of each prime index, firsts A355732.

Programs

  • Maple
    b:= proc(n,i) option remember; `if`(n=0, {[]},
          `if`(i=1, {[n]}, {seq(map(x-> `if`(j=0, x,
           sort([x[], i*j])), b(n-i*j, i-1))[], j=0..n/i)}))
        end:
    a:= n-> nops(b(n$2)):
    seq(a(n), n=0..50);  # Alois P. Heinz, Jul 01 2019
  • Mathematica
    u[n_, k_] := u[n, k] = Map[Total, Split[IntegerPartitions[n][[k]]]]; t[n_] := t[n] = DeleteDuplicates[Table[Sort[u[n, k]], {k, 1, PartitionsP[n]}]]; Table[Length[t[n]], {n, 0,   30}]
    Table[Length[Select[IntegerPartitions[n], Length[Select[Tuples[Divisors/@#],UnsameQ@@#&]]>0&]], {n,0,30}] (* Gus Wiseman, Mar 12 2024 *)

Extensions

Typo in definition corrected by Manfred Scheucher, May 29 2015
Name edited by Gus Wiseman, Mar 13 2024

A353846 Triangle read by rows where T(n,k) is the number of integer partitions of n with partition run-sum trajectory of length k.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 2, 1, 0, 0, 2, 2, 1, 0, 0, 3, 4, 0, 0, 0, 0, 4, 6, 1, 0, 0, 0, 0, 5, 9, 1, 0, 0, 0, 0, 0, 6, 11, 4, 1, 0, 0, 0, 0, 0, 8, 20, 2, 0, 0, 0, 0, 0, 0, 0, 10, 25, 7, 0, 0, 0, 0, 0, 0, 0, 0, 12, 37, 6, 1, 0, 0, 0, 0, 0, 0, 0
Offset: 0

Views

Author

Gus Wiseman, May 26 2022

Keywords

Comments

Every sequence can be uniquely split into a sequence of non-overlapping runs. For example, the runs of (2,2,1,1,1,3,2,2) are ((2,2),(1,1,1),(3),(2,2)), with sums (4,3,3,4). The run-sum trajectory is obtained by repeatedly taking run-sums (or condensations) until a strict partition is reached. For example, the trajectory of (2,1,1) is (2,1,1) -> (2,2) -> (4).
Also the number of integer partitions of n with Kimberling's depth statistic (see A237685, A237750) equal to k-1.

Examples

			Triangle begins:
   1
   0   1
   0   1   1
   0   2   1   0
   0   2   2   1   0
   0   3   4   0   0   0
   0   4   6   1   0   0   0
   0   5   9   1   0   0   0   0
   0   6  11   4   1   0   0   0   0
   0   8  20   2   0   0   0   0   0   0
   0  10  25   7   0   0   0   0   0   0   0
   0  12  37   6   1   0   0   0   0   0   0   0
   0  15  47  13   2   0   0   0   0   0   0   0   0
   0  18  67  15   1   0   0   0   0   0   0   0   0   0
   0  22  85  25   3   0   0   0   0   0   0   0   0   0   0
   0  27 122  26   1   0   0   0   0   0   0   0   0   0   0   0
For example, row n = 8 counts the following partitions (empty columns indicated by dots):
.  (8)    (44)        (422)     (4211)  .  .  .  .
   (53)   (332)       (32111)
   (62)   (611)       (41111)
   (71)   (2222)      (221111)
   (431)  (3221)
   (521)  (3311)
          (5111)
          (22211)
          (311111)
          (2111111)
          (11111111)
		

Crossrefs

Row-sums are A000041.
Column k = 1 is A000009.
Column k = 2 is A237685.
Column k = 3 is A237750.
The version for run-lengths instead of run-sums is A225485 or A325280.
This statistic (trajectory length) is ranked by A353841 and A326371.
The version for compositions is A353859, see also A353847-A353858.
A005811 counts runs in binary expansion.
A275870 counts collapsible partitions, ranked by A300273.
A304442 counts partitions with all equal run-sums, ranked by A353833.
A353832 represents the operation of taking run-sums of a partition
A353836 counts partitions by number of distinct run-sums.
A353838 ranks partitions with all distinct run-sums, counted by A353837.
A353840-A353846 pertain to partition run-sum trajectory.
A353845 counts partitions whose run-sum trajectory ends in a singleton.

Programs

  • Mathematica
    rsn[y_]:=If[y=={},{},NestWhileList[Reverse[Sort[Total/@ Split[Sort[#]]]]&,y,!UnsameQ@@#&]];
    Table[Length[Select[IntegerPartitions[n],Length[rsn[#]]==k&]],{n,0,15},{k,0,n}]

A353853 Trajectory of the composition run-sum transformation (or condensation) of n, using standard composition numbers.

Original entry on oeis.org

0, 1, 2, 3, 2, 4, 5, 6, 7, 4, 8, 9, 10, 8, 11, 10, 8, 12, 13, 14, 10, 8, 15, 8, 16, 17, 18, 19, 18, 20, 21, 17, 22, 23, 20, 24, 25, 26, 24, 27, 26, 24, 28, 20, 29, 21, 17, 30, 18, 31, 16, 32, 33, 34, 35, 34, 36, 32, 37, 38, 39, 36, 32, 40, 41, 42, 32
Offset: 0

Views

Author

Gus Wiseman, Jun 01 2022

Keywords

Comments

Every sequence can be uniquely split into a sequence of non-overlapping runs. For example, the runs of (2,2,1,1,1,3,2,2) are ((2,2),(1,1,1),(3),(2,2)), with sums (4,3,3,4).
The k-th composition in standard order (graded reverse-lexicographic, A066099) is obtained by taking the set of positions of 1's in the reversed binary expansion of k, prepending 0, taking first differences, and reversing again. This gives a bijective correspondence between nonnegative integers and integer compositions.
The run-sum trajectory is obtained by repeatedly taking the run-sum transformation (A353847) until the rank of an anti-run is reached. For example, the trajectory 11 -> 10 -> 8 given in row 11 corresponds to the trajectory (2,1,1) -> (2,2) -> (4).

Examples

			Triangle begins:
   0
   1
   2
   3  2
   4
   5
   6
   7  4
   8
   9
  10  8
  11 10  8
  12
  13
  14 10  8
For example, the trajectory of 29 is 29 -> 21 -> 17, corresponding to the compositions (1,1,2,1) -> (2,2,1) -> (4,1).
		

Crossrefs

These sequences for partitions are A353840-A353846.
This is the iteration of A353847, with partition version A353832.
Row-lengths are A353854, counted by A353859.
Final terms are A353855.
Counting rows by weight of final term gives A353856.
Rows ending in a power of 2 are A353857, counted by A353858.
A003242 counts anti-run compositions, ranked by A333489, complement A261983.
A005811 counts runs in binary expansion.
A011782 counts compositions.
A066099 lists compositions in standard order.
A318928 gives runs-resistance of binary expansion.
A329739 counts compositions with all distinct run-lengths.
A333627 ranks the run-lengths of standard compositions.
A351014 counts distinct runs in standard compositions, firsts A351015.
A353853-A353859 pertain to composition run-sum trajectory.
A353929 counts distinct runs in binary expansion, firsts A353930.
A353932 lists run-sums of standard compositions.

Programs

  • Mathematica
    stc[n_]:=Differences[Prepend[Join@@ Position[Reverse[IntegerDigits[n,2]],1],0]]//Reverse;
    stcinv[q_]:=Total[2^(Accumulate[Reverse[q]])]/2;
    Table[NestWhileList[stcinv[Total/@Split[stc[#]]]&,n,MatchQ[stc[#],{_,x_,x_,_}]&],{n,0,50}]

A353840 Trajectory of the partition run-sum transformation of n, using Heinz numbers.

Original entry on oeis.org

1, 2, 3, 4, 3, 5, 6, 7, 8, 5, 9, 7, 10, 11, 12, 9, 7, 13, 14, 15, 16, 7, 17, 18, 14, 19, 20, 15, 21, 22, 23, 24, 15, 25, 13, 26, 27, 13, 28, 21, 29, 30, 31, 32, 11, 33, 34, 35, 36, 21, 37, 38, 39, 40, 25, 13, 41, 42, 43, 44, 33, 45, 35, 46, 47, 48, 21, 49, 19
Offset: 1

Views

Author

Gus Wiseman, May 25 2022

Keywords

Comments

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 run-sum trajectory is obtained by repeatedly taking the run-sum transformation (A353832) until a squarefree number is reached. For example, the trajectory 12 -> 9 -> 7 given in row 12 corresponds to the partitions (2,1,1) -> (2,2) -> (4).
This is the iteration of the transformation f described by Kimberling at A237685.

Examples

			Triangle begins:
   1
   2
   3
   4  3
   5
   6
   7
   8  5
   9  7
  10
  11
  12  9  7
Row 87780 is the following trajectory (left column), with prime indices shown on the right:
  87780: {1,1,2,3,4,5,8}
  65835: {2,2,3,4,5,8}
  51205: {3,4,4,5,8}
  19855: {3,5,8,8}
   2915: {3,5,16}
		

Crossrefs

The version for run-lengths instead of sums is A325239 or A325277.
This is the iteration of A353832, with composition version A353847.
Row-lengths are A353841, counted by A353846.
Final terms are A353842.
Counting rows by final omega gives A353843.
Rows ending in a prime number are A353844, counted by A353845.
These sequences for compositions are A353853-A353859.
A001222 counts prime factors, distinct A001221.
A056239 adds up prime indices, row sums of A112798 and A296150.
A124010 gives prime signature, sorted A118914.
A182850 or A323014 gives frequency depth.
A300273 ranks collapsible partitions, counted by A275870.
A353833 ranks partitions with all equal run-sums, counted by A304442.
A353835 counts distinct run-sums of prime indices, weak A353861.
A353838 ranks partitions with all distinct run-sums, counted by A353837.
A353862 gives greatest run-sum of prime indices, least A353931.

Programs

  • Mathematica
    Table[NestWhileList[Times@@Prime/@Cases[FactorInteger[#],{p_,k_}:>PrimePi[p]*k]&,n,Not@*SquareFreeQ],{n,30}]

A353859 Triangle read by rows where T(n,k) is the number of integer compositions of n with composition run-sum trajectory of length k.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 3, 1, 0, 0, 4, 2, 2, 0, 0, 7, 7, 2, 0, 0, 0, 14, 14, 4, 0, 0, 0, 0, 23, 29, 12, 0, 0, 0, 0, 0, 39, 56, 25, 8, 0, 0, 0, 0, 0, 71, 122, 53, 10, 0, 0, 0, 0, 0, 0, 124, 246, 126, 16, 0, 0, 0, 0, 0, 0, 0, 214, 498, 264, 48, 0, 0, 0, 0, 0, 0, 0
Offset: 0

Views

Author

Gus Wiseman, Jun 02 2022

Keywords

Comments

Every sequence can be uniquely split into a sequence of non-overlapping runs. For example, the runs of (2,2,1,1,1,3,2,2) are ((2,2),(1,1,1),(3),(2,2)), with sums (4,3,3,4). The run-sum trajectory is obtained by repeatedly taking the run-sums transformation (or condensation, represented by A353847) until an anti-run is reached. For example, the trajectory (2,4,2,1,1) -> (2,4,2,2) -> (2,4,4) -> (2,8) is counted under T(10,4).

Examples

			Triangle begins:
   1
   0   1
   0   1   1
   0   3   1   0
   0   4   2   2   0
   0   7   7   2   0   0
   0  14  14   4   0   0   0
   0  23  29  12   0   0   0   0
   0  39  56  25   8   0   0   0   0
   0  71 122  53  10   0   0   0   0   0
   0 124 246 126  16   0   0   0   0   0   0
   0 214 498 264  48   0   0   0   0   0   0   0
For example, row n = 5 counts the following compositions:
  (5)    (113)    (1121)
  (14)   (122)    (1211)
  (23)   (221)
  (32)   (311)
  (41)   (1112)
  (131)  (2111)
  (212)  (11111)
		

Crossrefs

Column k = 1 is A003242, ranked by A333489, complement A261983.
Row sums are A011782.
Positive row-lengths are A070939.
The version for partitions is A353846, ranked by A353841.
This statistic (trajectory length) is ranked by A353854, firsts A072639.
Counting by length of last part instead of number of parts gives A353856.
A333627 ranks the run-lengths of standard compositions.
A353847 represents the run-sums of a composition, partitions A353832.
A353853-A353859 pertain to composition run-sum trajectory.
A353932 lists run-sums of standard compositions.

Programs

  • Mathematica
    rsc[y_]:=If[y=={},{},NestWhileList[Total/@Split[#]&,y,MatchQ[#,{_,x_,x_,_}]&]];
    Table[Length[Select[Join@@Permutations/@IntegerPartitions[n],Length[rsc[#]]==k&]],{n,0,10},{k,0,n}]

A353860 Number of collapsible integer compositions of n.

Original entry on oeis.org

0, 1, 2, 2, 5, 2, 12, 2, 26, 9, 36, 2, 206, 2, 132, 40, 677, 2, 1746, 2, 3398, 136, 2052, 2, 44388, 33, 8196, 730, 79166, 2, 263234, 2, 458330, 2056, 131076, 160, 8804349, 2, 524292, 8200, 13662156, 2, 36036674, 2, 48844526, 90282, 8388612, 2, 1971667502, 129
Offset: 0

Views

Author

Gus Wiseman, Jun 04 2022

Keywords

Comments

If a collapse is a joining of some number of adjacent equal parts of an integer composition, we call a composition collapsible iff by some sequence of collapses it can be reduced to a single part. An example of such a sequence of collapses is (1,1,1,3,2,1,1,2) -> (3,3,2,1,1,2) -> (3,3,2,2,2) -> (6,2,2,2) -> (6,6) -> (12), which shows that (1,1,1,3,2,1,1,2) is a collapsible composition of 12.

Examples

			The a(0) = 0 through a(6) = 12 compositions:
  .  (1)  (2)   (3)    (4)     (5)      (6)
          (11)  (111)  (22)    (11111)  (33)
                       (112)            (222)
                       (211)            (1113)
                       (1111)           (1122)
                                        (2112)
                                        (2211)
                                        (3111)
                                        (11112)
                                        (11211)
                                        (21111)
                                        (111111)
		

Crossrefs

The version for partitions is A275870, ranked by A300273.
A003242 counts anti-run compositions, ranked by A333489, complement A261983.
A011782 counts compositions.
A353847 represents the run-sums of a composition, partitions A353832.
A353853-A353859 pertain to composition run-sum trajectory.
A353932 lists run-sums of standard compositions.

Programs

  • Mathematica
    repcams[q_List]:=repcams[q]=Union[{q},If[UnsameQ@@q,{},Union@@repcams/@ Union[Insert[Drop[q,#],Plus@@Take[q,#],First[#]]&/@ Select[Tuples[Range[Length[q]],2],And[Less@@#,SameQ@@Take[q,#]]&]]]];
    Table[Length[Select[Join@@Permutations/@ IntegerPartitions[n],MemberQ[repcams[#],{n}]&]],{n,0,15}]
  • PARI
    a(n) = if(n==0, 0, 1 - sumdiv(n, d, if(d>1, moebius(d)*a(n/d)^d ))) \\ Andrew Howroyd, Feb 04 2023

Formula

Sum_{d|n} mu(d)*a(n/d)^d = 1 for n > 0. - Andrew Howroyd, Feb 04 2023

Extensions

Terms a(16) and beyond from Andrew Howroyd, Feb 04 2023

A353863 Number of integer partitions of n whose weak run-sums cover an initial interval of nonnegative integers.

Original entry on oeis.org

1, 1, 1, 2, 2, 3, 4, 6, 7, 10, 11, 16, 20, 24, 30, 43, 47, 62, 79, 94, 113, 143, 170, 211, 256, 307, 372, 449, 531, 648, 779, 926, 1100, 1323, 1562, 1864, 2190, 2595, 3053, 3611, 4242, 4977, 5834, 6825, 7973, 9344, 10844, 12641, 14699, 17072, 19822
Offset: 0

Views

Author

Gus Wiseman, Jun 04 2022

Keywords

Comments

A weak run-sum of a sequence is the sum of any consecutive constant subsequence. For example, the weak run-sums of (3,2,2,1) are {1,2,3,4}.
This is a kind of completeness property, cf. A126796.

Examples

			The a(1) = 1 through a(8) = 7 partitions:
  (1)  (11)  (21)   (211)   (311)    (321)     (3211)     (3221)
             (111)  (1111)  (2111)   (3111)    (4111)     (32111)
                            (11111)  (21111)   (22111)    (41111)
                                     (111111)  (31111)    (221111)
                                               (211111)   (311111)
                                               (1111111)  (2111111)
                                                          (11111111)
		

Crossrefs

For parts instead of weak run-sums we have A000009.
For multiplicities instead of weak run-sums we have A317081.
If weak run-sums are distinct we have A353865, the completion of A353864.
A003242 counts anti-run compositions, ranked by A333489, complement A261983.
A005811 counts runs in binary expansion.
A165413 counts distinct run-lengths in binary expansion, sums A353929.
A300273 ranks collapsible partitions, counted by A275870, comps A353860.
A353832 represents taking run-sums of a partition, compositions A353847.
A353833 ranks partitions with all equal run-sums, counted by A304442.
A353835 counts distinct run-sums of prime indices.
A353837 counts partitions with distinct run-sums, ranked by A353838.
A353840-A353846 pertain to partition run-sum trajectory.
A353861 counts distinct weak run-sums of prime indices.
A353932 lists run-sums of standard compositions.

Programs

  • Mathematica
    normQ[m_]:=m=={}||Union[m]==Range[Max[m]];
    msubs[s_]:=Join@@@Tuples[Table[Take[t,i],{t,Split[s]},{i,0,Length[t]}]];
    wkrs[y_]:=Union[Total/@Select[msubs[y],SameQ@@#&]];
    Table[Length[Select[IntegerPartitions[n],normQ[Rest[wkrs[#]]]&]],{n,0,15}]
  • PARI
    \\ isok(p) tests the partition.
    isok(p)={my(b=0, s=0, t=0); for(i=1, #p, if(p[i]<>t, t=p[i]; s=0); s += t; b = bitor(b, 1<<(s-1))); bitand(b,b+1)==0}
    a(n) = {my(r=0); forpart(p=n, r+=isok(p)); r} \\ Andrew Howroyd, Jan 15 2024

Extensions

a(31) onwards from Andrew Howroyd, Jan 15 2024

A353841 Length of the trajectory of the partition run-sum transformation of n, using Heinz numbers; a(1) = 0.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, May 25 2022

Keywords

Comments

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.
Starting with n, this is one plus the number of times one must apply A353832 to reach a squarefree number.
Also Kimberling's depth statistic (defined in A237685 and A237750) plus one.

Examples

			The trajectory for a(1080) = 4 is the following, with prime indices shown on the right:
  1080: {1,1,1,2,2,2,3}
   325: {3,3,6}
   169: {6,6}
    37: {12}
The trajectory for a(87780) = 5 is the following, with prime indices shown on the right:
  87780: {1,1,2,3,4,5,8}
  65835: {2,2,3,4,5,8}
  51205: {3,4,4,5,8}
  19855: {3,5,8,8}
   2915: {3,5,16}
The trajectory for a(39960) = 5 is the following, with prime indices shown on the right:
  39960: {1,1,1,2,2,2,3,12}
  12025: {3,3,6,12}
   6253: {6,6,12}
   1369: {12,12}
     89: {24}
		

Crossrefs

Positions of 1's are A005117.
The version for run-lengths instead of sums is A182850 or A323014.
Positions of first appearances are A353743.
These are the row-lengths of A353840.
Other sequences pertaining to this trajectory are A353842-A353845.
Counting partitions by this statistic gives A353846.
The version for compositions is A353854, run-lengths of A353853.
A001222 counts prime factors, distinct A001221.
A005811 counts runs in binary expansion.
A056239 adds up prime indices, row sums of A112798 and A296150.
A300273 ranks collapsible partitions, counted by A275870.
A318928 gives runs-resistance of binary expansion.
A353832 represents the operation of taking run-sums of a partition.
A353833 ranks partitions with all equal run-sums, counted by A304442.
A353835 counts distinct run-sums of prime indices, weak A353861.
A353838 ranks partitions with all distinct run-sums, counted by A353837.
A353866 ranks rucksack partitions, counted by A353864.

Programs

  • Mathematica
    Table[If[n==1,0,Length[NestWhileList[Times@@Prime/@Cases[If[#==1,{},FactorInteger[#]],{p_,k_}:>PrimePi[p]*k]&,n,!SquareFreeQ[#]&]]],{n,100}]
  • PARI
    pis_to_runs(n) = { my(runs=List([]), f=factor(n)); for(i=1,#f~,while(f[i,2], listput(runs,primepi(f[i,1])); f[i,2]--)); (runs); };
    A353832(n) = if(1==n,n,my(pruns = pis_to_runs(n), m=1, runsum=pruns[1]); for(i=2,#pruns,if(pruns[i] == pruns[i-1], runsum += pruns[i], m *= prime(runsum); runsum = pruns[i])); (m*prime(runsum)));
    A353841(n) = if(1==n,0,for(i=1,oo,if(issquarefree(n), return(i), n = A353832(n)))); \\ Antti Karttunen, Jan 20 2025

Formula

a(1) = 0, and for n > 1, if A008966(n) = 1 [n is in A005117], a(n) = 1, otherwise a(n) = 1+a(A353832(n)). [See comments] - Antti Karttunen, Jan 20 2025

Extensions

More terms from Antti Karttunen, Jan 20 2025

A237750 Number of partitions of n having depth 2; see Comments.

Original entry on oeis.org

0, 0, 0, 1, 0, 1, 1, 4, 2, 7, 6, 13, 15, 25, 26, 46, 53, 74, 92, 136, 157, 218, 274, 356, 443, 583, 703, 899, 1125, 1447, 1746, 2182, 2661, 3331, 4077, 4997, 6066, 7432, 8984, 10904, 13212, 15845, 19161, 22932, 27526, 32968, 39351, 46778, 55791, 66272, 78480
Offset: 1

Views

Author

Clark Kimberling, Feb 19 2014

Keywords

Comments

Suppose that P is a partition of n. Let x(1), x(2),...,x(k) be the distinct parts of P, and let m(i) be the multiplicity of x(i) in P. Let f(P) be the partition [m(1)*x(1), m(2)*x(2),...,x(k)*m(k)] of n. Define c(0,P) = P, c(1,P) = f(P),..., c(n,P) = f(c(n-1,P)), and define d(P) = least n such that c(n,P) has no repeated parts; d(P) is introduced here as the depth of P. Clearly d(P) = 0 if and only if P is a strict partition, as in A000009. Conjecture: if d >= 0, then 2^d is the least n that has a partition of depth d.

Examples

			The 11 partitions of 6 are partitioned by depth as follows:
  depth 0: 6, 51, 42, 321;
  depth 1: 411, 33, 222, 2211, 21111, 11111;
  depth 2: 3111.
Thus, a(6) = 6, A000009(6) = 4, A237750(6) = 1, A237978(6) = 0.
		

Crossrefs

Programs

  • Mathematica
    z = 60; c[n_] := c[n] = Map[Length[FixedPointList[Sort[Map[Total, Split[#]], Greater] &, #]] - 2 &, IntegerPartitions[n]]
    Table[Count[c[n], 1], {n, 1, z}] (* A237685 *)
    Table[Count[c[n], 2], {n, 1, z}] (* this sequence *)
    Table[Count[c[n], 3], {n, 1, z}] (* A237978 *)
    (* Peter J. C. Moses, Feb 19 2014 *)
Showing 1-10 of 20 results. Next