OFFSET
1,1
COMMENTS
If a(n) differs from 1, then a(n) is the minimal prime divisor of A064079(n).
LINKS
Max Alekseyev, Table of n, a(n) for n = 1..730 (first 153 terms from Robert G. Wilson v)
MAPLE
a:= proc(n) local f, p;
f:= numtheory:-factorset(3^n - 1);
for p in f do
if numtheory:-order(3, p) = n then return p fi
od:
1
end proc:
seq(a(n), n=1..100); # Robert Israel, Oct 13 2014
MATHEMATICA
p = 2; t = Table[0, {100}]; While[p < 100000001, a = MultiplicativeOrder[3, p]; If[0 < a < 101 && t[[a]] == 0, t[[a]] = p; Print[{a, p}]]; p = NextPrime@ p]; t (* Robert G. Wilson v, Oct 13 2014 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Shevelev, Aug 28 2008
EXTENSIONS
More terms from Robert G. Wilson v, Dec 11 2013
STATUS
approved