OFFSET
1,1
COMMENTS
These numbers are rare.
PROG
(PARI) cfracnum(m, f) = { cfr = vector(10000); 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(isprime(numer) & isprime(denom), print1(numer", ")); ) }
CROSSREFS
KEYWORD
nonn,frac,more
AUTHOR
Cino Hilliard, Sep 06 2003
STATUS
approved