OFFSET
0,2
COMMENTS
LINKS
Hugo Pfoertner, Table of n, a(n) for n = 0..10000
PROG
(PARI) /* It is assumed that the files a305575 and a305576 contain the second column of the corresponding b-files omitting the initial 0 */
a305575=readvec(a305575); a305576=readvec(a305576);
a174344=vector(10000); L=0; d=1; n=0;
for(r=1, 100, d=-d; k=floor(r/2)*d; for(j=1, L++, a174344[n++]=k); forstep(j=k-d, -floor((r+1)/2)*d+d, -d, a174344[n++]=j));
a274923=vector(10100); L=1; d=1; n=0;
for(r=1, 100, d=-d; k=floor(r/2)*d; for(j=1, L++, a274923[n++]=k); forstep(j=k-d, -floor((r+1)/2)*d+d, -d, a274923[n++]=j));
findinspiral(i, j)={my(size=(2*max(abs(i), abs(j))+1)^2); forstep(k=size, 1, -1, if(i==a174344[k]&&j==a274923[k], return(k)))};
print1(findinspiral(0, 0), ", "); for(n=1, 67, print1(findinspiral(a305575[n], a305576[n]), ", "));
CROSSREFS
KEYWORD
nonn
AUTHOR
Hugo Pfoertner, May 11 2019
STATUS
approved