A109674 a(n)^(n/a(n)) = A092975(n) and a(n) is a prime.

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

Offset

1,2

Links

Table of n, a(n) for n=1..96.

Maple

A092975 := proc(n) local d, a ; a :=0 ; for d in numtheory[divisors](n) do a := max(a, (n/d)^d) ; od: RETURN(a) ; end: A109674 := proc(n) local d, a092975, p ; if n = 1 then RETURN(1) ; fi ; a092975 := A092975(n) ; for d in numtheory[divisors](n) do p := n/d ; if isprime(p) and p^(n/p) = a092975 then RETURN(p) ; fi ; od: end: seq( A109674(n), n=1..120) ; # R. J. Mathar, Feb 12 2008

Crossrefs

Sequence in context: A039634 A078833 A210437 * A284260 A218993 A067629

Adjacent sequences: A109671 A109672 A109673 * A109675 A109676 A109677

Keyword

easy,nonn

Author

Vladeta Jovovic, Aug 06 2005

Extensions

More terms from R. J. Mathar, Feb 12 2008

Status

approved