A036132 a(n) = 7^n mod 71.

1, 7, 49, 59, 58, 51, 2, 14, 27, 47, 45, 31, 4, 28, 54, 23, 19, 62, 8, 56, 37, 46, 38, 53, 16, 41, 3, 21, 5, 35, 32, 11, 6, 42, 10, 70, 64, 22, 12, 13, 20, 69, 57, 44, 24, 26, 40, 67, 43, 17, 48, 52, 9, 63, 15, 34, 25, 33

Offset: 0

Author

N. J. A. Sloane

Keywords

• nonn
• easy

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, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 1).

Crossrefs

Cf. A000420 (7^n).

Programs

GAP

List([0..60],n->PowerMod(7,n,71)); # Muniru A Asiru, Oct 17 2018

Magma

[Modexp(7, n, 71): n in [0..100]]; // G. C. Greubel, Oct 17 2018

Maple

[ seq(primroot(ithprime(i))^j mod ithprime(i),j=0..100) ];

Mathematica

PowerMod[7, Range[0, 100], 71] (* G. C. Greubel, Oct 17 2018 *)

PARI

a(n)=lift(Mod(7,71)^n) \\ Charles R Greathouse IV, Mar 22 2016

Python

for n in range(0, 100): print(int(pow(7, n, 71)), end=' ') # Stefano Spezia, Oct 17 2018

Formula

a(n) = a(n+70). - R. J. Mathar, Jun 04 2016

a(n) = a(n-1) - a(n-35) + a(n-36). - G. C. Greubel, Oct 17 2018