login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A033951 Write 1,2,... in clockwise spiral; sequence gives numbers on positive x axis. 48
1, 8, 23, 46, 77, 116, 163, 218, 281, 352, 431, 518, 613, 716, 827, 946, 1073, 1208, 1351, 1502, 1661, 1828, 2003, 2186, 2377, 2576, 2783, 2998, 3221, 3452, 3691, 3938, 4193, 4456, 4727, 5006, 5293, 5588, 5891, 6202, 6521, 6848, 7183, 7526, 7877, 8236 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,2

COMMENTS

Ulam's spiral (S spoke of A054552). - Robert G. Wilson v, Oct 31 2011

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

LINKS

Ivan Panchenko, Table of n, a(n) for n = 0..1000

Robert G. Wilson v, Cover of the March 1964 issue of Scientific American

Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

FORMULA

a(n) = 4*n^2 + 3*n + 1.

G.f.: (1 + 5*x + 2*x^2)/(1-x)^3.

A014848(2n+1) = a(n).

Equals A132774 * [1, 2, 3, ...]; = binomial transform of [1, 7, 8, 0, 0, 0, ...]. - Gary W. Adamson, Aug 28 2007

a(n) = A016754(n) - n. - Reinhard Zumkeller, May 17 2009

a(n) = a(n-1) + 8*n-1 (with a(0)=1). - Vincenzo Librandi, Nov 17 2010

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

EXAMPLE

Spiral begins:

.

  65--66--67--68--69--70--71--72--73

   |                               |

  64  37--38--39--40--41--42--43  74

   |   |                       |   |

  63  36  17--18--19--20--21  44  75

   |   |   |               |   |   |

  62  35  16   5---6---7  22  45  76

   |   |   |   |       |   |   |   |

  61  34  15   4   1   8  23  46  77

   |   |   |   |   |   |   |   |

  60  33  14   3---2   9  24  47

   |   |   |           |   |   |

  59  32  13--12--11--10  25  48

   |   |                   |   |

  58  31--30--29--28--27--26  49

   |                           |

  57--56--55--54--53--52--51--50

MAPLE

A033951:=n->4*n^2 + 3*n + 1: seq(A033951(n), n=0..100); # Wesley Ivan Hurt, Feb 11 2017

MATHEMATICA

lst={}; Do[p=4*n^2+3*n+1; AppendTo[lst, p], {n, 1, 6!}]; lst (* Vladimir Joseph Stephan Orlovsky, Sep 01 2008 *)

LinearRecurrence[{3, -3, 1}, {1, 8, 23}, 60] (* Harvey P. Dale, Feb 07 2015 *)

CoefficientList[Series[(1 + 5 x + 2 x^2)/(1 - x)^3, {x, 0, 45}], x] (* Michael De Vlieger, Feb 12 2017 *)

PROG

(PARI) a(n)=4*n^2+3*n+1

CROSSREFS

Sequences from spirals: A001107, A002939, A007742, A033951, A033952, A033953, A033954, A033989, A033990, A033991, A002943, A033996, A033988.

Sequences on the four axes of the square spiral: Starting at 0: A001107, A033991, A007742, A033954; starting at 1: A054552, A054556, A054567, 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.

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.

Cf. A132774.

Sequence in context: A226600 A178072 A185257 * A209992 A212118 A175346

Adjacent sequences:  A033948 A033949 A033950 * A033952 A033953 A033954

KEYWORD

nonn,easy,nice

AUTHOR

Olivier Gorin (gorin(AT)roazhon.inra.fr)

EXTENSIONS

Extended (with formula) by Erich Friedman

STATUS

approved

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified October 18 18:10 EDT 2018. Contains 316323 sequences. (Running on oeis4.)