OFFSET
1,1
COMMENTS
The old name was "Primes p that divide Lucas((p+1)/2) = A000032((p+1)/2)".
Note that F(p+1) = F((p+1)/2)*Lucas((p+1)/2), where F = A000045. Since gcd(F(n),Lucas(n)) = 1 or 2 (because Lucas(n)^2 - 5*F(n)^2 = 4*(-1)^n), this sequence (under the old definition above) lists primes p such that p divides F(p+1) but does not divides F((p+1)/2). By Propositions 1.1 and 1.2 (the k = 3 case) of my link below, this is primes p == 3, 7 (mod 20). - Jianing Song, Jun 20 2025
REFERENCES
David A. Cox, "Primes of the Form x^2 + n y^2", Wiley, 1989; see p. 33.
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
Jianing Song, Entry points and periods of Lucas sequences
Eric Weisstein's World of Mathematics, Lucas Number.
Eric Weisstein's World of Mathematics, Gaussian Prime.
MATHEMATICA
Select[Prime[Range[1000]], IntegerQ[(Fibonacci[(#1+1)/2-1]+Fibonacci[(#1+1)/2+1])/#1]&]
Select[Prime[Range[300]], MemberQ[{3, 7}, Mod[#, 20]]&] (* Vincenzo Librandi, Jan 06 2013 *)
PROG
(Magma) [p: p in PrimesUpTo(1500) | p mod 20 in [3, 7]]; // Vincenzo Librandi, Jan 06 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
Alexander Adamchuk, Sep 16 2006
EXTENSIONS
I merged A216816 into this entry at the suggestion of Jianing Song, Jun 20 2025. - N. J. A. Sloane, Jun 22 2025
STATUS
approved
