login
a(n) = (6*n+1)*(6*n+5).
2

%I #32 Feb 19 2023 03:33:24

%S 5,77,221,437,725,1085,1517,2021,2597,3245,3965,4757,5621,6557,7565,

%T 8645,9797,11021,12317,13685,15125,16637,18221,19877,21605,23405,

%U 25277,27221,29237,31325,33485,35717,38021,40397,42845,45365,47957,50621,53357,56165,59045

%N a(n) = (6*n+1)*(6*n+5).

%H T. D. Noe, <a href="/A001513/b001513.txt">Table of n, a(n) for n = 0..1000</a>

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

%F Sum_{k>=0} 1/a(k) = Pi/(8*sqrt(3)) = 0.22672492... - _Jaume Oliver Lafont_, May 30 2010

%F a(n) = 72*n + a(n-1) with a(0)=5. - _Vincenzo Librandi_, Nov 12 2010

%F G.f.: (-5 - 62*x - 5*x^2) / (x-1)^3. - _R. J. Mathar_, Jan 19 2013

%F From _Amiram Eldar_, Feb 19 2023: (Start)

%F a(n) = A016921(n)*A016969(n).

%F Sum_{n>=0} (-1)^n/a(n) = log(2+sqrt(3))/(4*sqrt(3)).

%F Product_{n>=0} (1 - 1/a(n)) = 2*cos(sqrt(5)*Pi/6).

%F Product_{n>=0} (1 + 1/a(n)) = 2*cos(sqrt(3)*Pi/6). (End)

%t a[n_] := (6*n + 1)*(6*n + 5); Array[a, 40, 0] (* _Amiram Eldar_, Feb 19 2023 *)

%o (PARI) a(n)=(6*n+1)*(6*n+5) \\ _Charles R Greathouse IV_, Jun 17 2017

%Y Cf. A016921, A016969.

%K nonn,easy

%O 0,1

%A _N. J. A. Sloane_