四平方和 c++

四平方和定理，又称为拉格朗日定理：每个正整数都可以表示为至多4个正整数的平方和。
如果把0包括进去，就正好可以表示为4个数的平方和。
比如：
5 = 0^2 + 0^2 + 1^2 + 2^2
7 = 1^2 + 1^2 + 1^2 + 22（符号表示乘方的意思）
对于一个给定的正整数N，可能存在多种平方和的表示法。
要求你对4个数排序：0 <= a <= b <= c <= d
并对所有的可能表示法按 a,b,c,d 为联合主键升序排列，最后输出第一个表示法
输入
输入存在多组测试数据，每组测试数据输入一行为一个正整数N (N<5000000)
输出
对于每组测试数据，要求输出4个非负整数，按从小到大排序，中间用空格分开
样例输入
5
样例输出
0 0 1 2
样例输入
12
样例输出
0 2 2 2

数据范围为5000000；若暴力枚举，即使开平方数值也过大，复杂度n4无法通过（1s大概运行10的8次方），需要优化。所有的数范围要小于等于sqrt（n）。代码如下

#include<iostream>
using namespace std;  

int main()
{
	int n; cin >> n;
	int q = int (sqrt(n));//获得开根的整数
	for (int i = 0; i <= q; i++)
	{
		for (int j = 0; j <= q; j++)
		{   
			if (i * i + j * j > n) break;//优化，超过的直接下一个；
			for (int k = 0; k <= q;k++)
			{
				if (i * i + j * j+k*k >n) break;//优化，超过的直接下一个；
				double q = sqrt(n - i * i - j * j - k * k);
				if (q ==int(q))//判断是否是整数
				{
					cout << i << " " << j << " " << k << " " << q;
					return 0;//跳出三重for循环
				}
			}
		}
	}
	return 0;
}