题目描述
四平方和定理,又称为拉格朗日定理:
每个正整数都可以表示为至多 4 个正整数的平方和。
如果把 0 包括进去,就正好可以表示为 4 个数的平方和。
比如:
5 = 0^2 + 0^2 + 1^2 + 2^25=0^2+0^2+1^2+2^2;
7 = 1^2 + 1^2 + 1^2 + 2^27=1^2+1^2+1^2+2^2;
对于一个给定的正整数,可能存在多种平方和的表示法。
要求你对 4 个数排序:
0≤a≤b≤c≤d
并对所有的可能表示法按 a,b,c,d 为联合主键升序排列,最后输出第一个表示法。
输入描述
程序输入为一个正整数 N (N<5 \times 10^6)N(N<5×106)。
输出描述
要求输出 4 个非负整数,按从小到大排序,中间用空格分开
输入输出样例
示例
输入
12
输出
0 2 2 2
运行限制
- 最大运行时间:3s
- 最大运行内存: 256M
法一:加多个break及时退出,令输出的是第一个解(这个较繁琐,不好)
用flag设置一个标签,得到第一个解后即退出
from math import sqrt
n = int(input())
flag = True
for a in range(int(sqrt(n)) + 1):
for b in range(a, int(sqrt(n - a*a)) + 1):
for c in range(b, int(sqrt(n - a*a - b*b)) + 1):
d = sqrt(n - a*a - b*b - c*c)
if d == int(d) and d >= c:
res = [a, b, c, int(d)]
flag = False
break
if not flag:
break
if not flag:
break
res.sort()
for i in range(len(res)):
print(res[i], end=" ")
法二:用sys.exit(0)退出
import sys
from math import sqrt
n = int(input())
for i in range(int(sqrt(n)) + 1):
for j in range(i, int(sqrt(n - i*i)) + 1):
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)