A324886 -id:A324886 - cp's OEIS frontend

A373983 Lexicographically earliest infinite sequence such that a(i) = a(j) = A246277(A324886(i)) = A246277(A324886(j)) and A278226(A328768(i)) = A278226(A328768(j)), for all i, j >= 1.

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

Offset: 1

Antti Karttunen, Jun 25 2024

Restricted growth sequence transform of the ordered pair [A246277(A276086(A108951(n))), A046523(A276086(A328768(n)))].
For all i, j >= 1:
  A305800(i) = A305800(j) => a(i) = a(j),
  a(i) = a(j) => A329345(i) = A329345(j) => A329045(i) = A329045(j),
  a(i) = a(j) => A373982(i) = A373982(j) => A328771(i) = A328771(j).
It is hard to say for sure which graphical features in the scatter plot have their provenance in A373982, and which ones in A329345.

Antti Karttunen, Table of n, a(n) for n = 1..100000
Index entries for sequences related to primorial numbers

Cf. A002110, A046523, A108951, A246277, A276086, A278226, A324886, A328768.
Cf. A305800, A329045, A329345, A328771, A373982.
Cf. also A329620.

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; };
A046523(n) = { my(f=vecsort(factor(n)[, 2], , 4), p); prod(i=1, #f, (p=nextprime(p+1))^f[i]); };
A108951(n) = { my(f=factor(n)); prod(i=1, #f~,  prod(i=1, primepi(f[i, 1]), prime(i))^f[i, 2]); };
A246277(n) = if(1==n, 0, my(f = factor(n), k = primepi(f[1,1])-1); for (i=1, #f~, f[i,1] = prime(primepi(f[i,1])-k)); factorback(f)/2);
A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
A002110(n) = prod(i=1,n,prime(i));
A328768(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]*A002110(primepi(f[i,1])-1)/f[i, 1]));
Aux373983(n) = [A246277(A276086(A108951(n))), A046523(A276086(A328768(n)))];
v373983 = rgs_transform(vector(up_to, n, Aux373983(n)));
A373983(n) = v373983[n];

A329346 a(n) = A322356(A324886(n)).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 5, 7, 1, 1, 1, 1, 1, 1, 5, 1, 7, 1, 1, 1, 1, 1, 7, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 13, 7, 1, 1, 1, 1, 1, 13, 1, 1, 1, 1, 1, 7, 13, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 13, 1, 19, 1, 1, 1, 1, 13, 1, 7, 1, 1, 13, 1, 1, 1, 1, 1, 7, 1, 1, 13, 1, 1, 1, 1, 1, 1, 19, 1, 1, 1, 1, 7, 1, 13, 1, 13, 1, 1, 1, 1, 13

Offset: 1

Antti Karttunen, Nov 11 2019

nonn

			For n = 128 = 2^7, A108951(128) = A034386(2)^7 = 128. As 128 = 4 * 30 + 1*6 + 1* 2, A276086(128) = 36015 = 7^4 * 5^1 * 3^1, and there are two such primes that both p and p-2 divide n, and p-2 is also prime, namely, 7 and 5, thus a(128) = 7*5 = 35. This is also the first occurrence of composite number in this sequence.

Cf. A002110, A034386, A108951, A276086, A324886, A329040, A322356, A329343, , A329344, A329347, A329348.

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));
A322356(n) = { my(f = factor(n), m=1); for(i=1, #f~, if(isprime(f[i,1]+2)&&!(n%(f[i,1]+2)), m *= (f[i,1]+2))); (m); };
A329346(n) = A322356(A324886(n));

a(n) = A322356(A324886(n)).

A276086 Primorial base exp-function: digits in primorial base representation of n become the exponents of successive prime factors whose product a(n) is.

Original entry on oeis.org

1, 2, 3, 6, 9, 18, 5, 10, 15, 30, 45, 90, 25, 50, 75, 150, 225, 450, 125, 250, 375, 750, 1125, 2250, 625, 1250, 1875, 3750, 5625, 11250, 7, 14, 21, 42, 63, 126, 35, 70, 105, 210, 315, 630, 175, 350, 525, 1050, 1575, 3150, 875, 1750, 2625, 5250, 7875, 15750, 4375, 8750, 13125, 26250, 39375, 78750, 49, 98, 147, 294, 441, 882, 245, 490, 735, 1470, 2205, 4410, 1225, 2450

Offset: 0

Antti Karttunen, Aug 21 2016

Prime product form of primorial base expansion of n.
Sequence is a permutation of A048103. It maps the smallest prime not dividing n to the smallest prime dividing n, that is, A020639(a(n)) = A053669(n) holds for all n >= 1.
The sequence satisfies the exponential function identity, a(x + y) = a(x) * a(y), whenever A329041(x,y) = 1, that is, when adding x and y together will not generate any carries in the primorial base. Examples of such pairs of x and y are A328841(n) & A328842(n), and also A328770(n) (when added with itself). - Antti Karttunen, Oct 31 2019
From Antti Karttunen, Feb 18 2022: (Start)
The conjecture given in A327969 asks whether applying this function together with the arithmetic derivative (A003415) in some combination or another can eventually transform every positive integer into zero.
Another related open question asks whether there are any other numbers than n=6 such that when starting from that n and by iterating with A003415, one eventually reaches a(n). See comments in A351088.
This sequence is used in A351255 to list the terms of A099308 in a different order, by the increasing exponents of the successive primes in their prime factorization. (End)
From Bill McEachen, Oct 15 2022: (Start)
From inspection, the least significant decimal digits of a(n) terms form continuous chains of 30 as follows. For n == i (mod 30), i=0..5, there are 6 ordered elements of these 8 {1,2,3,6,9,8,7,4}. Then for n == i (mod 30), i=6..29, there are 12 repeated pairs = {5,0}.
Moreover, when the individual elements of any of the possible groups of 6 are transformed via (7*digit) (mod 10), the result matches one of the other 7 groupings (not all 7 may be seen). As example, {1,2,3,6,9,8} transforms to {7,4,1,2,3,6}. (End)
The least significant digit of a(n) in base 4 is given by A353486, and in base 6 by A358840. - Antti Karttunen, Oct 25 2022, Feb 17 2024

			For n = 24, which has primorial base representation (see A049345) "400" as 24 = 4*A002110(2) + 0*A002110(1) + 0*A002110(0) = 4*6 + 0*2 + 0*1, thus a(24) = prime(3)^4 * prime(2)^0 * prime(1)^0 = 5^4 = 625.
For n = 35 = "1021" as 35 = 1*A002110(3) + 0*A002110(2) + 2*A002110(1) + 1*A002110(0) = 1*30 + 0*6 + 2*2 + 1*1, thus a(35) = prime(4)^1 * prime(2)^2 * prime(1) = 7 * 3*3 * 2 = 126.

Antti Karttunen, Table of n, a(n) for n = 0..2310
Antti Karttunen, Program in LODA-assembly
Antti Karttunen, Program in LODA-assembly [Cached copy]
Index entries for sequences related to primorial base

Cf. A276085 (a left inverse) and also A276087, A328403.
Cf. A000040, A001221, A001222, A002110, A020639, A049345, A053669, A055396, A057588, A071178, A143293, A257993, A267263, A276084, A276088, A276092, A276093, A276147, A276150, A276151, A276153, A276156, A283477, A324198 (= gcd(n, a(n))), A328584 (= lcm(n, a(n))), A324646, A324289, A328386, A328403, A328475, A328571, A328572, A328578, A328612, A328613, A328620, A328624, A328627, A328763, A328766, A328828, A328835, A328841, A328842, A328843, A328844, A329041, A324580 [= n*a(n)], A324895 (largest proper divisor of a(n)), A351252, A353486 (reduced modulo 4), A358840 (modulo 6), A353489, A353516.
Cf. A048103 (terms sorted into ascending order), A100716 (natural numbers not present in this sequence).
Cf. A278226 (associated filter-sequence), A286626 (and its rgs-version), A328477.
Cf. A328316 (iterates started from zero).
Cf. A327858, A327859, A327860, A327963, A328097, A328098, A328099, A328110, A328112, A328382 for various combinations with arithmetic derivative (A003415).
Cf. also A327167, A329037.
Cf. A019565 and A054842 for base-2 and base-10 analogs and A276076 for the analogous "factorial base exp-function", from which this differs for the first time at n=24, where a(24)=625 while A276076(24)=7.
Cf. A327969, A351088, A351458 for sequences with conjectures involving this sequence.

b = MixedRadix[Reverse@ Prime@ Range@ 12]; Table[Function[k, Times @@ Power @@@ # &@ Transpose@ {Prime@ Range@ Length@ k, Reverse@ k}]@ IntegerDigits[n, b], {n, 0, 51}] (* Michael De Vlieger, Aug 23 2016, Version 10.2 *)
f[n_] := Block[{a = {{0, n}}}, Do[AppendTo[a, {First@ #, Last@ #} &@ QuotientRemainder[a[[-1, -1]], Times @@ Prime@ Range[# - i]]], {i, 0, #}] &@ NestWhile[# + 1 &, 0, Times @@ Prime@ Range[# + 1] <= n &]; Rest[a][[All, 1]]]; Table[Times @@ Flatten@ MapIndexed[Prime[#2]^#1 &, Reverse@ f@ n], {n, 0, 73}] (* Michael De Vlieger, Aug 30 2016, Pre-Version 10 *)
a[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];
Table[a[n], {n, 0, 100}] (* Jean-François Alcover, Dec 01 2021, after Antti Karttunen's Sage code *)

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; }; \\ Antti Karttunen, May 12 2017

A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); }; \\ (Better than above one, avoids unnecessary construction of primorials). - Antti Karttunen, Oct 14 2019

from sympy import prime
def a(n):
    i=0
    m=pr=1
    while n>0:
        i+=1
        N=prime(i)*pr
        if n%N!=0:
            m*=(prime(i)**((n%N)/pr))
            n-=n%N
        pr=N
    return m # Indranil Ghosh, May 12 2017, after Antti Karttunen's PARI code

from sympy import nextprime
def a(n):
    m, p = 1, 2
    while n > 0:
        n, r = divmod(n, p)
        m *= p**r
        p = nextprime(p)
    return m
print([a(n) for n in range(74)])  # Peter Luschny, Apr 20 2024

def A276086(n):
    m=1
    i=1
    while n>0:
        p = sloane.A000040(i)
        m *= (p**(n%p))
        n = floor(n/p)
        i += 1
    return (m)
# Antti Karttunen, Oct 14 2019, after Indranil Ghosh's Python code above, and my own leaner PARI code from Oct 14 2019. This avoids unnecessary construction of primorials.

(define (A276086 n) (let loop ((n n) (t 1) (i 1)) (if (zero? n) t (let* ((p (A000040 i)) (d (modulo n p))) (loop (/ (- n d) p) (* t (expt p d)) (+ 1 i))))))

(definec (A276086 n) (if (zero? n) 1 (* (expt (A053669 n) (A276088 n)) (A276086 (A276093 n))))) ;; Needs macro definec from http://oeis.org/wiki/Memoization#Scheme

(definec (A276086 n) (if (zero? n) 1 (* (A053669 n) (A276086 (- n (A002110 (A276084 n))))))) ;; Needs macro definec from http://oeis.org/wiki/Memoization#Scheme

a(0) = 1; for n >= 1, a(n) = A053669(n) * a(A276151(n)) = A053669(n) * a(n-A002110(A276084(n))).
a(0) = 1; for n >= 1, a(n) = A053669(n)^A276088(n) * a(A276093(n)).
a(n) = A328841(a(n)) + A328842(a(n)) = A328843(n) + A328844(n).
a(n) = a(A328841(n)) * a(A328842(n)) = A328571(n) * A328572(n).
a(n) = A328475(n) * A328580(n) = A328476(n) + A328580(n).
a(A002110(n)) = A000040(n+1). [Maps primorials to primes]
a(A143293(n)) = A002110(n+1). [Maps partial sums of primorials to primorials]
a(A057588(n)) = A276092(n).
a(A276156(n)) = A019565(n).
a(A283477(n)) = A324289(n).
a(A003415(n)) = A327859(n).
Here the text in brackets shows how the right hand side sequence is a function of the primorial base expansion of n:
A001221(a(n)) = A267263(n). [Number of nonzero digits]
A001222(a(n)) = A276150(n). [Sum of digits]
A067029(a(n)) = A276088(n). [The least significant nonzero digit]
A071178(a(n)) = A276153(n). [The most significant digit]
A061395(a(n)) = A235224(n). [Number of significant digits]
A051903(a(n)) = A328114(n). [Largest digit]
A055396(a(n)) = A257993(n). [Number of trailing zeros + 1]
A257993(a(n)) = A328570(n). [Index of the least significant zero digit]
A079067(a(n)) = A328620(n). [Number of nonleading zeros]
A056169(a(n)) = A328614(n). [Number of 1-digits]
A056170(a(n)) = A328615(n). [Number of digits larger than 1]
A277885(a(n)) = A328828(n). [Index of the least significant digit > 1]
A134193(a(n)) = A329028(n). [The least missing nonzero digit]
A005361(a(n)) = A328581(n). [Product of nonzero digits]
A072411(a(n)) = A328582(n). [LCM of nonzero digits]
A001055(a(n)) = A317836(n). [Number of carry-free partitions of n in primorial base]
Various number theoretical functions applied:
A000005(a(n)) = A324655(n). [Number of divisors of a(n)]
A000203(a(n)) = A324653(n). [Sum of divisors of a(n)]
A000010(a(n)) = A324650(n). [Euler phi applied to a(n)]
A023900(a(n)) = A328583(n). [Dirichlet inverse of Euler phi applied to a(n)]
A069359(a(n)) = A329029(n). [Sum a(n)/p over primes p dividing a(n)]
A003415(a(n)) = A327860(n). [Arithmetic derivative of a(n)]
Other identities:
A276085(a(n)) = n.          [A276085 is a left inverse]
A020639(a(n)) = A053669(n). [The smallest prime not dividing n -> the smallest prime dividing n]
A046523(a(n)) = A278226(n). [Least number with the same prime signature as a(n)]
A246277(a(n)) = A329038(n).
A181819(a(n)) = A328835(n).
A053669(a(n)) = A326810(n), A326810(a(n)) = A328579(n).
A257993(a(n)) = A328570(n), A328570(a(n)) = A328578(n).
A328613(a(n)) = A328763(n), A328620(a(n)) = A328766(n).
A328828(a(n)) = A328829(n).
A053589(a(n)) = A328580(n). [Greatest primorial number which divides a(n)]
A276151(a(n)) = A328476(n). [... and that primorial subtracted from a(n)]
A111701(a(n)) = A328475(n).
A328114(a(n)) = A328389(n). [Greatest digit of primorial base expansion of a(n)]
A328389(a(n)) = A328394(n), A328394(a(n)) = A328398(n).
A235224(a(n)) = A328404(n), A328405(a(n)) = A328406(n).
a(A328625(n)) = A328624(n), a(A328626(n)) = A328627(n). ["Twisted" variants]
a(A108951(n)) = A324886(n).
a(n) mod n = A328386(n).
a(a(n)) = A276087(n), a(a(a(n))) = A328403(n). [2- and 3-fold applications]
a(2n+1) = 2 * a(2n). - Antti Karttunen, Feb 17 2022

Name edited and new link-formulas added by Antti Karttunen, Oct 29 2019

Name changed again by Antti Karttunen, Feb 05 2022

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

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

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

A329045 Lexicographically earliest infinite sequence such that a(i) = a(j) => A046523(A329044(i)) = A046523(A329044(j)) for all i, j.

Original entry on oeis.org

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

Offset: 1

Antti Karttunen, Nov 08 2019

nonn

Restricted growth sequence transform of function f(n) = A046523(A329044(n)).
For all i, j:
  A305800(i) = A305800(j) => a(i) = a(j),
  a(i) = a(j) => A324888(i) = A324888(j),
  a(i) = a(j) => A329046(i) = A329046(j).

Cf. A034386, A064989, A108951, A276086, A305800, A324886, A324888, A329044, A329046.

Cf. also A278226, A286626.

up_to = 65537;
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
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));
A064989(n) = {my(f); f = factor(n); if((n>1 && f[1,1]==2), f[1,2] = 0); for (i=1, #f~, f[i,1] = precprime(f[i,1]-1)); factorback(f)};
A329044(n) = A064989(A324886(n));
A046523(n) = { my(f=vecsort(factor(n)[, 2], , 4), p); prod(i=1, #f, (p=nextprime(p+1))^f[i]); }; \\ From A046523
v329045 = rgs_transform(vector(up_to, n, A046523(A329044(n))));
A329045(n) = v329045[n];

cp's OEIS Frontend

A329622 a(n) = A056239(n) - A324888(n) = A001222(A108951(n)) - A001222(A324886(n)).

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A346100 a(n) = A100995(gcd(n, A064989(A319626(A324886(n))))).

Views

Author

Keywords

Crossrefs

Programs

PARI

Formula

A373983 Lexicographically earliest infinite sequence such that a(i) = a(j) = A246277(A324886(i)) = A246277(A324886(j)) and A278226(A328768(i)) = A278226(A328768(j)), for all i, j >= 1.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

A329346 a(n) = A322356(A324886(n)).

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

PARI

Formula

A276086 Primorial base exp-function: digits in primorial base representation of n become the exponents of successive prime factors whose product a(n) is.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Python

Python

Sage

Scheme

Scheme

Scheme

Formula

Extensions

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

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

Views

Author

Keywords

Comments

Links

Crossrefs