|
%I #3 Oct 01 2013 17:58:06
%S 2,5,12,29,169,408,470832,6625109,44560482149,345869461223138161,
%T 13558774610046711780701,
%U 87275373599917999482560755526644279276078854297832438763690789
%N Denominator when the numerator of the continued fraction rational approximation of sqrt(2) is prime. Also the denominators of A086395(n).
%o (PARI) cfracnum(m,f) = { local(x,r,cfr,i,m1); default(realprecision,3000); numer2=0; denom2=0; cfr = vector(m+100); x=f; for(n=0,m, i=floor(x); x=1/(x-i); cfr[n+1] = i; ); for(m1=0,m, r=cfr[m1+1]; forstep(n=m1,1,-1, r = 1/r; r+=cfr[n]; ); numer=numerator(r); denom=denominator(r); if(ispseudoprime(numer),print1(denom",")); ); default(realprecision,28);
%Y Cf. A086395.
%K frac,nonn
%O 1,1
%A _Cino Hilliard_, Jan 15 2005
|
