【二分】四平方和

链接：四平方和

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

每个正整数都可以表示为至多 44 个正整数的平方和。

如果把 00 包括进去，就正好可以表示为 44 个数的平方和。

比如：

 

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

要求你对 4个数排序：

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

输入格式

输入一个正整数 N。

输出格式

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

数据范围

0<N<5∗10^6

输入样例：

5

输出样例：

0 0 1 2

 

 思路：

两两枚举，比如我先枚举c，d（也可以先枚举a,b）

for(c = 0; c*c <= N; c++)
	for(d = c; d*d + c * c <= N; d++)
        把c,d存起来

for(a = 0; a * a  <= N; a++)
	for(b = a; b * b + a * a <= N; b++){
		t = N - b * b - a * a;
        如果t在c*c+d*d中出现过，找到最小的输出

			

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<string>

using namespace std;

const int n = 5*1e6+5;

int a, b, c, d;

struct Sum{
	int c, d, s;
	
}sum[n];

bool compare(Sum x1, Sum x2){
	if(x1.s != x2.s)	return x1.s < x2.s;
	if(x1.c != x2.c)	return x1.c < x2.c;
	return x1.d < x2.d;
}


int main(){
	int m = 0;
	int N;
	cin >> N;
	
	for(c = 0; c*c <= N; c++){
		for(d = c; d*d + c * c <= N; d++){
			sum[m++] = {c,  d, c * c + d* d};
		}
	}
	
	sort(sum, sum+ m, compare);
	
	for(a = 0; a * a  <= N; a++){
		for(b = a; b * b + a * a <= N; b++){
			int t = N - b * b - a * a;
			int l = 0, r = m;
			while(l < r){
				int mid = (l + r) / 2;
				if(t <= sum[mid].s)	r = mid;
				else l = mid+ 1;
			}
			
			if(t == sum[l].s){
				printf("%d %d %d %d\n",a,b,sum[l].c,sum[l].d);
				return 0;
			}	
			
		}
	}
	
	
	
	
	return 0;
}