A329605 -id:A329605 - cp's OEIS frontend

A329606 Lexicographically earliest infinite sequence such that a(i) = a(j) => A329605(i) = A329605(j) for all i, j.

1, 2, 3, 4, 5, 6, 7, 3, 8, 9, 10, 5, 11, 12, 13, 14, 15, 9, 16, 7, 17, 18, 19, 20, 21, 22, 7, 10, 23, 12, 24, 6, 25, 26, 27, 28, 29, 30, 31, 32, 33, 18, 34, 11, 10, 35, 36, 9, 37, 17, 38, 15, 39, 32, 40, 41, 42, 43, 44, 45, 46, 47, 11, 48, 49, 22, 50, 16, 51, 25, 52, 13, 53, 54, 18, 19, 55, 26, 56, 12, 57, 58, 59, 60, 61, 62, 63, 64, 65, 41, 66, 23, 67, 68, 69, 70, 71, 40, 15, 72, 73, 30, 74, 75, 22

Offset: 1

Antti Karttunen, Nov 18 2019

nonn

Restricted growth sequence transform of A329605(n) = A000005(A108951(n)).
For all i, j:
  a(i) = a(j) => A329614(i) = A329614(j).

Cf. A000005, A108951, A329605, A329608, A329614.

up_to = 100000;
rgs_transform(invec) = { my(om = Map(), outvec = vector(length(invec)), u=1); for(i=1, length(invec), if(mapisdefined(om,invec[i]), my(pp = mapget(om, invec[i])); outvec[i] = outvec[pp] , mapput(om,invec[i],i); outvec[i] = u; u++ )); outvec; };
A034386(n) = prod(i=1, primepi(n), prime(i));
A108951(n) = { my(f=factor(n)); prod(i=1, #f~, A034386(f[i, 1])^f[i, 2]) };  \\ From A108951
A329605(n) = numdiv(A108951(n));
v329606 = rgs_transform(vector(up_to, n, A329605(n)));
A329606(n) = v329606[n];

A329608 Lexicographically earliest infinite sequence such that a(i) = a(j) => f(i) = f(j), where f(n) = A329605(n) for all other n, except for odd primes p, f(p) = 0.

Original entry on oeis.org

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

Offset: 1

Antti Karttunen, Nov 18 2019

nonn

For all i, j:

A305801(i) = A305801(j) => a(i) = a(j) => A329614(i) = A329614(j).

Cf. A000005, A108951, A305801, A329605, A329606, A329614.

up_to = 100000;
rgs_transform(invec) = { my(om = Map(), outvec = vector(length(invec)), u=1); for(i=1, length(invec), if(mapisdefined(om,invec[i]), my(pp = mapget(om, invec[i])); outvec[i] = outvec[pp] , mapput(om,invec[i],i); outvec[i] = u; u++ )); outvec; };
A034386(n) = prod(i=1, primepi(n), prime(i));
A108951(n) = { my(f=factor(n)); prod(i=1, #f~, A034386(f[i, 1])^f[i, 2]) };  \\ From A108951
A329605(n) = numdiv(A108951(n));
Aux329608(n) = if((n%2)&&isprime(n),0,A329605(n));
v329608 = rgs_transform(vector(up_to, n, Aux329608(n)));
A329608(n) = v329608[n];

A331284 Number of values of k, 1 <= k <= n, with A329605(k) = A329605(n), where A329605 is the number of divisors of primorial inflation of n (A108951).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 1, 1, 1, 3, 2, 1, 2, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 2, 3, 1, 1, 3, 1, 2, 1, 2, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 2, 1, 2, 1, 2, 1, 2, 1, 1, 3, 2, 1, 2, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 1, 1, 2, 3, 1, 1, 2, 1, 1, 3

Offset: 1

Antti Karttunen, Jan 14 2020

nonn

Ordinal transform of A329605, or equally, of A329606.

Cf. A000005, A108951, A329605, A329606, A331285 (positions of the first occurrences of each n, also positions of records).

Cf. also A067004.

c[n_] := c[n] = If[n == 1, 1, Module[{f = FactorInteger[n], p, e}, If[Length[f] > 1, Times @@ c /@ Power @@@ f, {{p, e}} = f; Times @@ (Prime[Range[PrimePi[p]]]^e)]]];
A329605[n_] := DivisorSigma[0, c[n]];
Module[{b}, b[_] = 0;
a[n_] := With[{t = A329605[n]}, b[t] = b[t] + 1]];
Array[a, 105] (* Jean-François Alcover, Jan 12 2022 *)

up_to = 65537;
ordinal_transform(invec) = { my(om = Map(), outvec = vector(length(invec)), pt); for(i=1, length(invec), if(mapisdefined(om,invec[i]), pt = mapget(om, invec[i]), pt = 0); outvec[i] = (1+pt); mapput(om,invec[i],(1+pt))); outvec; };
A329605(n) = if(1==n,1,my(f=factor(n),e=1,m=1); forstep(i=#f~,1,-1, e += f[i,2]; m *= e^(primepi(f[i,1])-if(1==i,0,primepi(f[i-1,1])))); (m));
v331284 = ordinal_transform(vector(up_to, n, A329605(n)));
A331284(n) = v331284[n];

a(A331285(n)) = n for all n.

A331286 Odd part of number of divisors of primorial inflation of n: a(n) = A000265(A329605(n)).

Original entry on oeis.org

1, 1, 1, 3, 1, 3, 1, 1, 9, 3, 1, 1, 1, 3, 9, 5, 1, 3, 1, 1, 9, 3, 1, 5, 27, 3, 1, 1, 1, 3, 1, 3, 9, 3, 27, 15, 1, 3, 9, 5, 1, 3, 1, 1, 1, 3, 1, 3, 81, 9, 9, 1, 1, 5, 27, 5, 9, 3, 1, 15, 1, 3, 1, 7, 27, 3, 1, 1, 9, 9, 1, 9, 1, 3, 3, 1, 81, 3, 1, 3, 25, 3, 1, 15, 27, 3, 9, 5, 1, 5, 81, 1, 9, 3, 27, 7, 1, 27, 1, 45, 1, 3, 1, 5, 3

Offset: 1

Antti Karttunen, Jan 14 2020

nonn

Antti Karttunen, Table of n, a(n) for n = 1..10000
Antti Karttunen, Data supplement: n, a(n) computed for n = 1..65537

Cf. A000265, A108951, A212181, A329605.

Cf. also A322865.

A331286(n) = if(1==n,1,my(f=factor(n),e=1,m=1); forstep(i=#f~,1,-1, e += f[i,2]; m *= (e>>valuation(e,2))^(primepi(f[i,1])-if(1==i,0,primepi(f[i-1,1])))); (m));

a(n) = A000265(A329605(n)).

a(n) = A212181(A108951(n)).

A331283 a(n) = gcd(n, A329605(n)), where A329605(n) gives the number of divisors of primorial inflation of n, A108951(n).

Original entry on oeis.org

1, 2, 1, 1, 1, 6, 1, 4, 9, 2, 1, 4, 1, 2, 3, 1, 1, 6, 1, 4, 3, 2, 1, 2, 1, 2, 1, 4, 1, 6, 1, 2, 3, 2, 1, 3, 1, 2, 3, 20, 1, 6, 1, 4, 1, 2, 1, 12, 1, 2, 3, 4, 1, 2, 1, 8, 3, 2, 1, 30, 1, 2, 1, 1, 1, 6, 1, 4, 3, 2, 1, 18, 1, 2, 3, 4, 1, 6, 1, 8, 1, 2, 1, 12, 1, 2, 3, 8, 1, 10, 1, 4, 3, 2, 1, 2, 1, 2, 1, 5, 1, 6, 1, 8, 3

Offset: 1

Antti Karttunen, Jan 14 2020

nonn

Antti Karttunen, Table of n, a(n) for n = 1..11025
Antti Karttunen, Data supplement: n, a(n) computed for n = 1..65537
Index entries for sequences computed from indices in prime factorization
Index entries for sequences related to primorial numbers

Cf. A000142, A002110, A034386, A108951, A329605.

Cf. also A009191, A329612.

A331283(n) = if(1==n,1,my(f=factor(n),e=1,m=1); forstep(i=#f~,1,-1, e += f[i,2]; m *= e^(primepi(f[i,1])-if(1==i,0,primepi(f[i-1,1])))); gcd(n,m));

a(n) = gcd(n, A329605(n)).

a(A002110(n)) = gcd(A002110(n), A000142(1+n)) = A034386(1+n), for n >= 0.

A056239 If n = Product_{k >= 1} (p_k)^(c_k) where p_k is k-th prime and c_k >= 0 then a(n) = Sum_{k >= 1} k*c_k.

Original entry on oeis.org

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

Offset: 1

Leroy Quet, Aug 19 2000

A pseudo-logarithmic function in the sense that a(b*c) = a(b)+a(c) and so a(b^c) = c*a(b) and f(n) = k^a(n) is a multiplicative function. [Cf. A248692 for example.] Essentially a function from the positive integers onto the partitions of the nonnegative integers (1->0, 2->1, 3->2, 4->1+1, 5->3, 6->1+2, etc.) so each value a(n) appears A000041(a(n)) times, first with the a(n)-th prime and last with the a(n)-th power of 2. Produces triangular numbers from primorials. - Henry Bottomley, Nov 22 2001
Michael Nyvang writes (May 08 2006) that the Danish composer Karl Aage Rasmussen discovered this sequence in the 1990's: it has excellent musical properties.
All A000041(a(n)) different n's with the same value a(n) are listed in row a(n) of triangle A215366. - Alois P. Heinz, Aug 09 2012
a(n) is the sum of the parts of the partition having Heinz number n. We define the Heinz number of a partition p = [p_1, p_2, ..., p_r] as Product_{j=1..r} (p_j-th prime) (concept used by Alois P. Heinz in A215366 as an "encoding" of a partition). For example, for the partition [1, 1, 2, 4, 10] we get 2*2*3*7*29 = 2436. Example: a(33) = 7 because the partition with Heinz number 33 = 3 * 11 is [2,5]. - Emeric Deutsch, May 19 2015

			a(12) = 1*2 + 2*1 = 4, since 12 = 2^2 *3^1 = (p_1)^2 *(p_2)^1.

N. J. A. Sloane, Table of n, a(n) for n = 1..10000
Index entries for sequences computed from indices in prime factorization

Row sums of A112798.
Cf. A003963 (gives the corresponding products).
Cf. also A000041, A000217, A000720, A001221, A001222, A002110, A005940, A008475, A027748, A049084, A060437, A081401, A082090, A088314, A088318, A088850, A088880, A088881, A088887, A088902, A108951, A122111, A124010, A127668, A141128, A153734, A154351, A156552, A163517, A178503, A215366, A215369, A215501, A242422, A242423, A242424, A243070, A243353, A243354, A243503, A248692, A249336, A249337, A329605, A329902.

a056239 n = sum $ zipWith (*) (map a049084 $ a027748_row n) (a124010_row n)
-- Reinhard Zumkeller, Apr 27 2013

# To get 10000 terms. First make prime table: M:=10000; pl:=array(1..M); for i from 1 to M do pl[i]:=0; od: for i from 1 to M do if ithprime(i) > M then break; fi; pl[ithprime(i)]:=i; od:
# Decode Maple's amazing syntax for factoring integers: g:=proc(n) local e,p,t1,t2,t3,i,j,k; global pl; t1:=ifactor(n); t2:=nops(t1); if t2 = 2 and whattype(t1) <> `*` then p:=op(1,op(1,t1)); e:=op(2,t1); t3:=pl[p]*e; else
t3:=0; for i from 1 to t2 do j:=op(i,t1); if nops(j) = 1 then e:=1; p:=op(1,j); else e:=op(2,j); p:=op(1,op(1,j)); fi; t3:=t3+pl[p]*e; od: fi; t3; end; # N. J. A. Sloane, May 10 2006
A056239 := proc(n) add( numtheory[pi](op(1,p))*op(2,p), p = ifactors(n)[2]) ; end proc: # R. J. Mathar, Apr 20 2010
# alternative:
with(numtheory): a := proc (n) local B: B := proc (n) local nn, j, m: nn := op(2, ifactors(n)): for j to nops(nn) do m[j] := op(j, nn) end do: [seq(seq(pi(op(1, m[i])), q = 1 .. op(2, m[i])), i = 1 .. nops(nn))] end proc: add(B(n)[i], i = 1 .. nops(B(n))) end proc: seq(a(n), n = 1 .. 130); # Emeric Deutsch, May 19 2015

a[1] = 0; a[2] = 1; a[p_?PrimeQ] := a[p] = PrimePi[p];
a[n_] := a[n] = Total[#[[2]]*a[#[[1]]] & /@ FactorInteger[n]]; a /@ Range[91] (* Jean-François Alcover, May 19 2011 *)
Table[Total[FactorInteger[n] /. {p_, c_} /; p > 0 :> PrimePi[p] c], {n, 91}] (* Michael De Vlieger, Jul 12 2017 *)

A056239(n) = if(1==n,0,my(f=factor(n)); sum(i=1, #f~, f[i,2] * primepi(f[i,1]))); \\ Antti Karttunen, Oct 26 2014, edited Jan 13 2020

from sympy import primepi, factorint
def A056239(n): return sum(primepi(p)*e for p, e in factorint(n).items()) # Chai Wah Wu, Jan 01 2023

(require 'factor) ;; Uses the function factor available in Aubrey Jaffer's SLIB Scheme library.
(define (A056239 n) (apply + (map A049084 (factor n))))
;; Antti Karttunen, Oct 26 2014

Totally additive with a(p) = PrimePi(p), where PrimePi(n) = A000720(n).
a(n) = Sum_{k=1..A001221(n)} A049084(A027748(k))*A124010(k). - Reinhard Zumkeller, Apr 27 2013
From Antti Karttunen, Oct 11 2014: (Start)
a(n) = n - A178503(n).
a(n) = A161511(A156552(n)).
a(n) = A227183(A243354(n)).
For all n >= 0:
a(A002110(n)) = A000217(n). [Cf. Henry Bottomley's comment above.]
a(A005940(n+1)) = A161511(n).
a(A243353(n)) = A227183(n).
Also, for all n >= 1:
a(A241909(n)) = A243503(n).
a(A122111(n)) = a(n).
a(A242424(n)) = a(n).
A248692(n) = 2^a(n). (End)
a(n) < A329605(n), a(n) = A001222(A108951(n)), a(A329902(n)) = A112778(n). - Antti Karttunen, Jan 14 2020

A108951 Primorial inflation of n: Fully multiplicative with a(p) = p# for prime p, where x# is the primorial A034386(x).

Original entry on oeis.org

1, 2, 6, 4, 30, 12, 210, 8, 36, 60, 2310, 24, 30030, 420, 180, 16, 510510, 72, 9699690, 120, 1260, 4620, 223092870, 48, 900, 60060, 216, 840, 6469693230, 360, 200560490130, 32, 13860, 1021020, 6300, 144, 7420738134810, 19399380, 180180, 240, 304250263527210, 2520

Offset: 1

Paul Boddington, Jul 21 2005

This sequence is a permutation of A025487.
And thus also a permutation of A181812, see the formula section. - Antti Karttunen, Jul 21 2014
A previous description of this sequence was: "Multiplicative with a(p^e) equal to the product of the e-th powers of all primes at most p" (see extensions), Giuseppe Coppoletta, Feb 28 2015

			a(12) = a(2^2) * a(3) = (2#)^2 * (3#) = 2^2 * 6 = 24
a(45) = (3#)^2 * (5#) = (2*3)^2 * (2*3*5) = 1080 (as 45 = 3^2 * 5).

Amiram Eldar, Table of n, a(n) for n = 1..2370 (terms 1..256 from Antti Karttunen)
Index to divisibility sequences.
Index entries for sequences computed from indices in prime factorization.
Index entries for sequences related to primorial numbers.

Cf. A319626, A329900 (left inverses).

Cf. A034386, A002110, A025487, A048673, A064216, A064989, A085082, A122111, A124859, A161360, A181811, A181812, A181814, A181815, A181817, A181819, A181822, A238690, A283477, A283478, A307035, A324886, A324887, A324888, A324896, A325226, A329040, A329046, A329047, A329344, A329348, A329349, A329378, A329382, A329600, A329602, A329605, A329607, A329615, A329616, A329617, A329619, A329622, A319627, A329647, A331292, A337474, A346108, A346109, A344698, A344699.

a[n_] := a[n] = Module[{f = FactorInteger[n], p, e}, If[Length[f]>1, Times @@ a /@ Power @@@ f, {{p, e}} = f; Times @@ (Prime[Range[PrimePi[p]]]^e)]]; a[1] = 1; Table[a[n], {n, 1, 42}] (* Jean-François Alcover, Feb 24 2015 *)
Table[Times @@ Map[#1^#2 & @@ # &, FactorInteger[n] /. {p_, e_} /; e > 0 :> {Times @@ Prime@ Range@ PrimePi@ p, e}], {n, 42}] (* Michael De Vlieger, Mar 18 2017 *)

primorial(n)=prod(i=1,primepi(n),prime(i))
a(n)=my(f=factor(n)); prod(i=1,#f~, primorial(f[i,1])^f[i,2]) \\ Charles R Greathouse IV, Jun 28 2015

from sympy import primerange, factorint
from operator import mul
def P(n): return reduce(mul, [i for i in primerange(2, n + 1)])
def a(n):
    f = factorint(n)
    return 1 if n==1 else reduce(mul, [P(i)**f[i] for i in f])
print([a(n) for n in range(1, 101)]) # Indranil Ghosh, May 14 2017

def sharp_primorial(n): return sloane.A002110(prime_pi(n))
def p(f):
    return sharp_primorial(f[0])^f[1]
[prod(p(f) for f in factor(n)) for n in range (1,51)]
# Giuseppe Coppoletta, Feb 07 2015

Dirichlet g.f.: 1/(1-2*2^(-s))/(1-6*3^(-s))/(1-30*5^(-s))...
Completely multiplicative with a(p_i) = A002110(i) = prime(i)#. [Franklin T. Adams-Watters, Jun 24 2009; typos corrected by Antti Karttunen, Jul 21 2014]
From Antti Karttunen, Jul 21 2014: (Start)
a(1) = 1, and for n > 1, a(n) = n * a(A064989(n)).
a(n) = n * A181811(n).
a(n) = A002110(A061395(n)) * A331188(n). - [added Jan 14 2020]
a(n) = A181812(A048673(n)).
Other identities:
A006530(a(n)) = A006530(n). [Preserves the largest prime factor of n.]
A071178(a(n)) = A071178(n). [And also its exponent.]
a(2^n) = 2^n. [Fixes the powers of two.]
A067029(a(n)) = A007814(a(n)) = A001222(n). [The exponent of the least prime of a(n), that prime always being 2 for n>1, is equal to the total number of prime factors in n.]
(End)
From Antti Karttunen, Nov 19 2019: (Start)
Further identities:
a(A307035(n)) = A000142(n).
a(A003418(n)) = A181814(n).
a(A025487(n)) = A181817(n).
a(A181820(n)) = A181822(n).
a(A019565(n)) = A283477(n).
A001221(a(n)) = A061395(n).
A001222(a(n)) = A056239(n).
A181819(a(n)) = A122111(n).
A124859(a(n)) = A181821(n).
A085082(a(n)) = A238690(n).
A328400(a(n)) = A329600(n). (smallest number with the same set of distinct prime exponents)
A000188(a(n)) = A329602(n). (square root of the greatest square divisor)
A072411(a(n)) = A329378(n). (LCM of exponents of prime factors)
A005361(a(n)) = A329382(n). (product of exponents of prime factors)
A290107(a(n)) = A329617(n). (product of distinct exponents of prime factors)
A000005(a(n)) = A329605(n). (number of divisors)
A071187(a(n)) = A329614(n). (smallest prime factor of number of divisors)
A267115(a(n)) = A329615(n). (bitwise-AND of exponents of prime factors)
A267116(a(n)) = A329616(n). (bitwise-OR of exponents of prime factors)
A268387(a(n)) = A329647(n). (bitwise-XOR of exponents of prime factors)
A276086(a(n)) = A324886(n). (prime product form of primorial base expansion)
A324580(a(n)) = A324887(n).
A276150(a(n)) = A324888(n). (digit sum in primorial base)
A267263(a(n)) = A329040(n). (number of distinct nonzero digits in primorial base)
A243055(a(n)) = A329343(n).
A276088(a(n)) = A329348(n). (least significant nonzero digit in primorial base)
A276153(a(n)) = A329349(n). (most significant nonzero digit in primorial base)
A328114(a(n)) = A329344(n). (maximal digit in primorial base)
A062977(a(n)) = A325226(n).
A097248(a(n)) = A283478(n).
A324895(a(n)) = A324896(n).
A324655(a(n)) = A329046(n).
A327860(a(n)) = A329047(n).
A329601(a(n)) = A329607(n).
(End)
a(A181815(n)) = A025487(n), and A319626(a(n)) = A329900(a(n)) = n. - Antti Karttunen, Dec 29 2019
From Antti Karttunen, Jul 09 2021: (Start)
a(n) = A346092(n) + A346093(n).
a(n) = A346108(n) - A346109(n).
a(A342012(n)) = A004490(n).
a(A337478(n)) = A336389(n).
A336835(a(n)) = A337474(n).
A342002(a(n)) = A342920(n).
A328571(a(n)) = A346091(n).
A328572(a(n)) = A344592(n).
(End)
Sum_{n>=1} 1/a(n) = A161360. - Amiram Eldar, Aug 04 2022

More terms computed by Antti Karttunen, Jul 21 2014
The name of the sequence was changed for more clarity, in accordance with the above remark of Franklin T. Adams-Watters (dated Jun 24 2009). It is implicitly understood that a(n) is then uniquely defined by completely multiplicative extension. - Giuseppe Coppoletta, Feb 28 2015
Name "Primorial inflation" (coined by Matthew Vandermast in A181815) prefixed to the name by Antti Karttunen, Jan 14 2020

A002183 Number of divisors of n-th highly composite number.

Original entry on oeis.org

1, 2, 3, 4, 6, 8, 9, 10, 12, 16, 18, 20, 24, 30, 32, 36, 40, 48, 60, 64, 72, 80, 84, 90, 96, 100, 108, 120, 128, 144, 160, 168, 180, 192, 200, 216, 224, 240, 256, 288, 320, 336, 360, 384, 400, 432, 448, 480, 504, 512, 576, 600, 640, 672, 720, 768, 800, 864, 896

Offset: 1

N. J. A. Sloane

Record values of tau(n).
RECORDS transform of A000005.
All powers of 2 are present through 2^17. No power of 2 above that is present at least through 2^51. - Comment from Robert G. Wilson v, modified by Ray Chandler, Nov 10 2005
No power of 2 above 2^17 is contained in this sequence - see McRae link for proof. - Graeme McRae, Apr 27 2006
All numbers of the form 9*2^n are present for n=0 through n=30. - Richard Peterson, Sep 07 2024

S. Ramanujan, Collected Papers, Ed. G. H. Hardy et al., Cambridge 1927; Chelsea, NY, 1962, p. 87.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Ray Chandler, Table of n, a(n) for n = 1..10000 (first 1000 terms from T. D. Noe)
A. Flammenkamp, First 1200 highly composite numbers.
Graeme McRae, Highly Composite Numbers.
S. Ramanujan, Table of First 103 Highly Composite Numbers.
N. J. A. Sloane, Transforms.
Eric Weisstein's World of Mathematics, Highly Composite Number.

Cf. A000005, A002182, A002201, A006218, A061799, A070319, A243220, A261100, A329605, A329902.

import Data.List (nub)
a002183 n = a002183_list !! (n-1)
a002183_list = nub $ map (a000005 . a061799) [1..]
-- Reinhard Zumkeller, Apr 01 2011

Reap[ For[ record = 0; n = 1, n <= 10^9, n = If[n < 60, n+1, n+60], tau = DivisorSigma[0, n]; If[tau > record, record = tau; Print[tau]; Sow[tau]]]][[2, 1]] (* Jean-François Alcover, Aug 13 2013 *)
DeleteDuplicates[DivisorSigma[0,Range[3*10^6]],GreaterEqual] (* The program generates the first 42 terms of the sequence. *) (* Harvey P. Dale, Aug 12 2025 *)

a(n) = A000005(A002182(n)).
Also record values of differences A006218(p)-A006218(p-1). These record values occur for any p = A002182(q) where q>=2. - Philippe LALLOUET (philip.lallouet(AT)wanadoo.fr), Jun 23 2007
a(A261100(n)) = A070319(n). - Antti Karttunen, Jun 06 2017
a(n) = A329605(A329902(n)). - Antti Karttunen, Jan 14 2020

More terms from Robert G. Wilson v, Jul 24 2002

A329902 Primorial deflation of the n-th highly composite number: the unique integer k such that A108951(k) = A002182(n).

Original entry on oeis.org

1, 2, 4, 3, 6, 12, 9, 24, 10, 20, 15, 40, 30, 60, 28, 21, 56, 42, 84, 63, 168, 126, 336, 140, 66, 189, 280, 132, 99, 264, 198, 528, 220, 396, 297, 440, 792, 156, 117, 312, 234, 624, 260, 468, 351, 520, 936, 390, 1040, 1872, 780, 585, 306, 1560, 340, 612, 459, 680, 1224, 510, 1360, 2448, 1020, 765, 342, 2040, 1530, 684, 513

Offset: 1

Antti Karttunen, Dec 22 2019

nonn

Antti Karttunen, Table of n, a(n) for n = 1..10000 (computed from the b-file of A002182)
Michael De Vlieger, List of a(n) for n = 1..779674 (gzipped text file, prepared from Flammenkamp's data)
Michael De Vlieger, List of a(n) for n = 1..500000 (noncompressed text file, an abridged version of above)

Cf. A002182, A002183, A056239, A108951, A112778, A112779, A306802, A319626, A324381, A324382, A324581, A324582, A324888, A329040, A329605, A329900, A330743, A330748.

Subsequence of A181815.

Map[Times @@ Prime@(TakeWhile[Reap[FixedPointList[Block[{k = 1}, While[Mod[#, Prime@ k] == 0, k++]; Sow[k - 1]; #/Product[Prime@ i, {i, k - 1}]] &, #]][[-1, 1]], # > 0 &]) &, Take[Import["https://oeis.org/b002182.txt", "Data"][[All, -1]], 69] ] (* Michael De Vlieger, Jan 13 2020, imports b-file at A002182 *)

a(n) = A329900(A002182(n)) = A319626(A002182(n)).
a(n) = A181815(A306802(n)).
A108951(a(n)) = A002182(n). [Highly composite numbers (undeflated)]
A056239(a(n)) = A112778(n). [Number of prime factors, counted with multiplicity]
A001222(a(n)) = A112779(n). [Largest exponent in the prime factorization]
A329605(a(n)) = A002183(n). [Number of divisors]
A329040(a(n)) = A324381(n).
A324888(a(n)) = A324382(n).
a(A330748(n)) = A330743(n).

More linking formulas added by Antti Karttunen, Jan 13 2020

A329382 Product of exponents of prime factors of A108951(n), where A108951 is fully multiplicative with a(prime(i)) = prime(i)# = Product_{i=1..i} A000040(i).

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 3, 4, 2, 1, 3, 1, 2, 4, 4, 1, 6, 1, 3, 4, 2, 1, 4, 8, 2, 9, 3, 1, 6, 1, 5, 4, 2, 8, 8, 1, 2, 4, 4, 1, 6, 1, 3, 9, 2, 1, 5, 16, 12, 4, 3, 1, 12, 8, 4, 4, 2, 1, 8, 1, 2, 9, 6, 8, 6, 1, 3, 4, 12, 1, 10, 1, 2, 18, 3, 16, 6, 1, 5, 16, 2, 1, 8, 8, 2, 4, 4, 1, 12, 16, 3, 4, 2, 8, 6, 1, 24, 9, 16, 1, 6, 1, 4, 18

Offset: 1

Antti Karttunen, Nov 17 2019

nonn

Also the product of parts of the conjugate of the integer partition with Heinz number n, where the Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). For example, the partition (3,2) with Heinz number 15 has conjugate (2,2,1) with product a(15) = 4. - Gus Wiseman, Mar 27 2022

Antti Karttunen, Table of n, a(n) for n = 1..10201
Antti Karttunen, Data supplement: n, a(n) computed for n = 1..65537
Index entries for sequences computed from indices in prime factorization
Index entries for sequences related to primorial numbers

Cf. A000005, A005361, A019565, A034386, A108951, A284001, A331188, A329378, A329605, A329617.
This is the conjugate version of A003963 (product of prime indices).
The solutions to a(n) = A003963(n) are A325040, counted by A325039.
The Heinz number of the conjugate partition is given by A122111.
These are the row products of A321649 and of A321650.
A000700 counts self-conj partitions, ranked by A088902, complement A330644.
A008480 counts permutations of prime indices, conjugate A321648.
A056239 adds up prime indices, row sums of A112798 and of A296150.
A124010 gives prime signature, sorted A118914, sum A001222.
A238744 gives the conjugate of prime signature, rank A238745.
Cf. A000701, A000720, A001221, A046682, A175508, A290822, A303975, A316524, A324850, A352486-A352491, A353570.

Table[Times @@ FactorInteger[Times @@ Map[#1^#2 & @@ # &, FactorInteger[n] /. {p_, e_} /; e > 0 :> {Times @@ Prime@ Range@ PrimePi@ p, e}]][[All, -1]], {n, 105}] (* Michael De Vlieger, Jan 21 2020 *)

A005361(n) = factorback(factor(n)[, 2]); \\ from A005361
A034386(n) = prod(i=1, primepi(n), prime(i));
A108951(n) = { my(f=factor(n)); prod(i=1, #f~, A034386(f[i, 1])^f[i, 2]) };  \\ From A108951
A329382(n) = A005361(A108951(n));

A329382(n) = if(1==n,1,my(f=factor(n),e=0,m=1); forstep(i=#f~,1,-1, e += f[i,2]; m *= e^(primepi(f[i,1])-if(1==i,0,primepi(f[i-1,1])))); (m)); \\ Antti Karttunen, Jan 14 2020

a(n) = A005361(A108951(n)).
A329605(n) >= a(n) >= A329617(n) >= A329378(n).
a(A019565(n)) = A284001(n).
From Antti Karttunen, Jan 14 2020: (Start)
If n = p(k1)^e(k1) * p(k2)^e(k2) * p(k3)^e(k3) * ... * p(kx)^e(kx), with p(n) = A000040(n) and k1 > k2 > k3 > ... > kx, then a(n) = e(k1)^(k1-k2) * (e(k1)+e(k2))^(k2-k3) * (e(k1)+e(k2)+e(k3))^(k3-k4) * ... * (e(k1)+e(k2)+...+e(kx))^kx.
a(n) = A000005(A331188(n)) = A329605(A052126(n)).
(End)
a(n) = A003963(A122111(n)). - Gus Wiseman, Mar 27 2022

cp's OEIS Frontend

A329606 Lexicographically earliest infinite sequence such that a(i) = a(j) => A329605(i) = A329605(j) for all i, j.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

A329608 Lexicographically earliest infinite sequence such that a(i) = a(j) => f(i) = f(j), where f(n) = A329605(n) for all other n, except for odd primes p, f(p) = 0.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

A331284 Number of values of k, 1 <= k <= n, with A329605(k) = A329605(n), where A329605 is the number of divisors of primorial inflation of n (A108951).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A331286 Odd part of number of divisors of primorial inflation of n: a(n) = A000265(A329605(n)).

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

Formula

A331283 a(n) = gcd(n, A329605(n)), where A329605(n) gives the number of divisors of primorial inflation of n, A108951(n).

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

Formula

A056239 If n = Product_{k >= 1} (p_k)^(c_k) where p_k is k-th prime and c_k >= 0 then a(n) = Sum_{k >= 1} k*c_k.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Python

Scheme

Formula

A108951 Primorial inflation of n: Fully multiplicative with a(p) = p# for prime p, where x# is the primorial A034386(x).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Python

Sage

Formula

Extensions

A002183 Number of divisors of n-th highly composite number.