A308948 a(n) = A006190(A322907(n)+1) mod n.
0, 1, 1, 1, 3, 1, 6, 5, 1, 3, 10, 1, 8, 13, 4, 9, 16, 1, 18, 9, 13, 21, 1, 13, 18, 5, 1, 13, 12, 19, 30, 17, 10, 33, 6, 1, 31, 37, 25, 29, 32, 13, 1, 21, 19, 1, 46, 25, 48, 43, 16, 25, 1, 1, 21, 41, 37, 17, 58, 49, 1, 61, 55, 33, 18, 43, 66, 33, 1, 41, 70, 37
Offset: 1
Keywords
Examples
For n = 7, {A006190(n) mod 7 : n > 0} = {1, 3, 3, 5, 4, 3, 6, 0, 6, ...}, so a(7) = 6. Also, A322907(7) = 8, and M^8 mod 7 = [{6, 0}, {0, 6}], so a(7) = 6.
Links
- Amiram Eldar, Table of n, a(n) for n = 1..10000
Programs
-
Mathematica
a[n_] := For[k = 1, True, k++, If[Divisible[Fibonacci[k, 3], n], Return[ Mod[Fibonacci[k + 1, 3], n]]]]; Array[a, 100] (* Jean-François Alcover, Jul 05 2019 *)
-
PARI
a(n) = my(M=[3, 1; 1, 0]); for(k=1, 12*n/7, if((Mod(M,n)^k)[2,1]==0, return(lift((Mod(M,n)^k)[1,1]))))
Comments