A331971 -id:A331971 - cp's OEIS frontend

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

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, 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, 2, 6, 1, 8, 2, 5, 4, 6, 4, 9, 0, 6, 1, 5, 3, 10, 2, 8, 2

Offset: 2

Amiram Eldar, Feb 03 2020

nonn

The infinitary version of A048138.

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.

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

Amiram Eldar, Table of n, a(n) for n = 2..30000

Cf. A048138, A049417, A077609, A126168, A324938, A331971.

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

A372739 a(n) is the number of possible values of k such that the sum of aliquot coreful divisors of k (A336563) is n.

Original entry on oeis.org

0, 1, 1, 0, 1, 3, 1, 0, 0, 2, 1, 1, 1, 3, 2, 0, 1, 1, 1, 0, 2, 2, 1, 1, 0, 2, 0, 0, 1, 6, 1, 0, 2, 2, 2, 1, 1, 2, 3, 0, 1, 5, 1, 0, 0, 2, 1, 0, 0, 0, 2, 0, 1, 0, 2, 1, 2, 2, 1, 2, 1, 3, 0, 0, 2, 4, 1, 0, 2, 4, 1, 0, 1, 2, 0, 0, 2, 5, 1, 1, 0, 2, 1, 1, 2, 2, 2

Offset: 1

Amiram Eldar, May 12 2024

nonn

A coreful divisor d of n is a divisor that is divisible by every prime that divides n (see also A307958).

			a(2) = 1 since there is 1 possible value of k, k = 4, such that A336563(k) = 2.
a(6) = 3 since there are 3 possible values of k, k = 8, 12 and 18, such that A336563(k) = 6.

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

Cf. A307958, A336563, A372740, A372742.

Similar sequences: A048138, A324938, A331971, A331973.

f[p_, e_] := (p^(e + 1) - 1)/(p - 1) - 1; s[1] = 0; s[n_] := Times @@ f @@@ FactorInteger[n] - n; seq[max_] := Module[{v = Table[0, {max}], i}, Do[i = s[k]; If[0 < i <= max, v[[i]]++], {k, 1, max^2}]; v]; seq[100]

s(n) = {my(f = factor(n)); prod(i = 1, #f~, (f[i, 1]^(f[i, 2] + 1) - 1)/(f[i, 1] - 1) - 1) - n;}
lista(nmax) = {my(v = vector(nmax), i); for(k = 1, nmax^2, i = s(k); if(i > 0 && i <= nmax, v[i]++)); v;}

a(n) = 0 if and only if n is in A372740.

a(n) = 1 if and only if n is in A372742.

cp's OEIS Frontend

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

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica