This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A033951 #84 Mar 15 2025 04:17:39
%S A033951 1,8,23,46,77,116,163,218,281,352,431,518,613,716,827,946,1073,1208,
%T A033951 1351,1502,1661,1828,2003,2186,2377,2576,2783,2998,3221,3452,3691,
%U A033951 3938,4193,4456,4727,5006,5293,5588,5891,6202,6521,6848,7183,7526,7877,8236,8603,8978
%N A033951 Write 1,2,... in a clockwise spiral; sequence gives numbers on positive x axis.
%C A033951 Ulam's spiral (S spoke of A054552). - _Robert G. Wilson v_, Oct 31 2011
%C A033951 a(n) is the first term in a sum of 2*n + 1 consecutive integers that equals (2*n + 1)^3. - _Patrick J. McNab_, Dec 24 2016
%H A033951 Ivan Panchenko, <a href="/A033951/b033951.txt">Table of n, a(n) for n = 0..1000</a>
%H A033951 Robert G. Wilson v, <a href="/A244677/a244677.jpg">Cover of the March 1964 issue of Scientific American</a>
%H A033951 <a href="/index/Rec#order_03">Index entries for linear recurrences with constant coefficients</a>, signature (3,-3,1).
%F A033951 a(n) = 4*n^2 + 3*n + 1.
%F A033951 G.f.: (1 + 5*x + 2*x^2)/(1-x)^3.
%F A033951 A014848(2n+1) = a(n).
%F A033951 Equals A132774 * [1, 2, 3, ...]; = binomial transform of [1, 7, 8, 0, 0, 0, ...]. - _Gary W. Adamson_, Aug 28 2007
%F A033951 a(n) = A016754(n) - n. - _Reinhard Zumkeller_, May 17 2009
%F A033951 a(n) = a(n-1) + 8*n-1 (with a(0)=1). - _Vincenzo Librandi_, Nov 17 2010
%F A033951 a(0)=1, a(1)=8, a(2)=23, a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). - _Harvey P. Dale_, Feb 07 2015
%F A033951 E.g.f.: exp(x)*(1 + 7*x + 4*x^2). - _Stefano Spezia_, Apr 24 2024
%e A033951 Spiral begins:
%e A033951 .
%e A033951 65--66--67--68--69--70--71--72--73
%e A033951 | |
%e A033951 64 37--38--39--40--41--42--43 74
%e A033951 | | | |
%e A033951 63 36 17--18--19--20--21 44 75
%e A033951 | | | | | |
%e A033951 62 35 16 5---6---7 22 45 76
%e A033951 | | | | | | | |
%e A033951 61 34 15 4 1 8 23 46 77
%e A033951 | | | | | | | |
%e A033951 60 33 14 3---2 9 24 47
%e A033951 | | | | | |
%e A033951 59 32 13--12--11--10 25 48
%e A033951 | | | |
%e A033951 58 31--30--29--28--27--26 49
%e A033951 | |
%e A033951 57--56--55--54--53--52--51--50
%e A033951 From _Aaron David Fairbanks_, Mar 06 2025: (Start)
%e A033951 Illustration of initial terms:
%e A033951 o o o o
%e A033951 o o o o o o o o o
%e A033951 o o o o o o o o o o o o o o o
%e A033951 o o o o o o o o o o o o o o o o o o o o o o
%e A033951 o o o o o o o o o o o o o o o
%e A033951 o o o o o o o o o
%e A033951 o o o o
%e A033951 1 8 23 46
%e A033951 (End)
%p A033951 A033951:=n->4*n^2 + 3*n + 1: seq(A033951(n), n=0..100); # _Wesley Ivan Hurt_, Feb 11 2017
%t A033951 lst={};Do[p=4*n^2+3*n+1;AppendTo[lst, p], {n, 1, 6!}];lst (* _Vladimir Joseph Stephan Orlovsky_, Sep 01 2008 *)
%t A033951 LinearRecurrence[{3,-3,1},{1,8,23},60] (* _Harvey P. Dale_, Feb 07 2015 *)
%t A033951 CoefficientList[Series[(1 + 5 x + 2 x^2)/(1 - x)^3, {x, 0, 45}], x] (* _Michael De Vlieger_, Feb 12 2017 *)
%o A033951 (PARI) a(n)=4*n^2+3*n+1
%o A033951 (Python)
%o A033951 [4*n**2 + 3*n + 1 for n in range(46)] # _Michael S. Branicky_, Jan 08 2021
%Y A033951 Sequences from spirals: A001107, A002939, A007742, A033951, A033952, A033953, A033954, A033989, A033990, A033991, A002943, A033996, A033988.
%Y A033951 Sequences on the four axes of the square spiral: Starting at 0: A001107, A033991, A007742, A033954; starting at 1: A054552, A054556, A054567, A033951.
%Y A033951 Sequences on the four diagonals of the square spiral: Starting at 0: A002939 = 2*A000384, A016742 = 4*A000290, A002943 = 2*A014105, A033996 = 8*A000217; starting at 1: A054554, A053755, A054569, A016754.
%Y A033951 Sequences obtained by reading alternate terms on the X and Y axes and the two main diagonals of the square spiral: Starting at 0: A035608, A156859, A002378 = 2*A000217, A137932 = 4*A002620; starting at 1: A317186, A267682, A002061, A080335.
%Y A033951 Cf. A014848, A132774.
%K A033951 nonn,easy,nice
%O A033951 0,2
%A A033951 Olivier Gorin (gorin(AT)roazhon.inra.fr)
%E A033951 Extended (with formula) by _Erich Friedman_