|
|
A086383
|
|
Prime terms in the sequence of Pell numbers, A000129.
|
|
9
|
|
|
2, 5, 29, 5741, 33461, 44560482149, 1746860020068409, 68480406462161287469, 13558774610046711780701, 4125636888562548868221559797461449, 4760981394323203445293052612223893281
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
Previous Name: Primes found among the denominators of the continued fraction rational approximations to sqrt(2).
|
|
LINKS
|
|
|
FORMULA
|
|
|
EXAMPLE
|
|
|
MATHEMATICA
|
Select[Table[ChebyshevU[k, 3]-ChebyshevU[k-1, 3], {k, 0, 50}], PrimeQ] (* Ed Pegg Jr, May 10 2007 *)
Select[Denominator[Convergents[Sqrt[2], 150]], PrimeQ] (* Harvey P. Dale, Dec 19 2012 *)
Select[LinearRecurrence[{2, 1}, {0, 1}, 16], PrimeQ] (* Zak Seidov, Oct 21 2013 *)
|
|
PROG
|
(PARI) \\ Continued fraction rational approximation of numeric constants f. m=steps.
cfracdenomprime(m, f) = { default(realprecision, 3000); cf = vector(m+10); x=f; for(n=0, m, i=floor(x); x=1/(x-i); cf[n+1] = i; ); for(m1=0, m, r=cf[m1+1]; forstep(n=m1, 1, -1, r = 1/r; r+=cf[n]; ); numer=numerator(r); denom=denominator(r); if(ispseudoprime(denom), print1(denom, ", ")); ) }
(GAP) f:=[0, 1];; for n in [3..100] do f[n]:=2*f[n-1]+f[n-2]; od; a:=Filtered(f, IsPrime);; Print(a); # Muniru A Asiru, Jan 03 2019
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|