四平方和（枚举）

1.问题描述：

四平方和
四平方和定理，又称为拉格朗日定理：
每个正整数都可以表示为至多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

2.问题分析：

既然a,b,c,d是非递减序，也就是a<=b<=c<=d，并且呢，a^2+b^2+c^2+d^2=n。

不妨设a^2>n/4，那么a^2+b^2+c^2+d^2>n，矛盾。

对于b的初始值为a，不妨设b^2>n/3, 那么b^2+c^2+d^2>n，矛盾。

同理，得到c,d的上界。

这样我们进行了稍微的程序优化。

3.程序如下：

 

#include <stdio.h>
void fun(int n)
{
	for(int a = 0; a*a <= (n/4); a++)
	{
		for(int b = a; b*b <= (n/3); b++)
		{
			for(int c = b; c*c <= (n/2); c++)
			{
				for(int d = c; d*d <= n; d++)
				{
					if((a*a+b*b+c*c+d*d)== n)
					{
						printf("%d %d %d %d\n", a, b, c, d);
						return; 
					}
				}
			}
		}
	}
}
int main()
{
	int n;
	scanf("%d", &n);
	fun(n);
	return 0;
}