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

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

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

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

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

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

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

%o（鼠尾草）

%o定义a（n）：

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

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

%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）+1）^2，n，8*s（n）^2-n+2），“，”））\\_G.C.Greubel_，2019年2月5日

%Y参见A020703、A090861、A090925、A09092、A0909、A090930。

%K容易，不是

%O 1,2号机组

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

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