# Greetings from The On-Line Encyclopedia of Integer Sequences! http://oeis.org/ Search: id:a341544 Showing 1-1 of 1 %I A341544 #13 Feb 14 2021 05:53:11 %S A341544 36,256,2916,38416,527076,7311616,101727396,1416468496,19727326116, %T A341544 274760478976,3826898412516,53301739046416,742397156205156, %U A341544 10340257357947136,144021201787572516,2005956552488017936,27939370476391960356,389145229905568604416,5420093847412497929316 %N A341544 a(n) = sqrt( Product_{j=1..n} Product_{k=1..4} (4*sin((2*j-1)*Pi/(2*n))^2 + 4*sin((2*k-1)*Pi/4)^2) ). %H A341544 Index entries for linear recurrences with constant coefficients, signature (19, -76, 76, -19, 1) %F A341544 a(n) = 19*a(n-1) - 76*a(n-2) + 76*a(n-3) - 19*a(n-4) + a(n-5). %F A341544 a(n) = 18*a(n-1) - 58*a(n-2) + 18*a(n-3) - a(n-4) + 144. %F A341544 From _Vaclav Kotesovec_, Feb 14 2021: (Start) %F A341544 G.f.: 4*(4 - 67*x + 197*x^2 - 107*x^3 + 9*x^4) / ((1 - x)*(1 - 14*x + x^2)*(1 - 4*x + x^2)). %F A341544 a(n) = 6 + 4*(2 + sqrt(3))^n + 4*(2 - sqrt(3))^n + (7 + 4*sqrt(3))^n + (7 - 4*sqrt(3))^n. (End) %t A341544 Table[6 + 4 (2 + Sqrt[3])^n + 4 (2 - Sqrt[3])^n + (7 + 4 Sqrt[3])^n + (7 - 4 Sqrt[3])^n, {n, 1, 20}] // FullSimplify (* _Vaclav Kotesovec_, Feb 14 2021 *) %o A341544 (PARI) default(realprecision, 120); %o A341544 a(n) = round(sqrt(prod(j=1, n, prod(k=1, 4, 4*sin((2*j-1)*Pi/(2*n))^2+4*sin((2*k-1)*Pi/4)^2)))); %Y A341544 Column k=4 of A341533. %K A341544 nonn %O A341544 1,1 %A A341544 _Seiichi Manyama_, Feb 14 2021 # Content is available under The OEIS End-User License Agreement: http://oeis.org/LICENSE