# Представление натурального числа
# в виде суммы квадратов
# натуральных чисел.
k=int(input())
for i in range(1, k+1):
if k==i**2:
print(i)
for i in range(1, k+1):
for j in range(i, k+1):
if k==i**2+j**2:
print(i, j)
for i in range(1, k+1):
for j in range(i, k+1):
for l in range(j, k+1):
if k==i**2+j**2+l**2:
print(i, j, l)
for i in range(1, k+1):
for j in range(i, k+1):
for l in range(j, k+1):
for m in range(l, k+1):
if k==i**2+j**2+l**2+m**2:
print(i, j, l, m)
Разумеется, такой нерациональный код выполняется крайне медленно, например, на моём домашнем компьютере при вводе числа 2021 программа работала около минуты, тогда как Wolfram Alpha делает то же самое за секунду.
Пожалуйста, подскажите, как изменить код, дабы радикально ускорить данное вычисление.
N = 1^2 + 1^2 + ... + 1^2
(N раз) ?