四平方和问题

四平方和问题

四平方和定理，又称为拉格朗日定理：

每个正整数都可以表示为至多 4 个正整数的平方和。

如果把 0 包括进去，就正好可以表示为 4 个数的平方和。

比如：

5 = 0^2 +0^2 +1^2 +2^2;
7= 1^2 + 1^2 + 1^2 + 2^2;

对于一个给定的正整数，可能存在多种平方和的表示法。

要求你对 4 个数排序：

0≤a≤b≤c≤d

并对所有的可能表示法按 a,b,c,d 为联合主键升序排列，最后输出第一个表示法。

输入描述

程序输入为一个正整数 N (N < 5*10^6)。

输出描述

要求输出 4 个非负整数，按从小到大排序，中间用空格分开

输入输出样例

示例
输入

12

输出

0 2 2 2

运行限制
最大运行时间：3s
最大运行内存: 256M

解答思路：
这道题第一反应是可以暴力破解，即使用四层循环直接求解。

例如

import time
n = int(input()) 
for i in range(n+1): # 构建四重循环
    for j in range(n+1):
        for k in range(n+1):
          for l in range(n+1):
              result = [i,j,k,l]
result.sort() # sort()函数默认使列表里面的元素呈升序排列，但可以添加reverse=True使得它呈降序排列
print(result) #打印最终排列好的结果

for m in result:
    print(m,end=' ')

end = time.time()
print('\n运行时间为：', end='')
print(end - start)

但是这样做，我们可以发现，根据题目要求，输入数值最大为10的6次方，每重循环都需要循环10的6次方，通过四层循环，需要总共循环10的24次方，这样做，运行完整个程序，是不可能在题目规定的时间内通过更多的样例的，简单测试了一下：
样例只有50的时候都已经需要至少2秒，而题目样例范围是10的6次方，题目要求是运行时间不超过三秒，可能只能通过一个样例，甚至不能通过

因此，我们需要进行优化

对于求平方和问题，其实可以转化为根号问题，N = a^2, a的定义域是(0,根号a),所以不需要实际range到n,只需要range到(根号n) +1，如果能被开根号，就说明是符合条件的，依此类推，四平方和问题也是一样的思路

最终解答如下：

#代码
import sys
from math import sqrt
n = int(input()) 
for i in range(int(sqrt(n))+1): # range到根号n+1,这样才能包括根号n
    for j in range(i,int(sqrt(n-i*i))+1): #range直接从i开始，这样直接从内部就排序好了，数量级也变小了
        for k in range(j,int(sqrt(n-i*i-j*j))+1):
          l = sqrt(n-i*i-j*j-k*k)
          if l<k: 
              break
          if int(l) == l:
              print(i,j,k,int(l))
              sys.exit(0)