A071625 -id:A071625 - cp's OEIS frontend

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

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

Matthew Vandermast, Dec 07 2010

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

			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.

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

Column 1 of A353510 (see also A325239 and A325277).
Cf. A000040, A000265, A001511, A001222, A003963, A005361, A007814, A008578, A028234, A046523, A056239, A064553, A064989, A067029, A101296 (restricted growth sequence transform), A108951, A122111, A124010, A124859, A156552, A181820, A181821, A182850, A182855, A182857 (also A323014), A115621, A101296, A238690, A238745, A238747, A238748, A246029, A304465, A304647, A305732, A305733, A320118, A323022, A325501, A325502, A325507, A325508, A325755 (A353566), A325756, A328830 [= a(a(n))], A328835, A351564 (characteristic function of A130091), A351944, A351946, A353379.
Left inverse of A304660.

a181819 = product . map a000040 . a124010_row
-- Reinhard Zumkeller, Mar 26 2012

A181819 := proc(n)
    local a;
    a := 1;
    for pf in ifactors(n)[2] do
        a := a*ithprime(pf[2]) ;
    end do:
    a ;
end proc:
seq(A181819(n),n=1..80) ; # R. J. Mathar, Jan 09 2019
# second Maple program:
a:= n-> mul(ithprime(i[2]), i=ifactors(n)[2]):
seq(a(n), n=1..105);  # Alois P. Heinz, Nov 26 2024

{1}~Join~Table[Times @@ Prime@ Map[Last, FactorInteger@ n], {n, 2, 120}] (* Michael De Vlieger, Feb 07 2016 *)

