A338668 a(n) is the rightmost prime number in prime tower factorization of n; a(1) = 1.

1, 2, 3, 2, 5, 3, 7, 3, 2, 5, 11, 3, 13, 7, 5, 2, 17, 2, 19, 5, 7, 11, 23, 3, 2, 13, 3, 7, 29, 5, 31, 5, 11, 17, 7, 2, 37, 19, 13, 5, 41, 7, 43, 11, 5, 23, 47, 3, 2, 2, 17, 13, 53, 3, 11, 7, 19, 29, 59, 5, 61, 31, 7, 3, 13, 11, 67, 17, 23, 7, 71, 2, 73, 37, 2

Offset

1,2

Comments

The prime tower factorization of a number is defined in A182318.

Links

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

Rémy Sigrist, Illustration of initial terms

Formula

a(n) <= n with equality iff n = 1 or n is a prime number.

a(n) = a(A053585(n)).

Example

See Links section.

Prog

(PARI) a(n) = if (n==1, 1, my (f=factor(n), w=#f~); if (f[w, 2]==1, f[w, 1], a(f[w, 2])))

Crossrefs

Cf. A006530, A182318, A338669.

Sequence in context: A327076 A215041 A347241 * A284695 A081812 A356427

Adjacent sequences: A338665 A338666 A338667 * A338669 A338670 A338671

Keyword

nonn

Author

Rémy Sigrist, Apr 23 2021

Status

approved