|
%I #16 May 10 2023 09:46:42
%S 0,1,0,-1,-2,-1,0,1,2,2,1,1,0,0,-1,-2,-3,-2,-1,-1,0,0,1,2,3,4,3,2,2,1,
%T 0,-1,-1,-2,-3,-3,-4,-4,-3,-2,-2,-1,0,1,1,2,3,3,4,5,4,3,3,2,1,0,0,-1,
%U -2,-2,-3,-4,-4,-5,-5,-5,-4,-3,-3,-2,-1,0,0,1,2,2,3,4,4,5,5,6
%N a(n) is the x-coordinate of the n-th point in a square spiral mapped to a square grid rotated by Pi/4 using the distance-limited strip bijection described in A307110.
%H Hugo Pfoertner, <a href="/A362955/b362955.txt">Table of n, a(n) for n = 0..3000</a>
%H Hugo Pfoertner, <a href="https://oeis.org/plot2a?name1=A362955&name2=A362956&tform1=untransformed&tform2=untransformed&shift=0&radiop1=xy&drawlines=true">Plot of mapped spiral</a>, using Plot 2.
%o (PARI) \\ ax(n), ay(n) after Kevin Ryde's functions in A174344 and A274923,
%o \\ p(i,j) given in A307110
%o ax(n) = {my(m=sqrtint(n), k=ceil(m/2)); n -= 4*k^2; if(n<0, if(n<-m, k, -k-n), if(n<m, -k, n-3*k))};
%o ay(n) = {my(m=sqrtint(n), k=ceil(m/2)); n -= 4*k^2; if(n<0, if(n<-m, 3*k+n, k), if(n<m, k-n, -k))};
%o p(i, j) = {my(C=cos(Pi/8), S=sin(Pi/8), T=S/C, gx=i*C-j*S, gy=i*S+j*C, k, xm, ym, v=[0, 0]); k=round(gy/C); ym=C*k; xm=gx+(gy-ym)*T; v[1]=round((xm-ym*T)*C); v[2]=round((ym+v[1]*S)/C); v};
%o for (k=0, 81, print1 (p(ax(k),ay(k))[1]", "))
%Y A362956 gives the corresponding y-coordinates.
%Y Cf. A174344, A274923, A307110, A307731.
%K sign
%O 0,5
%A _Hugo Pfoertner_, May 10 2023
|
