2016-8-四平方和

博客围绕四平方和定理展开,该定理指出每个正整数可表示为至多4个正整数平方和,包含0时正好为4个数平方和。对于给定正整数,需按特定规则输出第一个表示法。采用暴力枚举并优化,如缩小变量取值范围、减少变量个数,避免超时。


四平方和

四平方和定理,又称为拉格朗日定理: 每个正整数都可以表示为至多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                                        12                                          773535
则程序应该输出:
0 0 1 2                               0 2 2 2                                   1 1 267 838

思路:暴力枚举+优化   

优化① 缩小变量的取值范围    double tmp = sqrt(n);

       ② 减少变量个数   计算出d的值    int d  =  (int) sqrt( n-a*a-b*b-c*c-d*d );

重点在优化,否则会超时

 

代码

#include <cstdio>
#include <cmath>
#include<iostream>
#define MAXN 2300     //(int)sqrt(5000000)=2236
using namespace std;
//注意范围,题目给的范围枚举到2300就够了
/*
直接枚举从小的开始,枚举到答案就可以return了,
不过注意的是只枚举三个数,最后一个数可以倒着求出来,然后正过来验证,
*/

int main(){
    int n;
    scanf("%d",&n);
    double tmp=sqrt(n);  //tmp<MAXN  这是关键,不要让时间复杂度太大(超时问题)
    
    for(int a=0; a<tmp; a++){
        
        for(int b=a; b<tmp; b++){
        
            for(int c=b; c<tmp; c++){
            	
                int d = (int)sqrt(n-a*a-b*b-c*c);     // sqrt()返回值是double,强制转换为整形  算出第4个数
                if( a*a+b*b+c*c+d*d==n ){
                    printf("%d %d %d %d\n",a,b,c,d);
                    return 0;              
                }
            }
        }
    }
    return 0;
}

 

### 实现平方和定理的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]。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值