Matlab找出10万以内的亲和数对

最新推荐文章于 2023-09-19 16:36:24 发布

Burannnn

最新推荐文章于 2023-09-19 16:36:24 发布

阅读量430

点赞数

分类专栏： matlab 文章标签：算法 matlab

本文链接：https://blog.csdn.net/Burannnn/article/details/129716525

版权

matlab 专栏收录该内容

1 篇文章 0 订阅

订阅专栏

% 找出10万以内的所有亲和数对

N = 100000;

sum_divisors = zeros(1,N);
for ii = 1:N-1
    divisors = 1;
    for jj = 2:sqrt(ii)
        if mod(ii,jj) == 0
            divisors = [divisors jj ii/jj];
        end
    end
    sum_divisors(ii) = sum(divisors);
end

amicable_nums = [];
for ii = 1:N-1
    if sum_divisors(ii) < N && sum_divisors(sum_divisors(ii)) == ii && sum_divisors(ii) ~= ii
        amicable_nums = [amicable_nums; ii sum_divisors(ii)];
    end
end

disp(amicable_nums)

结果