2016蓝桥杯C++： 四平方和

四平方和定理，又称为拉格朗日定理：

每个正整数都可以表示为至多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

思路

最开始直接暴力枚举，明显超时

#include <bits/stdc++.h>
using namespace std;
int main()
{
    double N;
    scanf("%lf",&N);
    double n=ceil (sqrt(N));
    for (int i = 0; i < n/4; ++i) {
        for (int j = i; j < n/3; ++j) {
            for (int k = j; k <n/2 ; ++k) {
                for (int l = k; l <n ; ++l) {
                    if(i*i+j*j+k*k+l*l==N)
                    {
                        printf("%d %d %d %d",i,j,k,l);
                        exit(0);
                    }
                }
            }
        }
    }
}

后来想到最后一层循环大可不必

#include <bits/stdc++.h>
using namespace std;
int main()
{
    double N;
    scanf("%lf",&N);
    double n=ceil (sqrt(N));
    for (int i = 0; i < n/4; ++i) {
        for (int j = i; j < n/3; ++j) {
            for (int k = j; k <n/2 ; ++k) {
               double l=sqrt(N-i*i-j*j-k*k);
                    if(l==(int)l)
                    if(i*i+j*j+k*k+l*l==N)
                    {
                        printf("%d %d %d %d",i,j,k,int(l));
                        exit(0);
                    }

            }
        }
    }
}

但还是超时，后来借鉴大佬的方法，用map哈希表缓存kk+ll,再计算N-ii-jj是否等于kk+ll,即可得到k,l

#include <bits/stdc++.h>
using namespace std;
map <int,int> fun;
int main()
{
    double N;
    scanf("%lf",&N);
    double n=ceil (sqrt(N));
    for (int k = 0; k <n/2 ; ++k) {
        for (int l = k; l <n ; ++l) {
            if(fun.find(k*k+l*l)==fun.end())
                fun[k*k+l*l]=k;
        }
    }
    for (int i = 0; i < n/4; ++i) {
        for (int j = i; j < n/3; ++j) {
            if(fun.find(N-i*i-j*j)!=fun.end())
            {
                int k=fun[N-i*i-j*j];
                int l=int(sqrt(N-i*i-j*j-k*k));
                printf("%d %d %d %d",i,j,k,l);
                exit(0);
            }
            }
        }
}