%I #39 Jan 18 2022 05:54:26
%S 0,0,2,3,8,5,10,7,18,12,14,11,21,13,18,16,32,17,24,19,27,20,26,23,36,
%T 20,30,27,33,29,30,31,50,28,38,24,40,37,42,32,44,41,36,43,45,33,50,47,
%U 55,28,36,40,51,53,44,32,52,44,62,59,48,61,66,39,72,36,48,67,63,52,42,71,60,73
%N a(0) = 0; for n >= 1, a(n) = (number of primes, counted with repetition) * (sum of primes, counted with repetition).
%C More precisely, a(n) = (number of prime factors of n, counted with repetition) * (sum of primes factors of n, counted with repetition): a(n) = A001414(n) * A001222(n).
%C Suggested by Mike Klein in an email, Dec 31 2021.
%C Conjecture (Mike Klein): Iterating n -> a(n) eventually leads to one of the fixed points {primes union 0, 27, 30} or the loop (28, 33).
%C It appears that with the exception of {1,4,16,27,30}, n divides a(n) iff n is prime. - _Gary Detlefs_, Jan 11 2022
%C From _Gary Detlefs_, Jan 14 2022: (Start)
%C The loop (28,33) referenced above would more correctly be denoted by (33,28). The only value of n which reaches 28 before 33 is 49.
%C Values of n for which the trajectory terminates at 27 are {9,12,20,21,25}. (End)
%H Rémy Sigrist, <a href="/A349824/b349824.txt">Table of n, a(n) for n = 0..10000</a>
%e If n = 27 = 3^3, a(n) = 3*(3+3+3) = 27.
%e If we start with n = 4, iterating this map produces the trajectory 4, 8, 18, 24, 36, 40, 44, 45, 33, 28, 33, 28, 33, 28, ...
%e If we start with n = 6, iterating this map produces the trajectory 6, 10, 14, 18, 24, 36, 40, 44, 45, 33, 28, 33, 28, 33, 28, ...
%t {0, 0}~Join~Array[Total[#]*Length[#] &@ Flatten[ConstantArray[#1, #2] & @@@ FactorInteger[#]] &, 72, 2] (* _Michael De Vlieger_, Jan 02 2022 *)
%o (PARI) a(n) = { if (n==0, 0, my (f=factor(n)); bigomega(f)*sum(k=1, #f~, f[k,1]*f[k,2])) } \\ _Rémy Sigrist_, Jan 01 2022
%o (Python)
%o from sympy import factorint
%o def a(n):
%o if n == 0: return 0
%o f = factorint(n)
%o return sum(f.values()) * sum(p*e for p, e in f.items())
%o print([a(n) for n in range(74)]) # _Michael S. Branicky_, Jan 02 2022
%Y Cf. A001222, A001414.
%K nonn
%O 0,3
%A _N. J. A. Sloane_, Jan 01 2022
|