蓝桥 四平方和

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

法一：枚举

#include <cstdio>
#include <cmath>
using namespace std;
const int mx=2250;
int a[mx];
int n,s,x,y;
bool f;
int main()
{
    for (int i=0;i<mx;i++) a[i]=i*i;
    while (~scanf ("%d",&n)) {
        f=false;
        s=0;
        for (int i=0;i<mx;i++) {
            s=a[i];
            for (int j=i;j<mx;j++) {
                s+=a[j];
                if (s>n) {
                    s-=a[j];
                    continue;
                }
                for (int k=j;k<mx;k++) {
                    s+=a[k];
                    if (s>n) {
                        s-=a[k];
                        continue;
                    }
                    x=n-s;
                    y=sqrt(x);
                    if (a[y]==x) {
                        printf ("%d %d %d %d\n",i,j,k,y);
                        f=true;
                        break;
                    }
                    else s-=a[k];
                }
                s-=a[j];
                if (f) break;
            }
            s-=a[i];
            if (f) break;
        }
    }
    return 0;
}

法二：学长写的dfs

#include <bits/stdc++.h>
using namespace std;
int n, limit;
vector<int> vec;
bool dfs(int lst, int x, int cnt){
	if(cnt == 3){
		int tmp = (int)sqrt(x);
		if(x == tmp*tmp){
			for(int i = 0; i < 3; i++)
				printf("%d ", vec[i]);
			printf("%d\n", tmp);
			return true;
		}
		return false;
	}
	for(int i = lst; i <= limit; i++){
		if(x - i * i < 0) break;
		vec.push_back(i);
		if(dfs(i, x - i * i, cnt + 1)) return true;
		vec.pop_back();
	}
	return false;
}
int main(){
	cin >> n;
	limit = (int)sqrt(n*1.0);
	dfs(0, n, 0);
	return 0;
}