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

A181819 Prime shadow of n: a(1) = 1; for n>1, if n = Product prime(i)^e(i), then a(n) = Product prime(e(i)).

Original entry on oeis.org

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

Views

Author

Matthew Vandermast, Dec 07 2010

Keywords

Comments

a(n) depends only on prime signature of n (cf. A025487). a(m) = a(n) iff m and n have the same prime signature, i.e., iff A046523(m) = A046523(n).
Because A046523 (the smallest representative of prime signature of n) and this sequence are functions of each other as A046523(n) = A181821(a(n)) and a(n) = a(A046523(n)), it implies that for all i, j: a(i) = a(j) <=> A046523(i) = A046523(j) <=> A101296(i) = A101296(j), i.e., that equivalence-class-wise this is equal to A101296, and furthermore, applying any function f on this sequence gives us a sequence b(n) = f(a(n)) whose equivalence class partitioning is equal to or coarser than that of A101296, i.e., b is then a sequence that depends only on the prime signature of n (the multiset of exponents of its prime factors), although not necessarily in a very intuitive way. - Antti Karttunen, Apr 28 2022

Examples

			20 = 2^2*5 has the exponents (2,1) in its prime factorization. Accordingly, a(20) = prime(2)*prime(1) = A000040(2)*A000040(1) = 3*2 = 6.
		

Crossrefs

Programs

Formula

From Antti Karttunen, Feb 07 2016: (Start)
a(1) = 1; for n > 1, a(n) = A000040(A067029(n)) * a(A028234(n)).
a(1) = 1; for n > 1, a(n) = A008578(A001511(n)) * a(A064989(n)).
Other identities. For all n >= 1:
a(A124859(n)) = A122111(a(n)) = A238745(n). - from Matthew Vandermast's formulas for the latter sequence.
(End)
a(n) = A246029(A156552(n)). - Antti Karttunen, Oct 15 2016
From Antti Karttunen, Apr 28 & Apr 30 2022: (Start)
A181821(a(n)) = A046523(n) and a(A046523(n)) = a(n). [See comments]
a(n) = A329900(A124859(n)) = A319626(A124859(n)).
a(n) = A246029(A156552(n)).
a(a(n)) = A328830(n).
a(A304660(n)) = n.
a(A108951(n)) = A122111(n).
a(A185633(n)) = A322312(n).
a(A025487(n)) = A181820(n).
a(A276076(n)) = A275735(n) and a(A276086(n)) = A328835(n).
As the sequence converts prime exponents to prime indices, it effects the following mappings:
A001221(a(n)) = A071625(n). [Number of distinct indices --> Number of distinct exponents]
A001222(a(n)) = A001221(n). [Number of indices (i.e., the number of prime factors with multiplicity) --> Number of exponents (i.e., the number of distinct prime factors)]
A056239(a(n)) = A001222(n). [Sum of indices --> Sum of exponents]
A066328(a(n)) = A136565(n). [Sum of distinct indices --> Sum of distinct exponents]
A003963(a(n)) = A005361(n). [Product of indices --> Product of exponents]
A290103(a(n)) = A072411(n). [LCM of indices --> LCM of exponents]
A156061(a(n)) = A290107(n). [Product of distinct indices --> Product of distinct exponents]
A257993(a(n)) = A134193(n). [Index of the least prime not dividing n --> The least number not among the exponents]
A055396(a(n)) = A051904(n). [Index of the least prime dividing n --> Minimal exponent]
A061395(a(n)) = A051903(n). [Index of the greatest prime dividing n --> Maximal exponent]
A008966(a(n)) = A351564(n). [All indices are distinct (i.e., n is squarefree) --> All exponents are distinct]
A007814(a(n)) = A056169(n). [Number of occurrences of index 1 (i.e., the 2-adic valuation of n) --> Number of occurrences of exponent 1]
A056169(a(n)) = A136567(n). [Number of unitary prime divisors --> Number of exponents occurring only once]
A064989(a(n)) = a(A003557(n)) = A295879(n). [Indices decremented after <--> Exponents decremented before]
Other mappings:
A007947(a(n)) = a(A328400(n)) = A329601(n).
A181821(A007947(a(n))) = A328400(n).
A064553(a(n)) = A000005(n) and A000005(a(n)) = A182860(n).
A051903(a(n)) = A351946(n).
A003557(a(n)) = A351944(n).
A258851(a(n)) = A353379(n).
A008480(a(n)) = A309004(n).
a(A325501(n)) = A325507(n) and a(A325502(n)) = A038754(n+1).
a(n!) = A325508(n).
(End)

