四平方和
四平方和定理,又称为拉格朗日定理:
每个正整数都可以表示为至多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<5000000)
要求输出4个非负整数,按从小到大排序,中间用空格分开
例如,输入:
5
则程序应该输出:
0 0 1 2
再例如,输入:
12
则程序应该输出:
0 2 2 2
再例如,输入:
773535
则程序应该输出:
1 1 267 838
解题思路:
定义四个变量,用循环求出,但由于数比较大,需进一步优化。
可以从开方数最大值开始循环,如果值大于输入值跳出循环等等……还需进一步优化
部分结果:
#include<iostream>
#include<math.h>
#include<algorithm>
using namespace std;
int main()
{
int a, b, c, d;
int s[4];
long n;
cin >> n;
for (a = 0; a < sqrt(n)/2; a++)
for (b = 0; b < sqrt(n); b++)
for (c = sqrt(n); c >= 0; c--)
{
int d = n - c*c - b*b - a*a;
if (a*a + b*b + c*c + d*d == n)
{
s[0] = a;
s[1] = b;
s[2] = c;
s[3] = d;
sort(s, s + 4);
for (int i = 0; i < 4; i++)
{
if (i == 3)
cout << s[i] << endl;
else cout << s[i] << " ";
}
}
if (a*a + b*b + c*c + d*d > n)
continue;
}
}