login
a(n) = Product_{k=1..floor((n-1)/2)} (1 + 4*cos(2*Pi*k/n)^2).
0

%I #33 Oct 17 2021 13:28:54

%S 1,1,1,2,1,5,4,13,9,34,25,89,64,233,169,610,441,1597,1156,4181,3025,

%T 10946,7921,28657,20736,75025,54289,196418,142129,514229,372100,

%U 1346269,974169,3524578,2550409,9227465,6677056,24157817,17480761,63245986,45765225

%N a(n) = Product_{k=1..floor((n-1)/2)} (1 + 4*cos(2*Pi*k/n)^2).

%H <a href="/index/Rec#order_06">Index entries for linear recurrences with constant coefficients</a>, signature (0,2,0,2,0,-1).

%F Lim_{n->infinity} sqrt(a(n+2)/a(n)) = (sqrt(5) + 1)/2.

%F G.f.: (1+x^6-x^5-3*x^4-x^2+x)/((x^2+1)*(x^2+x-1)*(x^2-x-1)). - Maksym Voznyy (voznyy(AT)mail.ru), Jul 26 2009

%F For n > 0, a(n) = Fibonacci(n) for n odd, and Fibonacci(n/2)^2 for n even. - _Greg Dresden_, Oct 16 2021

%t a[n_] := Product[(1 + 4*Cos[2*Pi*k/n]^2), {k, 1, Floor[(n - 1)/2]}]; a = Table[N[a[n]], {n, 0, 30}]

%t Join[{1}, Table[If[EvenQ[n], Fibonacci[(n)/2]^2, Fibonacci[n]], {n, 1, 30}]] (* _Greg Dresden_, Oct 16 2021 *)

%o (PARI) a(n) = round(prod(k=1, floor((n-1)/2), (1+4*cos(2*Pi*k/n)^2))) \\ _Colin Barker_, Apr 11 2014

%Y Cf. A152189, A152191.

%K nonn,easy

%O 0,4

%A _Roger L. Bagula_ and _Gary W. Adamson_, Nov 28 2008

%E More terms and edited by _Colin Barker_, _Michel Marcus_, and _Joerg Arndt_, Apr 11 2014