login
A081466
Consider the mapping f(a/b) = (a^2+b^2)/(a^2-b^2) from rationals to rationals. Starting with 2/1 (a=2, b=1) and applying the mapping to each new (reduced) rational number gives 2/1, 5/3, 17/8, 353/225, ... Sequence gives values of the denominators.
2
1, 3, 8, 225, 36992, 6308330625, 21009822254496776192, 3255818067933293622186199316985612890625, 3264008661830516310447364816658205121507617681188862393654856638929469798612992
OFFSET
1,2
PROG
(PARI) {r=2; for(n=1, 9, a=numerator(r); b=denominator(r); print1(b, ", "); r=(a^2+b^2)/(a^2-b^2))}
CROSSREFS
KEYWORD
nonn
AUTHOR
Amarnath Murthy, Mar 22 2003
EXTENSIONS
Edited and extended by Klaus Brockhaus, Mar 24 2003
STATUS
approved