AcWing 1221. 四平方和

题目

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

每个正整数都可以表示为至多 4

个正整数的平方和。

如果把 0
包括进去，就正好可以表示为 4

个数的平方和。

比如：

5=02+02+12+22

7=12+12+12+22

对于一个给定的正整数，可能存在多种平方和的表示法。

要求你对 4

个数排序：

0≤a≤b≤c≤d

并对所有的可能表示法按 a,b,c,d
为联合主键升序排列，最后输出第一个表示法。

思路

二分做法

#include 
#include 
#include 
using namespace std;
const int N = 2500010;
struct number{
    int s,c,d;
    bool operator < (const number& t){
        if (s != t.s) return s < t.s;
        if (c != t.c) return c < t.c;
        return d < t.d;
    }
}sum[N];
int main(){
    int n;
    cin >> n;
    int idx=0;
    for(int c=0;c*c<=n;c++){//先把所有的c，d枚举一遍
        for(int d=c;d*d+c*c<=n;d++){
            sum[idx++]={c*c+d*d,c,d};
        }
    }
    sort(sum,sum+idx);
            
    for (int a = 0; a * a <= n; a ++ )
        for (int b = 0; a * a + b * b <= n; b ++ )
        {
            int t = n - a * a - b * b;
            int l = 0, r = idx - 1;
            while (l < r)//用二分取查找最左边的c，d，就是最小解
            {
                int mid = l + r >> 1;
                if (sum[mid].s >= t) r = mid;
                else l = mid + 1;
            }
            if (sum[l].s == t)
            {
                printf("%d %d %d %d\n", a, b, sum[l].c, sum[l].d);
                return 0;
            }
        }
}

哈希做法

#include 
#include 
#include 
#include 

#define x first
#define y second
using namespace std;
const int N = 2500010;
typedef pair<int, int> PII;
int n, m;
unordered_map<int, PII> S;
int main(){
    cin >> n;
    int idx=0;
    for(int c=0;c*c<=n;c++){
        for(int d=c;d*d+c*c<=n;d++){
            if (S.count(d*d+c*c) == 0) S[d*d+c*c] = {c, d};//保证是最小解
        }
    }
    for (int a = 0; a * a <= n; a ++ )
        for (int b = 0; a * a + b * b <= n; b ++ )
        {
            int t = n - a * a - b * b;
            int l = 0, r = idx - 1;
            if(S.count(t)){
                printf("%d %d %d %d\n", a, b,S[t].x, S[t].y);
                return 0;
                
            }
        }
}