四平方和定理,又称为拉格朗日定理:
每个正整数都可以表示为至多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
资源约定:
峰值内存消耗 < 256M
CPU消耗 < 3000ms
请严格按要求输出,不要画蛇添足地打印类似:“请您输入...” 的多余内容。
所有代码放在同一个源文件中,调试通过后,拷贝提交该源码。
注意: main函数需要返回0
注意: 只使用ANSI C/ANSI C++ 标准,不要调用依赖于编译环境或操作系统的特殊函数。
注意: 所有依赖的函数必须明确地在源文件中 #include <xxx>, 不能通过工程设置而省略常用头文件。
提交时,注意选择所期望的编译器类型。
分析:先利用map建立一个索引将c与c*c+d*d建立映射,之后枚举a与b,在索引中寻找N-a*a-b*b是否在索引中,如果找到,即找到a,b,c
#include <iostream>
#include <map>
#include <cmath>
using namespace std;
map <int,int>x;
int main()
{
int n;
cin>>n;
for(int c=0;c*c<=n/2;c++)
for(int d=c;c*c+d*d<=n;d++)
{
if(x.find(c*c+d*d)==x.end())
x[c*c+d*d]=c;
}
for(int a=0;a*a<=n/4;a++)
for(int b=a;a*a+b*b<=n/2;b++)
{
if(x.find(n-a*a-b*b)!=x.end())
{
int c=x[n-a*a-b*b];
int d=(int)sqrt(n-a*a-b*b-c*c);
printf("%d %d %d %d\n",a,b,c,d);
return 0;
}
}
return 0;
}