login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A349824 a(0) = 0; for n >= 1, a(n) = (number of primes, counted with repetition) * (sum of primes, counted with repetition). 5

%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

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 18:04 EDT 2024. Contains 371254 sequences. (Running on oeis4.)