Большие числа

n=int(input())
a=0
d=[]
c=[0]
f=[]
o=0
for i in range(1,n+1):
    for x in range(-n,n+1):
        for y in range(-n,n+1):
            if x<=n:
                if y<=n:
                    if x**2 + y**2==i**2:
                        a+=1 
                        y+=1 
                    else:
                        y+=1 
                else:
                    x+=1 
            else:
                i+=1 
    for j in range(1):
        d.append(a)
    for h in range(1):
        c.append(a)
    f=[d-c for d,c in zip(d,c)]
    o=f.index(max(f))
print(o+1)

from math import sqrt
 
def isSquare(n):
    q=int(sqrt(n))
    if q*q==n or (q+1)*(q+1)==n:
        return True
    else:
        return False
 
def task(n):
    ma=2
    rr=1
    for r in range(2,n+1):
        s=2
        for x in range(1,r+1):
            k=r**2-x**2
            if isSquare(k):
                s+=1
        if s>ma:
            ma=s
            rr=r
    return (rr,ma)
    
print(task(2000))

n=int(input())
res = dict()
quat = [i*i for i in range(1,n+1)]
m, amax = 0, 0
for i in range(1,n+1):
    a = 0
    for x in range(1,int(i / 1.4)+1):
        if i*i - x*x in quat :
            a+=1 
            #if a == 4 :
                #print(i,x,quat.index(i*i-x*x)+1)
    if a > amax :
        amax = a
        m = i
print(amax*8 + 4,m)
print("end")

pifpr = [5, 13, 17, 29, 37, 41]
n = int(input())
s = n
i, a = 0, 1
while True :
    if pifpr[i] <= n :
        a *= pifpr[i]
        n //= pifpr[i]
    elif pifpr[i - 1] <= n :
        a *= pifpr[i - 1]
        n //= pifpr[i - 1]
    else :
        break
    i += 1
# a - искомый радиус
print(a)
# Для проверки координат при "а"
s = 0
for x in range(1,int(a / 1.4)+1):
    x2 = round((a * a - x * x) ** 0.5)
    if  x2 * x2 + x * x == a * a :
        s += 1 
        print(x, x2, a, ' -> № ', s)
# умножение на 8 - эта точка по модулю 
# во всех четвертях (4) умножено на 2 -
# "х" и "у" можно менять местами
# "+ 4" - точки на осях
print(s * 8 + 4)
print("end")

pifpr = [5, 13, 17, 29, 37, 41, 53, 61, 73, 89, 97]
n = int(input())
s = n
i, a = 0, 1
while True :
    if pifpr[i] <= n :
        a *= pifpr[i]
        n //= pifpr[i]
    elif pifpr[i - 1] <= n :
        a *= pifpr[i - 1]
        n //= pifpr[i - 1]
    else :
        break
    i += 1
# a - искомый радиус
print(a)

pifpr = [5, 13, 17, 29, 37, 41, 53, 61, 73, 89, 97]
n = int(input())
s = n
i, a = 0, 1
while n - 1 :
    if pifpr[i] <= n :
        a *= pifpr[i]
        n //= pifpr[i]
        i += 1
    else :
        while i :
            if pifpr[i - 1] <= n :
                a *= pifpr[i - 1]
                n //= pifpr[i - 1]
            i -= 1
# a - искомый радиус
print(a)