OFFSET
0,2
REFERENCES
I. M. Vinogradov, Elements of Number Theory, pp. 220 ff.
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..10000
Index entries for linear recurrences with constant coefficients, signature (1,0,0,0,0,0,0,0,0,0,0,0,0,0,-1,1).
FORMULA
a(n) = a(n+30). - R. J. Mathar, Jun 04 2016
a(n) = a(n-1) - a(n-15) + a(n-16). - G. C. Greubel, Oct 16 2018
a(n) = 31 - a(n+15) for all n in Z. - Michael Somos, Oct 17 2018
MAPLE
[ seq(primroot(ithprime(i))^j mod ithprime(i), j=0..100) ];
MATHEMATICA
PowerMod[3, Range[0, 100], 31] (* G. C. Greubel, Oct 16 2018 *)
PROG
(PARI) a(n)=lift(Mod(3, 31)^n) \\ Charles R Greathouse IV, Mar 22 2016
(Magma) [Modexp(3, n, 31): n in [0..100]]; // G. C. Greubel, Oct 18 2018
(GAP) List([0..65], n->PowerMod(3, n, 31)); # Muniru A Asiru, Oct 18 2018
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
STATUS
approved