|
%I #22 Feb 22 2022 23:07:25
%S 1,2,4,11,19,22,29,31,38,44,58,59,62,71,76,79,101,116,118,121,124,131,
%T 139,142,151,158,179,181,191,199,202,209,211,229,236,239,242,251,262,
%U 271,278,284,302,311,316,319,331,341,349,358,359,361,362,379,382,398
%N Numbers with 1 zero in Fibonacci numbers mod m.
%C Conjecture: m is on this list iff m is an odd number all of whose factors are on this list or m is 2 or 4 times such an odd number.
%C A001176(a(n)) = A128924(a(n),1) = 1. - _Reinhard Zumkeller_, Jan 16 2014
%C Also numbers n such that A001175(n) = A001177(n). - _Daniel Suteu_, Aug 08 2018
%H Reinhard Zumkeller, <a href="/A053031/b053031.txt">Table of n, a(n) for n = 1..1000</a>
%H Brennan Benfield and Michelle Manes, <a href="https://arxiv.org/abs/2202.08986">The Fibonacci Sequence is Normal Base 10</a>, arXiv:2202.08986 [math.NT], 2022.
%H M. Renault, <a href="http://webspace.ship.edu/msrenault/fibonacci/fib.htm">Fibonacci sequence modulo m</a>
%t With[{s = {1}~Join~Table[Count[Drop[NestWhile[Append[#, Mod[Total@ Take[#, -2], n]] &, {1, 1}, If[Length@ # < 3, True, Take[#, -2] != {1, 1}] &], -2], 0], {n, 2, 400}]}, Position[s, 1][[All, 1]] ] (* _Michael De Vlieger_, Aug 08 2018 *)
%o (Haskell)
%o a053031 n = a053031_list !! (n-1)
%o a053031_list = filter ((== 1) . a001176) [1..]
%o -- _Reinhard Zumkeller_, Jan 16 2014
%o (PARI) entryp(p)=my(k=p+[0, -1, 1, 1, -1][p%5+1], f=factor(k)); for(i=1, #f[, 1],for(j=1, f[i, 2], if((Mod([1, 1; 1, 0], p)^(k/f[i, 1]))[1, 2], break); k/=f[i, 1])); k
%o entry(n)=if(n==1, return(1)); my(f=factor(n), v); v=vector(#f~, i, if(f[i, 1]>1e14, entryp(f[i, 1]^f[i, 2]), entryp(f[i, 1])*f[i, 1]^(f[i, 2]-1))); if(f[1, 1]==2&&f[1, 2]>1, v[1]=3<<max(f[1, 2]-2, 1)); lcm(v)
%o fibmod(n, m)=((Mod([1, 1; 1, 0], m))^n)[1, 2]
%o is(n)=fibmod(entry(n)+1,n)==1 \\ _Charles R Greathouse IV_, Dec 14 2016
%Y Cf. A001176, A053032.
%Y Cf. A053029 (with 4 zeros), A053030 (with 2 zeros).
%K nonn
%O 1,2
%A _Henry Bottomley_, Feb 23 2000
|
