OFFSET
1,2
COMMENTS
In other words, a(n) is constructed by incrementally shifting the elements of the ordered exponent multiset {e_1,e_2,...,e_k} to the right one place at a time, a total of k = A001221(n) times, and adding the results of these exponents applied to the distinct size ordered primes p_1,p_2,...,p_k which divide n. The first shift gives {e_k,e_1,e_2,...,e_(k-1)}, then {e_(k-1),e_k,...,e_(k-2)} ... and the last yields {e_1,e_2,...,e_k} (which is n itself). At each shift the last exponent of the previous shift becomes the first of the next. There are A001221(n) summands in the computation of a(n), each having the same rad (A007947) and omega (A001221) values as n.
a(n) >= n with equality for n > 1 in A000961.
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..10000
Michael De Vlieger, Log log scatterplot of a(n), n = 1..2^20.
FORMULA
a(n) = n for n > 1 in A000961.
Let b(0) = n and let b(n+1) = A105119(b(n)) for n >= 0 and let omega(n) be the number of distinct prime factors of n. Then a(n) = Sum_{i = 1..omega(n)} b(i). - David A. Corneth, Nov 07 2023
EXAMPLE
a(1) = 0, the empty sum.
a(6) = a(2^1*3^1) = 6 + 2^1*3^1 = 6 + 6 = 12.
a(12) = a(2^2*3^1) = 12 + 2^1*3^2 = 12 + 18 = 30.
a(2250) = a(2^1*3^2*5^3) = 2^3*3^1*5^2 + 2^2*3^3*5^1 + 2^1*3^2*5^3 = 3390.
MATHEMATICA
{0}~Join~Array[Total@ Table[Times @@ Power @@@ Transpose@ {#1, RotateRight[#2, k]}, {k, PrimeNu[#]}] & @@ Transpose@ FactorInteger[#] &, 64, 2] (* Michael De Vlieger, Nov 09 2023 *)
PROG
(PARI)
a(n) = {
my(f = factor(n), res = 0);
for(i = 1, #f~,
res+=prod(j = 1, #f~, f[j, 1]^f[(i+j-1)%#f~ + 1, 2])
);
res
} \\ David A. Corneth, Nov 07 2023
CROSSREFS
KEYWORD
nonn
AUTHOR
David James Sycamore and David A. Corneth Nov 07 2023
EXTENSIONS
More terms and revised name from David A. Corneth, Nov 07 2023
STATUS
approved