A322906 -id:A322906 - cp's OEIS frontend

A172236 Array A(n,k) = n*A(n,k-1) + A(n,k-2) read by upward antidiagonals, starting A(n,0) = 0, A(n,1) = 1.

0, 0, 1, 0, 1, 1, 0, 1, 2, 2, 0, 1, 3, 5, 3, 0, 1, 4, 10, 12, 5, 0, 1, 5, 17, 33, 29, 8, 0, 1, 6, 26, 72, 109, 70, 13, 0, 1, 7, 37, 135, 305, 360, 169, 21, 0, 1, 8, 50, 228, 701, 1292, 1189, 408, 34, 0, 1, 9, 65, 357, 1405, 3640, 5473, 3927, 985, 55, 0, 1, 10, 82, 528, 2549, 8658, 18901, 23184, 12970, 2378, 89, 0, 1, 11, 101, 747, 4289, 18200, 53353, 98145, 98209, 42837, 5741, 144

Offset: 1

Roger L. Bagula, Jan 29 2010

Equals A073133 with an additional column A(.,0).
If the first column and top row are deleted, antidiagonal reading yields A118243.
Adding a top row of 1's and antidiagonal reading downwards yields A157103.
Antidiagonal sums are 0, 1, 2, 5, 12, 32, 93, 297, 1035, 3911, 15917, 69350, ....
From Jianing Song, Jul 14 2018: (Start)
All rows have strong divisibility, that is, gcd(A(n,k_1), A(n,k_2)) = A(n,gcd(k_1,k_2)) holds for all k_1, k_2 >= 0.
Let E(n,m) be the smallest number l such that m divides A(n,l), we have: for odd primes p that are not divisible by n^2 + 4, E(n,p) divides p - ((n^2+4)/p) if p == 3 (mod 4) and (p - ((n^2+4)/p))/2 if p == 1 (mod 4). E(n,p) = p for odd primes p that are divisible by n^2 + 4. E(n,2) = 2 for even n and 3 for odd n. Here ((n^2+4)/p) is the Legendre symbol. A prime p such that p^2 divides T(n,E(n,p)) is called an n-Wall-Sun-Sun prime.
E(n,p^e) <= p^(e-1)*E(n,p) for all primes p. If p^2 does not divide A(n, E(n,p)), then E(n,p^e) = p^(e-1)*E(n,p) for all exponent e. Specially, for primes p >= 5 that are divisible by n^2 + 4, p^2 is never divisible by A(n,p), so E(n,p^e) = p^e as described above. E(n,m_1*m_2) = lcm(E(n,m_1),E(n,m_2)) if gcd(m_1,m_2) = 1.
Let pi(n,m) be the Pisano period of A(n, k) modulo m, i.e, the smallest number l such that A(n, k+1) == T(n,k) (mod m) holds for all k >= 0, we have: for odd primes p that are not divisible by n^2 - 4, pi(n,p) divides p - 1 if ((n^2+4)/p) = 1 and 2(p+1) if ((n^2+4)/p) = -1. pi(n,p) = 4p for odd primes p that are divisible by n^2 + 4. pi(n,2) = 2 even n and 3 for odd n.
pi(n,p^e) <= p^(e-1)*pi(n,p) for all primes p. If p^2 does not divide A(n, E(n,p)), then pi(n,p^e) = p^(e-1)*pi(n,p) for all exponent e. Specially, for primes p >= 5 that are divisible by n^2 + 4, p^2 is never divisible by A(n, p), so pi(n,p^e) = 4p^e as described above. pi(n,m_1*m_2) = lcm(pi(n,m_1), pi(n,m_2)) if gcd(m_1,m_2) = 1.
If n != 2, the largest possible value of pi(n,m)/m is 4 for even n and 6 for odd n. For even n, pi(n,p^e) = 4p^e; for odd n, pi(n,2p^e) = 12p^e, where p is any odd prime factor of n^2 + 4. For n = 2 it is 8/3, obtained by m = 3^e.
Let z(n,m) be the number of zeros in a period of A(n, k) modulo m, i.e., z(n,m) = pi(n,m)/E(n,m), then we have: z(n,p) = 4 for odd primes p that are divisible by n^2 + 4. For other odd primes p, z(n,p) = 4 if E(n,p) is odd; 1 if E(n,p) is even but not divisible by 4; 2 if E(n,p) is divisible by 4; see the table below. z(n,2) = z(n,4) = 1.
Among all values of z(n,p) when p runs through all odd primes that are not divisible by n^2 + 4, we have:
((n^2+4)/p)...p mod 8....proportion of 1.....proportion of 2.....proportion of 4
......1..........1......1/6 (conjectured)...2/3 (conjectured)...1/6 (conjectured)*
......1..........5......1/2 (conjectured)...........0...........1/2 (conjectured)*
......1.........3,7.............1...................0...................0
.....-1.........1,5.............0...................0...................1
.....-1.........3,7.............0...................1...................0
* The result is that among all odd primes that are not divisible by n^2 + 4, 7/24 of them are with z(n,p) = 1, 5/12 are with z(n,p) = 2 and 7/24 are with z(n,p) = 4 if n^2 + 4 is a twice a square; 1/3 of them are with z(n,p) = 1, 1/3 are with z(n,p) = 2 and 1/3 are with z(n,p) = 4 otherwise. [Corrected by Jianing Song, Jul 06 2019]
z(n,p^e) = z(n,p) for all odd primes p; z(n,2^e) = 1 for even n and 2 for odd n, e >= 3.
(End)
From Michael A. Allen, Mar 06 2023: (Start)
Removing the first (n=0) row of A352361 gives this sequence.
Row n is the n-metallonacci sequence.
A(n,k) is (for k>0) the number of tilings of a (k-1)-board (a board with dimensions (k-1) X 1) using unit squares and dominoes (with dimensions 2 X 1) if there are n kinds of squares available. (End)

			The array, A(n, k), starts in row n = 1 with columns k >= 0 as
  0      1      1      2      3      5      8
  0      1      2      5     12     29     70
  0      1      3     10     33    109    360
  0      1      4     17     72    305   1292
  0      1      5     26    135    701   3640
  0      1      6     37    228   1405   8658
  0      1      7     50    357   2549  18200
  0      1      8     65    528   4289  34840
  0      1      9     82    747   6805  61992
  0      1     10    101   1020  10301 104030
  0      1     11    122   1353  15005 166408
