|
%I #45 Sep 27 2019 06:30:49
%S 5,33,869,48449,1504375,124787549,119224713221,10003972882859,
%T 9610660732622149,3811875515172994001,40741092389313992153,
%U 1050927826884051298685761,707754011714996709527574437,184405400463251288725766546203,496687160874729261988243149308101
%N a(n) = numerator(1 + Sum_{k=1..n} n^2 / Product_{j=1..k} 4*j^2).
%H Charles R Greathouse IV, <a href="/A212296/b212296.txt">Table of n, a(n) for n = 1..225</a>
%e r(n) = 5/4, 33/16, 869/256, 48449/9216, 1504375/196608, 124787549/11796480, ....
%e From _Petros Hadjicostas_, Sep 26 2019: (Start)
%e a(3) = numerator(1 + 3^2/(4*1^2) + 3^2/(4*1^2 * 4*2^2) + 3^2/(4*1^2 * 4*2^2 * 4*3^2)) = numerator(1 + 9/4 + 9/64 + 9/2304) = numerator(869/256) = 869.
%e a(4) = numerator(1 + 4^2/(4*1^2) + 4^2/(4*1^2 * 4*2^2) + 4^2/(4*1^2 * 4*2^2 * 4*3^2) + 4^2/(4*1^2 * 4*2^2 * 4*3^2 * 4*4^2)) = numerator(1 + 16/4 + 16/64 + 16/2304 + 16/147456) = denominator(48449/9216) = 48449.
%e (End)
%p a := n -> numer(1 + add(n^2 / mul(4*j^2, j=1..k), k=1..n)):
%p seq(a(n), n=1..15); # _Peter Luschny_, Sep 26 2019
%t G[n_] := Module[{N=1, D=1}, Sum[N*=2*k-1; D*=2*k; (n/D)^2, {k, 1, n}] + 1]; a[n_] := Numerator[G[n]]; Array[a, 15] (* _Jean-François Alcover_, Sep 05 2015, translated from PARI *)
%o (PARI) G(n)=my(N=1,D=1); sum(k=1,n, N*=2*k-1; D*=2*k; (n/D)^2)+1
%o a(n)=numerator(G(n))
%o vector(15, n, a(n))
%Y Denominators are A212297.
%K nonn,frac
%O 1,1
%A _Charles R Greathouse IV_, Jul 02 2012
%E Redefinition according to the data by _Peter Luschny_, Sep 26 2019
|
