A327633 -id:A327633 - cp's OEIS frontend

A348271 a(n) is the sum of noninfinitary divisors of n.

0, 0, 0, 2, 0, 0, 0, 0, 3, 0, 0, 8, 0, 0, 0, 14, 0, 9, 0, 12, 0, 0, 0, 0, 5, 0, 0, 16, 0, 0, 0, 12, 0, 0, 0, 41, 0, 0, 0, 0, 0, 0, 0, 24, 18, 0, 0, 56, 7, 15, 0, 28, 0, 0, 0, 0, 0, 0, 0, 48, 0, 0, 24, 42, 0, 0, 0, 36, 0, 0, 0, 45, 0, 0, 20, 40, 0, 0, 0, 84, 39

Offset: 1

Amiram Eldar, Oct 09 2021

nonn

			a(12) = 8 since 12 has 2 noninfinitary divisors, 2 and 6, and 2 + 6 = 8.

Amiram Eldar, Table of n, a(n) for n = 1..10000
Eric Weisstein's World of Mathematics, Infinitary Divisor.

Cf. A000203, A036537, A049417, A077609, A162643, A327633.

Similar sequences: A034448, A048146, A051377, A188999.

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]; a[n_]:= DivisorSigma[1,n] - isigma[n]; Array[a, 100]

a(n) = A000203(n) - A049417(n).
a(n) = 0 if and only if the number of divisors of n is a power of 2, (i.e., n is in A036537).
a(n) > 0 if and only if the number of divisors of n is not a power of 2, (i.e., n is in A162643).

A348274 Noninfinitary abundant numbers: numbers k such that A348271(k) > k.

Original entry on oeis.org

36, 48, 80, 144, 180, 240, 252, 288, 300, 324, 336, 396, 400, 432, 468, 528, 560, 576, 588, 612, 624, 684, 720, 768, 784, 816, 828, 880, 900, 912, 960, 1008, 1040, 1044, 1104, 1116, 1200, 1232, 1260, 1280, 1296, 1332, 1360, 1392, 1440, 1456, 1476, 1488, 1520, 1548

Offset: 1

Amiram Eldar, Oct 09 2021

nonn

The first odd term is a(3577) = 99225.

The number of terms not exceeding 10^k, for k = 1, 2, ..., are 0, 3, 31, 360, 3605, 36160, 360840, 3618980, 36144059, ... Apparently this sequence has an asymptotic density 0.0361...

			36 is a term since A348271(36) = 41 > 36.

Amiram Eldar, Table of n, a(n) for n = 1..10000

Cf. A348271, A327633.
Subsequence of A005101.
Similar sequences: A034683, A064597, A129575, A129656, A292982.

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]; Select[Range[1500], s[#] > # &]

A322609 Numbers k such that s(k) = 2*k, where s(k) is the sum of divisors of k that have a square factor (A162296).

Original entry on oeis.org

24, 54, 112, 150, 294, 726, 1014, 1734, 1984, 2166, 3174, 5046, 5766, 8214, 10086, 11094, 13254, 16854, 19900, 20886, 22326, 26934, 30246, 31974, 32512, 37446, 41334, 47526, 56454, 61206, 63654, 68694, 71286, 76614, 96774, 102966, 112614, 115926, 133206

Offset: 1

Amiram Eldar, Dec 20 2018

nonn

This sequence is infinite since 6*p^2 is included for all primes p. Terms that are not of the form 6*p^2: 112, 1984, 19900, 32512, 134201344, ...
Includes 4*k if k is an even perfect number: see A000396. - Robert Israel, Jan 06 2019
From Amiram Eldar, Oct 01 2022: (Start)
24 = 6*prime(1)^2 = 4*A000396(1) is the only term that is common to the two forms that are mentioned above.
19900 is the only term below 10^11 which is not of any of these two forms. Are there any other such terms?
All the known nonunitary perfect numbers (A064591) are also of the form 4*k, where k is an even perfect number.
Equivalently, numbers k such that A325314(k) = -k. (End)

			24 is a term since A162296(24) = 48 = 2*24.

Amiram Eldar, Table of n, a(n) for n = 1..12091 (terms below 10^11; terms 1..300 from Robert Israel)

Cf. A162296, A325314.
Subsequence of A005101 and A013929.
Numbers k such that A162296(k) = m*k: A005117 (m=0), A001248 (m=1), this sequence (m=2), A357493 (m=3), A357494 (m=4).
Similar sequences: A000396, A002827, A007357, A054979, A064591, A322486, A323757, A324707, A327633.

filter:= proc(n) convert(remove(numtheory:-issqrfree,numtheory:-divisors(n)),`+`)=2*n end proc:
select(filter, [$1..200000]); # Robert Israel, Jan 06 2019

s[1]=0; s[n_] := DivisorSigma[1,n] - Times@@(1+FactorInteger[n][[;;,1]]); Select[Range[10000], s[#] == 2# &]

s(n) = sumdiv(n, d, d*(1-moebius(d)^2)); \\ A162296
isok(n) = s(n) == 2*n; \\ Michel Marcus, Dec 20 2018

from sympy import divisors, factorint
A322609_list = [k for k in range(1,10**3) if sum(d for d in divisors(k,generator=True) if max(factorint(d).values(),default=1) >= 2) == 2*k] # Chai Wah Wu, Sep 19 2021

A348343 Smaller 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).

Original entry on oeis.org

336, 1792, 5376, 6096, 21504, 32004, 97536, 34062336, 64512000, 118008576, 30064771072

Offset: 1

Amiram Eldar, Oct 13 2021

The larger counterparts are in A348344.

			336 is a term since A348271(336) = 448 and A348271(448) = 336.

Cf. A327633, A348271, A348344.

Similar sequences: A002025, A002952, A126165, A126169, A259038, A292980, A322541, A324708.

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, n]], {n,1,10^4}]; seq

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).

Original entry on oeis.org

448, 2032, 8128, 7168, 24384, 41984, 130048, 41940480, 102222432, 221316608, 34359738352

Offset: 1

Amiram Eldar, Oct 13 2021

The terms are ordered according to their smaller counterparts (A348343).

			448 is a term since A348271(448) = 336 and A348271(336) = 448.

Cf. A327633, A348271, A348343.

Similar sequences: A002046, A002953, A126166, A126170, A259039, A292981, A322542, A324709.

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

cp's OEIS Frontend

A348271 a(n) is the sum of noninfinitary divisors of n.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A348274 Noninfinitary abundant numbers: numbers k such that A348271(k) > k.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A322609 Numbers k such that s(k) = 2*k, where s(k) is the sum of divisors of k that have a square factor (A162296).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

A348343 Smaller 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).

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

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).

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica