#来自在线整数序列百科全书的问候!http://oeis.org/搜索:id:a090928显示第1-1页,共1页%一A090928%S A090928 1,6,7,8,9,2,3,4,5,18,19,20,21,22,23,24,25,10,11,12,13,14,15,16,17,38,%电话A090928 39,40,41,42,43,44,45,46,47,48,49,26,27,28,29,30,31,32,33,34,35,36,37,%U A090928 66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,50,51,52,53,54,55,56%N A090928由平方螺旋产生的自然数的排列。%C A090928用正方形逆时针螺旋线写出自然数:%C A090928。%川A090928 17--16--15--14--13%川A090928 ||%川A090928 18 5---4---3 12%C A090928 | | ||%川A090928 19 6 1---2 11%川A090928 | ||%川A090928 20 7---8---9--10%川A090928|%川A090928 21--22--23--24--25%C A090928。%cA090928现在以逆时针螺旋的方式读出数字:1->6->7->8->9->2->3->4->5->18->等。%H A090928埃里克M。施密特,n=1..1000的n,a(n)表%t A090928 With[{x=Floor[(Floor[Sqrt[n-1]]+1)/2]},表[If[n+4*x<=(2*x+1)^2,n+4*x,n-4*x],{n,1,75}]](*\u G。C。格雷贝尔,2019年2月5日*)%o A090928(鼠尾草)%o A090928定义a(n):%o A090928 x=(isqrt(n-1)+1)//2%o A090928如果n+4*x<=(2*x+1)^2否则n-4*x返回n+4*x%o A090928[a(n)表示(1..75)中的n。施密特,2016年5月18日%o A090928(PARI){s(n)=((sqrtint(n-1)+1)/2)\1};%o A090928 for(n=1,75,打印1(如果(n+4*s(n)<=(2*s(n)+1)^2,n+4*s(n),n-4*s(n)),“,”)\\\ u G。C。格雷贝尔,2019年2月5日%Y A090928,参见A020703、A090861、A090915、A090925、A090929、A090930。%K A090928轻松,不%O A090928 1,2%A A090928 Felix Tubiana(fat2(AT)columbia.edu),2004年2月26日%E A090928偏移量由∗Eric M.校正。施密特,2016年5月18日#根据OEIS最终用户许可协议提供内容:http://oeis.org/LICENSE