%I #34 Feb 11 2024 17:27:15
%S 1,33,97,193,321,481,673,897,1153,1441,1761,2113,2497,2913,3361,3841,
%T 4353,4897,5473,6081,6721,7393,8097,8833,9601,10401,11233,12097,12993,
%U 13921,14881,15873,16897,17953,19041,20161,21313,22497,23713,24961,26241
%N Centered 32-gonal numbers.
%C Sequence found by reading the line from 1, in the direction 1, 33, ..., in the square spiral whose vertices are the generalized decagonal numbers A074377. Semi-axis opposite to A016802 in the same spiral.
%H Vincenzo Librandi, <a href="/A195315/b195315.txt">Table of n, a(n) for n = 1..10000</a>
%H <a href="/index/Rec#order_03">Index entries for linear recurrences with constant coefficients</a>, signature (3,-3,1).
%F a(n) = 16*n^2 - 16*n + 1.
%F G.f.: -x*(1 + 30*x + x^2) / (x-1)^3. - _R. J. Mathar_, Sep 18 2011
%F Sum_{n>=1} 1/a(n) = Pi*tan(sqrt(3)*Pi/4)/(8*sqrt(3)). - _Amiram Eldar_, Feb 11 2022
%t Table[16*n^2 - 16*n + 1, {n, 1, 41}] (* _Amiram Eldar_, Feb 11 2022 *)
%t LinearRecurrence[{3,-3,1},{1,33,97},50] (* _Harvey P. Dale_, Feb 11 2024 *)
%o (Magma) [(16*n^2-16*n+1): n in [1..50]]; // _Vincenzo Librandi_, Sep 19 2011
%o (PARI) a(n)=16*n^2-16*n+1 \\ _Charles R Greathouse IV_, Oct 07 2015
%Y Bisection of A195146.
%Y Cf. A003154, A069129, A069133, A069190, A195314, A195316, A195317, A195318.
%K nonn,easy
%O 1,2
%A _Omar E. Pol_, Sep 16 2011
|