|
|
A070386
|
|
a(n) = 5^n mod 39.
|
|
1
|
|
|
1, 5, 25, 8, 1, 5, 25, 8, 1, 5, 25, 8, 1, 5, 25, 8, 1, 5, 25, 8, 1, 5, 25, 8, 1, 5, 25, 8, 1, 5, 25, 8, 1, 5, 25, 8, 1, 5, 25, 8, 1, 5, 25, 8, 1, 5, 25, 8, 1, 5, 25, 8, 1, 5, 25, 8, 1, 5, 25, 8, 1, 5, 25, 8, 1, 5, 25, 8, 1, 5, 25, 8, 1, 5, 25, 8, 1, 5, 25, 8, 1, 5, 25, 8, 1, 5, 25, 8, 1, 5, 25, 8, 1
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
LINKS
|
|
|
FORMULA
|
a(n) = a(n-4).
G.f.: ( -1-5*x-25*x^2-8*x^3 ) / ( (x-1)*(1+x)*(x^2+1) ). (End)
a(n) = (13*(-1)^n-3*((8+i)*(-i)^n+(8-i)*i^n)+39)/4, where i=sqrt(-1). - Bruno Berselli, Feb 07 2011
E.g.f.: (1/2)*(26*cosh(x) + 13*sinh(x) - 24*cos(x) - 3*sin(x)). - G. C. Greubel, Mar 16 2016
|
|
MATHEMATICA
|
PowerMod[5, Range[0, 120], 39] (* or *) PadRight[{}, 120, {1, 5, 25, 8}] (* Harvey P. Dale, Nov 19 2014 *)
|
|
PROG
|
(Sage) [power_mod(5, n, 39) for n in range(0, 93)] # Zerinvary Lajos, Nov 26 2009
(PARI) a(n) = lift(Mod(5, 39)^n); \\ Altug Alkan, Mar 16 2016
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|