This site is supported by donations to The OEIS Foundation. Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A069129 Centered 16-gonal numbers. 32

%I

%S 1,17,49,97,161,241,337,449,577,721,881,1057,1249,1457,1681,1921,2177,

%T 2449,2737,3041,3361,3697,4049,4417,4801,5201,5617,6049,6497,6961,

%U 7441,7937,8449,8977,9521,10081,10657,11249,11857,12481,13121,13777,14449

%N Centered 16-gonal numbers.

%C Also, sequence found by reading the line from 1, in the direction 1, 17,..., in the square spiral whose vertices are the triangular numbers A000217. Opposite numbers to the members of A139098 in the same spiral. - _Omar E. Pol_, Apr 26 2008

%C The subsequence of primes begins: 17, 97, 241, 337, 449, 577, 881, 1249, 3041, 3361, 3697, 4049, 4801, 6961, 7937, 9521, 10657, 13121, 14449. See A184899: n such that the n-th centered 12-gonal number is prime. Indices of prime star numbers. - _Jonathan Vos Post_, Feb 27, 2011

%C Binomial transform of [1, 16, 16, 0, 0, 0,...] and Narayana transform (A001263) of [1, 16, 0, 0, 0,...]. - _Gary W. Adamson_, Jul 28 2011

%C Centered hexadecagonal numbers or centered hexakaidecagonal numbers. - _Omar E. Pol_, Oct 03 2011

%C a(n) = m(n,n) for an array constructed by using the terms in A016813 as the antidiagonals; the first few antidiagonals are 1; 5,9; 13,17,21; 25,29,33,37. - _J. M. Bergot_, Jul 05 2013

%C [The first five rows begin: 1,9,21,37,57; 5,17,33,53,77; 13,29,49,73,101; 25,45,69,97,129; 41,65,93,125,161.]

%H Vincenzo Librandi, <a href="/A069129/b069129.txt">Table of n, a(n) for n = 1..1000</a>

%H Omar E. Pol, <a href="http://polprimos.com">Determinacion geometrica de los numeros primos y perfectos</a>.

%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/CenteredPolygonalNumber.html">Centered Polygonal Numbers</a>

%H <a href="/index/Ce#CENTRALCUBE">Index entries for sequences related to centered polygonal numbers</a>

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

%F a(n) = 8*n^2 - 8*n + 1.

%F a(n) = A035008(n-1) + 1. - _Omar E. Pol_, Apr 26 2008

%F a(n) = 16*n + a(n-1) - 16 with n>1, a(1)=1. - _Vincenzo Librandi_, Aug 08 2010

%F G.f.: -x*(1+14*x+x^2) / (x-1)^3. - _R. J. Mathar_, Feb 04 2011

%F E.g.f.: (8*x^2 + 1)*exp(x). - _G. C. Greubel_, Jul 18 2017

%F a(n) = A056220(2n-1). - _Bruce J. Nicholson_, Aug 31 2017

%F Sum_{n>=1} 1/a(n) = Pi * tan(Pi/(2*sqrt(2))) / (4*sqrt(2)). - _Vaclav Kotesovec_, Jul 23 2019

%e a(5) = 161 because 8*5^2 - 8*5 + 1 = 200 - 40 + 1 = 161.

%t FoldList[#1 + #2 &, 1, 16 Range@ 45] (* _Robert G. Wilson v_, Feb 02 2011 *)

%t Rest[CoefficientList[Series[-x(1+14x+x^2)/(x-1)^3,{x,0,50}],x]] (* _Harvey P. Dale_, Apr 22 2011 *)

%o (MAGMA) [8*n^2-8*n+1: n in [0..50]]; // _Vincenzo Librandi_, Feb 05 2013

%o (PARI) a(n)=8*n^2-8*n+1 \\ _Charles R Greathouse IV_, Sep 24 2015

%Y Cf. A005448, A001844, A005891, A003215, A069099, A000217, A035008, A139098.

%Y Bisection of A077221.

%K easy,nice,nonn

%O 1,2

%A _Terrel Trotter, Jr._, Apr 07 2002

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified September 17 18:50 EDT 2019. Contains 327136 sequences. (Running on oeis4.)