OFFSET
1,3
COMMENTS
90625^6 = 553972386755049228668212890625 hence 90625 is in the sequence.
LINKS
MATHEMATICA
okQ[n_]:=Module[{idn=IntegerDigits[n], id6n=IntegerDigits[n^6]}, idn==Take[id6n, -Length[idn]]]
Select[Range[120000], okQ] (* Harvey P. Dale, Jan 16 2011 *)
PROG
(Sage)
def automorphic(maxdigits, pow, base=10) :
morphs = [[0]]
for i in range(maxdigits):
T=[d*base^i+x for x in morphs[-1] for d in range(base)]
morphs.append([x for x in T if x^pow % base^(i+1) == x])
res = list(set(sum(morphs, []))); res.sort()
return res
# (call with pow=6 for this sequence), Eric M. Schmidt, Jul 29 2013
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Benoit Cloitre, Mar 08 2002
EXTENSIONS
More terms from Eric M. Schmidt, Jul 29 2013
STATUS
approved