Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Write 1,2,... in a clockwise spiral on a square lattice, writing each digit at a separate lattice point, starting with 1 at the origin and 2 at x=0, y=-1; sequence gives the numbers on the positive x-axis.
26

%I #16 Aug 12 2018 21:26:49

%S 1,8,6,2,3,6,6,0,3,1,8,0,2,7,1,3,9,3,4,3,9,6,0,6,8,9,6,1,2,4,2,1,5,9,

%T 4,1,0,7,7,1,7,8,0,2,6,3,4,2,7,1,8,2,0,2,2,3,5,6,6,3,2,4,1,4,1,5,6,4,

%U 2,9,2,5,5,6,8,5,0,6,4,6,7,9,0,7,6,6,7,7,7,6,4,8,0,9,1,9,5,5,9,1,0,4,0,0,9

%N Write 1,2,... in a clockwise spiral on a square lattice, writing each digit at a separate lattice point, starting with 1 at the origin and 2 at x=0, y=-1; sequence gives the numbers on the positive x-axis.

%C Same as the South spoke of the Champernowne spiral (A244677).

%H Andrey Zabolotskiy, <a href="/A033952/b033952.txt">Table of n, a(n) for n = 1..10000</a>

%e The spiral begins

%e .

%e 3---1---4---1---5

%e | |

%e 1 5---6---7 1

%e | | | |

%e 2 4 1 8 6

%e | | | | |

%e 1 3---2 9 1

%e | | |

%e 1---1---0---1 7

%e .

%t almostNatural[n_, b_] := Block[{m = 0, d = n, i = 1, l, p}, While[m <= d, l = m; m = (b - 1) i*b^(i - 1) + l; i++]; i--; p = Mod[d - l, i]; q = Floor[(d - l)/i] + b^(i - 1); If[p != 0, IntegerDigits[q, b][[p]], Mod[q - 1, b]]];

%t f[n_] := 4n^2 - 5n + 2; Array[ almostNatural[ f@#, 10] &, 105] (* _Robert G. Wilson v_, Aug 08 2014 *)

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

%K nonn,base,easy,nice

%O 1,2

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

%E More terms from Andrew J. Gacek (andrew(AT)dgi.net)

%E Edited by _Charles R Greathouse IV_, Nov 01 2009