OFFSET
1,2
COMMENTS
See A227070 for more details and for the numbers n such that n = a(n).
The entries in the b-file have been tentatively obtained by comparing the terms < 10^30 in the sets s(n). - Giovanni Resta, Jul 30 2013
LINKS
Giovanni Resta, Table of n, a(n) for n = 1..10000 (warning: contains tentative terms)
FORMULA
Conjecture: a(n+1) = A132741(n) + 1. - Eric M. Schmidt, Jul 30 2013
MATHEMATICA
ts = {{}}; t2 = {1}; te = {1}; Do[s = Select[Range[0, 10^7], PowerMod[#, n, 10^IntegerLength[#]] == # &]; If[MemberQ[ts, s], AppendTo[t2, te[[Position[ts, s, 1, 1][[1, 1]]]]], AppendTo[ts, s]; AppendTo[te, n]; AppendTo[t2, n]], {n, 2, 82}]; t2
CROSSREFS
KEYWORD
nonn,base
AUTHOR
T. D. Noe, Jul 29 2013
EXTENSIONS
Mathematica program and some entries corrected by Giovanni Resta, Jul 30 2013
STATUS
approved