A328841 -id:A328841 - 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

A276085 Primorial base log-function: fully additive with a(p) = p#/p, where p# = A034386(p).

Original entry on oeis.org

0, 1, 2, 2, 6, 3, 30, 3, 4, 7, 210, 4, 2310, 31, 8, 4, 30030, 5, 510510, 8, 32, 211, 9699690, 5, 12, 2311, 6, 32, 223092870, 9, 6469693230, 5, 212, 30031, 36, 6, 200560490130, 510511, 2312, 9, 7420738134810, 33, 304250263527210, 212, 10, 9699691, 13082761331670030, 6, 60, 13, 30032, 2312, 614889782588491410, 7, 216, 33, 510512, 223092871, 32589158477190044730, 10

Offset: 1

Antti Karttunen, Aug 21 2016

nonn

Completely additive with a(p^e) = e * A002110(A000720(p)-1).
This is a left inverse of A276086 ("primorial base exp-function"), hence the name "primorial base log-function". When the domain is restricted to the terms of A048103, this works also as a right inverse, as A276086(a(A048103(n))) = A048103(n) for all n >= 1. - Antti Karttunen, Apr 24 2022
On average, every third term is a multiple of 4. See A369001. - Antti Karttunen, May 26 2024

A left inverse of A276086.
Cf. A000040, A000720, A002110, A028234, A034386, A048103, A049345, A055396, A067029, A108951, A143293, A276154, A328316, A328624, A328625, A328768, A328832, A346105, A351576, A376398 (partial sums).
Positions of multiples of k in this sequence, for k=2, 3, 4, 5, 8, 27, 3125: A003159, A339746, A369002, A373140, A373138, A377872, A377878.
Cf. A036554 (positions of odd terms), A035263, A096268 (parity of terms).
Cf. A372575 (rgs-transform), A372576 [a(n) mod 360], A373842 [= A003415(a(n))].
Cf. A373145 [= gcd(A003415(n), a(n))], A373361 [= gcd(n, a(n))], A373362 [= gcd(A001414(n), a(n))], A373485 [= gcd(A083345(n), a(n))], A373835 [= gcd(bigomega(n), a(n))], and also A373367 and A373147 [= A003415(n) mod a(n)], A373148 [= a(n) mod A003415(n)].
Other completely additive sequences with primes p mapped to a function of p include: A001222 (with a(p)=1), A001414 (with a(p)=p), A059975 (with a(p)=p-1), A341885 (with a(p)=p*(p+1)/2), A373149 (with a(p)=prevprime(p)), A373158 (with a(p)=p#).
Cf. also A276075 for factorial base and A048675, A054841 for base-2 and base-10 analogs.

nn = 60; b = MixedRadix[Reverse@ Prime@ Range@ PrimePi[nn + 1]]; Table[FromDigits[#, b] &@ Reverse@ If[n == 1, {0}, Function[k, ReplacePart[Table[0, {PrimePi[k[[-1, 1]]]}], #] &@ Map[PrimePi@ First@ # -> Last@ # &, k]]@ FactorInteger@ n], {n, nn}] (* Version 10.2, or *)
f[w_List] := Total[Times @@@ Transpose@ {Map[Times @@ # &, Prime@ Range@ Range[0, Length@ w - 1]], Reverse@ w}]; Table[f@ Reverse@ If[n == 1, {0}, Function[k, ReplacePart[Table[0, {PrimePi[k[[-1, 1]]]}], #] &@ Map[PrimePi@ First@ # -> Last@ # &, k]]@ FactorInteger@ n], {n, 60}] (* Michael De Vlieger, Aug 30 2016 *)

A276085(n) = { my(f = factor(n), pr=1, i=1, s=0); for(k=1, #f~, while(i <= primepi(f[k, 1])-1, pr *= prime(i); i++); s += f[k, 2]*pr); (s); }; \\ Antti Karttunen, Nov 11 2024

from sympy import primorial, primepi, factorint
def a002110(n):
    return 1 if n<1 else primorial(n)
def a(n):
    f=factorint(n)
    return sum(f[i]*a002110(primepi(i) - 1) for i in f)
print([a(n) for n in range(1, 101)]) # Indranil Ghosh, Jun 22 2017

a(1) = 0; for n > 1, a(n) = a(A028234(n)) + (A067029(n) * A002110(A055396(n)-1)).
a(1) = 0, a(n) = (e1*A002110(i1-1) + ... + ez*A002110(iz-1)) when n = prime(i1)^e1 * ... * prime(iz)^ez.
Other identities.
For all n >= 0:
a(A276086(n)) = n.
a(A000040(1+n)) = A002110(n).
a(A002110(1+n)) = A143293(n).
From Antti Karttunen, Apr 24 & Apr 29 2022: (Start)
a(A283477(n)) = A283985(n).
a(A108951(n)) = A346105(n). [The latter has a similar additive formula as this sequence, but instead of primorials, uses their partial sums]
When applied to sequences where a certain subset of the divisors of n has been multiplicatively encoded with the help of A276086, this yields a corresponding number-theoretical sequence, i.e. completes their computation:
a(A319708(n)) = A001065(n) and a(A353564(n)) = A051953(n).
a(A329350(n)) = A069359(n) and a(A329380(n)) = A323599(n).
In the following group, the sum of the rhs-sequences is n [on each row, as say, A328841(n)+A328842(n)=n], because the pointwise product of the corresponding lhs-sequences is A276086:
a(A053669(n)) = A053589(n) and a(A324895(n)) = A276151(n).
a(A328571(n)) = A328841(n) and a(A328572(n)) = A328842(n).
a(A351231(n)) = A351233(n) and a(A327858(n)) = A351234(n).
a(A351251(n)) = A351253(n) and a(A324198(n)) = A351254(n).
The sum or difference of the rhs-sequences is A108951:
a(A344592(n)) = A346092(n) and a(A346091(n)) = A346093(n).
a(A346106(n)) = A346108(n) and a(A346107(n)) = A346109(n).
Here the two sequences are inverse permutations of each other:
a(A328624(n)) = A328625(n) and a(A328627(n)) = A328626(n).
a(A346102(n)) = A328622(n) and a(A346233(n)) = A328623(n).
a(A346101(n)) = A289234(n). [Self-inverse]
Other correspondences:
a(A324350(x,y)) = A324351(x,y).
a(A003961(A276086(n))) = A276154(n). [The primorial base left shift]
a(A276076(n)) = A351576(n). [Sequence reinterpreting factorial base representation as a primorial base representation]
(End)

Name amended by Antti Karttunen, Apr 24 2022

Name simplified, the old name moved to the comments - Antti Karttunen, Jun 23 2024

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

A276156 Numbers obtained by reinterpreting base-2 representation of n in primorial base: a(0) = 0, a(2n) = A276154(a(n)), a(2n+1) = 1 + A276154(a(n)).

Original entry on oeis.org

0, 1, 2, 3, 6, 7, 8, 9, 30, 31, 32, 33, 36, 37, 38, 39, 210, 211, 212, 213, 216, 217, 218, 219, 240, 241, 242, 243, 246, 247, 248, 249, 2310, 2311, 2312, 2313, 2316, 2317, 2318, 2319, 2340, 2341, 2342, 2343, 2346, 2347, 2348, 2349, 2520, 2521, 2522, 2523, 2526, 2527, 2528, 2529, 2550, 2551, 2552, 2553, 2556, 2557, 2558, 2559, 30030, 30031

Offset: 0

Antti Karttunen, Aug 24 2016

Numbers that are sums of distinct primorial numbers, A002110.

Numbers with no digits larger than one in primorial base, A049345.

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

Cf. A000040, A001511, A002110, A007088, A007814, A019565, A049345, A257993, A276084, A276085, A276154, A351073, A328461, A328473, A328474, A328571, A328831, A328836, A341518, A344591.
Complement of A177711.
Subsequences: A328233, A328832, A328462 (odd bisection).
Conjectured subsequences: A328110, A380527.
Fixed points of A328841, positions of zeros in A328828, A328842, and A329032, positions of ones in A328581, A328582, and A381032.
Positions of terms < 2 in A328114.
Indices where A327860 and A329029 coincide.
Cf. also table A328464 (and its rows).
Cf. also A059590, A283985, A290249, A342921.

nn = 65; b = MixedRadix[Reverse@ Prime@ Range[IntegerLength[nn, 2] - 1]]; Table[FromDigits[IntegerDigits[n, 2], b], {n, 0, 65}] (* Version 10.2, or *)
Table[Total[Times @@@ Transpose@ {Map[Times @@ # &, Prime@ Range@ Range[0, Length@ # - 1]], Reverse@ #}] &@ IntegerDigits[n, 2], {n, 0, 65}] (* Michael De Vlieger, Aug 26 2016 *)

A276156(n) = { my(s=0, p=1, r=1); while(n, if(n%2, s += r); n>>=1; p = nextprime(1+p); r *= p); (s); }; \\ Antti Karttunen, Feb 03 2022

from sympy import prime, primorial, primepi, factorint
from operator import mul
def a002110(n): return 1 if n<1 else primorial(n)
def a276085(n):
    f=factorint(n)
    return sum([f[i]*a002110(primepi(i) - 1) for i in f])
def a019565(n): return reduce(mul, (prime(i+1) for i, v in enumerate(bin(n)[:1:-1]) if v == '1')) # after Chai Wah Wu
def a(n): return 0 if n==0 else a276085(a019565(n))
print([a(n) for n in range(101)]) # Indranil Ghosh, Jun 23 2017

a(0) = 0, a(2n) = A276154(a(n)), a(2n+1) = 1+A276154(a(n)).
Other identities. For all n >= 0:
a(n) = A276085(A019565(n)).
A049345(a(n)) = A007088(n).
A257993(a(n)) = A001511(n).
A276084(a(n)) = A007814(n).
A051903(a(n)) = A351073(n).

A329041 Square array read by antidiagonals: A(n, k) = A327936(A276086(n) * A276086(k)).

Original entry on oeis.org

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

Offset: 0

Antti Karttunen, Nov 03 2019

Array A(n, k) is symmetric, and is read as (n,k) = (0, 0), (1, 0), (0, 1), (2, 0), (1, 1), (0, 2), (3, 0), (2, 1), (1, 2), (0, 3), ...

If A(n, k) is 1, it tells that adding of n and k do not generate any carries, when done in primorial base (A049345). If A(n, k) is larger than one, then its prime factors indicate in which specific moduli (digit positions) the sum was larger than allowed for that position.

			The top left corner of the array:
        0  1  2  3  4  5  6  7  8  9 10 11 12
      +--------------------------------------
   0: | 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, ...
   1: | 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, ...
   2: | 1, 1, 1, 1, 3, 3, 1, 1, 1, 1, 3, 3, 1, ...
   3: | 1, 2, 1, 2, 3, 6, 1, 2, 1, 2, 3, 6, 1, ...
   4: | 1, 1, 3, 3, 3, 3, 1, 1, 3, 3, 3, 3, 1, ...
   5: | 1, 2, 3, 6, 3, 6, 1, 2, 3, 6, 3, 6, 1, ...
   6: | 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, ...
   7: | 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, ...
   8: | 1, 1, 1, 1, 3, 3, 1, 1, 1, 1, 3, 3, 1, ...
   9: | 1, 2, 1, 2, 3, 6, 1, 2, 1, 2, 3, 6, 1, ...
  10: | 1, 1, 3, 3, 3, 3, 1, 1, 3, 3, 3, 3, 1, ...
  11: | 1, 2, 3, 6, 3, 6, 1, 2, 3, 6, 3, 6, 1, ...
  12: | 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, ...
...
A(11,25) = A(25,11) = 10 because 11 is written in primorial base representation (A049345) as "121" and 25 as "401", and when these are added together digit by digit, we see that the maximal allowed digits "421" for the rightmost three positions are exceeded in positions 1 and 3, with the 1st and 3rd primes 2 and 5 as their moduli, thus A(11,25) = 2*5 = 10.

Cf. A049345, A276086, A327936.

Cf. also A317836, A324351, A328770, A328841, A328842.

up_to = 105;
A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
A327936(n) = { my(f = factor(n)); for(k=1, #f~, f[k,2] = (f[k,2]>=f[k,1])); factorback(f); };
A329041sq(row,col) = A327936(A276086(row)*A276086(col));
A329041list(up_to) = { my(v = vector(up_to), i=0); for(a=0,oo, for(col=0,a, if(i++ > up_to, return(v)); v[i] = A329041sq(a-col,col))); (v); };
v329041 = A329041list(up_to);
A329041(n) = v329041[1+n];

A(n, k) = A327936(A276086(n) * A276086(k)).

For all n, A(A328841(n), A328842(n)) = 1 and A(A328770(n), A328770(n)) = 1.

A328842 Decrement each nonzero digit by one in primorial base representation of n, then convert back to decimal.

Original entry on oeis.org

0, 0, 0, 0, 2, 2, 0, 0, 0, 0, 2, 2, 6, 6, 6, 6, 8, 8, 12, 12, 12, 12, 14, 14, 18, 18, 18, 18, 20, 20, 0, 0, 0, 0, 2, 2, 0, 0, 0, 0, 2, 2, 6, 6, 6, 6, 8, 8, 12, 12, 12, 12, 14, 14, 18, 18, 18, 18, 20, 20, 30, 30, 30, 30, 32, 32, 30, 30, 30, 30, 32, 32, 36, 36, 36, 36, 38, 38, 42, 42, 42, 42, 44, 44, 48, 48, 48, 48, 50, 50, 60, 60

Offset: 0

Antti Karttunen, Oct 30 2019

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

Cf. A002110, A003557, A276085, A276086, A328114, A328572, A328841, A328844.
Cf. A276156 (positions of zeros).
Cf. A276009 for analogous sequence.

A328842(n) = { my(p=2, r=1, s=0); while(n, if((n%p)>0, s += ((n%p)-1)*r); r *= p; n = n\p; p = nextprime(1+p)); (s); };

a(n) = n - A328841(n).
a(n) = A276085(A328572(n)) = A276085(A003557(A276086(n))).
For all n >= 0, a(A276086(n)) = A328844(n).
For all n >= 1, A328114(a(n)) = A328114(n) - 1.

A343046 Array T(n, k), n, k >= 0, read by antidiagonals; lunar multiplication table for the primorial base.

Original entry on oeis.org

0, 0, 0, 0, 1, 0, 0, 2, 2, 0, 0, 3, 6, 3, 0, 0, 2, 8, 8, 2, 0, 0, 3, 6, 9, 6, 3, 0, 0, 6, 8, 8, 8, 8, 6, 0, 0, 7, 30, 9, 12, 9, 30, 7, 0, 0, 8, 32, 36, 14, 14, 36, 32, 8, 0, 0, 9, 36, 39, 30, 15, 30, 39, 36, 9, 0, 0, 8, 38, 38, 32, 36, 36, 32, 38, 38, 8, 0

Offset: 0

Rémy Sigrist, Apr 05 2021

To compute T(n, k):
- write the primorial base representations of n and of k on two lines, right aligned,
- to "multiply" two digits: take the smallest,
- to "add" two digits: take the largest,
- for example, for T(9, 10):
      9   ->   1 1 1
     10   -> x 1 2 0
             -------
               0 0 0
             1 1 1
         + 1 1 1
         -----------
           1 1 1 1 0  ->  248 = T(9, 10)
See A343044 for the corresponding addition table.

			Array T(n, k) begins:
  n\k|  0  1   2   3   4   5    6    7    8    9   10   11   12
  ---+---------------------------------------------------------
    0|  0  0   0   0   0   0    0    0    0    0    0    0    0
    1|  0  1   2   3   2   3    6    7    8    9    8    9    6  ->  A328841
    2|  0  2   6   8   6   8   30   32   36   38   36   38   30
    3|  0  3   8   9   8   9   36   39   38   39   38   39   36
    4|  0  2   6   8  12  14   30   32   36   38   42   44   60
    5|  0  3   8   9  14  15   36   39   38   39   44   45   66
    6|  0  6  30  36  30  36  210  216  240  246  240  246  210
    7|  0  7  32  39  32  39  216  217  248  249  248  249  216
    8|  0  8  36  38  36  38  240  248  246  248  246  248  240
    9|  0  9  38  39  38  39  246  249  248  249  248  249  246
   10|  0  8  36  38  42  44  240  248  246  248  252  254  270
   11|  0  9  38  39  44  45  246  249  248  249  254  255  276
   12|  0  6  30  36  60  66  210  216  240  246  270  276  420

Rémy Sigrist, Table of n, a(n) for n = 0..10010
Rémy Sigrist, PARI program for A343046
Index entries for sequences related to dismal (or lunar) arithmetic
Index entries for sequences related to primorial base

Cf. A087062, A235168, A328841, A343042, A343044, A343047.

PARI
```
See Links section.
```

T(n, k) = T(k, n).
T(m, T(n, k)) = T(T(m, n), k).
T(n, 0) = 0.
T(n, 1) = A328841(n).
T(n, n) = A343047(n).

cp's OEIS Frontend

A329031 a(n) = A327860(A328841(n)).

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

Formula

A328843 a(n) = A328841(A276086(n)).

Views

Author

Keywords

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

A276085 Primorial base log-function: fully additive with a(p) = p#/p, where p# = A034386(p).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

Extensions

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

A276156 Numbers obtained by reinterpreting base-2 representation of n in primorial base: a(0) = 0, a(2n) = A276154(a(n)), a(2n+1) = 1 + A276154(a(n)).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

A329041 Square array read by antidiagonals: A(n, k) = A327936(A276086(n) * A276086(k)).

Views