Antidiagonal triangle, T(n, k), begins as:
  0;
  0, 1;
  0, 1, 1;
  0, 1, 2,  2;
  0, 1, 3,  5,   3;
  0, 1, 4, 10,  12,   5;
  0, 1, 5, 17,  33,  29,    8;
  0, 1, 6, 26,  72, 109,   70,   13;
  0, 1, 7, 37, 135, 305,  360,  169,  21;
  0, 1, 8, 50, 228, 701, 1292, 1189, 408, 34;

Jianing Song, Antidiagonals n = 1..100, flattened (the first 30 antidiagonals from Vincenzo Librandi)
Michael A. Allen and Kenneth Edwards, Fence tiling derived identities involving the metallonacci numbers squared or cubed, Fib. Q. 60:5 (2022) 5-17.

Rows n include: A000045 (n=1), A000129 (n=2), A006190 (n=3), A001076 (n=4), A052918 (n=5), A005668 (n=6), A054413 (n=7), A041025 (n=8), A099371 (n=9), A041041 (n=10), A049666 (n=11), A041061 (n=12), A140455 (n=13), A041085 (n=14), A154597 (n=15), A041113 (n=16), A178765 (n=17), A041145 (n=18), A243399 (n=19), A041181 (n=20). (Note that there are offset shifts for rows n = 5, 7, 8, 10, 12, 14, 16..20.)
Columns k include: A000004 (k=0), A000012 (k=1), A000027 (k=2), A002522 (k=3), A054602 (k=4), A057721 (k=5), A124152 (k=6).
Entry points for A(n,k) modulo m: A001177 (n=1), A214028 (n=2), A322907 (n=3).
Pisano period for A(n,k) modulo m: A001175 (n=1), A175181 (n=2), A175182 (n=3), A175183 (n=4), A175184 (n=5), A175185 (n=6).
Number of zeros in a period for A(n,k) modulo m: A001176 (n=1), A214027 (n=2), A322906 (n=3).
Cf. A084844, A084845, A118243, A157103, A271782, A316269.
Sums include: A304357, A304359.
Similar to: A073133.

A172236:= func< n,k | k le 1 select k else Evaluate(DicksonSecond(k-1,-1), n-k) >;
[A172236(n,k): k in [0..n-1], n in [1..13]]; // G. C. Greubel, Sep 29 2024

A172236[n_,k_]:=Fibonacci[k, n-k];
Table[A172236[n, k], {n,15}, {k,0,n-1}]//Flatten

A(n, k) = if (k==0, 0, if (k==1, 1, n*A(n, k-1) + A(n, k-2)));
tabl(nn) = for(n=1, nn, for (k=0, nn, print1(A(n, k), ", ")); print); \\ Jianing Song, Jul 14 2018 (program from Michel Marcus; see also A316269)

A(n, k) = ([n, 1; 1, 0]^k)[2, 1] \\ Jianing Song, Nov 10 2018

