#acwing四平方和

acwing四平方和

题目
四平方和定理，又称为拉格朗日定理：
每个正整数都可以表示为至多 4 个正整数的平方和。
如果把 0 包括进去，就正好可以表示为 4 个数的平方和。
比如：
5=02+02+12+22
7=12+12+12+22
对于一个给定的正整数，可能存在多种平方和的表示法。
要求你对 4 个数排序：
0≤a≤b≤c≤d
并对所有的可能表示法按 a,b,c,d 为联合主键升序排列，最后输出第一个表示法。

输入格式
输入一个正整数 N。

输出格式
输出4个非负整数，按从小到大排序，中间用空格分开。

数据范围
0<N<5∗106

输入样例：

5

输出样例：

0 0 1 2

分析
数据的取值范围是：0<N<5∗10⁶，每个数都<= √ N，因此，每个数的取值范围大约是0~2236，两层循环大约是10⁶，三层循环大约是10⁹,会超时，所以最多只能枚举两层（最坏情况下）
1.枚举三个数a,b,c。d=√ (N-(a*a+b*b+c*c)
2.先枚举两个数，并把这两个数的和（a²+b²）,a,b用结构体数组存起来存起来，排序，再枚举另两个数（c^2^+d^2^）,求一个差值t=n-c*c-d*d，再利用二分法在存储的结构体中的和中找到我们一个值等与t，找了话，就从c,d,m[l].a,m[l]b;依次输出，那么输出的就是顺序就是题目所谓的字典序。

暴力法

#include<cstdio>
#include<iostream>
#include<cmath>
#include<algorithm>

using namespace std;

int main()
{
    int N;
    scanf("%d",&N);
    //有a*a+b*b+c*c<=N做限制条件，3层循环实际的时间复杂度是小于10^6,因此不会超出时间限制
    for(int a=0;a*a<=N;a++)
    for(int b=a;a*a+b*b<=N;b++)//0≤a≤b
    for(int c=b;a*a+b*b+c*c<=N;c++)//0≤a≤b≤c
    {
        int t=N-a*a-b*b-c*c;
        int d=sqrt(t);
        if(d*d==t)
        {
            printf("%d %d %d %d\n",a,b,c,d);
            return 0;
        }
    }
}

二分法

#include<cstdio>
#include<iostream>
#include<cmath>
#include<algorithm>

using namespace std;
const int N=5000010;
int n,mm;
struct sum//定义一个数组，记录枚举两个数的值和各自分别的值
{
    int s,a,b;
    bool operator< (const sum &t)const//重载小于符号，为了后面排序
    {
        if(s!=t.s)return s<t.s;
        if(a!=t.a)return a<t.a;
        return b<t.b;
    }
}m[N];

int main()
{
    int n,a,b,c,d;
    scanf("%d",&n);
    for(a=0;a*a<=n;a++)
    for(b=a;a*a+b*b<=n;b++)
    {
        m[mm++]={a*a+b*b,a,b};
    }
    
    sort(m,m+mm);//结构数组m排序
    
    for(c=0;c*c<=n;c++)
    for(d=c;d*d<=n;d++)//从小到大枚举c和d
    {
        int t=n-c*c-d*d;//找出c与d的差值
        int l=0,r=mm-1,mid;
        while(l<r)//在前面m数组中用二分法找出到差值
        {
         mid=l+r>>1;
         if(t<=m[mid].s)r=mid;
         else l=mid+1;
        }
        if(m[r].s==t)//找到了就输出结果
        {
            cout<<c<<' '<<d<<' '<<m[l].a<<' '<<m[l].b<<endl;//c与d是从小到大枚举，因此，c与d是最小的。
            return 0;
        }
    }
    return 0;
}