cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A348344 Larger member of a noninfinitary amicable pair: numbers (k, m) such that nisigma(k) = m and nisigma(m) = k, where nisigma(k) is the sum of the noninfinitary divisors of k (A348271).

This page as a plain text file.
%I A348344 #4 Oct 14 2021 21:51:20
%S A348344 448,2032,8128,7168,24384,41984,130048,41940480,102222432,221316608,
%T A348344 34359738352
%N A348344 Larger member of a noninfinitary amicable pair: numbers (k, m) such that nisigma(k) = m and nisigma(m) = k, where nisigma(k) is the sum of the noninfinitary divisors of k (A348271).
%C A348344 The terms are ordered according to their smaller counterparts (A348343).
%e A348344 448 is a term since A348271(448) = 336 and A348271(336) = 448.
%t A348344 f[p_, e_] := Module[{b = IntegerDigits[e, 2], m}, m = Length[b]; Product[If[b[[j]] > 0, 1 + p^(2^(m - j)), 1], {j, 1, m}]]; isigma[1] = 1; isigma[n_] := Times @@ f @@@ FactorInteger[n]; s[n_] := DivisorSigma[1,n] - isigma[n]; seq={}; Do[m=s[n]; If[m>n && s[m]==n, AppendTo[seq, m]], {n,1,10^4}]; seq
%Y A348344 Cf. A327633, A348271, A348343.
%Y A348344 Similar sequences: A002046, A002953, A126166, A126170, A259039, A292981, A322542, A324709.
%K A348344 nonn,more
%O A348344 1,1
%A A348344 _Amiram Eldar_, Oct 13 2021