OFFSET
1,2
COMMENTS
a(16) > 10^45 if it exists. - Michael S. Branicky, Jun 19 2022
MATHEMATICA
perfectPowerQ[n_] := n==1||GCD @@ FactorInteger[n][[All, 2]] > 1; (* A001597 *) Select[Range[10^5], perfectPowerQ[#] && Max[Differences[d=IntegerDigits[#]]]<1 && Count[d, 0]==0&] (* Stefano Spezia, Jul 01 2025 *)
PROG
(Python)
from sympy import perfect_power as pp
from itertools import count, islice, combinations_with_replacement as mc
def agen():
yield 1
for d in count(1):
nd = (int("".join(m)) for m in mc("987654321", d))
yield from sorted(filter(pp, nd))
print(list(islice(agen(), 14))) # Michael S. Branicky, Jun 16 2022
(PARI) isok(m) = if (ispower(m), my(d=digits(m)); vecmin(d) && (d == vecsort(d, , 4))); \\ Michel Marcus, Jun 17 2022
CROSSREFS
KEYWORD
nonn,base,more
AUTHOR
Jon E. Schoenfield, Jun 16 2022
STATUS
approved
