%我#2022年7月5日17日11:13:49

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

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

%U 74,75,76,77,78,79,81,50,52,53,54,55,56,57,58,59,60,61,62,63,64

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

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

%C、。

%C 17-16-15-14-13

%C ||

%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->8->9->2->3->4->5->6->7->22->等等。

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

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

%o（鼠尾草）

%o定义a（n）：

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

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

%o[a（n）代表n in（1..75）]

%o#_Eric M.Schmidt，2016年5月18日

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

%o表示（n=1,75，打印1（如果（n+6*s（n）<=（2*s（n+1）^2，n+6*s（n），n-2*s

%Y参见A020703、A090861、A090915、A090925、A090%928、A090930。

%K容易，不是

%O 1,2号机组

%2004年2月26日，A _图比亚纳邮件

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