Extensions

Name "Prime shadow" (coined by Gus Wiseman in A325755) prefixed to the definition by Antti Karttunen, Apr 27 2022

A353864 Number of rucksack partitions of n: every consecutive constant subsequence has a different sum.

Original entry on oeis.org

1, 1, 2, 3, 4, 6, 8, 11, 14, 19, 25, 33, 39, 51, 65, 82, 101, 126, 154, 191, 232, 284, 343, 416, 496, 600, 716, 855, 1018, 1209, 1430, 1691, 1991, 2345, 2747, 3224, 3762, 4393, 5116, 5946, 6897, 7998, 9257, 10696, 12336, 14213, 16343, 18781, 21538, 24687, 28253, 32291, 36876, 42057
Offset: 0

Views

Author

Gus Wiseman, May 23 2022

Keywords

Comments

In a knapsack partition (A108917), every submultiset has a different sum, so these are run-knapsack partitions or rucksack partitions for short. Another variation of knapsack partitions is A325862.

Examples

			The a(0) = 1 through a(7) = 11 partitions:
  ()  (1)  (2)   (3)    (4)     (5)      (6)       (7)
           (11)  (21)   (22)    (32)     (33)      (43)
                 (111)  (31)    (41)     (42)      (52)
                        (1111)  (221)    (51)      (61)
                                (311)    (222)     (322)
                                (11111)  (321)     (331)
                                         (411)     (421)
                                         (111111)  (511)
                                                   (2221)
                                                   (4111)
                                                   (1111111)
		

Crossrefs

Knapsack partitions are counted by A108917, ranked by A299702.
The strong case is A353838, counted by A353837, complement A353839.
The perfect case is A353865, ranked by A353867.
These partitions are ranked by A353866.
A000041 counts partitions, strict A000009.
A300273 ranks collapsible partitions, counted by A275870.
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.
A353840-A353846 pertain to partition run-sum trajectory.
A353852 ranks compositions with all distinct run-sums, counted by A353850.
A353863 counts partitions whose weak run-sums cover an initial interval.

