A201001 Primes in the Lucas U(3,-1) sequence.
3, 109, 141481, 183642229, 12435170468929559540592573684601
Offset: 1
Keywords
Links
- Wikipedia, Lucas sequence
This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
seq=RecurrenceTable[{a[n]==3*a[n-1]+a[n-2],a[0]==1,a[1]==5},a,{n,1,1000}]; Select[Range[1000],PrimeQ[seq[[#]]]&]
b(n,k)=([n,1;1,0]^k)[1,2] a(n)=for(k=1,2^12,if(ispseudoprime(b(n,k)),return(k)))
Comments