题目
四平方和定理,又称为拉格朗日定理:
每个正整数都可以表示为至多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)
解题思路
使用暴力循环可以解开这题,但因为数据规模不能得全分,需要进行优化。从两方面入手。
第一,引入散列表用来存入键值;
第二,将四重嵌套循环分开成两个两重循环,将第一个两重循环得结果存入散列表,再在第二个两重循环中,查看是否有键值相等的散列表。
代码
暴力循环:要注意变量的范围,加上等号
#include<iostream>
using namespace std;
int main(){
int n = 0;
scanf("%d",&n);
for(int a=0;a*a<=n/4;a++){ //a*a<=n/4,不能写<n/4
for(int b=a;b*b<=n/3;b++){
for(int c=b;c*c<=n/2;c++){
for(int d=c;d*d<=n;d++){
if(a*a+b*b+c*c+d*d==n){
printf("%d %d %d %d",a,b,c,d);
return 0;
}
}
}
}
}
return 0;
}
优化:
#include<iostream>
#include<cmath>
#include<map>
using namespace std;
int n = 0;
map<int,int> cache;
int main(){
scanf("%d",&n);
for(int c=0;c*c<=n/2;c++){
for(int d=c;d*d<=n;d++){
cache[c*c+d*d] = c;
}
}
for(int a=0;a*a<=n/4;a++){
for(int b=a;b*b<=n/3;b++){
if(cache.find(n-a*a-b*b)!=cache.end()){
int c = cache[n-a*a-b*b];
printf("%d %d %d %d",a,b,c,(int)sqrt(n-a*a-c*c-b*b)); //这里需要类型转换
return 0;
}
}
}
return 0;
}
注意点:
- 两重循环的顺序不能颠倒
- map的find方法返回元素的位置,如果没有找到则和end方法返回的位置相同。