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.

A331973 a(n) is the number of values of m such that the sum of proper infinitary divisors of m (A126168) is n.

This page as a plain text file.
%I A331973 #10 Feb 04 2020 03:56:33
%S A331973 0,0,0,0,1,1,2,1,2,1,2,2,2,2,2,3,1,3,2,4,3,3,2,3,2,3,3,4,2,4,1,4,3,4,
%T A331973 3,5,0,3,2,4,3,5,1,4,3,4,2,6,2,5,2,5,3,7,1,6,2,4,2,7,1,5,4,5,3,8,0,5,
%U A331973 2,6,1,8,2,5,4,6,4,9,0,6,1,5,3,10,2,8,2
%N A331973 a(n) is the number of values of m such that the sum of proper infinitary divisors of m (A126168) is n.
%C A331973 The infinitary version of A048138.
%C A331973 The offset is 2 as in A048138 since there are infinitely many numbers k (the primes and squares of primes) for which A126168(k) = 1.
%H A331973 Amiram Eldar, <a href="/A331973/b331973.txt">Table of n, a(n) for n = 2..30000</a>
%e A331973 a(8) = 2 since 8 is the sum of the proper infinitary divisors of 2 numbers: 10 (1 + 2 + 5) and 12 (1 + 3 + 4).
%t A331973 fun[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 @@ (fun @@@ FactorInteger[n]); is[n_] := isigma[n] - n; m = 300; v = Table[0, {m}]; Do[i = is[k]; If[2 <= i <= m, v[[i]]++], {k, 1, m^2}]; Rest@v
%Y A331973 Cf. A048138, A049417, A077609, A126168, A324938, A331971.
%K A331973 nonn
%O A331973 2,7
%A A331973 _Amiram Eldar_, Feb 03 2020