首先我们需要确定我们在哪个周期(距中心的距离)和扇区(北、东、南或西)。然后我们可以确定数字的确切位置。
每个周期的第一个数字如下:1, 9, 25
这是一个二次序列:first(n) = (2n-1)^2 = 4n^2 - 4n + 1
与此相反的是循环数:cycle(i) = floor((sqrt(i) + 1) / 2)
一个周期的长度是:length(n) = first(n+1) - first(n) = 8n
该部门将是:
sector(i) = floor(4 * (i - first(cycle(i))) / length(cycle(i)))
最后,为了得到位置,我们需要从循环和扇区中第一个数字的位置进行推断。
把它们放在一起:
def first(cycle):
x = 2 * cycle - 1
return x * x
def cycle(index):
return (isqrt(index) + 1)//2
def length(cycle):
return 8 * cycle
def sector(index):
c = cycle(index)
offset = index - first(c)
n = length(c)
return 4 * offset / n
def position(index):
c = cycle(index)
s = sector(index)
offset = index - first(c) - s * length(c) // 4
if s == 0: #north
return -c, -c + offset + 1
if s == 1: #east
return -c + offset + 1, c
if s == 2: #south
return c, c - offset - 1
# else, west
return c - offset - 1, -c
def isqrt(x):
"""Calculates the integer square root of a number"""
if x < 0:
raise ValueError('square root not defined for negative numbers')
n = int(x)
if n == 0:
return 0
a, b = divmod(n.bit_length(), 2)
x = 2**(a+b)
while True:
y = (x + n//x)//2
if y >= x:
return x
x = y
例子:
>>> position(9)
(-2, -1)
>>> position(4)
(1, 1)
>>> position(123456)
(-176, 80)