login
a(n) is the least power of a prime that n reaches under the map x -> A070229(x), or -1 if we never reach a power of a prime.
2

%I #13 Jan 12 2019 10:29:23

%S 1,2,3,4,5,9,7,8,9,25,11,25,13,49,25,16,17,49,19,25,49,121,23,27,25,

%T 169,27,49,29,49,31,32,121,289,49,169,37,361,169,121,41,49,43,121,121,

%U 529,47,289,49,121,289,169,53,361,121,121,361,841,59,169,61,961

%N a(n) is the least power of a prime that n reaches under the map x -> A070229(x), or -1 if we never reach a power of a prime.

%C The powers of primes correspond to A000961.

%C Apparently, a(n) > 0 for all n > 0.

%H Rémy Sigrist, <a href="/A323136/b323136.txt">Table of n, a(n) for n = 1..10000</a>

%F a(n) = n iff n belongs to A000961.

%F A006530(a(n)) >= A006530(n) when a(n) > 0.

%e For n = 12:

%e - 12 = 2^2 * 3,

%e - A070229(12) = 12 + 3 = 15 = 3 * 5,

%e - A070229(15) = 15 + 5 = 20 = 2^2 * 5,

%e - A070229(20) = 20 + 5 = 25 = 5^5,

%e - hence a(12) = 25.

%o (PARI) a(n) = for (k=0, oo, if (omega(n) <= 1, return (n), my (f=factor(n)); n += f[#f~, 1]))

%Y See A323135 for the corresponding number of iterations.

%Y Cf. A000961, A006530, A070229.

%K nonn

%O 1,2

%A _Rémy Sigrist_, Jan 05 2019