A167166 a(n) = n^7 mod 16.

0, 1, 0, 11, 0, 13, 0, 7, 0, 9, 0, 3, 0, 5, 0, 15, 0, 1, 0, 11, 0, 13, 0, 7, 0, 9, 0, 3, 0, 5, 0, 15, 0, 1, 0, 11, 0, 13, 0, 7, 0, 9, 0, 3, 0, 5, 0, 15, 0, 1, 0, 11, 0, 13, 0, 7, 0, 9, 0, 3, 0, 5, 0, 15, 0, 1, 0, 11, 0, 13, 0, 7, 0, 9, 0, 3, 0, 5, 0, 15, 0, 1, 0, 11, 0, 13, 0, 7, 0, 9, 0, 3, 0

Offset

0,4

Comments

Equivalently: n^(4*m+7) mod 16. - G. C. Greubel, Jun 04 2016

Links

Table of n, a(n) for n=0..92.

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

Formula

From R. J. Mathar, Sep 30 2013: (Start)

a(n) = a(n-16).

G.f. -x*(1 +11*x^2 +13*x^4 +7*x^6 +9*x^8 +3*x^10 +5*x^12 +15*x^14) / ( (x-1)*(1+x)*(1+x^2)*(1+x^4)*(1+x^8) ). (End)

a(n) = A130909(A001015(n)). - Michel Marcus, Jun 04 2016

Mathematica

Table[Mod[n^7, 16], {n, 0, 10}] (* G. C. Greubel, Jun 04 2016 *)
PowerMod[Range[0, 100], 7, 16] (* or *) PadRight[{}, 100, {0, 1, 0, 11, 0, 13, 0, 7, 0, 9, 0, 3, 0, 5, 0, 15}] (* Harvey P. Dale, Jul 29 2018 *)

Prog

(Sage) [power_mod(n, 7, 16)for n in range(0, 93)] #
(PARI) a(n)=n^7%16 \\ Charles R Greathouse IV, Apr 06 2016

Crossrefs

Sequence in context: A256480 A088621 A088623 * A271572 A180059 A147765

Adjacent sequences: A167163 A167164 A167165 * A167167 A167168 A167169

Keyword

nonn,easy

Author

Zerinvary Lajos, Oct 29 2009

Status

approved