%I #66 Mar 15 2024 02:22:40
%S 1,2,4,8,7,5,1,2,4,8,7,5,1,2,4,8,7,5,1,2,4,8,7,5,1,2,4,8,7,5,1,2,4,8,
%T 7,5,1,2,4,8,7,5,1,2,4,8,7,5,1,2,4,8,7,5,1,2,4,8,7,5,1,2,4,8,7,5,1,2,
%U 4,8,7,5,1,2,4,8,7,5,1,2,4,8,7,5,1,2,4,8,7,5,1,2,4,8,7,5,1,2,4,8,7,5,1,2,4,8,7,5
%N Period 6: repeat [1, 2, 4, 8, 7, 5].
%C Digital root of 2^n.
%C A regular version of Pitoun's sequence: a(n) = A029898(n+1).
%C Also obtained from permutations of A141425, A020806, A070366, A153110, A153990, A154127, A154687, or A154815.
%C This sequence and its (again period 6) repeated differences produce the table:
%C 1, 2, 4, 8, 7, 5, 1, 2, 4, 8, 7, ...
%C 1, 2, 4, -1, -2, -4, 1, 2, 4, -1, -2, ...
%C 1, 2, -5, -1, -2, 5, 1, 2, -5, -1, -2, ...
%C 1, -7, 4, -1, 7, -4, 1, -7, 4, -1, 7, ...
%C -8, 11, -5, 8,-11, 5, -8, 11, -5, 8,-11, ...
%C 19,-16, 13,-19, 16,-13, 19,-16, 13,-19, 16, ...
%C -35, 29,-32, 35,-29, 32,-35, 29,-32, 35,-29, ...
%C 64,-61, 67,-64, 61,-67, 64,-61, 67,-64, 61, ...
%C If each entry of this table is read modulo 9 we obtain the very regular table:
%C 1, 2, 4, 8, 7, 5, 1, 2, 4, 8, 7, ...
%C 1, 2, 4, 8, 7, 5, 1, 2, 4, 8, 7, ...
%C 1, 2, 4, 8, 7, 5, 1, 2, 4, 8, 7, ...
%C 1, 2, 4, 8, 7, 5, 1, 2, 4, 8, 7, ...
%C 1, 2, 4, 8, 7, 5, 1, 2, 4, 8, 7, ...
%C 1, 2, 4, 8, 7, 5, 1, 2, 4, 8, 7, ...
%C Also the decimal expansion of the constant 125/1001. - _R. J. Mathar_, Jan 23 2009
%C Digital root of the powers of any number congruent to 2 mod 9. - _Alonso del Arte_, Jan 26 2014
%D Cecil Balmond, Number 9: The Search for the Sigma Code. Munich, New York: Prestel (1998): 203.
%H <a href="/index/Rec#order_04">Index entries for linear recurrences with constant coefficients</a>, signature (1,0,-1,1).
%F a(n) + a(n+3) = 9 = A010734(n).
%F G.f.: (1+x+2x^2+5x^3)/((1-x)(1+x)(1-x+x^2)). - _R. J. Mathar_, Jan 23 2009
%F a(n) = A082365(n) mod 9. - _Paul Curtz_, Mar 31 2009
%F a(n) = -1/2*cos(Pi*n) - 3*cos(1/3*Pi*n) - 3^(1/2)*sin(1/3*Pi*n) + 9/2. - _Leonid Bedratyuk_, May 13 2012
%F a(n) = A010888(A004000(n+1)). - _Ivan N. Ianakiev_, Nov 27 2014
%F From _Wesley Ivan Hurt_, Apr 20 2015: (Start)
%F a(n) = a(n-6) for n>5.
%F a(n) = a(n-1) - a(n-3) + a(n-4) for n>3.
%F a(n) = (2+3*(n-1 mod 3))*(n mod 2) + (1+3*(-n mod 3))*(n-1 mod 2). (End)
%F a(n) = 2^n mod 9. - _Nikita Sadkov_, Oct 06 2018
%p seq(op([1, 2, 4, 8, 7, 5]), n=0..40); # _Wesley Ivan Hurt_, Jul 05 2016
%t Flatten[Table[{1, 2, 4, 8, 7, 5}, {20}]] (* _Paul Curtz_, Dec 19 2008 *)
%t Table[Mod[2^n, 9], {n, 0, 99}] (* _Alonso del Arte_, Jan 26 2014 *)
%o (PARI) a(n)=lift(Mod(2,9)^n) \\ _Charles R Greathouse IV_, Apr 21 2015
%o (Magma) &cat [[1, 2, 4, 8, 7, 5]^^30]; // _Wesley Ivan Hurt_, Jul 05 2016
%Y Cf. A030132, A145389, A189510.
%Y Cf. digital roots of powers of c mod 9: c = 4, A100402; c = 5, A070366; c = 7, A070403; c = 8, A010689.
%K nonn,easy
%O 0,2
%A _Paul Curtz_, Dec 19 2008
%E Edited by _R. J. Mathar_, Apr 09 2009