欧拉计划 58

从1开始逆时针旋转，可以得到一个边长为7的螺旋正方形。

37 36 35 34 33 32 31
38 17 16 15 14 13 30
39 18 5 4 3 12 29
40 19 6 1 2 11 28
41 20 7 8 9 10 27
42 21 22 23 24 25 26
43 44 45 46 47 48 49

有趣的是奇数的平方数都处于右下对角线上。更有趣的是，对角线上的13个数字中有8个质数，其百分比是8/13 ≈ 62%。

如果在上面的螺旋正方形上再加一层，可以得到一个边长为9的螺旋正方形。

如果这个过程继续，到螺旋正方形的边长为多少时，对角线上的质数百分比第一次降到10%以下？

def is_prime(x):
    """ x是否质数 """
    x_sqrt = int(pow(x, 0.5))
    for i in range(3, x_sqrt + 1, 2):
        if x % i == 0:
            return False
    return True


i = 1
diagonals_num = 1
prime_num = 0
while 1:
    i += 2
    diagonals_num += 4
    for v in [i ** 2 - (i - 1), i ** 2 - 2 * (i - 1), i ** 2 - 3 * (i - 1)]:
        if is_prime(v):
            prime_num += 1
    if prime_num / diagonals_num < 0.1:
        break
print(i)