求圆上和谐点数

问题描述：

假设一个圆是X**2+Y**2=N，定义在圆上的整数坐标的点为和谐点，例如

N为25时，(±3,±4),(±4,±3), (±5,0), (0,±5) 这12个点为和谐点。

解法：（Python3）

import sys
import math
def find_pointnum(n):
    count = 0
    n=float(n)
    l = int(math.sqrt(n))
    for i in range(l+1):  
        for j in range(i,l+1):  #对1/8象限内做遍历
            if i**2+j**2 == n:
                if i==0 or j==0 or i==j:
                    count+=1  #边界加1
                else:
                    count+=2  #中心加2
    return count*4  #四个象限
  
if __name__ == '__main__':
    n = sys.stdin.readline().strip()
    print (find_pointnum(n))