OFFSET
1,2
LINKS
Robert Price, Table of n, a(n) for n = 1..2000
EXAMPLE
For a(12), 12, 14, and 15 are not prime powers, and 13 is a prime power but it is prime. Since 16 = 2^4 is a prime power, a(12) = 16.
MATHEMATICA
Module[{ppwrs=Join[{1}, Sort[Flatten[Table[Prime[Range[5]]^p, {p, 2, 10}]]]]}, Flatten[ Table[Select[ppwrs, #>=n&, 1], {n, 90}]]] (* Harvey P. Dale, Oct 06 2014 *)
PROG
(PARI) isA025475(n) = (omega(n) == 1 & !isprime(n)) || (n == 1)
A167184(n) = {local(m); m=n; while(!isA025475(m), m++); m}
(Python)
from itertools import count
from sympy import factorint
def A167184(n): return next(filter(lambda m:len(f:=factorint(m))<=1 and max(f.values(), default=2)>1, count(n))) # Chai Wah Wu, Oct 25 2024
CROSSREFS
KEYWORD
easy,nonn,changed
AUTHOR
Michael B. Porter, Oct 29 2009
STATUS
approved