OFFSET
1,4
COMMENTS
An equivalent definition of this sequence is: a(1) = 0; and for n > 1, a(n) = n*w(n), where if p1, ..., pk are the prime divisors of n (with repetition) and p1 <= p2 <= ... <= pk, then w(n) = 1/pk + 1/(pk-1*pk) + ... + 1/(p2*p3*...*pk) + ... + 1/(p1*p2*...*pk). Since 2 is smallest prime w(n) <= 1/2 + 1/(2^2) + ... + 1/(2^k), a partial sum of a series which -->1 as n-->oo. Therefore w(n) < 1 and n-a(n) is a sequence of positive numbers (1,1,2,1,4,3,6,...). For n=p^k, p a prime and n >= 1, a(n) = a(p^k) = p^(n-1) + p^(n-2) + ... + p^2 + p + 1 = (p^k-1)/(p-1); e.g., a(2^k) = 2^k - 1.
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..16383
Antti Karttunen, Data supplement: n, a(n) computed for n = 1..65537
FORMULA
From Antti Karttunen, Feb 28 2019 , Mar 04 2019: (Start)
If n is of the form p^k, k >= 1, then a(n) = A000203(A003557(n)). [Based on author's comments above] (End)
a(n) = Sum_{k=1..bigomega(n)} F^k(n), where F^k(n) is the k-th iterate of F(n) = A052126(n). - Ridouane Oudra, Aug 17 2024
EXAMPLE
For any prime p, a(p) = p*a(p/p)+1 = p*a(1)+1 = 1, because a(1) = 0.
For n = 6, the least prime divisor is 2, so a(6) = 2*a(6/2)+1 = 2*a(3)+1 = 3.
Using the equivalent definition we get w(6) = 1/3 + 1/6 = 1/2, so a(6) = 6*w(6) = 6*1/2 = 3. For n=32, a(32) = a(2^5) = 2^5 - 1 = 32 - 1 = 31.
PROG
(PARI) a(n) = if (n==1, 0, my(p = vecmin(factor(n)[, 1])); p*a(n/p)+1); \\ Michel Marcus, Jan 25 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
David James Sycamore, Dec 05 2018
EXTENSIONS
More terms from Michel Marcus, Jan 25 2019
STATUS
approved