%I#15 2022年7月5日11:12:46

%S 1,4,5,6,7,8,9,2,3,14,15,16,17,18,19,20,21,22,23,24,25,10,11,12,13,32，

%电话：33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,26,27,28,29,30,31，

%U 58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80

%N由正方形螺旋线产生的自然数的置换。

%C用正方形逆时针螺旋写出自然数：

%C、。

%C 17-16-15-14-13

%抄送||

%C 18 5---4---3 12

%C|||（C）|

%C 19 6 1----2 11

%抄送|||

%C 20 7--8-9-10

%C类|

%C 21-22-23-24-25

%C、。

%C现在读出一个正方形逆时针螺旋中的数字：1->4->5->6->7->8->9->2->3->14->等等。

%H Eric M.Schmidt，n的表，n=1..1000的a（n）</a>

%t带有[{x=楼层[（楼层[Sqrt[n-1]]+1）/2]}，表[如果[n+2*x<=（2*x+1）^2，n+2*x，n-6*x]，{n，1，75}]（*_G.C.格鲁贝尔，2019年2月5日*）

%o（鼠尾草）

%o定义a（n）：

%o x=（isqrt（n-1）+1）//2

%o如果n+2*x<=（2*x+1）^2，则返回n+2*x，否则返回n-6*x

%o【a（n）代表n in（1..75）】#_Eric M.Schmidt，2016年5月18日

%o（PARI）{s（n）=（（平方（n-1）+1）/2）\1}；

%o表示（n=1,75，打印1（如果（n+2*s（n）<=（2*s（n+1）^2，n+2*s（n），n-6*s（n）），“，”））\\_G.C.Greubel_，2019年2月5日

%Y参见A020703、A090861、A090915、A090928、A09092、A090930。

%K简单，无

%O 1,2号机组

%阿费利克斯·图比亚纳，2004年2月26日

%E偏移由_Eric M.Schmidt修正，2016年5月18日