1,2
Table of n, a(n) for n=1..5.
a(3) = 125 = 5^3 = 31*a(2) + 1 where a(2) = 4=2^2 = 3*1 + 1.
(PARI) v=vector(15):v[2]=4:print1("1, 4, "):for(n=3, 15, for(k=2, 10^8, if(((k^n)%(v[n-1]))==1, print1(k^n", "):v[n]=k^n:break)))
Sequence in context: A073351 A195797 A286017 * A116944 A124260 A253099
Adjacent sequences: A068891 A068892 A068893 * A068895 A068896 A068897
nonn
Amarnath Murthy, Mar 20 2002
Corrected and extended by Ralf Stephan, Mar 21 2003
approved
