A001578 -id:A001578 - cp's OEIS frontend

A001605 Indices of prime Fibonacci numbers.

3, 4, 5, 7, 11, 13, 17, 23, 29, 43, 47, 83, 131, 137, 359, 431, 433, 449, 509, 569, 571, 2971, 4723, 5387, 9311, 9677, 14431, 25561, 30757, 35999, 37511, 50833, 81839, 104911, 130021, 148091, 201107, 397379, 433781, 590041, 593689, 604711, 931517, 1049897, 1285607, 1636007, 1803059, 1968721, 2904353, 3244369, 3340367

Offset: 1

N. J. A. Sloane

Some of the larger entries may only correspond to probable primes.
Since F(n) divides F(mn) (cf. A001578, A086597), all terms of this sequence are primes except for a(2) = 4 = 2 * 2 but F(2) = 1. - M. F. Hasler, Dec 12 2007
What is the next larger twin prime after F(4) = 3, F(5) = 5, F(7) = 13? The next candidates seem to be F(104911) or F(1968721) (greater of a pair), or F(397379), F(931517) (lesser of a pair). - M. F. Hasler, Jan 30 2013, edited Dec 24 2016, edited Sep 23 2017 by Bobby Jacobs
_Henri Lifchitz_ confirms that the data section gives the full list (49 terms) as far as we know it today of indices of prime Fibonacci numbers (including proven primes and PRPs). - N. J. A. Sloane, Jul 09 2016
Terms n such that n-2 is also a term are listed in A279795. - M. F. Hasler, Dec 24 2016
There are no Fibonacci numbers that are twin primes after F(7) = 13. Every Fibonacci prime greater than F(4) = 3 is of the form F(2*n+1). Since F(2*n+1)+2 and F(2*n+1)-2 are F(n+2)*L(n-1) and F(n-1)*L(n+2) in some order, and F(n+2) > 1, L(n-1) > 1, F(n-1) > 1, and L(n+2) > 1 for n > 3. - Bobby Jacobs, Sep 23 2017
These primes are occurring with about the same normalized frequency as Repunit primes (see Generalized Repunit Conjecture Ref).  Assuming a base=1.618 (ratio of sequential terms), then the best fit coefficient is 0.60324 for the first 56 terms, which is already approaching Euler's constant 0.56145948. - Paul Bourdelais, Aug 23 2024

Clifford A. Pickover, Mazes for the Mind, St. Martin's Press, NY, 1992, p. 350.
Clifford A. Pickover, A Passion for Mathematics, Wiley, 2005; see p. 54.
Paulo Ribenboim, The Little Book of Big Primes, Springer-Verlag, NY, 1991, p. 178.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Paul Bourdelais, Table of n, a(n) for n = 1..56 (first 51 terms from Henri Lifchitz)
Paul Bourdelais, A Generalized Repunit Conjecture
John Brillhart, Peter L. Montgomery and Robert D. Silverman, Tables of Fibonacci and Lucas factorizations, Math. Comp. 50 (1988), 251-260.
David Broadhurst, Fibonacci Numbers
David Broadhurst, Proof that F(81839) is prime, NMBRTHRY Mailing List, Apr 22 2001.
Chris K. Caldwell, The Prime Glossary, Fibonacci prime
Rosina Campbell, Duc Van Huynh, Tyler Melton, and Andrew Percival, Elliptic Curves of Fibonacci order over F_p, arXiv:1710.05687 [math.NT], 2017.
Harvey Dubner and Wilfrid Keller, New Fibonacci and Lucas Primes, Math. Comp. 68 (1999) 417-427.
Dudley Fox, Search for Possible Fibonacci Primes
Shyam Sunder Gupta, Fabulous Fibonacci Numbers, Lucas Numbers, and Golden Ratio, Exploring the Beauty of Fascinating Numbers, Springer (2025) Ch. 8, 223-274.
Dov Jarden, Recurring Sequences, Riveon Lematematika, Jerusalem, 1966. [Annotated scanned copy] See p. 36.
Ron Knott, Mathematics of the Fibonacci Series
Alex Kontorovich and Jeff Lagarias, On Toric Orbits in the Affine Sieve, arXiv:1808.03235 [math.NT], 2018.
Henri & Renaud Lifchitz, PRP Records.
Tony D. Noe and Jonathan Vos Post, Primes in Fibonacci n-step and Lucas n-step Sequences, J. of Integer Sequences, Vol. 8 (2005), Article 05.4.4.
Rudolf Ondrejka, The Top Ten: a Catalogue of Primal Configurations
PRP Top Records, Search for: F(n)
Lawrence Somer and Michal Křížek, On Primes in Lucas Sequences, Fibonacci Quart. 53 (2015), no. 1, 2-23.
Eric Weisstein's World of Mathematics, Fibonacci Prime
Eric Weisstein's World of Mathematics, Integer Sequence Primes

