蓝桥杯省赛真题-四平方和

四平方和

四平方和定理，又称为拉格朗日定理：
每个正整数都可以表示为至多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>， 不能通过工程设置而省略常用头文件。

提交时，注意选择所期望的编译器类型。

1.思路

看到这个题目的时候，就很容易想到暴力解法：如四重for循环求和，但时间复杂度太高了，只能舍弃这个。然后还有三重for循环
下面来看看吧

三重for循环

#include<cstdio>
#include<iostream>
#include<cmath>
const int maxn=1e5+5;
typedef long long ll;
using namespace std;
int main(){
	ll n;
	cin>>n;
	for(int i=0;i<1005;i++){
		for(int j=0;j<2005;j++){
			for(int q=0;q<2005;q++){
				int tmp=sqrt(n-i*i-j*j-q*q);
				if(n==i*i+j*j+q*q+tmp*tmp){
					cout<<i<<" "<<j<<" "<<q<<" "<<tmp<<endl;
					return 0;
				}
			}
		}
	}
}

肯定时间复杂度更低的方法，怎么实现呢？
aa+bb+cc+dd=n;
用一个数组存c*c+d*d 就用到了双重for循环
然后再用一个双重for循环来计算a和b

双重for循环map

#include<iostream>
#include<cstdio>
#include<map>
#include<cmath>
using namespace std;

int n;
map<int,int>cache;
int main(){
	scanf("%d",&n);
	for(int c=0;c*c<=n/2;c++){
		for(int d=c;c*c+d*d<=n;d++){
			if(cache.find(c*c+d*d)==cache.end())     //找不到就存进去 
			cache[c*c+d*d]=c;
		}
	}
	for(int a=0;a*a<=n/4;a++){
		for(int b=a;a<=n/2;b++){
			if(cache.find(n-a*a-b*b)!=cache.end()){
			int c=cache[n-a*a-b*b];
			int d=int(sqrt(n-a*a-b*b-c*c));
			printf("%d %d %d %d",a,b,c,d);
			return 0; 
		}
	}
}
}