login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A215907 Odd numbers n such that the Lucas number L(n) is the sum of two squares. 3
1, 3, 7, 13, 19, 31, 37, 43, 49, 61, 67, 73, 79, 91, 111, 127, 163, 169, 183, 199, 223, 307, 313, 349, 361, 397, 433, 511, 523, 541, 613, 619, 709, 823, 907, 1087, 1123, 1129, 1147, 1213, 1279, 1434 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,2
COMMENTS
These Lucas numbers L(n) have no prime factor congruent to 3 mod 4 to an odd power.
Also, numbers n such that L(n) can be written in the form a^2 + 5*b^2.
Subsequence of A124132.
Is this A124132 without the 6? - Joerg Arndt, Sep 07 2012
Any prime factor of Lucas(n) for the prime values of n is always of the form 1 (mod 10) or 9 (mod 10).
A number n can be written in the form a^2 + 5*b^2 if and only if n is 0, or of the form 2^(2i) 5^j Product_{p==1 or 9 mod 20} p^k Product_{q==3 or 7 mod 20) q^(2m) or of the form 2^(2i+1) 5^j Product_{p==1 or 9 mod 20} p^k Product_{q==3 or 7 mod 20) q^(2m+1), for integers i,j,k,m, for primes p,q.
1501 <= a(42) <= 1531. 1531, 1651, 1747, 1849, 1951, 2053, 2413, 2449, 2467, 4069, 5107, 5419, 5851, 7243, 7741, 8467, 13963, 14449, 14887, 15511, 15907, 35449, 51169, 193201, 344293, 387433, 574219, 901657, 1051849 are terms. - Chai Wah Wu, Jul 22 2020
LINKS
Eric W. Weisstein, Sum of squares function
EXAMPLE
Lucas(19) = 9349 = 95^2 + 18^2.
Lucas(19) = 9349 = 23^2 + 5*42^2.
PROG
(PARI) for(i=2, 500, a=factorint(fibonacci(i-1)+fibonacci(i+1))~; has=0; for(j=1, #a, if(a[1, j]%4==3&&a[2, j]%2==1, has=1; break)); if(has==0&&i%2==1, print(i", "))) \\ a^2 + b^2 form.
(PARI) for(i=2, 500, a=factorint(fibonacci(i-1)+fibonacci(i+1))~; flag=0; flip=0; for(j=1, #a, if(((a[1, j]%20>10))&&a[2, j]%2==1, flag=1); if(((a[1, j]%20==2)||(a[1, j]%20==3)||(a[1, j]%20==7))&&a[2, j]%2==1, flip=flip+1)); if(flag==0&&flip%2==0, print(i", "))) \\ a^2 + 5*b^2 form.
CROSSREFS
Cf. A180363.
Cf. A020669, A033205 (numbers and primes of the form x^2 + 5*y^2).
Sequence in context: A086148 A262086 A205956 * A007645 A144919 A215801
KEYWORD
nonn,more
AUTHOR
V. Raman, Aug 26 2012
EXTENSIONS
17 more terms from V. Raman, Aug 28 2012
A215940 merged into this sequence by T. D. Noe, Sep 21 2012
a(38)-a(41) from Chai Wah Wu, Jul 22 2020
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 25 12:53 EDT 2024. Contains 371969 sequences. (Running on oeis4.)