%I #45 May 08 2021 08:33:31
%S 1,18,960,76440,37437400,157024707840,2777798704721040,
%T 4169982785629476816,331259342780844858796416,
%U 743322803326470921519628462800,163037651356772148158514292729628880,187555796967791569325602741834073910082560,3838658658324493911932517275499048601188128008800
%N a(n) is the denominator of Sum_{i > 0} 1/(Fibonacci(i)*Fibonacci(i+2n)).
%C The numerators are given in A333088.
%C Sum_{i > 0} 1/(Fibonacci(i)*Fibonacci(i+2n)) is a fraction for n > 0.
%C Sum_{i > 0} 1/Fibonacci(i)^2, i.e., the n = 0 case, is known to be transcendental. See A105393.
%C Sum_{i > 0} 1/(Fibonacci(i)*Fibonacci(i+2n-1)) is an irrational number for n > 0. See for instance A290565 (n = 1).
%H A.H.M. Smeets, <a href="/A333089/b333089.txt">Table of n, a(n) for n = 1..62</a>
%H Brother Alfred Brousseau, <a href="https://www.mathstat.dal.ca/FQ/Scanned/7-2/brousseau1.pdf">Summation of Infinite Fibonacci Series</a>, The Fibonacci Quarterly, Vol. 7, No. 2 (1969), pp. 143-168. See (5) and (6) p. 148.
%H Stanley Rabinowitz, <a href="https://www.mathstat.dal.ca/FQ/Scanned/37-2/rabinowitz1.pdf">Algorithmic summation of reciprocals of products of Fibonacci numbers</a>, The Fibonacci Quarterly, Vol. 37 (1999), pp. 122-127. See (23) and (25) p. 5.
%F a(n) = denominator of (1/Fibonacci(2n)) * Sum_{0 < i <= n} 1/(Fibonacci(2i-1)*Fibonacci(2i)).
%e These infinite sums begin: 1, 7/18, 143/960, ...
%t a[n_] := Denominator[Sum[1/(Fibonacci[2i-1]*Fibonacci[2i]),{i,1,n}] / Fibonacci[2n]]; Array[a, 13] (* _Amiram Eldar_, Mar 10 2020 *)
%o (PARI) a(n) = denominator(sum(i=1, n, 1/(fibonacci(2*i-1)*fibonacci(2*i)))/ fibonacci(2*n)); \\ _Michel Marcus_, Mar 10 2020
%o (Python)
%o from math import gcd
%o f0, f1, snum, sden, n = 1, 1, 0, 1, 0
%o while n < 13:
%o snum, sden, n = f0*f1*snum+sden, sden*f0*f1, n+1
%o d = gcd(snum, sden*f0)
%o print(n, sden*f0//d)
%o f0, f1 = 2*f0+f1, f0+f1 # _A.H.M. Smeets_, May 16 2020
%Y Cf. A105393, A290565, A333088 (numerator).
%K nonn,frac
%O 1,2
%A _A.H.M. Smeets_, Mar 07 2020
|