Cf. A000045, A001578, A005478, A080345, A086597, A117595.
Subsequence of A046022.
Column k=1 of A303215.

Select[Range[10^4], PrimeQ[Fibonacci[#]] &] (* Harvey P. Dale, Nov 20 2012 *)
(* Start ~ 1.8x faster than the above *)
Select[Range[10^4], PrimeQ[#] && PrimeQ[Fibonacci[#]] &] (* Eric W. Weisstein, Nov 07 2017 *)
Select[Prime[Range[PrimePi[10^4]]], PrimeQ[Fibonacci[#]] &] (* Eric W. Weisstein, Nov 07 2017 *)
(* End *)

v=[3,4]; forprime(p=5,1e5, if(ispseudoprime(fibonacci(p)), v=concat(v,p))); v \\ Charles R Greathouse IV, Feb 14 2011

is_A001605(n)={n==4 || isprime(n) & ispseudoprime(fibonacci(n))}  \\ M. F. Hasler, Sep 29 2012

Prime(i) = a(n) for some n <=> A080345(i) <= 1. - M. F. Hasler, Dec 12 2007

Additional comments from Robert G. Wilson v, Aug 18 2000
More terms from David Broadhurst, Nov 08 2001
Two more terms (148091 and 201107) from T. D. Noe, Feb 12 2003 and Mar 04 2003
397379 from T. D. Noe, Aug 18 2003
433781, 590041, 593689 from Henri Lifchitz submitted by Ray Chandler, Feb 11 2005
604711 from Henri Lifchitz communicated by Eric W. Weisstein, Nov 29 2005
931517, 1049897, 1285607 found by Henri Lifchitz circa Nov 01 2008 and submitted by Alexander Adamchuk, Nov 28 2008
1636007 from Henri Lifchitz March 2009, communicated by Eric W. Weisstein, Apr 24 2009
1803059 and 1968721 from Henri Lifchitz, November 2009, submitted by Alex Ratushnyak, Aug 08 2012
a(49)=2904353 from Henri Lifchitz, Jul 15 2014
a(50)=3244369 from Henri Lifchitz, Nov 04 2017
a(51)=3340367 from Henri Lifchitz, Apr 25 2018
a(52)-a(56) from Ryan Propper added by Paul Bourdelais, Aug 23 2024

A058036 Smallest primitive prime factor of the n-th Lucas number (A000032); i.e., L(n), L(0) = 2, L(1) = 1 and L(n) = L(n-1) + L(n-2).

Original entry on oeis.org

2, 1, 3, 1, 7, 11, 1, 29, 47, 19, 41, 199, 23, 521, 281, 31, 2207, 3571, 107, 9349, 2161, 211, 43, 139, 1103, 101, 90481, 5779, 14503, 59, 2521, 3010349, 1087, 9901, 67, 71, 103681, 54018521, 29134601, 79, 1601, 370248451, 83, 6709, 263, 181, 4969

Offset: 0

Robert G. Wilson v, Nov 16 2000

nonn

A Lucas number can have more than one primitive factor; the primitive factors of L(22) are 43 and 307.

T. D. Noe, Table of n, a(n) for n = 0..1000 (using Blair Kelly's data).
Mansur S. Boase, A Result About the Primes Dividing Fibonacci Numbers, The Fibonacci Quarterly, 39.5 (2001) 386.
J. Brillhart, P. L. Montgomery and R. D. Silverman, Tables of Fibonacci and Lucas factorizations, Math. Comp. 50 (1988), 251-260, S1-S15. Math. Rev. 89h:11002.
Blair Kelly, Fibonacci and Lucas Factorizations

Cf. A000032, A086600 (number of primitive prime factors in L(n)).

Cf. A001578 (analog for Fibonacci).

a=3; b=-1; prms={}; Table[c=a+b; a=b; b=c; f=First/@FactorInteger[c]; p=Complement[f, prms]; prms=Join[prms, p]; If[p=={}, 1, First[p]], {47}]

lucas(n) = fibonacci(n+1)+fibonacci(n-1); \\ A000032
a(n) = {n++; my(v = vector(n, k, k--; lucas(k))); my(vf = vector(n, k, factor(v[k])[,1]~)); for (k=1, n-1, vf[n] = setminus(vf[n], vf[k]);); if (#vf[n], vecmin(vf[n]), 1);} \\ Michel Marcus, May 11 2021

A152012 Indices of Fibonacci numbers having exactly one primitive prime factor.

Original entry on oeis.org

3, 4, 5, 7, 8, 9, 10, 11, 13, 14, 15, 16, 17, 18, 20, 21, 22, 23, 24, 25, 26, 28, 29, 30, 32, 33, 34, 35, 36, 38, 39, 40, 42, 43, 45, 47, 48, 51, 52, 54, 56, 60, 62, 63, 65, 66, 72, 74, 75, 76, 82, 83, 93, 94, 98, 105, 106, 108, 111, 112, 119, 121, 122, 123, 124, 125, 131

Offset: 1

Max Alekseyev, Nov 19 2008

nonn

It is known that Fibonacci number A000045(n) has a primitive prime factor for all n, except n=0, 1, 2, 6 and 12. This sequence lists such indices n that A000045(n) has exactly one primitive prime factor (equal A001578(n)). Sister sequence A152013 provides indices of Fibonacci numbers with at least 2 prime factors. The current sequence A152012 and its sister sequence A152013 along with the finite set {0,1,2,6,12} form a partition of the natural numbers.
Numbers k such that A086597(k) = 1.
For prime p, all prime factors of Fibonacci(p) are primitive. Hence, the only primes in this sequence are the prime numbers in A001605, which gives the indices of prime Fibonacci numbers.

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

Cf. A000045, A001578, A152013.

primitivePrimeFactors[n_] := Cases[FactorInteger[Fibonacci[n]][[All, 1]], p_ /; And @@ (GCD[p, #] == 1 & /@ Array[Fibonacci, n-1])]; Reap[For[n=3, n <= 200, n++, If[Length[primitivePrimeFactors[n]] == 1, Print[n]; Sow[n]]]][[2, 1]] (* Jean-François Alcover, Dec 12 2014 *)

isok(pf, vp) = sum(i=1, #pf, vecsearch(vp, pf[i]) == 0) == 1;
lista(nn) = {vp = []; for (n=3, nn, pf = factor(fibonacci(n))[,1]; if (isok(pf, vp), print1(n, ", ")); vp = vecsort(concat(vp, pf),, 8););} \\ Michel Marcus, Nov 29 2014

A178763 Product of primitive prime factors of Fibonacci(n).

Original entry on oeis.org

1, 1, 2, 3, 5, 1, 13, 7, 17, 11, 89, 1, 233, 29, 61, 47, 1597, 19, 4181, 41, 421, 199, 28657, 23, 3001, 521, 5777, 281, 514229, 31, 1346269, 2207, 19801, 3571, 141961, 107, 24157817, 9349, 135721, 2161, 165580141, 211, 433494437, 13201, 109441, 64079

Offset: 1

T. D. Noe, Jun 10 2010

nonn

Same as A001578 for the first 18 terms.

Let b(n) be the greatest divisor of Fibonacci(n) that is coprime to Fibonacci(m) for all positive integers m < n, then a(n) = b(n) for all n, provided that no Wall-Sun-Sun prime exists. Otherwise, if p is a Wall-Sun-Sun prime and A001177(p) = k (then A001177(p^2) = k), then p^2 divides b(k), but by definition a(k) is squarefree. - Jianing Song, Jul 02 2019

T. D. Noe, Table of n, a(n) for n = 1..1000
Dov Jarden, Recurring Sequences, Riveon Lematematika, Jerusalem, 1966. [Annotated scanned copy] See pp. 60-64 for a table of the first 385 terms.
Florian Luca, Carl Pomerance, and Stephen Wagner, Fibonacci Integers (preprint)

Cf. A061446, A086597, A152012 (Indices of prime terms).

a(n)=my(d=divisors(n), f=fibonacci(n), t); t=lcm(apply(fibonacci,d[1..#d-1])); while((t=gcd(t,f))>1, f/=t); f \\ Charles R Greathouse IV, Nov 30 2016

a(n) = A061446(n) / A178764(n).

a(n) = A061446(n) / gcd(A061446(n), n) if n != 5, 6, provided that no Wall-Sun-Sun prime exists. - Jianing Song, Jul 02 2019

A246556 a(n) = smallest prime which divides Pell(n) = A000129(n) but does not divide any Pell(k) for k

Original entry on oeis.org

2, 5, 3, 29, 7, 13, 17, 197, 41, 5741, 11, 33461, 239, 269, 577, 137, 199, 37, 19, 45697, 23, 229, 1153, 1549, 79, 53, 113, 44560482149, 31, 61, 665857, 52734529, 103, 1800193921, 73, 593, 9369319, 389, 241, 1746860020068409, 4663, 11437, 43, 6481, 47, 3761, 97, 293, 45245801, 101, 22307, 68480406462161287469, 7761799, 109, 1535466241

Offset: 2

Eric Chen, Nov 15 2014

nonn

First differs from A264137 (Largest prime factor of the n-th Pell number) at n=17; see Example section. - Jon E. Schoenfield, Dec 10 2016

			a(2) = 2 because Pell(2) = 2 and Pell(k) < 2 for k < 2.
a(4) = 3 because Pell(4) = 12 = 2^2 * 3, but 2 is not a primitive prime factor since Pell(2) = 2, so therefore 3 is the primitive prime factor.
a(5) = 29 because Pell(5) = 29, which is prime.
a(6) = 7 because Pell(6) = 70 = 2 * 5 * 7, but neither 2 nor 5 is a primitive prime factor, so therefore 7 is the primitive prime factor.
a(17) = 137 because Pell(17) = 1136689 = 137 * 8297, and both of them are primitive factors, we choose the smallest. (Pell(17) is the smallest Pell number with more than one primitive prime factor.)

Table of n, a(n) for n = 2..630
R. D. Carmichael, On the numerical factors of the arithmetic forms α^n ± β^n, Annals of Math., 15 (1/4) (1913), 30-70.

Cf. A001578 (for Fibonacci(n)), A000129 (Pell numbers), A008555, A086383, A096650, A120947, A175181, A214028, A264137.

prms={}; Table[f=First/@FactorInteger[Pell[n]]; p=Complement[f, prms]; prms=Join[prms, p]; If[p=={}, 1, First[p]], {n, 36}]

a(n) >= 2 for all n >= 2, by Carmichael's theorem. - Jonathan Sondow, Dec 08 2017

Edited by N. J. A. Sloane, Nov 29 2014
Terms up to a(612) in b-file added by Sean A. Irvine, Sep 23 2019
Terms a(613)-a(630) in b-file added by Max Alekseyev, Aug 26 2021

A061488 Factorize the Fibonacci numbers in order, skipping F(0)-F(2), F(6)=8 and F(12)=144; at each step at least one new prime will occur; append to the sequence the smallest such new prime.

Original entry on oeis.org

2, 3, 5, 13, 7, 17, 11, 89, 233, 29, 61, 47, 1597, 19, 37, 41, 421, 199, 28657, 23, 3001, 521, 53, 281, 514229, 31, 557, 2207, 19801, 3571, 141961, 107, 73, 9349, 135721, 2161, 2789, 211, 433494437, 43, 109441, 139, 2971215073, 1103, 97, 101, 6376021

Offset: 3

N. J. A. Sloane, Nov 08 2001

Carmichael showed that the sequence is well defined.
Same as A001578 without the "1" terms.
Given the definition, in particular omission of F(6) and F(12), setting offset=1 would be more adequate; offset=5 (= number of omitted terms) would give A001578 for n > 12 on. - M. F. Hasler, Oct 21 2012

T. D. Noe, Table of n, a(n) for n=3..998
Ron Knott, Mathematics of the Fibonacci Series
Hisanori Mishima, Factorizations of many number sequences
Hisanori Mishima, Factorizations of many number sequences
Hisanori Mishima, Factorizations of many number sequences
Hisanori Mishima, Factorizations of many number sequences
Hisanori Mishima, Factorizations of many number sequences

Cf. A000045, A061446.

f[n_] := Block[{p = First /@ FactorInteger[ Fibonacci[ n]]}, k = 1; lmt = 1 + Length@ p; While[k < lmt && MemberQ[lst, p[[k]]], k++]; If[k < lmt, AppendTo[lst, p[[k]]]]]; lst = {}; Do[ f[n], {n, 3, 51}]; lst (* Robert G. Wilson v, Oct 23 2012 *)

a(n) = A001578(n+2) from n=11 on. - M. F. Hasler, Oct 21 2012

More terms from Vladeta Jovovic and Lior Manor, Nov 09 2001

Corrected by T. D. Noe, Feb 10 2007

A152013 Indices of Fibonacci numbers having at least two distinct primitive prime factors.

Original entry on oeis.org

19, 27, 31, 37, 41, 44, 46, 49, 50, 53, 55, 57, 58, 59, 61, 64, 67, 68, 69, 70, 71, 73, 77, 78, 79, 80, 81, 84, 85, 86, 87, 88, 89, 90, 91, 92, 95, 96, 97, 99, 100, 101, 102, 103, 104, 107, 109, 110, 113, 114, 115, 116, 117, 118, 120, 126, 127, 128, 129, 130, 133, 134

Offset: 1

Max Alekseyev, Nov 19 2008

nonn

It is known that Fibonacci number A000045(n) has a primitive prime factor for all n, except n=0, 1, 2, 6 and 12. This sequence lists such indices n that A000045(n) has at least two primitive prime factors. Its sister sequence A152012 provides indices of Fibonacci numbers with exactly one primitive prime factor. The current sequence A152013 and its sister sequence A152012 along with the finite set {0,1,2,6,12} form a partition of the natural numbers.

Cf. A000045, A001578, A152012.

A119997 Sum of all matrix elements of n X n matrix M[i,j] = (-1)^(i+j)*Fibonacci[i+j-1].

Original entry on oeis.org

1, 1, 4, 5, 17, 32, 97, 225, 628, 1573, 4225, 10880, 28769, 74849, 196708, 513765, 1347025, 3523360, 9229441, 24154625, 63251156, 165571781, 433507969, 1134881280, 2971250497, 7778684737, 20365103812, 53316141125, 139584105233, 365434903328, 956722661665

Offset: 1

Alexander Adamchuk, Aug 03 2006

Prime p divides a(p-1) for p={5,11,19,29,31,41,59,61,71,...} = A038872[n] Primes congruent to {0, 1, 4} mod 5. Also odd primes where 5 is a square mod p. p^2 divides a(p-1) for prime p={11,19,29,31,41,59,61,71,...} = A045468[n] Primes congruent to {1, 4} mod 5. Square prime divisors of a(n) up to n=50 are{2,3,5,7,11,13,19,23,29,31,41,47,89,101,139,151,199,211,461,521,3571,9349}, It appears that square prime divisors of a(n) belong to A061446[n] Primitive part of Fibonacci(n), A001578[n] Smallest primitive prime factor of Fibonacci number F(n) and A072183[n] Sequence arising from factorization of the Fibonacci numbers. Sum[Sum[Fibonacci[i+j-1],{i,1,n}],{j,1,n}] = A120297[n]. Sum[Sum[i+j-1,{i,1,n}],{j,1,n}] = n^3. Sum[Sum[(-1)^(i+j)*(i+j-1),{i,1,n}],{j,1,n}] = n for odd n and = 0 for even n.

			Matrix begins:
1 -1 2 -3 5
-1 2 -3 5 -8
2 -3 5 -8 13
-3 5 -8 13 -21
5 -8 13 -21 34

Colin Barker, Table of n, a(n) for n = 1..1000 [Terms up to n=200 from _Vincenzo Librandi_]
Index entries for linear recurrences with constant coefficients, signature (3,1,-7,5,-1).

Cf. A120297, A000045, A038872, A001924, A062381, A038872, A045468, A061446, A001578, A072183.

Table[Sum[Sum[(-1)^(i+j)*Fibonacci[i+j-1],{i,1,n}],{j,1,n}],{n,1,50}]

a(n) = sum(i=1, n, sum(j=1, n, (-1)^(i+j)*fibonacci(i+j-1))) \\ Colin Barker, Mar 26 2015

Vec(-x*(x^3+2*x-1)/((x-1)*(x^2-3*x+1)*(x^2-x-1)) + O(x^100)) \\ Colin Barker, Mar 26 2015

a(n) = Sum[Sum[(-1)^(i+j)*Fibonacci[i+j-1],{i,1,n}],{j,1,n}].
a(n) = 3*a(n-1)+a(n-2)-7*a(n-3)+5*a(n-4)-a(n-5) for n>5. - Colin Barker, Mar 26 2015
G.f.: -x*(x^3+2*x-1) / ((x-1)*(x^2-3*x+1)*(x^2-x-1)). - Colin Barker, Mar 26 2015

A121708 Numerator of Sum/Product of first n Fibonacci numbers A000045[n].

Original entry on oeis.org

1, 2, 2, 7, 2, 1, 11, 3, 11, 1, 29, 47, 29, 1, 19, 41, 19, 1, 199, 23, 199, 1, 521, 281, 521, 1, 31, 2207, 31, 1, 3571, 107, 3571, 1, 9349, 2161, 9349, 1, 211, 13201, 211, 1, 64079, 1103, 64079, 1, 15251, 90481, 15251, 1, 5779, 14503, 5779, 1, 1149851, 2521

Offset: 1

Alexander Adamchuk, Aug 16 2006, Sep 21 2006

a(1) = 1 and a(4k+2) = 1 for k>0.

For k >1 a(4k-1) = a(4k+1) = A072183(2k+1) = A061447(2k+1) Primitive part of Lucas(n).

Cf. A000045, A000032, A079451, A001578, A096362, A058036, A121709, A090585.

Cf. A061447, A072183.

Table[Numerator[Sum[Fibonacci[k],{k,1,n}]/Product[Fibonacci[k],{k,1,n}]],{n,1,100}]
With[{fibs=Fibonacci[Range[60]]},Numerator[Accumulate[fibs]/Rest[ FoldList[ Times,1,fibs]]]] (* This is significantly faster than the first program above *) (* Harvey P. Dale, Aug 19 2012 *)

a(n) = numerator( sum(k=1..n, Fibonacci(k)) / prod(k=1..n, Fibonacci(k)) ).

cp's OEIS Frontend

A001605 Indices of prime Fibonacci numbers.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Formula

Extensions

A058036 Smallest primitive prime factor of the n-th Lucas number (A000032); i.e., L(n), L(0) = 2, L(1) = 1 and L(n) = L(n-1) + L(n-2).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

A152012 Indices of Fibonacci numbers having exactly one primitive prime factor.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

A178763 Product of primitive prime factors of Fibonacci(n).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

Formula

A246556 a(n) = smallest prime which divides Pell(n) = A000129(n) but does not divide any Pell(k) for k

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

Extensions

A061488 Factorize the Fibonacci numbers in order, skipping F(0)-F(2), F(6)=8 and F(12)=144; at each step at least one new prime will occur; append to the sequence the smallest such new prime.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Formula

Extensions

A152013 Indices of Fibonacci numbers having at least two distinct primitive prime factors.

Views

Author

Keywords

Comments

Crossrefs

A119997 Sum of all matrix elements of n X n matrix M[i,j] = (-1)^(i+j)*Fibonacci[i+j-1].

Views

Author

Keywords

Comments

Examples