login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A323136 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
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, 169, 27, 49, 29, 49, 31, 32, 121, 289, 49, 169, 37, 361, 169, 121, 41, 49, 43, 121, 121, 529, 47, 289, 49, 121, 289, 169, 53, 361, 121, 121, 361, 841, 59, 169, 61, 961 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

COMMENTS

The powers of primes correspond to A000961.

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

LINKS

Rémy Sigrist, Table of n, a(n) for n = 1..10000

FORMULA

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

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

EXAMPLE

For n = 12:

- 12 = 2^2 * 3,

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

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

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

- hence a(12) = 25.

PROG

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

CROSSREFS

See A323135 for the corresponding number of iterations.

Cf. A000961, A006530, A070229.

Sequence in context: A079871 A273291 A130064 * A273285 A068795 A222257

Adjacent sequences:  A323133 A323134 A323135 * A323137 A323138 A323139

KEYWORD

nonn

AUTHOR

Rémy Sigrist, Jan 05 2019

STATUS

approved

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified June 25 10:29 EDT 2019. Contains 324351 sequences. (Running on oeis4.)