def Divisors(N): N=abs(N) L,U=[],[] k=1 while k*k <=N: if N%k== 0: L.append(k) if k*k!=N: U.append(N//k) k+=1 return L+U[::-1] def Floor_Root(a,k): """floor(a^(1/k)) を求める. a:非負整数 k:正の整数 """ assert 0<=a and 0a: x-=1 return x #================================================== def f(c,x): if S%(c*c-x*x)==0: M=S//(c*c-x*x)+c*c if M>=0: L=Floor_Root(M,2) if L*L==M: return L else: return -1 else: return -1 else: return -1 def sq_area(a,b,c): s=(a+b+c)//2 return 16*s*(s-a)*(s-b)*(s-c) def is_valid(a,b,c): return min(a,b,c)>=1 and sq_area(a,b,c)==S and abs(b-c)=0: a=( x+y)//2 b=(-x+y)//2 if is_valid(a,b,c): a,b,c=sorted([a,b,c]) E.add((a,b,c)) print(len(E))