|
|
A036247
|
|
Numerator of fraction equal to the continued fraction [ 2, 3, 5, ...prime(n) ].
|
|
5
|
|
|
2, 7, 37, 266, 2963, 38785, 662308, 12622637, 290982959, 8451128448, 262275964847, 9712661827787, 398481410904114, 17144413330704689, 806185907954024497, 42744997534894003030, 2522761040466700203267
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
LINKS
|
|
|
FORMULA
|
a(1)=2; a(2)=7; for n>2, a(n)=prime(n)*a(n-1)+a(n-2) [From Zak Seidov, Nov 07 2008]
|
|
EXAMPLE
|
a(2)=7 because 2+1/3 = 7/3. a(3)=37 because 2+1/(3+1/5) = 37/16. a(4)=266 because 2+1/(3+1/(5+1/7)) = 266/115.
|
|
MATHEMATICA
|
a[1]=2; a[2]=7; a[n_]:=Prime[n]*a[n-1]+a[n-2]; Table[a[n], {n, 20}] [From Zak Seidov, Nov 07 2008]
With[{prs=Prime[Range[20]]}, Numerator[Table[FromContinuedFraction[ Take[ prs, n]], {n, 20}]]] (* Harvey P. Dale, Dec 02 2011 *)
|
|
PROG
|
(PARI) a(n)=if(n<0, 0, contfracpnqn(vector(n, i, prime(i)))[1, 1])
|
|
CROSSREFS
|
|
|
KEYWORD
|
frac,nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|