OFFSET
1,3
COMMENTS
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
LINKS
FORMULA
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.
From Antti Karttunen, Apr 24 & Apr 29 2022: (Start)
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:
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:
The sum or difference of the rhs-sequences is A108951:
Here the two sequences are inverse permutations of each other:
Other correspondences:
a(A276076(n)) = A351576(n). [Sequence reinterpreting factorial base representation as a primorial base representation]
(End)
MATHEMATICA
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 *)
PROG
(Scheme, with memoization-macro definec)
(definec (A276085 n) (cond ((= 1 n) (- n 1)) (else (+ (* (A067029 n) (A002110 (+ -1 (A055396 n)))) (A276085 (A028234 n))))))
(PARI) A276085(n) = { my(f=factor(n)); sum(k=1, #f~, f[k, 2]*prod(i=1, primepi(f[k, 1]-1), prime(i))); }; \\ Antti Karttunen, Mar 15 2021, Jun 23 2024
(Python)
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
CROSSREFS
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.
Positions of multiples of k in this sequence, for k=2, 3, 4, 5, 8: A003159, A339746, A369002, A373140, A373138.
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)].
KEYWORD
nonn
AUTHOR
Antti Karttunen, Aug 21 2016
EXTENSIONS
Name amended by Antti Karttunen, Apr 24 2022
Name simplified, the old name moved to the comments - Antti Karttunen, Jun 23 2024
STATUS
approved