Programs

  • Mathematica
    msubs[s_]:=Join@@@Tuples[Table[Take[t,i],{t,Split[s]},{i,0,Length[t]}]];
    Table[Length[Select[IntegerPartitions[n],UnsameQ@@Total/@Select[msubs[#],SameQ@@#&]&]],{n,0,30}]

Extensions

a(50)-a(53) from Robert Price, Apr 03 2025

A225485 Number of partitions of n that have frequency depth k, an array read by rows.

Original entry on oeis.org

0, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 2, 3, 1, 3, 4, 3, 1, 1, 4, 8, 1, 1, 3, 6, 9, 3, 1, 2, 8, 12, 7, 1, 3, 11, 17, 10, 1, 1, 11, 26, 17, 1, 5, 19, 25, 27, 1, 1, 17, 44, 38, 1, 3, 25, 53, 52, 1, 1, 3, 29, 63, 76, 4
Offset: 1

Views

Author

Clark Kimberling, May 08 2013

Keywords

Comments

Let S = {x(1),...,x(k)} be a multiset whose distinct elements are y(1),...,y(h). Let f(i) be the frequency of y(i) in S. Define F(S) = {f(1),..,f(h)}, F(1,S) = F(S), and F(m,S) = F(F(m-1),S) for m>1. Then lim(F(m,S)) = {1} for every S, so that there is a least positive integer i for which F(i,S) = {1}, which we call the frequency depth of S.
Equivalently, the frequency depth of an integer partition is the number of times one must take the multiset of multiplicities to reach (1). For example, the partition (32211) has frequency depth 5 because we have: (32211) -> (221) -> (21) -> (11) -> (2) -> (1). - Gus Wiseman, Apr 19 2019
From Clark Kimberling, Sep 26 2023: (Start)
Below, m^n abbreviates the sum m+...+m of n terms. In the following list, the numbers p_1,...,p_k are distinct, m >= 1, and k >= 1. The forms of the partitions being counted are as follows:
column 1: [n],
column 2: [m^k],
column 3: [p_1^m,...,p_k^m],
column 4: [(p_1^m_1)^m,..., (p_k^m_k)^m], distinct numbers m_i.
Column 3 is of special interest. Assume first that m = 1, so that the form of partition being counted is p = [p_1,...,p_k], with conjugate given by [q_1,...,q_m] where q_i is the number of parts of p that are >= i. Since the p_i are distinct, the distinct parts of q are the integers 1,2,...,k. For the general case that m >= 1, the distinct parts of q are the integers m,...,km. Let S(n) denote the set of partitions of n counted by column 3. Then if a and b are in the set S*(n) of conjugates of partitions in S(n), and if a > b, then a - b is also in S*(n). Call this the subtraction property. Conversely, if a partition q has the subtraction property, then q must consist of a set of numbers m,..,km for some m. Thus, column 3 counts the partitions of n that have the subtraction property. (End)

Examples

			The first 9 rows:
  n = 1 .... 0
  n = 2 .... 1..1
  n = 3 .... 1..1..1
  n = 4 .... 1..2..1..1
  n = 5 .... 1..1..2..3
  n = 6 .... 1..3..4..3
  n = 7 .... 1..1..4..8..1
  n = 8 .... 1..3..6..9..3
  n = 9 .... 1..2..8.12..7
For the 7 partitions of 5, successive frequencies are shown here:
  5 -> 1 (depth 1)
  41 -> 11 -> 2 -> 1 (depth 3)
  32 -> 11 -> 2 -> 1 (depth 3)
  311 -> 12 -> 11 -> 2 -> 1 (depth 4)
  221 -> 12 -> 11 -> 2 -> 1 (depth 4)
  2111 -> 13 -> 11 -> 2 -> 1 (depth 4)
  11111 -> 5 -> 1 (depth 2)
Summary: 1 partition has depth 1; 1 has depth 2; 2 have 3; and 3 have 4, so that the row for n = 5 is 1..1..2..3 .
		

Crossrefs

Row sums are A000041.
Column k = 2 is A032741.
Column k = 3 is A325245.
a(n!) = A325272(n).
Integer partition triangles: A008284 (first omega), A116608 (second omega), A325242 (third omega), A325268 (second-to-last omega), A225485 or A325280 (length/frequency depth).

Programs

  • Mathematica
    c[s_] := c[s] = Select[Table[Count[s, i], {i, 1, Max[s]}], # > 0 &]
    f[s_] := f[s] = Drop[FixedPointList[c, s], -2]
    t[s_] := t[s] = Length[f[s]]
    u[n_] := u[n] = Table[t[Part[IntegerPartitions[n], i]],
      {i, 1, Length[IntegerPartitions[n]]}];
    Flatten[Table[Count[u[n], k], {n, 2, 25}, {k, 1, Max[u[n]]}]]

A325277 Irregular triangle read by rows where row 1 is {1} and row n is the sequence starting with n and repeatedly applying A181819 until a prime number is reached.

Original entry on oeis.org

1, 2, 3, 4, 3, 5, 6, 4, 3, 7, 8, 5, 9, 3, 10, 4, 3, 11, 12, 6, 4, 3, 13, 14, 4, 3, 15, 4, 3, 16, 7, 17, 18, 6, 4, 3, 19, 20, 6, 4, 3, 21, 4, 3, 22, 4, 3, 23, 24, 10, 4, 3, 25, 3, 26, 4, 3, 27, 5, 28, 6, 4, 3, 29, 30, 8, 5, 31, 32, 11, 33, 4, 3
Offset: 1

Views

Author

Gus Wiseman, Apr 15 2019

Keywords

Comments

The function A181819 maps p^i*...*q^j to prime(i)*...*prime(j) where p through q are distinct primes.

Examples

			Triangle begins:
   1            26 4 3        51 4 3          76 6 4 3
   2            27 5          52 6 4 3        77 4 3
   3            28 6 4 3      53              78 8 5
   4 3          29            54 10 4 3       79
   5            30 8 5        55 4 3          80 14 4 3
   6 4 3        31            56 10 4 3       81 7
   7            32 11         57 4 3          82 4 3
   8 5          33 4 3        58 4 3          83
   9 3          34 4 3        59              84 12 6 4 3
  10 4 3        35 4 3        60 12 6 4 3     85 4 3
  11            36 9 3        61              86 4 3
  12 6 4 3      37            62 4 3          87 4 3
  13            38 4 3        63 6 4 3        88 10 4 3
  14 4 3        39 4 3        64 13           89
  15 4 3        40 10 4 3     65 4 3          90 12 6 4 3
  16 7          41            66 8 5          91 4 3
  17            42 8 5        67              92 6 4 3
  18 6 4 3      43            68 6 4 3        93 4 3
  19            44 6 4 3      69 4 3          94 4 3
  20 6 4 3      45 6 4 3      70 8 5          95 4 3
  21 4 3        46 4 3        71              96 22 4 3
  22 4 3        47            72 15 4 3       97
  23            48 14 4 3     73              98 6 4 3
  24 10 4 3     49 3          74 4 3          99 6 4 3
  25 3          50 6 4 3      75 6 4 3       100 9 3
		

Crossrefs

Row lengths are 1 for n = 1 and A323014(n) for n > 1.

Programs

  • Mathematica
    red[n_]:=Times@@Prime/@Last/@If[n==1,{},FactorInteger[n]];
    Table[NestWhileList[red,n,#>1&&!PrimeQ[#]&],{n,30}]

Formula

T(n,k) = A325239(n,k) for k <= A323014(n).
A001222(T(n,k)) = A323023(n,k) for n > 1.

A325238 First positive integer with each omega-sequence.

Original entry on oeis.org

1, 2, 4, 6, 8, 12, 16, 24, 30, 32, 36, 48, 60, 64, 96, 120, 128, 192, 210, 216, 240, 256, 360, 384, 420, 480, 512, 720, 768, 840, 900, 960, 1024, 1260, 1296, 1440, 1536, 1680, 1920, 2048, 2310, 2520, 2880, 3072, 3360, 3840, 4096, 4620, 5040, 5760, 6144, 6720
Offset: 1

Views

Author

Gus Wiseman, Apr 14 2019

Keywords

Comments

We define the omega-sequence of n (row n of A323023) to have length A323014(n) = frequency depth of n, and the k-th part is Omega(red^{k-1}(n)), where Omega = A001222 and red^{k} is the k-th functional iteration of red = A181819, given by red(n = p^i*...*q^j) = prime(i)*...*prime(j), i.e., the product of primes indexed by the prime exponents of n.

Examples

			The sequence of terms together with their omega-sequences begins:
    1:
    2: 1
    4: 2 1
    6: 2 2 1
    8: 3 1
   12: 3 2 2 1
   16: 4 1
   24: 4 2 2 1
   30: 3 3 1
   32: 5 1
   36: 4 2 1
   48: 5 2 2 1
   60: 4 3 2 2 1
   64: 6 1
   96: 6 2 2 1
  120: 5 3 2 2 1
  128: 7 1
  192: 7 2 2 1
  210: 4 4 1
  216: 6 2 1
  240: 6 3 2 2 1
  256: 8 1
  360: 6 3 3 1
  384: 8 2 2 1
  420: 5 4 2 2 1
		

Crossrefs

Programs

  • Mathematica
    tomseq[n_]:=If[n<=1,{},Most[FixedPointList[Sort[Length/@Split[#]]&,Sort[Last/@FactorInteger[n]]]]];
    omseqs=Table[Total/@tomseq[n],{n,1000}];
    Sort[Table[Position[omseqs,x][[1,1]],{x,Union[omseqs]}]]

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

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

A225486 Maximal frequency depth for the partitions of n.

Original entry on oeis.org

0, 2, 3, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7
Offset: 1

Views

Author

Clark Kimberling, May 08 2013

Keywords

Comments

See A225485 for the definition of frequency depth.
The frequency depth of an integer partition is the number of times one must take the multiset of multiplicities to reach (1). For example, the partition (32211) has frequency depth 5 because we have: (32211) -> (221) -> (21) -> (11) -> (2) -> (1). Differs from A325282 at a(0) and a(1). - Gus Wiseman, Apr 19 2019

Examples

			(See A225485.)
		

Crossrefs

Run lengths are A325258, i.e., first differences of Levine's sequence A011784.

Programs

  • Mathematica
    c[s_] := c[s] = Select[Table[Count[s, i], {i, 1, Max[s]}], # > 0 &]
    f[s_] := f[s] = Drop[FixedPointList[c, s], -2]
    t[s_] := t[s] = Length[f[s]]
    u[n_] := u[n] = Table[t[Part[IntegerPartitions[n], k]],
        {k, 1, Length[IntegerPartitions[n]]}];
    Prepend[Table[Max[u[n]], {n, 2, 10}], 0]
    (* second program *)
    grw[q_]:=Join@@Table[ConstantArray[i,q[[Length[q]-i+1]]],{i,Length[q]}];
    Join@@MapIndexed[ConstantArray[#2[[1]]-1,#1]&,Length[#]-Last[#]&/@NestList[grw,{1,1},6]] (* Gus Wiseman, Apr 19 2019 *)

Formula

a(n) = number of terms in row n of the array in A225485, for n > 0.

Extensions

More terms from Gus Wiseman, Apr 19 2019

A353865 Number of complete rucksack partitions of n. Partitions whose weak run-sums are distinct and cover an initial interval of nonnegative integers.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Jun 04 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). A weak run-sum is the sum of any consecutive constant subsequence.
Do all positive integers appear only finitely many times in this sequence?

Examples

			The a(n) compositions for n = 1, 3, 9, 15, 18:
  (1)  (21)   (4311)       (54321)            (543321)
       (111)  (51111)      (532221)           (654111)
              (111111111)  (651111)           (7611111)
                           (81111111)         (111111111111111111)
                           (111111111111111)
For example, the weak runs of y = {7,5,4,4,3,3,3,1,1} are {}, {1}, {1,1}, {3}, {4}, {5}, {3,3}, {7}, {4,4}, {3,3,3}, with sums 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, which are all distinct and cover an initial interval, so y is counted under a(31).
		

Crossrefs

Perfect partitions are counted by A002033, ranked by A325780.
Knapsack partitions are counted by A108917, ranked by A299702.
This is the complete case of A353864, ranked by A353866.
These partitions are ranked by A353867.
A000041 counts partitions, strict A000009.
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.
A353837 counts partitions with distinct run-sums, ranked by A353838.
A353840-A353846 pertain to partition run-sum trajectory.
A353850 counts compositions with all distinct run-sums, ranked by A353852.
A353863 counts partitions whose weak run-sums cover an initial interval.

Programs

  • Mathematica
    norqQ[m_]:=Sort[m]==Range[0,Max[m]];
    msubs[s_]:=Join@@@Tuples[Table[Take[t,i],{t,Split[s]},{i,0,Length[t]}]];
    Table[Length[Select[IntegerPartitions[n],norqQ[Total/@Select[msubs[#],SameQ@@#&]]&]],{n,0,15}]
  • PARI
    a(n) = my(c=0, s, v); if(n, forpart(p=n, if(p[1]==1, v=List([s=1]); for(i=2, #p, if(p[i]==p[i-1], listput(v, s+=p[i]), listput(v, s=p[i]))); s=#v; listsort(v, 1); if(s==#v&&s==v[s], c++))); c, 1); \\ Jinyuan Wang, Feb 21 2025

Extensions

More terms from Jinyuan Wang, Feb 21 2025

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