(Problem 21)Amicable numbers

Let d(n) be defined as the sum of proper divisors of n (numbers less than n which divide evenly into n).
If d(a) = b and d(b) = a, where a ≠ b, then a and b are an amicable pair and each of a and b are called amicable numbers.

For example, the proper divisors of 220 are 1, 2, 4, 5, 10, 11, 20, 22, 44, 55 and 110; therefore d(220) = 284. The proper divisors of 284 are 1, 2, 4, 71 and 142; so d(284) = 220.

Evaluate the sum of all the amicable numbers under 10000.

#include<stdio.h>
#include<math.h>
#include<string.h>
#include<ctype.h>
#include<stdlib.h>
#include<stdbool.h>

int FactorSum(int n)  //计算n的所有小于n的因素和
{
	int i;
	int sum=1;
	for(i=2; i<=n/2; i++)
	{
		if(n%i==0)
			sum+=i;
	}
	return sum;
}

int main()
{
	int t,i=2;
	int sum=0;
	while(i<10000)
	{
		t=FactorSum(i);
		if(t!=i && FactorSum(t)==i) 
			sum+=i;
		i++;
	}
	printf("%d\n",sum);
	return 0;
}

Answer:
31626

转载于:https://www.cnblogs.com/cpoint/p/3367356.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
这段代码存在一个问题,就是没有考虑到亲密数对的概念。实际上,亲密数对是指两个数,其一个数的所有因子之和等于另一个数,且另一个数的所有因子之和等于第一个数。因此,在判断一个数是否为亲密数对时,需要先计算该数的所有因子之和,再判断是否与另一个数相等。另外,代码输出的是完全平方数和其紧随其后的质数,这与亲密数对的定义不符。 以下是修改后的代码,可以找到指定区间内的亲密数对: ``` import math # 判断一个数是否为素数 def is_prime(n): if n <= 1: return False for i in range(2, int(math.sqrt(n))+1): if n % i == 0: return False return True # 判断一个数是否为完全平方数 def is_square(n): sqrt_n = int(math.sqrt(n)) return sqrt_n**2 == n # 计算一个数的所有因子之和 def divisor_sum(n): result = 1 for i in range(2, int(math.sqrt(n))+1): if n % i == 0: result += i if n // i != i: result += n // i return result # 找出指定区间内的亲密数对 def find_amicable_numbers(s, t): result = [] for i in range(s, t+1): j = divisor_sum(i) if i != j and j <= t and divisor_sum(j) == i: result.append((i, j)) return result # 测试 s, t = map(int, input().split()) pairs = find_amicable_numbers(s, t) for pair in pairs: print(pair[0], pair[1]) ``` 这里新增了一个 `divisor_sum` 函数来计算一个数的所有因子之和,同时修改了 `find_amicable_numbers` 函数以找到所有亲密数对。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值