1,1

(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);

Cf. A086395.

Sequence in context: A036671 A152171 A132807 this_sequence A042789 A101911 A052109

Adjacent sequences: A101408 A101409 A101410 this_sequence A101412 A101413 A101414

frac,nonn

Cino Hilliard (hillcino368(AT)gmail.com), Jan 15 2005