四平方和（前缀和+二分查找）

最新推荐文章于 2024-09-12 19:02:20 发布

余晴尽欢

最新推荐文章于 2024-09-12 19:02:20 发布

阅读量190

点赞数 2

文章标签：算法数据结构

本文链接：https://blog.csdn.net/lsfff666/article/details/136769898

版权

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

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

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

比如：

$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。

输出格式

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

数据范围

0<N<5e6

输入样例：

输出样例：

#include<iostream>
#include<algorithm>
using namespace std;

const int N = 5e6 + 3;
int n, m;

struct Sum
{
	int s, c, d;
	bool operator<( Sum& t)
	{
		if (s != t.s) return s < t.s;
		if (c != t.c)return c < t.c;
		return d < t.c;
	}
}sum[N];

int find(int t)
{
	int l = 0, r = m - 1;
	while (l < r)
	{
		int mid = l + r >> 1;
		if (sum[mid].s >= t) r = mid;
		else l = mid + 1;
	}
	return l;
}
int main()
{
	scanf("%d", &n);

	for (int a = 0; a * a <= n; a++)
		for (int b = a; b * b <= n; b++)
			sum[m++] = { a * a + b * b,a,b };

	sort(sum, sum + m);

	for(int a=0;a*a<=n;a++)
		for (int b = a; b * b <= n; b++)
		{
			int t = n - a * a - b * b;
			int u = find(t);
			if (sum[u].s == t)
			{
				printf("%d %d %d %d\n", a, b, sum[u].c, sum[u].d);
				return 0;
			}
		}
	return 0;
}