A324653 -id:A324653 - cp's OEIS frontend

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

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

A323243 a(1) = 0; for n > 1, a(n) = A000203(A156552(n)).

Original entry on oeis.org

0, 1, 3, 4, 7, 6, 15, 8, 12, 13, 31, 12, 63, 18, 18, 24, 127, 14, 255, 20, 39, 48, 511, 24, 28, 84, 24, 48, 1023, 32, 2047, 32, 54, 176, 42, 40, 4095, 258, 144, 56, 8191, 38, 16383, 68, 36, 800, 32767, 48, 60, 31, 252, 132, 65535, 30, 91, 72, 528, 1302, 131071, 44, 262143, 2736, 60, 104, 126, 96, 524287, 304, 774, 42, 1048575, 72, 2097151, 4356, 42

Offset: 1

Antti Karttunen, Jan 10 2019

nonn

Antti Karttunen, Table of n, a(n) for n = 1..10000 (based on Hans Havermann's factorization of A156552)
Index entries for sequences related to binary expansion of n
Index entries for sequences computed from indices in prime factorization
Index entries for sequences related to sigma(n)

Cf. A000203, A156552, A323244, A323247, A323248, A324118, A324543 (Möbius transform), A324396, A324823.

Cf. A323173, A324054, A324184, A324545 for other permutations of sigma, and also A324573, A324653.

Array[If[# == 0, 0, DivisorSigma[1, #]] &@ Floor@ Total@ Flatten@ MapIndexed[#1 2^(#2 - 1) &, Flatten[Table[2^(PrimePi@ #1 - 1), {#2}] & @@@ FactorInteger@ #]] &, 75] (* Michael De Vlieger, Apr 21 2019 *)

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)};
A156552(n) = if(1==n, 0, if(!(n%2), 1+(2*A156552(n/2)), 2*A156552(A064989(n))));
A323243(n) = if(1==n, 0, sigma(A156552(n)));

\\ For computing terms a(n), with n > ~4000 use Hans Havermann's factorization file https://oeis.org/A156552/a156552.txt
v156552sigs = readvec("a156552.txt"); \\ First read it in as a PARI-vector.
A323243(n) = if(n<=2,n-1,my(prsig=v156552sigs[n],ps=prsig[1],es=prsig[2]); prod(i=1,#ps,((ps[i]^(1+es[i]))-1)/(ps[i]-1))); \\ Then play sigma
\\ Antti Karttunen, Mar 15 2019

from sympy import divisor_sigma, primepi, factorint
def A323243(n): return divisor_sigma(sum((1< 1 else 0 # Chai Wah Wu, Mar 10 2023

a(1) = 0; for n > 1, a(n) = A000203(A156552(n)).
a(n) = 2*A156552(n) - A323244(n).
a(n) = A323247(n) - A323248(n).
From Antti Karttunen, Mar 12 2019: (Start)
a(A000040(n)) = A000225(n).
a(n) = Sum_{d|n} A324543(d).
For n > 1, a(2*A246277(n)) = A324118(n).
gcd(a(n), A156552(n)) = A324396(n).
A000035(a(n)) = A324823(n).
(End)

A327860 Arithmetic derivative of the primorial base exp-function: a(n) = A003415(A276086(n)).

Original entry on oeis.org

0, 1, 1, 5, 6, 21, 1, 7, 8, 31, 39, 123, 10, 45, 55, 185, 240, 705, 75, 275, 350, 1075, 1425, 3975, 500, 1625, 2125, 6125, 8250, 22125, 1, 9, 10, 41, 51, 165, 12, 59, 71, 247, 318, 951, 95, 365, 460, 1445, 1905, 5385, 650, 2175, 2825, 8275, 11100, 30075, 4125, 12625, 16750, 46625, 63375, 166125, 14, 77, 91, 329, 420

Offset: 0

Antti Karttunen, Sep 30 2019

Are there any other fixed points after 0, 1, 7, 8 and 2556? (A328110, see also A351087 and A351088).
Out of the 30030 initial terms, 19220 are multiples of 5. (See A327865).
Proof that a(n) is even if and only if n is a multiple of 4: Consider Charlie Neder's Feb 25 2019 comment in A235992. As A276086 is never a multiple of 4, and as it toggles the parity, we only need to know when A001222(A276086(n)) = A276150(n) is even. The condition for that is given in the latter sequence by David A. Corneth's Feb 27 2019 comment. From this it also follows that A166486 gives similarly the parity of terms of A342002, A351083 and A345000. See also comment in A327858. - Antti Karttunen, May 01 2022

			2556 has primorial base expansion [1,1,1,1,0,0] as 1*A002110(5) + 1*A002110(4) + 1*A002110(3) + 1*A002110(2) = 2310 + 210 + 30 + 6 = 2556. That in turn is converted by A276086 to 13^1 * 11^1 * 7^1 * 5^1 = 5005, whose arithmetic derivative is 5' * 1001 + 1001' * 5 = 1*1001 + 311*5 = 2556, thus 2556 is one of the rare fixed points (A328110) of this sequence.

Antti Karttunen, Table of n, a(n) for n = 0..2310
Antti Karttunen, Data supplement: n, a(n) computed for n = 0..30030
Index entries for sequences related to primorial base

Cf. A002110 (positions of 1's), A003415, A048103, A276086, A327858, A327859, A327865, A328110 (fixed points), A328233 (positions of primes), A328242 (positions of squarefree terms), A328388, A328392, A328571, A328572, A329031, A329032, A329041, A342002.
Cf. A345000, A351074, A351075, A351076, A351077, A351080, A351083, A351084, A351087 (numbers k such that a(k) is a multiple of k), A351088.
Coincides with A329029 on positions given by A276156.
Cf. A166486 (a(n) mod 2), A353630 (a(n) mod 4).
Cf. A267263, A276150, A324650, A324653, A324655 for omega, bigomega, phi, sigma and tau applied to A276086(n).
Cf. also A351950 (analogous sequence).

Block[{b = MixedRadix[Reverse@ Prime@ Range@ 12]}, Array[Function[k, If[# < 2, 0, # Total[#2/#1 & @@@ FactorInteger[#]] ] &@ Abs[Times @@ Power @@@ # &@ Transpose@{Prime@ Range@ Length@ k, Reverse@ k}]]@ IntegerDigits[#, b] &, 65, 0]] (* Michael De Vlieger, Mar 12 2021 *)

A003415(n) = {my(fac); if(n<1, 0, fac=factor(n); sum(i=1, matsize(fac)[1], n*fac[i, 2]/fac[i, 1]))}; \\ From A003415
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; };
A327860(n) = A003415(A276086(n));

A327860(n) = { my(s=0, m=1, p=2, e); while(n, e = (n%p); m *= (p^e); s += (e/p); n = n\p; p = nextprime(1+p)); (s*m); }; \\ (Standalone version) - Antti Karttunen, Nov 07 2019

a(n) = A003415(A276086(n)).
a(A002110(n)) = 1 for all n >= 0.
From Antti Karttunen, Nov 03 2019: (Start)
Whenever A329041(x,y) = 1, a(x + y) = A003415(A276086(x)*A276086(y)) = a(x)*A276086(y) + a(y)*A276086(x). For example, we have:
a(n) = a(A328841(n)+A328842(n)) = A329031(n)*A328572(n) + A329032(n)*A328571(n).
A051903(a(n)) = A328391(n).
A328114(a(n)) = A328392(n).
(End)
From Antti Karttunen, May 01 2022: (Start)
a(n) = A328572(n) * A342002(n).
For all n >= 0, A000035(a(n)) = A166486(n). [See comments]
(End)

Verbal description added to the definition by Antti Karttunen, May 01 2022

A324583 Numbers k such that k and A276086(k) are coprime, where A276086 is the primorial base exp-function.

Original entry on oeis.org

0, 1, 2, 4, 5, 6, 7, 8, 11, 12, 13, 14, 16, 17, 18, 19, 22, 23, 24, 26, 28, 29, 30, 31, 32, 34, 36, 37, 38, 41, 43, 44, 46, 47, 48, 52, 53, 54, 58, 59, 60, 61, 62, 64, 65, 66, 67, 68, 71, 72, 73, 74, 76, 78, 79, 82, 83, 86, 88, 89, 90, 92, 94, 95, 96, 97, 101, 102, 103, 104, 106, 107, 108, 109, 113, 114, 116, 118, 120, 121

Offset: 1

Antti Karttunen, Mar 10 2019

nonn

Numbers k for which A324198(k) = 1.
For terms k > 0 it holds that:
  A000005(A324580(k)) = A000005(k) * A324655(k),
  A000010(A324580(k)) = A000010(k) * A324650(k),
  A000203(A324580(k)) = A000203(k) * A324653(k),
and similarly for any multiplicative function.

Cf. A324584 (complement), A356162 (characteristic function).
Cf. A276086, A324580, A324650, A324653, A324655.
Some subsequences are: A055932 ⊃ A025487 ⊃ A002182, and also A002110.
Subsequence of A356316.
Positions of 1's in A324198, positions 0's in A351254, A356302 and A356303, positions of fixed points in A351250 and in A356309.
Cf. also A355821, A356311.

A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
A324198(n) = gcd(n,A276086(n));
for(n=0,oo,if(1==A324198(n),print1(n,", ")));

Initial 0 prepended by Antti Karttunen, Nov 03 2022

A324650 a(n) = A000010(A276086(n)).

Original entry on oeis.org

1, 1, 2, 2, 6, 6, 4, 4, 8, 8, 24, 24, 20, 20, 40, 40, 120, 120, 100, 100, 200, 200, 600, 600, 500, 500, 1000, 1000, 3000, 3000, 6, 6, 12, 12, 36, 36, 24, 24, 48, 48, 144, 144, 120, 120, 240, 240, 720, 720, 600, 600, 1200, 1200, 3600, 3600, 3000, 3000, 6000, 6000, 18000, 18000, 42, 42, 84, 84, 252, 252, 168, 168, 336, 336, 1008, 1008

Offset: 0

Antti Karttunen, Mar 10 2019

nonn

Terms are duplicated because phi(2*(2n+1)) = phi(2n+1) for all n >= 0.

Cf. A000010, A002110, A276086, A324651 (bisection).
Cf. also A267263, A276150, A324653, A324655 for omega, bigomega, sigma and tau analogs.
Cf. also A290077.

A324650(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)-1)*(prime(i)^(((n%nextpr)/pr)-1)); n-=(n%nextpr));pr=nextpr); (m); };

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; };
A324650(n) = eulerphi(A276086(n));

a(n) = A000010(A276086(n)).
a(2n+1) = a(2n) for all n >= 0.
For n >= 1, a(A002110(n-1)) = A000040(n)-1.

A324655 a(n) = A000005(A276086(n)).

Original entry on oeis.org

1, 2, 2, 4, 3, 6, 2, 4, 4, 8, 6, 12, 3, 6, 6, 12, 9, 18, 4, 8, 8, 16, 12, 24, 5, 10, 10, 20, 15, 30, 2, 4, 4, 8, 6, 12, 4, 8, 8, 16, 12, 24, 6, 12, 12, 24, 18, 36, 8, 16, 16, 32, 24, 48, 10, 20, 20, 40, 30, 60, 3, 6, 6, 12, 9, 18, 6, 12, 12, 24, 18, 36, 9, 18, 18, 36, 27, 54, 12, 24, 24, 48, 36, 72, 15, 30, 30, 60, 45, 90, 4, 8, 8

Offset: 0

Antti Karttunen, Mar 10 2019

nonn

Alternative construction: write n down in primorial base (as in A049345, taking care of not mangling digits larger than 9), increment all the digits by one, and multiply together to get a(n). a(0) = 1 either as an empty product, or as a product of any number of 1's. See examples.

			For n = 11, its primorial base representation is "121" as 11 = 1*A002110(2) + 2*A002110(1) + 1*A002110(0) = 1*6 + 2*2 + 1*1, thus a(11) = (1+1)*(2+1)*(1+1) = 12.
For n = 13, its primorial base representation is "201" as 13 = 2*6 + 0*2 + 1*1, thus a(13) = (2+1)*(0+1)*(1+1) = 6.

Cf. A000005, A002110 (positions of 2's), A049345, A276086.

Cf. also A267263, A276150, A324650, A324653 for omega, bigomega, phi and sigma analogs.

A324655(n) = { my(t=1,m); forprime(p=2, , if(!n, return(t)); m = n%p; t *= (1+m); n = (n-m)/p); };

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; };
A324655(n) = numdiv(A276086(n));

a(n) = A000005(A276086(n)).

a(A002110(n)) = 2.

A346471 a(n) = A344695(A276086(n)), where A344695(x) = gcd(psi(x), sigma(x)), and A276086 gives the prime product form of primorial base expansion of n.

Original entry on oeis.org

1, 3, 4, 12, 1, 3, 6, 18, 24, 72, 6, 18, 1, 3, 4, 12, 1, 3, 6, 18, 24, 72, 12, 36, 1, 3, 4, 12, 1, 3, 8, 24, 32, 96, 8, 24, 48, 144, 192, 576, 48, 144, 8, 24, 32, 96, 8, 24, 48, 144, 192, 576, 96, 288, 8, 24, 32, 96, 8, 24, 1, 3, 4, 12, 3, 9, 6, 18, 24, 72, 18, 54, 3, 9, 12, 36, 3, 9, 12, 36, 48, 144, 36, 108, 3, 9, 12, 36

Offset: 0

Antti Karttunen, Jul 21 2021

Cf. A000203, A001615, A276086, A324653, A344695, A346470.

Cf. also A346469.

A001615(n) = if(1==n,n, my(f=factor(n)); prod(i=1, #f~, f[i, 1]^f[i, 2] + f[i, 1]^(f[i, 2]-1))); \\ After code in A001615
A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
A344695(n) = gcd(sigma(n), A001615(n));
A346471(n) = A344695(A276086(n));

a(n) = A344695(A276086(n)) = gcd(A324653(n), A346470(n)).

A331733 a(n) = sigma(A225546(n)), where sigma is the sum of divisors.

Original entry on oeis.org

1, 3, 7, 4, 31, 15, 511, 12, 13, 63, 131071, 28, 8589934591, 1023, 127, 6, 36893488147419103231, 39, 680564733841876926926749214863536422911, 124, 2047, 262143, 231584178474632390847141970017375815706539969331281128078915168015826259279871, 60, 121, 17179869183, 91, 2044

Offset: 1

Antti Karttunen, Feb 02 2020

nonn

Cf. A000203, A048675, A225546, A331309, A331734, A331735, A331741.

Cf. A323243, A323173, A324054, A324184, A324545 for other permutations of sigma, and also A324573, A324653.

Array[If[# == 1, 1, DivisorSigma[1, #] &@ Apply[Times, Flatten@ Map[Function[{p, e}, Map[Prime[Log2@ # + 1]^(2^(PrimePi@ p - 1)) &, DeleteCases[NumberExpand[e, 2], 0]]] @@ # &, FactorInteger[#]]]] &, 28] (* Michael De Vlieger, Feb 08 2020 *)

A048675(n) = { my(f = factor(n)); sum(k=1, #f~, f[k, 2]*2^primepi(f[k, 1]))/2; };
A331733(n) = if(1==n,1,my(f=factor(n),u=#binary(vecmax(f[, 2])),prods=vector(u,x,1),m=1,e); for(i=1,u,for(k=1,#f~, if(bitand(f[k,2],m),prods[i] *= f[k,1])); m<<=1); prod(i=1,u,(prime(i)^(1+A048675(prods[i]))-1)/(prime(i)-1)));

a(n) = A000203(A225546(n)).

For all n >= 1, A000035(a(A016754(n))) = 1. [Result is odd for all odd squares]

A324654 a(n) = A033879(A276086(n)).

Original entry on oeis.org

1, 1, 2, 0, 5, -3, 4, 2, 6, -12, 12, -54, 19, 7, 26, -72, 47, -309, 94, 32, 126, -372, 222, -1584, 469, 157, 626, -1872, 1097, -7959, 6, 4, 10, -12, 22, -60, 22, -4, 18, -156, 6, -612, 102, -44, 58, -876, -74, -3372, 502, -244, 258, -4476, -474, -17172, 2502, -1244, 1258, -22476, -2474, -86172, 41, 25, 66, -96, 141, -459, 148, -46, 102

Offset: 0

Antti Karttunen, Mar 10 2019

sign

Interestingly, this kind of sampling of deficiency (A033879; recall that the range of A276086 does not cover the whole N) biases it strongly towards negative values: of the first 2310 terms, 1565 are negative (~ 68%) and of the first 30030 terms, 22507 are negative (~ 75%). Compare also to A323174, which covers whole N.

Cf. A033879, A276086, A324653.

Cf. also A323174, A324055.

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; };
A033879(n) = (2*n-sigma(n));
A324654(n) = A033879(A276086(n));

a(n) = A033879(A276086(n)).

a(n) = 2*A276086(n) - A324653(n).

A346470 a(n) = psi(A276086(n)), where psi is Dedekind psi function A001615, and A276086 is the prime product form of primorial base expansion of n.

Original entry on oeis.org

1, 3, 4, 12, 12, 36, 6, 18, 24, 72, 72, 216, 30, 90, 120, 360, 360, 1080, 150, 450, 600, 1800, 1800, 5400, 750, 2250, 3000, 9000, 9000, 27000, 8, 24, 32, 96, 96, 288, 48, 144, 192, 576, 576, 1728, 240, 720, 960, 2880, 2880, 8640, 1200, 3600, 4800, 14400, 14400, 43200, 6000, 18000, 24000, 72000, 72000, 216000, 56, 168, 224, 672

Offset: 0

Antti Karttunen, Jul 21 2021

Antti Karttunen, Table of n, a(n) for n = 0..2310
Index entries for sequences related to primorial base

Cf. A001615.
Other number-theoretical functions similarly applied to A276086: A267263 (omega), A276150 (bigomega), A324650 (phi), A324653 (sigma), A324655 (tau), A327860 (arithmetic derivative).
Cf. also A346471, A346475.

A346470(n) = { my(m=1, p=2, e); while(n, e = (n%p); if(e, m *= (p+1)*(p^(e-1))); n = n\p; p = nextprime(1+p)); (m); };

a(n) = A001615(A276086(n)).

cp's OEIS Frontend

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

A323243 a(1) = 0; for n > 1, a(n) = A000203(A156552(n)).

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Python

Formula

A327860 Arithmetic derivative of the primorial base exp-function: a(n) = A003415(A276086(n)).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Formula

Extensions

A324583 Numbers k such that k and A276086(k) are coprime, where A276086 is the primorial base exp-function.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

Extensions

A324650 a(n) = A000010(A276086(n)).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

PARI

Formula

A324655 a(n) = A000005(A276086(n)).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

PARI

PARI

Formula