A108951 -id:A108951 - cp's OEIS frontend

A324886 a(n) = A276086(A108951(n)).

2, 3, 5, 9, 7, 25, 11, 15, 35, 49, 13, 625, 17, 121, 117649, 225, 19, 1225, 23, 2401, 1771561, 169, 29, 875, 717409, 289, 55, 14641, 31, 184877, 37, 21, 4826809, 361, 36226650889, 1500625, 41, 529, 24137569, 77, 43, 143, 47, 28561, 1127357, 841, 53, 1715, 902613283, 514675673281, 47045881, 83521, 59, 3025, 8254129, 214358881, 148035889, 961, 61

Offset: 1

Antti Karttunen, Mar 30 2019

nonn

Cf. A034386, A108951, A117366, A276086, A324887, A324888, A324896, A329047, A342920, A344592, A346106, A346107.
Permutation of A324576.
Cf. also A346105.

With[{b = MixedRadix[Reverse@ Prime@ Range@ 120]}, Array[Function[k, Times @@ Power @@@ # &@ Transpose@ {Prime@ Range@ Length@ k, Reverse@ k}]@ IntegerDigits[Apply[Times, Map[#1^#2 & @@ # &, FactorInteger[#] /. {p_, e_} /; e > 0 :> {Times @@ Prime@ Range@ PrimePi@ p, e}]], b] &, 58]] (* Michael De Vlieger, Nov 18 2019 *)
A276086[n0_] := Module[{m = 1, i = 1, n = n0, p}, While[n > 0, p = Prime[i]; m *= p^Mod[n, p]; n = Quotient[n, p]; i++]; m];
(* b is A108951 *)
b[n_] := b[n] = Module[{pe = FactorInteger[n], p, e}, If[Length[pe] > 1, Times @@ b /@ Power @@@ pe, {{p, e}} = pe; Times @@ (Prime[Range[ PrimePi[p]]]^e)]]; b[1] = 1;
a[n_] := A276086[b[n]];
Array[a, 100] (* Jean-François Alcover, Dec 01 2021, after _Antti Karttunen in A296086 *)

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
A276086(n) = { my(i=0,m=1,pr=1,nextpr); while((n>0),i=i+1; nextpr = prime(i)*pr; if((n%nextpr),m*=(prime(i)^((n%nextpr)/pr));n-=(n%nextpr));pr=nextpr); m; };
A324886(n) = A276086(A108951(n));

a(n) = A276086(A108951(n)).
a(n) = A117366(n) * A324896(n).
A001222(a(n)) = A324888(n).
A020639(a(n)) = A117366(n).
A032742(a(n)) = A324896(n).
a(A000040(n)) = A000040(1+n).
From Antti Karttunen, Jul 09 2021: (Start)
For n > 1, a(n) = A003961(A329044(n)).
a(n) = A346091(n) * A344592(n).
a(n) = A346106(n) /  A346107(n).
A003415(a(n)) = A329047(n).
A003557(a(n)) = A344592(n).
A342001(a(n)) = A342920(n) = A329047(n) / A344592(n).
(End)

A324888 Minimal number of primorials (A002110) that add to A108951(n).

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 2, 2, 2, 1, 4, 1, 2, 6, 4, 1, 4, 1, 4, 6, 2, 1, 4, 6, 2, 2, 4, 1, 6, 1, 2, 6, 2, 10, 8, 1, 2, 6, 2, 1, 2, 1, 4, 6, 2, 1, 4, 8, 12, 6, 4, 1, 4, 6, 8, 6, 2, 1, 6, 1, 2, 6, 4, 14, 12, 1, 4, 6, 10, 1, 6, 1, 2, 10, 4, 18, 12, 1, 4, 8, 2, 1, 4, 12, 2, 6, 8, 1, 12, 18, 4, 6, 2, 8, 8, 1, 16, 12, 8, 1, 12, 1, 8, 8

Offset: 1

Antti Karttunen, Mar 30 2019

nonn

Sum of digits when A108951(n) is written in primorial base (A049345).

Cf. A001222, A002110, A034386, A049345, A108951, A276086, A276150, A324886.

Cf. A324383, A324386, A324387 (permutations of this sequence).

With[{b = Reverse@ Prime@ Range@ 120}, Array[Total@ IntegerDigits[#, MixedRadix[b]] &@ Apply[Times, Map[#1^#2 & @@ # &, FactorInteger[#] /. {p_, e_} /; e > 0 :> {Times @@ Prime@ Range@ PrimePi@ p, e}]] &, 105] ] (* Michael De Vlieger, Nov 18 2019 *)

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
A276150(n) = { my(s=0,m); forprime(p=2, , if(!n, return(s)); m = n%p; s += m; n = (n-m)/p); };
A324888(n) = A276150(A108951(n));

a(n) = A276150(A108951(n)).

a(n) = A001222(A324886(n)).

A329605 Number of divisors 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, 2, 4, 3, 8, 6, 16, 4, 9, 12, 32, 8, 64, 24, 18, 5, 128, 12, 256, 16, 36, 48, 512, 10, 27, 96, 16, 32, 1024, 24, 2048, 6, 72, 192, 54, 15, 4096, 384, 144, 20, 8192, 48, 16384, 64, 32, 768, 32768, 12, 81, 36, 288, 128, 65536, 20, 108, 40, 576, 1536, 131072, 30, 262144, 3072, 64, 7, 216, 96, 524288, 256, 1152, 72, 1048576, 18, 2097152, 6144, 48

Offset: 1

Antti Karttunen, Nov 18 2019

nonn

Cf. A000005, A000040, A000142, A002110, A010052, A034386, A052126, A056239, A108951, A329902, A329378, A329382, A329614, A329617, A331283, A331286 (odd part).
Cf. A329606 (rgs-transform), A329608, A331284 (ordinal transform).
Cf. A331285 (the position where for the first time some term has occurred n times in this sequence).

Block[{a}, 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; Array[DivisorSigma[0, a@ #] &, 75]] (* Michael De Vlieger, Jan 24 2020, after Jean-François Alcover at A108951 *)

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

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)); \\ Antti Karttunen, Jan 14 2020

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

a(n) = A000005(A108951(n)).
a(n) >= A329382(n) >= A329617(n) >= A329378(n).
A020639(a(n)) = A329614(n).
From Antti Karttunen, Jan 14 2020: (Start)
a(A052126(n)) = A329382(n).
a(A002110(n)) = A000142(1+n), for all n >= 0.
a(n) > A056239(n).
a(A329902(n)) = A002183(n).
A000265(a(n)) = A331286(n).
gcd(n,a(n)) = A331283(n).
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 > ... > kx, then a(n) = (1+e(k1))^(k1-k2) * (1+e(k1)+e(k2))^(k2-k3) * ... * (1+e(k1)+e(k2)+...+e(kx))^kx.
A000035(a(n)) = A000035(A000005(n)) = A010052(n).
(End)

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

A329040 Number of distinct primorials in the greedy sum of primorials adding to A108951(n).

Original entry on oeis.org

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

Offset: 1

Antti Karttunen, Nov 11 2019

nonn

The greedy sum is also the sum with the minimal number of primorials used in the primorial base representation.

			For n = 18 = 2 * 3^2, A108951(18) = A034386(2) * A034386(3)^2 = 2 * 6^2 = 72 = 2*A002110(3) + 2*A002110(2) = 2*30 + 2*6, and because there occurs only two distinct primorials (30 and 6) in the sum, we have a(18) = 2.

Cf. A001221, A002110, A034386, A108951, A267263, A276086, A324886, A324888, A329051 (positions of records), A329343.

Cf. also A329045, A329046.

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
A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
A324886(n) = A276086(A108951(n));
A329040(n) = omega(A324886(n));

a(n) = A001221(A324886(n)).
a(n) = A267263(A108951(n)).
a(n) <= A324888(n).

A329348 The least significant nonzero digit in the primorial base expansion of primorial inflation of n, A108951(n).

Original entry on oeis.org

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

Offset: 1

Antti Karttunen, Nov 11 2019

nonn

Number of occurrences of the least primorial present in the greedy sum of primorials adding to A108951(n).

The greedy sum is also the sum with the minimal number of primorials, used for example in the primorial base representation.

			For n = 24 = 2^3 * 3, A108951(24) = A034386(2)^3 * A034386(3) = 2^3 * 6 = 48 = 1*30 + 3*6, and as the factor of the least primorial in the sum is 3, we have a(24) = 3.

Cf. A002110, A034386, A067029, A108951, A117366, A276086, A276088, A324886, A324888, A329040, A329345, A329343, A329344, A329349, A331188, A331289, A331290, A331291.

Cf. also A331292 (the next digit to left of this one), A331293.

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
A276088(n) = { my(e=0, p=2); while(n && !(e=(n%p)), n = n/p; p = nextprime(1+p)); (e); };
A329348(n) = A276088(A108951(n));

A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
A324886(n) = A276086(A108951(n));
A067029(n) = if(1==n, 0, factor(n)[1, 2]); \\ From A067029
A329348(n) = A067029(A324886(n));

A002110(n) = prod(i=1, n, prime(i));
A329348(n) = if(1==n, n, my(f=factor(n), p=nextprime(1+vecmax(f[, 1]))); prod(i=1, #f~, A002110(primepi(f[i, 1]))^(f[i, 2]-(#f~==i)))%p); \\ Antti Karttunen, Jan 15 2020

a(n) = A067029(A324886(n)) = A276088(A108951(n)).
a(n) <= A324888(n).
From Antti Karttunen, Jan 15-17 2020: (Start)
a(n) = A331188(n) mod A117366(n).
a(n) = A001511(A246277(A324886(n))).
(End)

Name changed by Antti Karttunen, Jan 17 2020

A344592 a(n) = A003557(A276086(A108951(n))).

Original entry on oeis.org

1, 1, 1, 3, 1, 5, 1, 1, 1, 7, 1, 125, 1, 11, 16807, 15, 1, 35, 1, 343, 161051, 13, 1, 25, 9317, 17, 1, 1331, 1, 2401, 1, 1, 371293, 19, 253333223, 42875, 1, 23, 1419857, 1, 1, 1, 1, 2197, 14641, 29, 1, 49, 371293, 6684099653, 2476099, 4913, 1, 55, 37349, 19487171, 6436343, 31, 1, 5929, 1, 37, 20449, 21, 582622237229761, 1792160394037

Offset: 1

Antti Karttunen, May 26 2021

nonn

Cf. A003557, A085731, A108951, A276086, A324886, A328572, A329047, A342002, A342920, A346091, A346092 [= A276085(a(n))].

Cf. A344591 (positions of ones), A344593 (rgs-transform).

Block[{b = MixedRadix[Reverse@ Prime@ Range@ 20]}, Array[#/(Times @@ FactorInteger[#][[All, 1]]) &@ Apply[Times, Power @@@ #] &@ Transpose@ {Prime@ Range@ Length@ #, Reverse@ #} &@ IntegerDigits[#, b] &@ Apply[Times, Map[(Times @@ Prime@ Range@ PrimePi@ #1)^#2 & @@ # &, FactorInteger[#]]] &, 66]] (* Michael De Vlieger, Jul 14 2021 *)

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
A328572(n) = { my(m=1, p=2); while(n, if(n%p, m *= p^((n%p)-1)); n = n\p; p = nextprime(1+p)); (m); };
A344592(n) = A328572(A108951(n));

A003557(n) = { my(f=factor(n)); for (i=1, #f~, f[i, 2] = f[i, 2]-1); factorback(f); }; \\ From A003557
A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
A344592(n) = A003557(A276086(A108951(n)));

a(n) = A328572(A108951(n)) = A003557(A324886(n)) = A003557(A276086(A108951(n))).
a(n) = A329047(n) / A342920(n).
a(n) = A085731(A324886(n)) = gcd(A324886(n), A329047(n)) = A324886(n) / A346091(n). - Antti Karttunen, Jul 09 2021

A329343 Difference between the indices of the smallest and the largest primorial in the greedy sum of primorials adding to A108951(n).

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 1, 0, 2, 0, 0, 1, 0, 2, 0, 0, 1, 1, 0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 1, 2, 1, 0, 0, 0, 2, 1, 0, 0, 0, 0, 1, 0, 0, 1, 2, 1, 0, 0, 0, 0, 1, 0, 2, 0, 0, 2, 0, 1, 0, 0, 1, 2, 0, 0, 1, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 1, 0, 2, 1, 2, 0, 0, 0, 0, 2

Offset: 1

Antti Karttunen, Nov 11 2019

nonn

The greedy sum is also the sum with the minimal number of primorials, used for example in the primorial base representation.

Positions of the records (and conjecturally, the positions of the first occurrences of each n) begin as 1, 8, 27, 162, 289, 529, 841, 1369, 1681, 2209, 2809, 3481, 4489, 5041, 5329, 6889, ..., that after 162 all seem to be squares of certain primes. See also A329051.

			For n = 18 = 2 * 3^2, A108951(18) = A034386(2) * A034386(3)^2 = 2 * 6^2 = 72 = 30 + 30 + 6 + 6, and as the largest primorial in the sum is 30 = A002110(3), and the least primorial is 6 = A002110(2), we have a(18) = 3-2 = 1.

Cf. A002110, A034386, A108951, A243055, A276086, A324886, A324888, A329040, A329051.

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
A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
A324886(n) = A276086(A108951(n));
A243055(n) = if(1==n,0,my(f = factor(n), lpf = f[1, 1], gpf = f[#f~, 1]); (primepi(gpf)-primepi(lpf)));
A329343(n) = A243055(A324886(n));

a(n) = A243055(A324886(n)).

A329344 Number of times most frequent primorial is present in the greedy sum of primorials adding to A108951(n).

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 1, 1, 2, 1, 4, 1, 2, 6, 2, 1, 2, 1, 4, 6, 2, 1, 3, 4, 2, 1, 4, 1, 5, 1, 1, 6, 2, 8, 4, 1, 2, 6, 1, 1, 1, 1, 4, 5, 2, 1, 3, 6, 8, 6, 4, 1, 2, 4, 8, 6, 2, 1, 3, 1, 2, 3, 2, 13, 12, 1, 4, 6, 5, 1, 3, 1, 2, 5, 4, 16, 12, 1, 2, 6, 2, 1, 2, 11, 2, 6, 8, 1, 10, 12, 4, 6, 2, 7, 6, 1, 12, 10, 6, 1, 12, 1, 8, 4

Offset: 1

Antti Karttunen, Nov 11 2019

nonn

The greedy sum is also the sum with the minimal number of primorials, used for example in the primorial base representation.

			For n = 24 = 2^3 * 3, A108951(24) = A034386(2)^3 * A034386(3) = 2^3 * 6 = 48 = 30 + 6 + 6 + 6, and as the most frequent primorial in the sum is 6 = A002110(2), we have a(24) = 3.

Cf. A002110, A034386, A051903, A108951, A276086, A324886, A324888, A328114, A329040, A329045, A329343, A329348, A329349.

With[{b = Reverse@ Prime@ Range@ 120}, Array[Max@ IntegerDigits[#, MixedRadix[b]] &@ Apply[Times, Map[#1^#2 & @@ # &, FactorInteger[#] /. {p_, e_} /; e > 0 :> {Times @@ Prime@ Range@ PrimePi@ p, e}]] &, 105] ] (* Michael De Vlieger, Nov 18 2019 *)

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
A328114(n) = { my(s=0, p=2); while(n, s = max(s,(n%p)); n = n\p; p = nextprime(1+p)); (s); };
A329344(n) = A328114(A108951(n));

A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
A324886(n) = A276086(A108951(n));
A051903(n) = if((1==n),0,vecmax(factor(n)[, 2]));
A329344(n) = A051903(A324886(n));

a(n) = A328114(A108951(n)) = A051903(A324886(n)).

cp's OEIS Frontend

A324886 a(n) = A276086(A108951(n)).

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A324888 Minimal number of primorials (A002110) that add to A108951(n).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A329605 Number of divisors of A108951(n), where A108951 is fully multiplicative with a(prime(i)) = prime(i)# = Product_{i=1..i} A000040(i).

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

PARI

PARI

Formula

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

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

Formula

Extensions

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

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Formula

A329040 Number of distinct primorials in the greedy sum of primorials adding to A108951(n).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

PARI

Formula

A329348 The least significant nonzero digit in the primorial base expansion of primorial inflation of n, A108951(n).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

PARI

PARI

PARI

Formula