from math import sqrt def fermat(n): N=round(sqrt(n))+1 for i in range(1,N): if i**2==n: return str(i)+"^2" for i in range(1,N): for j in range(1,N): if i**2+j**2==n: return str(i)+"^2 + "+str(j)+"^2" for i in range(1,N): for j in range(1,N): for k in range(1,N): if i**2+j**2+k**2==n: return str(i)+"^2 + "+str(j)+"^2 + "+str(k)+"^2" else: for l in range(1,N): if i**2+j**2+k**2+l**2==n: return str(i)+"^2 + "+str(j)+"^2 + "+str(k)+"^2 + "+str(l)+"^2"