A070417 a(n) = 7^n mod 33.

1, 7, 16, 13, 25, 10, 4, 28, 31, 19, 1, 7, 16, 13, 25, 10, 4, 28, 31, 19, 1, 7, 16, 13, 25, 10, 4, 28, 31, 19, 1, 7, 16, 13, 25, 10, 4, 28, 31, 19, 1, 7, 16, 13, 25, 10, 4, 28, 31, 19, 1, 7, 16, 13, 25, 10, 4, 28, 31, 19, 1, 7, 16, 13, 25, 10, 4, 28, 31, 19, 1, 7, 16, 13, 25, 10, 4

Offset

0,2

Links

G. C. Greubel, Table of n, a(n) for n = 0..1000

Index entries for linear recurrences with constant coefficients, signature (1,-1,1,-1,1,-1,1,-1,1).

Formula

From R. J. Mathar, Apr 20 2010: (Start)

a(n) = a(n-1) - a(n-2) + a(n-3) - a(n-4) + a(n-5) - a(n-6) + a(n-7) - a(n-8) + a(n-9).

G.f.: ( -1-6*x-10*x^2-3*x^3-22*x^4+12*x^5-16*x^6-12*x^7-19*x^8 ) / ( (x-1)*(1-x+x^2-x^3+x^4)*(1+x+x^2+x^3+x^4) ). (End)

a(n) = a(n-10). - G. C. Greubel, Mar 20 2016

Mathematica

PowerMod[7, Range[0, 80], 33] (* Harvey P. Dale, Mar 04 2015 *)

Prog

(Sage) [power_mod(7, n, 33) for n in range(0, 77)] # Zerinvary Lajos, Nov 27 2009
(PARI) a(n)=lift(Mod(7, 33)^n) \\ Charles R Greathouse IV, Mar 22 2016
(Magma) [Modexp(7, n, 33): n in [0..100]]; // Bruno Berselli, Mar 22 2016

Crossrefs

Sequence in context: A061039 A258771 A063593 * A101681 A067776 A193688

Adjacent sequences: A070414 A070415 A070416 * A070418 A070419 A070420

Keyword

nonn,easy

Author

N. J. A. Sloane, May 12 2002

Status

approved