a(n) = {my(f=factor(n)); prod(k=1, #f~, prime(f[k,2]));} \\ Michel Marcus, Nov 16 2015

;; With memoization-macro definec.
(definec (A181819 n) (cond ((= 1 n) 1) (else (* (A000040 (A067029 n)) (A181819 (A028234 n))))))
;; Antti Karttunen, Feb 05 2016

;; With memoization-macro definec.
(definec (A181819 n) (cond ((= 1 n) 1) ((even? n) (* (A000040 (A007814 n)) (A181819 (A000265 n)))) (else (A181819 (A064989 n)))))
;; Antti Karttunen, Feb 07 2016

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)

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

A181821 a(n) = smallest integer with factorization as Product p(i)^e(i) such that Product p(e(i)) = n.

Original entry on oeis.org

1, 2, 4, 6, 8, 12, 16, 30, 36, 24, 32, 60, 64, 48, 72, 210, 128, 180, 256, 120, 144, 96, 512, 420, 216, 192, 900, 240, 1024, 360, 2048, 2310, 288, 384, 432, 1260, 4096, 768, 576, 840, 8192, 720, 16384, 480, 1800, 1536, 32768, 4620, 1296, 1080, 1152, 960, 65536

Offset: 1

Matthew Vandermast, Dec 07 2010

nonn

A permutation of A025487. a(n) is the member m of A025487 such that A181819(m) = n. a(n) is also the member of A025487 whose prime signature is conjugate to the prime signature of A108951(n).

If n = Product_i prime(e(i)) with the e(i) weakly decreasing, then a(n) = Product_i prime(i)^e(i). For example, 90 = prime(3) * prime(2) * prime(2) * prime(1), so a(90) = prime(1)^3 * prime(2)^2 * prime(3)^2 * prime(4)^1 = 12600. - Gus Wiseman, Jan 02 2019

			The canonical factorization of 24 is 2^3*3^1. Therefore, p(e(i)) = prime(3)*prime(1)(i.e., A000040(3)*A000040(1)), which equals 5*2 = 10. Since 24 is the smallest integer for which p(e(i)) = 10, a(10) = 24.

Alois P. Heinz, Table of n, a(n) for n = 1..10000
Eric Weisstein's World of Mathematics, Conjugate Partition

Other rearrangements of A025487 include A036035, A059901, A063008, A077569, A085988, A086141, A087443, A108951, A181822.
Cf. A046523, A181819, A181820.
Cf. A000040, A001221, A001222, A002110, A056239, A071625, A112798, A118914, A122111, A124859, A182850, A305936, A361808.

a:= n-> (l-> mul(ithprime(i)^l[i], i=1..nops(l)))(sort(map(i->
             numtheory[pi](i[1])$i[2], ifactors(n)[2]), `>`)):
seq(a(n), n=1..70);  # Alois P. Heinz, Sep 05 2018

With[{s = Array[If[# == 1, 1, Times @@ Map[Prime@ Last@ # &, FactorInteger@ #]] &, 2^16]}, Array[First@ FirstPosition[s, #] &, LengthWhile[Differences@ Union@ s, # == 1 &]]] (* Michael De Vlieger, Dec 17 2018 *)
Table[Times@@MapIndexed[Prime[#2[[1]]]^#1&,Reverse[Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]]],{n,30}] (* Gus Wiseman, Jan 02 2019 *)

A181821(n) = { my(f=factor(n),p=0,m=1); forstep(i=#f~,1,-1,while(f[i,2], f[i,2]--; m *= (p=nextprime(p+1))^primepi(f[i,1]))); (m); }; \\ Antti Karttunen, Dec 10 2018

from math import prod
from sympy import prime, primepi, factorint
def A181821(n): return prod(prime(i)**e for i, e in enumerate(sorted(map(primepi,factorint(n,multiple=True)),reverse=True),1)) # Chai Wah Wu, Sep 15 2023

If A108951(n) = Product p(i)^e(i), then a(n) = Product A002110(e(i)). I.e., a(n) = A108951(A181819(A108951(n))).
a(A181819(n)) = A046523(n). - [See also A124859]. Antti Karttunen, Dec 10 2018
a(n) = A025487(A361808(n)). - Pontus von Brömssen, Mar 25 2023
a(n) = A108951(A122111(n)). - Antti Karttunen, Sep 15 2023

Definition corrected by Gus Wiseman, Jan 02 2019

A006939 Chernoff sequence: a(n) = Product_{k=1..n} prime(k)^(n-k+1).

Original entry on oeis.org

1, 2, 12, 360, 75600, 174636000, 5244319080000, 2677277333530800000, 25968760179275365452000000, 5793445238736255798985527240000000, 37481813439427687898244906452608585200000000, 7517370874372838151564668004911177464757864076000000000, 55784440720968513813368002533861454979548176771615744085560000000000

Offset: 0

N. J. A. Sloane

Product of first n primorials: a(n) = Product_{i=1..n} A002110(i).
Superprimorials, from primorials by analogy with superfactorials.
Smallest number k with n distinct exponents in its prime factorization, i.e., A071625(k) = n.
Subsequence of A130091. - Reinhard Zumkeller, May 06 2007
Hankel transform of A171448. - Paul Barry, Dec 09 2009
This might be a good place to explain the name "Chernoff sequence" since his name does not appear in the References or Links as of Mar 22 2014. - Jonathan Sondow, Mar 22 2014
Pickover (1992) named this sequence after Paul Chernoff of California, who contributed this sequence to his book. He was possibly referring to American mathematician Paul Robert Chernoff (1942 - 2017), a professor at the University of California. - Amiram Eldar, Jul 27 2020

			a(4) = 360 because 2^3 * 3^2 * 5 = 1 * 2 * 6 * 30 = 360.
a(5) = 75600 because 2^4 * 3^3 * 5^2 * 7 = 1 * 2 * 6 * 30 * 210 = 75600.

Clifford A. Pickover, Mazes for the Mind, St. Martin's Press, NY, 1992, p. 351.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
James K. Strayer, Elementary number theory, Waveland Press, Inc., Long Grove, IL, 1994. See p. 37.

T. D. Noe, Table of n, a(n) for n=0..25

Cf. A000178 (product of first n factorials), A007489 (sum of first n factorials), A060389 (sum of first n primorials).
Cf. A002110, A051357.
A000142 counts divisors of superprimorials.
A000325 counts uniform divisors of superprimorials.
A008302 counts divisors of superprimorials by bigomega.
A022915 counts permutations of prime indices of superprimorials.
A076954 is a sister-sequence.
A118914 has row a(n) equal to {1..n}.
A124010 has row a(n) equal to {n..1}.
A130091 lists numbers with distinct prime multiplicities.
A317829 counts factorizations of superprimorials.
A336417 counts perfect-power divisors of superprimorials.
A336426 gives non-products of superprimorials.
Cf. A001221, A001222, A005117, A022559, A071625, A181796, A181819, A336419, A336420, A336496.

a006939 n = a006939_list !! n
a006939_list = scanl1 (*) a002110_list -- Reinhard Zumkeller, Jul 21 2012

[1] cat [(&*[NthPrime(k)^(n-k+1): k in [1..n]]): n in [1..15]]; // G. C. Greubel, Oct 14 2018

a := []; printlevel := -1; for k from 0 to 20 do a := [op(a),product(ithprime(i)^(k-i+1),i=1..k)] od; print(a);

Rest[FoldList[Times,1,FoldList[Times,1,Prime[Range[15]]]]] (* Harvey P. Dale, Jul 07 2011 *)
Table[Times@@Table[Prime[i]^(n - i + 1), {i, n}], {n, 12}] (* Alonso del Arte, Sep 30 2011 *)

a(n)=prod(k=1,n,prime(k)^(n-k+1)) \\ Charles R Greathouse IV, Jul 25 2011

from math import prod
from sympy import prime
def A006939(n): return prod(prime(k)**(n-k+1) for k in range(1,n+1)) # Chai Wah Wu, Aug 12 2025

a(n) = m(1)*m(2)*m(3)*...*m(n), where m(n) = n-th primorial number. - N. J. A. Sloane, Feb 20 2005
a(0) = 1, a(n) = a(n - 1)p(n)#, where p(n)# is the n-th primorial A002110(n) (the product of the first n primes). - Alonso del Arte, Sep 30 2011
log a(n) = n^2(log n + log log n - 3/2 + o(1))/2. - Charles R Greathouse IV, Mar 14 2011
A181796(a(n)) = A000110(n+1). It would be interesting to have a bijective proof of this theorem, which is stated at A181796 without proof. See also A336420. - Gus Wiseman, Aug 03 2020

Corrected and extended by Labos Elemer, May 30 2001

A323014 a(1) = 0; a(prime) = 1; otherwise a(n) = 1 + a(A181819(n)).

Original entry on oeis.org

0, 1, 1, 2, 1, 3, 1, 2, 2, 3, 1, 4, 1, 3, 3, 2, 1, 4, 1, 4, 3, 3, 1, 4, 2, 3, 2, 4, 1, 3, 1, 2, 3, 3, 3, 3, 1, 3, 3, 4, 1, 3, 1, 4, 4, 3, 1, 4, 2, 4, 3, 4, 1, 4, 3, 4, 3, 3, 1, 5, 1, 3, 4, 2, 3, 3, 1, 4, 3, 3, 1, 4, 1, 3, 4, 4, 3, 3, 1, 4, 2, 3, 1, 5, 3, 3, 3, 4, 1, 5, 3, 4, 3, 3, 3, 4, 1, 4, 4, 3, 1, 3, 1, 4, 3

Offset: 1

Gus Wiseman, Jan 02 2019

nonn

Except for n = 2, same as A182850. Unlike A182850, the terms of this sequence depend only on the prime signature (A101296, A118914) of the index.

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

Positions of 1's are the prime numbers A000040.
Positions of 2's are the proper prime powers A246547.
Positions of 3's are A182853.
Row lengths of A323023.
Cf. A001221, A001222, A046523, A056239, A070175, A071625, A101296, A112798, A118914, A181819, A181821, A182850, A182857, A304465, A323022.

dep[n_]:=If[n==1,0,If[PrimeQ[n],1,1+dep[Times@@Prime/@Last/@FactorInteger[n]]]];
Array[dep,100]

A181819(n) = factorback(apply(e->prime(e),(factor(n)[,2])));
A323014(n) = if(1==n,0,if(isprime(n),1, 1+A323014(A181819(n)))); \\ Antti Karttunen, Jun 10 2022

For all n >= 1, a(n) = a(A046523(n)). [See comment] - Antti Karttunen, Jun 10 2022

Terms a(88) and beyond from Antti Karttunen, Jun 10 2022

A304818 If n = Product_i p(y_i) where p(i) is the i-th prime number and y_i <= y_j for i < j, then a(n) = Sum_i y_i*i.

Original entry on oeis.org

0, 1, 2, 3, 3, 5, 4, 6, 6, 7, 5, 9, 6, 9, 8, 10, 7, 11, 8, 12, 10, 11, 9, 14, 9, 13, 12, 15, 10, 14, 11, 15, 12, 15, 11, 17, 12, 17, 14, 18, 13, 17, 14, 18, 15, 19, 15, 20, 12, 16, 16, 21, 16, 19, 13, 22, 18, 21, 17, 21, 18, 23, 18, 21, 15, 20, 19, 24, 20, 19

Offset: 1

Gus Wiseman, May 18 2018

If n > 1 is not a prime number, we have a(n) >= A056239(n) >= Omega(n) >= omega(n) >= A071625(n) >= ... >= Omicron(n) >= omicron(n) > 1, where Omega = A001222, omega = A001221, Omicron = A304687 and omicron = A304465.

			The multiset of prime indices (see A112798) of 216 is {1,1,1,2,2,2}, which becomes {1,2,3,4,4,5,5,6,6} under A304660, so a(216) = 1+2+3+4+4+5+5+6+6 = 36.

Alois P. Heinz, Table of n, a(n) for n = 1..20000

Cf. A000720, A001221, A001222, A001358, A055932, A056239, A071625, A112798, A181819, A182850, A182857, A275870, A304465, A304660.

a:= n-> (l-> add(i*numtheory[pi](l[i]), i=1..nops(l)))(
             sort(map(i-> i[1]$i[2], ifactors(n)[2]))):
seq(a(n), n=1..100);  # Alois P. Heinz, May 20 2018

primeMS[n_]:=If[n===1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Table[With[{y=primeMS[n]},Sum[y[[i]]*i,{i,Length[y]}]],{n,20}]

a(n) = {my(f = factor(n), s = 0, i = 0); for (k=1, #f~, for (kk = 1, f[k, 2], i++; s += i*primepi(f[k,1]););); s;} \\ Michel Marcus, May 19 2018

vf(n) = {my(f=factor(n), nb = bigomega(n), g = vector(nb), i = 0); for (k=1, #f~, for (kk = 1, f[k, 2], i++; g[i] = primepi(f[k,1]););); return(g);} \\ A112798
a(n) = {my(g = vf(n)); sum(k=1, #g, k*g[k]);} \\ Michel Marcus, May 19 2018

a(n) = A056239(A304660(n)).

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

Gus Wiseman, May 23 2022

nonn

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.

			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.

Antti Karttunen, Table of n, a(n) for n = 1..20000
Mathematics Stack Exchange, What is a sequence run? (answered 2011-12-01)
Index entries for sequences related to prime indices in the factorization of n.

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

Table[Times@@Prime/@Cases[If[n==1,{},FactorInteger[n]],{p_,k_}:>PrimePi[p]*k],{n,100}]

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

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

More terms from Antti Karttunen, Jan 20 2025

A323023 Irregular triangle read by rows where row n is the omega-sequence of n.

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Jan 02 2019

We define the omega-sequence of n to have length A323014(n), and the k-th term is Omega(red^{k-1}(n)), where Omega = A001222 and red^{k} is the k-th functional iteration of A181819.

Except for n = 1, all rows end with 1. If n is not prime, the term in row n prior to the last is A304465(n).

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

Row lengths are A323014, or A182850 if we assume A182850(2) = 1.
First column is empty if n = 1 and otherwise A001222(n).
Second column is empty if n is 1 or prime and otherwise A001221(n).
Third column is empty if n is 1, prime, or a power of a prime and otherwise A071625(n).
Cf. A024619, A056239, A067340, A118914, A181819, A181821, A182857, A304464, A304465, A323022.

red[n_]:=Times@@Prime/@Last/@If[n==1,{},FactorInteger[n]];
omg[n_,k_]:=If[k==1,PrimeOmega[n],omg[red[n],k-1]];
dep[n_]:=If[n==1,0,If[PrimeQ[n],1,1+dep[Times@@Prime/@Last/@If[n==1,{},FactorInteger[n]]]]];
Table[omg[n,k],{n,100},{k,dep[n]}]

A353837 Number of integer partitions of n with all distinct run-sums.

Original entry on oeis.org

1, 1, 2, 3, 4, 7, 10, 14, 17, 28, 35, 49, 62, 85, 107, 149, 174, 238, 305, 384, 476, 614, 752, 950, 1148, 1451, 1763, 2205, 2654, 3259, 3966, 4807, 5773, 7039, 8404, 10129, 12140, 14528, 17288, 20668, 24505, 29062, 34437, 40704, 48059, 56748, 66577, 78228

Offset: 0

Gus Wiseman, May 26 2022

nonn

The run-sums of a sequence are the sums of its maximal consecutive constant subsequences (runs). For example, the run-sums of (2,2,1,1,1,3,2,2) are (4,3,3,4). The first partition whose run-sums are not all distinct is (2,1,1).

			The a(0) = 1 through a(6) = 10 partitions:
  ()  (1)  (2)   (3)    (4)     (5)      (6)
           (11)  (21)   (22)    (32)     (33)
                 (111)  (31)    (41)     (42)
                        (1111)  (221)    (51)
                                (311)    (222)
                                (2111)   (321)
                                (11111)  (411)
                                         (2211)
                                         (21111)
                                         (111111)

Max Alekseyev, Table of n, a(n) for n = 0..100

For multiplicities instead of run-sums we have A098859, ranked by A130091.
For equal run-sums we have A304442, ranked by A353833 (nonprime A353834).
These partitions are ranked by A353838, complement A353839.
The version for compositions is A353850, ranked by A353852.
The weak version (rucksack partitions) is A353864, ranked by A353866.
The weak perfect version is A353865, ranked by A353867.
A005811 counts runs in binary expansion.
A275870 counts collapsible partitions, ranked by A300273.
A351014 counts distinct runs in standard compositions.
A353832 represents the operation of taking run-sums of a partition.
A353840-A353846 pertain to partition run-sum trajectory.
A353849 counts distinct run-sums in standard compositions.
Cf. A000041, A008284, A047966, A071625, A073093, A116608, A175413, A181819, A333755, A353848, A353867.

Table[Length[Select[IntegerPartitions[n],UnsameQ@@Total/@Split[#]&]],{n,0,15}]

a353837 = lambda n: sum( abs(BipartiteGraph( Matrix(len(p), len(D:=list(set.union(*map(lambda t: set(divisors(t)),p)))), lambda i,j: p[i]%D[j]==0) ).matching_polynomial()[len(D)-len(p)]) for p in Partitions(n,max_slope=-1) ) # Max Alekseyev, Sep 11 2023

A182857 Smallest number that requires exactly n iterations to reach a fixed point under the x -> A181819(x) map.

Original entry on oeis.org

1, 3, 4, 6, 12, 60, 2520, 1286485200, 35933692027611398678865941374040400000

Offset: 0

Matthew Vandermast, Jan 05 2011

nonn

a(9) has 296 digits.
Related to Levine's sequence (A011784): A011784(n) = A001222(a(n)) = A001221(a(n+1)) = A051903(a(n+2)) = A071625(a(n+2)). Also see A182858.
Values of n where A182850(n) increases to a record.
The multiplicity of prime(k) in a(n+1) is the k-th largest prime index of a(n), which is A296150(a(n),k). - Gus Wiseman, May 13 2018

			From _Gus Wiseman_, May 13 2018: (Start)
Like A001462 the following sequence of multisets whose Heinz numbers belong to this sequence is a run-length describing sequence, as the number of k's in row n + 1 is equal to the k-th term of row n.
{2}
{1,1}
{1,2}
{1,1,2}
{1,1,2,3}
{1,1,1,2,2,3,4}
{1,1,1,1,2,2,2,3,3,4,4,5,6,7}
{1,1,1,1,1,1,1,2,2,2,2,2,2,3,3,3,3,3,4,4,4,4,5,5,5,5,6,6,6,7,7,7,8,8,9,9,10,10,11,12,13,14}
(End)

Gus Wiseman, Table of n, a(n) for n = 0..9

Cf. A001462, A007755, A007916, A009287, A012257, A112798, A181819, A182850-A182858, A296150, A304455, A304464, A304465.

Prepend[Function[m,Times@@Prime/@m]/@NestList[Join@@Table[Table[i,{Reverse[#][[i]]}],{i,Length[#]}]&,{2},8],1] (* Gus Wiseman, May 13 2018 *)

For n > 0, a(n) = A181819(a(n+1)). For n > 1, a(n) = A181821(a(n-1)).

A353833 Numbers whose multiset of prime indices has all equal run-sums.

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 8, 9, 11, 12, 13, 16, 17, 19, 23, 25, 27, 29, 31, 32, 37, 40, 41, 43, 47, 49, 53, 59, 61, 63, 64, 67, 71, 73, 79, 81, 83, 89, 97, 101, 103, 107, 109, 112, 113, 121, 125, 127, 128, 131, 137, 139, 144, 149, 151, 157, 163, 167, 169, 173, 179

Offset: 1

Gus Wiseman, May 23 2022

nonn

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.

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 prime indices of 12 are {1,1,2}, with run-sums (2,2), so 12 is in the sequence.

Mathematics Stack Exchange, What is a sequence run? (answered 2011-12-01)

For parts instead of run-sums we have A000961, counted by A000005.
For run-lengths instead of run-sums we have A072774, counted by A047966.
These partitions are counted by A304442.
These are the positions of powers of primes in A353832.
The restriction to nonprimes is A353834.
For distinct instead of equal run-sums we have A353838, counted by A353837.
The version for compositions is A353848, counted by A353851.
A001222 counts prime factors, distinct A001221.
A005811 counts runs in binary expansion, distinct run-lengths A165413.
A056239 adds up prime indices, row sums of A112798 and A296150.
A124010 gives prime signature, sorted A118914.
A300273 ranks collapsible partitions, counted by A275870.
A353835 counts distinct run-sums of prime indices, weak A353861.
A353840-A353846 deal with iterated run-sums for partitions.
A353862 gives greatest run-sum of prime indices, least A353931.
A353866 ranks rucksack partitions, counted by A353864.
Cf. A007947, A071625, A073093, A181819, A238279, A304660, A323014, A333755, A353839.

Select[Range[100],SameQ@@Cases[FactorInteger[#],{p_,k_}:>PrimePi[p]*k]&]

cp's OEIS Frontend

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

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Scheme

Scheme

Formula

Extensions

A181821 a(n) = smallest integer with factorization as Product p(i)^e(i) such that Product p(e(i)) = n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

Formula

Extensions

A006939 Chernoff sequence: a(n) = Product_{k=1..n} prime(k)^(n-k+1).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

Python

Formula

Extensions

A323014 a(1) = 0; a(prime) = 1; otherwise a(n) = 1 + a(A181819(n)).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A304818 If n = Product_i p(y_i) where p(i) is the i-th prime number and y_i <= y_j for i < j, then a(n) = Sum_i y_i*i.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

Formula

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

Views

Author