def A172236(n,k): return sum(binomial(k-j-1,j)*(n-k)^(k-2*j-1) for j in range(1+(k-1)//2))
flatten([[A172236(n,k) for k in range(n)] for n in range(1,14)]) # G. C. Greubel, Sep 29 2024

A(n,k) = (((n + sqrt(n^2 + 4))/2)^k - ((n-sqrt(n^2 + 4))/2)^k)/sqrt(n^2 + 4), n >= 1, k >= 0. - Jianing Song, Jun 27 2018
For n >= 1, Sum_{i=1..k} 1/A(n,2^i) = ((u^(2^k-1) + v^(2^k-1))/(u + v)) * (1/A(n,2^k)), where u = (n + sqrt(n^2 + 4))/2, v = (n - sqrt(n^2 + 4))/2 are the two roots of the polynomial x^2 - n*x - 1. As a result, Sum_{i>=1} 1/A(n,2^i) = (n^2 + 4 - n*sqrt(n^2 + 4))/(2*n). - Jianing Song, Apr 21 2019
From G. C. Greubel, Sep 29 2024: (Start)
A(n, k) = F_{k}(n) (Fibonacci polynomials F_{n}(x)) (array).
T(n, k) = F_{k}(n-k) (antidiagonal triangle).
Sum_{k=0..n-1} T(n, k) = A304357(n) - (1-(-1)^n)/2.
Sum_{k=0..n-1} (-1)^k*T(n, k) = (-1)*A304359(n) + (1-(-1)^n)/2.
T(2*n, n) = A084844(n).
T(2*n+1, n+1) = A084845(n). (End)

More terms from Jianing Song, Jul 14 2018

A053028 Odd primes p with 4 zeros in any period of the Fibonacci numbers mod p.

Original entry on oeis.org

5, 13, 17, 37, 53, 61, 73, 89, 97, 109, 113, 137, 149, 157, 173, 193, 197, 233, 257, 269, 277, 293, 313, 317, 337, 353, 373, 389, 397, 421, 433, 457, 557, 577, 593, 613, 617, 653, 661, 673, 677, 701, 733, 757, 761, 773, 797, 821, 829, 853, 857, 877, 937, 953

Offset: 1

Henry Bottomley, Feb 23 2000

nonn

Also, primes that do not divide any Lucas number. - T. D. Noe, Jul 25 2003
Although every prime divides some Fibonacci number, this is not true for the Lucas numbers. In fact, exactly 1/3 of all primes do not divide any Lucas number. See Lagarias and Moree for more details. The Lucas numbers separate the primes into three disjoint sets: (A053028) primes that do not divide any Lucas number, (A053027) primes that divide Lucas numbers of even index and (A053032) primes that divide Lucas numbers of odd index. - T. D. Noe, Jul 25 2003; revised by N. J. A. Sloane, Feb 21 2004
From Jianing Song, Jun 16 2024: (Start)
Primes p such that A001176(p) = 4.
For p > 2, p is in this sequence if and only if A001175(p) == 4 (mod 8), and if and only if A001177(p) is odd. For a proof of the equivalence between A001176(p) = 4 and A001177(p) being odd, see Section 2 of my link below.
This sequence contains all primes congruent to 13, 17 (mod 20). This corresponds to case (1) for k = 3 in the Conclusion of Section 1 of my link below. (End) [Comment rewritten by Jianing Song, Jun 16 2024 and Jun 25 2024]

T. D. Noe, Table of n, a(n) for n=1..1000
C. Ballot and M. Elia, Rank and period of primes in the Fibonacci sequence; a trichotomy, Fib. Quart., 45 (No. 1, 2007), 56-63 (The sequence B2).
Nicholas Bragman and Eric Rowland, Limiting density of the Fibonacci sequence modulo powers of p, arXiv:2202.00704 [math.NT], 2022.
J. C. Lagarias, The set of primes dividing the Lucas numbers has density 2/3, Pacific J. Math., 118. No. 2, (1985), 449-461.
J. C. Lagarias, Errata to: The set of primes dividing the Lucas numbers has density 2/3, Pacific J. Math., 162, No. 2, (1994), 393-396.
Diego Marques and Pavel Trojovsky, The order of appearance of the product of five consecutive Lucas numbers, Tatra Mountains Math. Publ. 59 (2014), 65-77.
Pieter Moree, Counting Divisors of Lucas Numbers, Pacific J. Math, Vol. 186, No. 2, 1998, pp. 267-284.
M. Renault, Fibonacci sequence modulo m
H. Sedaghat, Zero-Avoiding Solutions of the Fibonacci Recurrence Modulo A Prime, Fibonacci Quart. 52 (2014), no. 1, 39-45. See p. 45.
Jianing Song, Lucas sequences and entry point modulo p
Eric Weisstein's World of Mathematics, Lucas Number

Cf. A001175, A001177.
Cf. A000204 (Lucas numbers), A001602 (index of the smallest Fibonacci number divisible by prime(n)).
Let {x(n)} be a sequence defined by x(0) = 0, x(1) = 1, x(n+2) = m*x(n+1) + x(n). Let w(k) be the number of zeros in a fundamental period of {x(n)} modulo k.
                             |   m=1    |   m=2   |   m=3
-----------------------------+----------+---------+---------
The sequence {x(n)}          | A000045  | A000129 | A006190
The sequence {w(k)}          | A001176  | A214027 | A322906
Primes p such that w(p) = 1  | A112860* | A309580 | A309586
Primes p such that w(p) = 2  | A053027  | A309581 | A309587
Primes p such that w(p) = 4  | this seq | A261580 | A309588
Numbers k such that w(k) = 1 | A053031  | A309583 | A309591
Numbers k such that w(k) = 2 | A053030  | A309584 | A309592
Numbers k such that w(k) = 4 | A053029  | A309585 | A309593
* and also A053032 U {2}

Lucas[n_] := Fibonacci[n+1] + Fibonacci[n-1]; badP={}; Do[p=Prime[n]; k=1; While[k0, k++ ]; If[k==p, AppendTo[badP, p]], {n, 200}]; badP

A prime p = prime(i) is in this sequence if p > 2 and A001602(i) is odd. - T. D. Noe, Jul 25 2003

Edited: Name clarified. Moree and Renault link updated. Ballot and Elia reference linked. - Wolfdieter Lang, Jan 20 2015

A053032 Odd primes p with one zero in Fibonacci numbers mod p.

Original entry on oeis.org

11, 19, 29, 31, 59, 71, 79, 101, 131, 139, 151, 179, 181, 191, 199, 211, 229, 239, 251, 271, 311, 331, 349, 359, 379, 419, 431, 439, 461, 479, 491, 499, 509, 521, 541, 571, 599, 619, 631, 659, 691, 709, 719, 739, 751, 809, 811, 839, 859, 911, 919, 941, 971

Offset: 1

Henry Bottomley, Feb 23 2000

nonn

Also, odd primes that divide Lucas numbers of odd index. - T. D. Noe, Jul 25 2003
From Charles R Greathouse IV, Dec 14 2016: (Start)
It seems that this sequence contains about 1/3 of the primes. In particular, members of this sequence constitute:
        35 of the first 10^2 primes
       330 of the first 10^3 primes
      3328 of the first 10^4 primes
     33371 of the first 10^5 primes
    333329 of the first 10^6 primes
   3333720 of the first 10^7 primes
  33333463 of the first 10^8 primes
etc. (End)
Of the Fibonacci-like sequences modulo a prime p that are not A000004, one of them has a period length less than A001175(p) if and only if p = 5 or p is in this sequence. - Isaac Saffold, Dec 18 2018
Odd primes in A053031. - Jianing Song, Jun 19 2019

			From _Michael B. Porter_, Jan 25 2019: (Start)
The Fibonacci numbers (mod 7) repeat the pattern 0, 1, 1, 2, 3, 5, 1, 6, 0, 6, 6, 5, 4, 2, 6, 1. Since there are two zeros, 7 is not in the sequence.
The Fibonacci numbers (mod 11) repeat the pattern 0, 1, 1, 2, 3, 5, 8, 2, 10, 1 which has only one zero, so 11 is in the sequence.
(End)

T. D. Noe, Table of n, a(n) for n = 1..1000
C. Ballot and M. Elia, Rank and period of primes in the Fibonacci sequence; a trichotomy, Fib. Quart., 45 (No. 1, 2007), 56-63 (The sequence B1).
Nicholas Bragman and Eric Rowland, Limiting density of the Fibonacci sequence modulo powers of p, arXiv:2202.00704 [math.NT], 2022.
M. Renault, Fibonacci sequence modulo m

Cf. A001175, A001177. See A112860 for another version.
Cf. A000204 (Lucas numbers), A001602 (index of the smallest Fibonacci number divisible by prime(n)).
Let {x(n)} be a sequence defined by x(0) = 0, x(1) = 1, x(n+2) = m*x(n+1) + x(n). Let w(k) be the number of zeros in a fundamental period of {x(n)} modulo k.
                             |    m=1     |   m=2   |   m=3
-----------------------------+------------+---------+---------
The sequence {x(n)}          | A000045    | A000129 | A006190
The sequence {w(k)}          | A001176    | A214027 | A322906
Primes p such that w(p) = 1  | A112860*   | A309580 | A309586
Primes p such that w(p) = 2  | A053027**  | A309581 | A309587
Primes p such that w(p) = 4  | A053028*** | A261580 | A309588
Numbers k such that w(k) = 1 | A053031    | A309583 | A309591
Numbers k such that w(k) = 2 | A053030    | A309584 | A309592
Numbers k such that w(k) = 4 | A053029    | A309585 | A309593
* and also this sequence U {2}
** also primes dividing Lucas numbers of even index
*** also primes dividing no Lucas number

Prime@ Rest@ Position[Table[Count[Drop[NestWhile[Append[#, Mod[Total@ Take[#, -2], n]] &, {1, 1}, If[Length@ # < 3, True, Take[#, -2] != {1, 1}] &], -2], 0], {n, Prime@ Range@ 168}], 1][[All, 1]] (* Michael De Vlieger, Aug 08 2018 *)

fibmod(n,m)=(Mod([1, 1; 1, 0], m)^n)[1, 2]
is(n)=my(k=n+[0, -1, 1, 1, -1][n%5+1]); k>>=valuation(k,2)-1; fibmod(k,n)==0 && fibmod(k/2,n) && isprime(n) \\ Charles R Greathouse IV, Dec 14 2016

A prime p = prime(i) is in this sequence if p > 2 and A001602(i)/2 is odd. - T. D. Noe, Jul 25 2003

A053027 Odd primes p with 2 zeros in Fibonacci numbers mod p.

Original entry on oeis.org

3, 7, 23, 41, 43, 47, 67, 83, 103, 107, 127, 163, 167, 223, 227, 241, 263, 281, 283, 307, 347, 367, 383, 401, 409, 443, 449, 463, 467, 487, 503, 523, 547, 563, 569, 587, 601, 607, 641, 643, 647, 683, 727, 743, 769, 787, 823, 827, 863, 881, 883, 887, 907, 929

Offset: 1

Henry Bottomley, Feb 23 2000

nonn

Also, odd primes that divide Lucas numbers of even index. - T. D. Noe, Jul 25 2003
Primes in A053030. - Jianing Song, Jun 19 2019
From Jianing Song, Jun 16 2024: (Start)
Primes p such that A001176(p) = 2.
For p > 2, p is in this sequence if and only if 8 divides of A001175(p), and if and only if 4 divides A001177(p). For a proof of the equivalence between A001176(p) = 2 and 4 dividing A001177(p), see Section 2 of my link below.
This sequence contains all primes congruent to 3, 7 (mod 20). This corresponds to case (2) for k = 3 in the Conclusion of Section 1 of my link below.
Conjecturely, this sequence has density 1/3 in the primes. (End) [Comment rewritten by Jianing Song, Jun 16 2024 and Jun 25 2024]

T. D. Noe, Table of n, a(n) for n=1..1000
C. Ballot and M. Elia, Rank and period of primes in the Fibonacci sequence; a trichotomy, Fib. Quart., 45 (No. 1, 2007), 56-63 (The sequence B3).
Nicholas Bragman and Eric Rowland, Limiting density of the Fibonacci sequence modulo powers of p, arXiv:2202.00704 [math.NT], 2022.
M. Renault, Fibonacci sequence modulo m
Jianing Song, Lucas sequences and entry point modulo p

Cf. A001175, A001177.
Cf. A000204 (Lucas numbers), A001602 (index of the smallest Fibonacci number divisible by prime(n)).
Let {x(n)} be a sequence defined by x(0) = 0, x(1) = 1, x(n+2) = m*x(n+1) + x(n). Let w(k) be the number of zeros in a fundamental period of {x(n)} modulo k.
                             |    m=1    |   m=2   |   m=3
-----------------------------+-----------+---------+---------
The sequence {x(n)}          | A000045   | A000129 | A006190
The sequence {w(k)}          | A001176   | A214027 | A322906
Primes p such that w(p) = 1  | A112860*  | A309580 | A309586
Primes p such that w(p) = 2  | this seq  | A309581 | A309587
Primes p such that w(p) = 4  | A053028** | A261580 | A309588
Numbers k such that w(k) = 1 | A053031   | A309583 | A309591
Numbers k such that w(k) = 2 | A053030   | A309584 | A309592
Numbers k such that w(k) = 4 | A053029   | A309585 | A309593
* and also A053032 (primes dividing Lucas numbers of odd index) U {2}
** also primes dividing no Lucas number

A prime p = prime(i) is in this sequence if p > 2 and A001602(i)/2 is even. - T. D. Noe, Jul 25 2003

A214027 The number of zeros in the fundamental Pisano period of the sequence A000129 mod n.

Original entry on oeis.org

1, 1, 2, 1, 4, 2, 1, 1, 2, 2, 2, 2, 4, 1, 2, 1, 2, 2, 2, 1, 2, 2, 1, 1, 4, 2, 2, 1, 4, 2, 1, 1, 2, 2, 2, 2, 4, 2, 2, 1, 1, 2, 2, 2, 2, 1, 1, 1, 1, 2, 2, 1, 4, 2, 2, 1, 2, 2, 2, 2, 4, 1, 2, 1, 4, 2, 2, 2, 2, 2, 1, 1, 2, 2, 2, 2, 2, 2, 1, 1, 2, 1, 2, 2, 2, 2, 2, 1

Offset: 1

Art DuPre, Jul 04 2012

nonn

This is intimately connected with A175181 and A214028, much as A001176 is intimately connected with A001175 and A001177. In fact, A175181(n)/a(n) = A214028(n). This is the same divisibility relation that holds between A001175, A001176 and A001177.

T. D. Noe, Table of n, a(n) for n = 1..1000

Cf. A175181, A214028.

Similar sequences: A001176, A322906.

Join[{1}, Table[s = t = Mod[{0, 1}, n]; zeros = 0; While[tmp = Mod[2*t[[2]] + t[[1]], n]; t[[1]] = t[[2]]; t[[2]] = tmp; s != t, If[tmp == 0, zeros++]]; zeros, {n, 2, 100}]] (* T. D. Noe, Jul 09 2012 *)

A000129(m) = ([2, 1; 1, 0]^m)[2, 1]
a(n) = my(i=1); while(A000129(i)%n!=0, i++); znorder(Mod(A000129(i+1), n)) \\ Jianing Song, Aug 10 2019

From Jianing Song, Sep 12 2018: (Start)
For odd primes p, a(p^e) = 4 if A214028(p) is odd; 1 if A214028(p) is even but not divisible by 4; 2 if A214028(p) is divisible by 4.
a(n) = 2 for n == 3 (mod 8). For primes p, a(p^e) = 1 if p == 7 (mod 8), 4 if p == 5 (mod 8). Conjecture: 1/6 of the primes congruent to 1 mod 8 satisfy a(p^e) = 1, 2/3 of them satisfy a(p^e) = 2 and 1/6 of them satisfy a(p^e) = 4.
(End)

A053029 Numbers with 4 zeros in Fibonacci numbers mod m.

Original entry on oeis.org

5, 10, 13, 17, 25, 26, 34, 37, 50, 53, 61, 65, 73, 74, 85, 89, 97, 106, 109, 113, 122, 125, 130, 137, 146, 149, 157, 169, 170, 173, 178, 185, 193, 194, 197, 218, 221, 226, 233, 250, 257, 265, 269, 274, 277, 289, 293, 298, 305, 313, 314, 317, 325, 337, 338, 346

Offset: 1

Henry Bottomley, Feb 23 2000

nonn

Conjecture: m is on this list iff m is an odd number all of whose factors are on this list or m is twice such an odd number.

A001176(a(n)) = A128924(a(n),1) = 4. - Reinhard Zumkeller, Jan 17 2014

Reinhard Zumkeller, Table of n, a(n) for n = 1..1000
Brennan Benfield and Michelle Manes, The Fibonacci Sequence is Normal Base 10, arXiv:2202.08986 [math.NT], 2022.
Brennan Benfield and Oliver Lippard, Connecting Zeros in Pisano Periods to Prime Factors of K-Fibonacci Numbers, arXiv:2407.20048 [math.NT], 2024.
M. Renault, Fibonacci sequence modulo m

Let {x(n)} be a sequence defined by x(0) = 0, x(1) = 1, x(n+2) = m*x(n+1) + x(n). Let w(k) be the number of zeros in a fundamental period of {x(n)} modulo k.
                             |   m=1    |   m=2   |   m=3
-----------------------------+----------+---------+---------
The sequence {x(n)}          | A000045  | A000129 | A006190
The sequence {w(k)}          | A001176  | A214027 | A322906
Primes p such that w(p) = 1  | A112860* | A309580 | A309586
Primes p such that w(p) = 2  | A053027  | A309581 | A309587
Primes p such that w(p) = 4  | A053028  | A261580 | A309588
Numbers k such that w(k) = 1 | A053031  | A309583 | A309591
Numbers k such that w(k) = 2 | A053030  | A309584 | A309592
Numbers k such that w(k) = 4 | this seq | A309585 | A309593
* and also A053032 U {2}

a053029 n = a053029_list !! (n-1)
a053029_list = filter ((== 4) . a001176) [1..]
-- Reinhard Zumkeller, Jan 17 2014

A053030 Numbers with 2 zeros in Fibonacci numbers mod m.

Original entry on oeis.org

3, 6, 7, 8, 9, 12, 14, 15, 16, 18, 20, 21, 23, 24, 27, 28, 30, 32, 33, 35, 36, 39, 40, 41, 42, 43, 45, 46, 47, 48, 49, 51, 52, 54, 55, 56, 57, 60, 63, 64, 66, 67, 68, 69, 70, 72, 75, 77, 78, 80, 81, 82, 83, 84, 86, 87, 88, 90, 91, 92, 93, 94, 95, 96, 98, 99, 100, 102, 103, 104

Offset: 1

Henry Bottomley, Feb 23 2000

nonn

m is on this list iff m does not have 1 or 4 zeros in the Fibonacci sequence modulo m.

A001176(a(n)) = A128924(a(n),1) = 2. - Reinhard Zumkeller, Jan 17 2014

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
Brennan Benfield and Oliver Lippard, Connecting Zeros in Pisano Periods to Prime Factors of K-Fibonacci Numbers, arXiv:2407.20048 [math.NT], 2024. See p. 2.
Brennan Benfield and Michelle Manes, The Fibonacci Sequence is Normal Base 10, arXiv:2202.08986 [math.NT], 2022.
M. Renault, Fibonacci sequence modulo m

Let {x(n)} be a sequence defined by x(0) = 0, x(1) = 1, x(n+2) = m*x(n+1) + x(n). Let w(k) be the number of zeros in a fundamental period of {x(n)} modulo k.
                             |   m=1    |   m=2   |   m=3
-----------------------------+----------+---------+---------
The sequence {x(n)}          | A000045  | A000129 | A006190
The sequence {w(k)}          | A001176  | A214027 | A322906
Primes p such that w(p) = 1  | A112860* | A309580 | A309586
Primes p such that w(p) = 2  | A053027  | A309581 | A309587
Primes p such that w(p) = 4  | A053028  | A261580 | A309588
Numbers k such that w(k) = 1 | A053031  | A309583 | A309591
Numbers k such that w(k) = 2 | this seq | A309584 | A309592
Numbers k such that w(k) = 4 | A053029  | A309585 | A309593
* and also A053032 U {2}

a053030 n = a053030_list !! (n-1)
a053030_list = filter ((== 2) . a001176) [1..]
-- Reinhard Zumkeller, Jan 17 2014

A053031 Numbers with 1 zero in Fibonacci numbers mod m.

Original entry on oeis.org

1, 2, 4, 11, 19, 22, 29, 31, 38, 44, 58, 59, 62, 71, 76, 79, 101, 116, 118, 121, 124, 131, 139, 142, 151, 158, 179, 181, 191, 199, 202, 209, 211, 229, 236, 239, 242, 251, 262, 271, 278, 284, 302, 311, 316, 319, 331, 341, 349, 358, 359, 361, 362, 379, 382, 398

Offset: 1

Henry Bottomley, Feb 23 2000

nonn

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.
A001176(a(n)) = A128924(a(n),1) = 1. - Reinhard Zumkeller, Jan 16 2014
Also numbers n such that A001175(n) = A001177(n). - Daniel Suteu, Aug 08 2018

Reinhard Zumkeller, Table of n, a(n) for n = 1..1000
Brennan Benfield and Michelle Manes, The Fibonacci Sequence is Normal Base 10, arXiv:2202.08986 [math.NT], 2022.
Brennan Benfield and Oliver Lippard, Connecting Zeros in Pisano Periods to Prime Factors of K-Fibonacci Numbers, arXiv:2407.20048 [math.NT], 2024.
M. Renault, Fibonacci sequence modulo m

Let {x(n)} be a sequence defined by x(0) = 0, x(1) = 1, x(n+2) = m*x(n+1) + x(n). Let w(k) be the number of zeros in a fundamental period of {x(n)} modulo k.
                             |   m=1    |   m=2   |   m=3
-----------------------------+----------+---------+---------
The sequence {x(n)}          | A000045  | A000129 | A006190
The sequence {w(k)}          | A001176  | A214027 | A322906
Primes p such that w(p) = 1  | A112860* | A309580 | A309586
Primes p such that w(p) = 2  | A053027  | A309581 | A309587
Primes p such that w(p) = 4  | A053028  | A261580 | A309588
Numbers k such that w(k) = 1 | this seq | A309583 | A309591
Numbers k such that w(k) = 2 | A053030  | A309584 | A309592
Numbers k such that w(k) = 4 | A053029  | A309585 | A309593
* and also A053032 U {2}

a053031 n = a053031_list !! (n-1)
a053031_list = filter ((== 1) . a001176) [1..]
-- Reinhard Zumkeller, Jan 16 2014

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 *)

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
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<Charles R Greathouse IV, Dec 14 2016

A112860 2 together with A053032.

Original entry on oeis.org

2, 11, 19, 29, 31, 59, 71, 79, 101, 131, 139, 151, 179, 181, 191, 199, 211, 229, 239, 251, 271, 311, 331, 349, 359, 379, 419, 431, 439, 461, 479, 491, 499, 509, 521, 541, 571, 599, 619, 631, 659, 691, 709, 719, 739, 751, 809, 811, 839, 859, 911, 919, 941, 971

Offset: 1

N. J. A. Sloane, Nov 30 2007

nonn

Consists of the primes that are in neither A053027 nor A053028.
From Jianing Song, Jun 16 2024: (Start)
Primes p such that A001176(p) = 1.
For p > 2, p is in this sequence if and only if A001175(p) == 2 (mod 4), and if and only if A001177(p) == 2 (mod 4). For a proof of the equivalence between A001176(p) = 1 and A001177(p) == 2 (mod 4), see Section 2 of my link below.
This sequence contains all primes congruent to 11, 19 (mod 20). This corresponds to case (3) for k = 3 in the Conclusion of Section 1 of my link below.
Conjecturely, this sequence has density 1/3 in the primes. (End) [Comment rewritten by Jianing Song, Jun 16 2024 and Jun 25 2024]

T. D. Noe, Table of n, a(n) for n=1..1001
C. Ballot and M. Elia, Rank and period of primes in the Fibonacci sequence; a trichotomy, Fib. Quart., 45 (No. 1, 2007), 56-63 (The sequence B1).
Jianing Song, Lucas sequences and entry point modulo p

Cf. A001175, A001177.
Let {x(n)} be a sequence defined by x(0) = 0, x(1) = 1, x(n+2) = m*x(n+1) + x(n). Let w(k) be the number of zeros in a fundamental period of {x(n)} modulo k.
                             |   m=1     |   m=2   |   m=3
-----------------------------+-----------+---------+---------
The sequence {x(n)}          | A000045   | A000129 | A006190
The sequence {w(k)}          | A001176   | A214027 | A322906
Primes p such that w(p) = 1  | this seq* | A309580 | A309586
Primes p such that w(p) = 2  | A053027   | A309581 | A309587
Primes p such that w(p) = 4  | A053028   | A261580 | A309588
Numbers k such that w(k) = 1 | A053031   | A309583 | A309591
Numbers k such that w(k) = 2 | A053030   | A309584 | A309592
Numbers k such that w(k) = 4 | A053029   | A309585 | A309593
* and also A053032 U {2}

A261580 Primes p such that A214028(p) is odd.

Original entry on oeis.org

5, 13, 29, 37, 53, 61, 101, 109, 137, 149, 157, 173, 181, 197, 229, 269, 277, 293, 317, 349, 373, 389, 397, 421, 461, 509, 521, 541, 557, 569, 593, 613, 653, 661, 677, 701, 709, 733, 757, 773, 797, 821, 829, 853, 857, 877, 941, 953, 997, 1013, 1021, 1061, 1069

Offset: 1

Michel Marcus, Aug 25 2015

nonn

From Jianing Song, Aug 13 2019: (Start)
Primes p with 4 zeros in a fundamental period of A000129 mod p, that is, primes p such that A214027(p) = 4. For a proof of the equivalence between A214027(p) = 4 and A214028(p) being odd, see Section 2 of my link below.
For p > 2, p is in this sequence if and only if A175181(p) == 4 (mod 8).
This sequence contains all primes congruent to 5 modulo 8. This corresponds to case (1) for k = 6 in the Conclusion of Section 1 of my link below.
Conjecturely, since (k+2)/2 = 4 is a square, this sequence has density 7/24 in the primes; see the end of Section 1 of my link. (End) [Comment rewritten by Jianing Song, Jun 16 2024 and Jun 20 2024]
The conjecture above is an analog of Hasse's result that the set {p prime : ord(2,p) is odd} has density 7/24 in the primes, where ord(a,m) is the multiplicative order of a modulo m; see A014663. - Jianing Song, Jun 26 2025

			The smallest Pell number divisible by the prime 5 has index 3, which is odd, so 5 is in the sequence.

Jianing Song, Table of n, a(n) for n = 1..1280
Bernadette Faye and Florian Luca, Pell Numbers whose Euler Function is a Pell Number, arXiv:1508.05714 [math.NT], 2015.
Jianing Song, Lucas sequences and entry point modulo p

Cf. A214028, A261581.
Cf. also A175181.
Let {x(n)} be a sequence defined by x(0) = 0, x(1) = 1, x(n+2) = m*x(n+1) + x(n). Let w(k) be the number of zeros in a fundamental period of {x(n)} modulo k.
                             |   m=1    |   m=2    |   m=3
-----------------------------+----------+----------+---------
The sequence {x(n)}          | A000045  | A000129  | A006190
The sequence {w(k)}          | A001176  | A214027  | A322906
Primes p such that w(p) = 1  | A112860* | A309580  | A309586
Primes p such that w(p) = 2  | A053027  | A309581  | A309587
Primes p such that w(p) = 4  | A053028  | this seq | A309588
Numbers k such that w(k) = 1 | A053031  | A309583  | A309591
Numbers k such that w(k) = 2 | A053030  | A309584  | A309592
Numbers k such that w(k) = 4 | A053029  | A309585  | A309593
* and also A053032 U {2}

f[n_] := Block[{k = 1}, While[Mod[Simplify[((1 + Sqrt@ 2)^k - (1 - Sqrt@ 2)^k)/(2 Sqrt@ 2)], n] != 0, k++]; k]; Select[Prime@ Range@ 180, OddQ@ f@ # &] (* Michael De Vlieger, Aug 25 2015 *)

pell(n) = polcoeff(Vec(x/(1-2*x-x^2) + O(x^(n+1))), n);
z(n) = {k=1; while (pell(k) % n, k++); k;}
lista(nn) = {forprime(p=2, nn, if (z(p) % 2, print1(p, ", ")););}

forprime(p=2, 1100, if(A214027(p)==4, print1(p, ", "))) \\ Jianing Song, Aug 13 2019

cp's OEIS Frontend

A172236 Array A(n,k) = n*A(n,k-1) + A(n,k-2) read by upward antidiagonals, starting A(n,0) = 0, A(n,1) = 1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Mathematica

PARI

PARI

SageMath

Formula

Extensions

A053028 Odd primes p with 4 zeros in any period of the Fibonacci numbers mod p.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Formula

Extensions

A053032 Odd primes p with one zero in Fibonacci numbers mod p.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A053027 Odd primes p with 2 zeros in Fibonacci numbers mod p.

Views

Author

Keywords

Comments

Links

Crossrefs

Formula

A214027 The number of zeros in the fundamental Pisano period of the sequence A000129 mod n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A053029 Numbers with 4 zeros in Fibonacci numbers mod m.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

A053030 Numbers with 2 zeros in Fibonacci numbers mod m.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

A053031 Numbers with 1 zero in Fibonacci numbers mod m.

Views

Author