四平方和(前缀和+二分查找)

四平方和定理,又称为拉格朗日定理:

每个正整数都可以表示为至多 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

输入样例:
5
输出样例:
#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;
}
0 0 1 2
### 实现平方和定理的Python代码 为了验证平方和定理,即任意自然数可以表示成最多个整数的平方之和,下面提供了一种通过穷举法来寻找符合条件的个整数的方法[^1]。 ```python #!/usr/bin/python3 # -*- coding: utf-8 -*- # @desc: 验证平方和定理 def find_four_squares(number): """ 查找给定自然数由个整数平方组成的组合。 参数: number (int): 用户输入的一个自然数 返回: tuple or None: 如果找到合适的个整数组合则返回该元组;否则返回None """ for x1 in range(int(number**0.5), -1, -1): for x2 in range(int((number-x1*x1)**0.5), -1, -1): for x3 in range(int((number-x1*x1-x2*x2)**0.5), -1, -1): x4 = (number-x1*x1-x2*x2-x3*x3)**0.5 if x4.is_integer(): return (x1, x2, x3, int(x4)) return None if __name__ == "__main__": try: num_input = int(input("请输入一个正整数:")) result = find_four_squares(num_input) if result is not None: print(f"{num_input}={result[0]}²+{result[1]}²+{result[2]}²+{result[3]}²") else: print("未找到满足条件的结果") except ValueError as e: print("输入错误,请确认您输入的是有效的正整数。") ``` 这段代码定义了一个`find_four_squares()`函数用于查找能够组成目标自然数n的最大可能值作为第一个变量x1,并依次减少直到最小可能性(-1),接着对于剩余部分继续尝试其他三个变量x2,x3,x4的可能性直至成功匹配或遍历结束。当找到了一组解时立即停止搜索并输出结果[^3]。
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

余晴尽欢

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值