A327566 -id:A327566 - cp's OEIS frontend

A326488 Numbers m such that A327566(m) = Sum_{k=1..m} isigma(k) is divisible by m, where isigma(k) is the sum of infinitary divisors of k (A049417).

1, 2, 160, 285, 2340, 2614, 8903, 81231, 171710, 182712, 434887, 2651907, 56517068, 143714354, 922484770, 5162883263, 39421525873

Offset: 1

Amiram Eldar, Sep 20 2019

The infinitary version of A056550.

The corresponding quotients, A327566(a(n))/a(n), are 1, 2, 118, 209, 1711, 1910, 6506, 59357, 125473, 133513, 317781, 1937798, 41298052, 105014703, 674076450, 3772612983, 28806028088, ...

			2 is in the sequence since isigma(1) + isigma(2) = 1 + 3 = 4 is divisible by 2.

Cf. A049417 (isigma), A327566 (sums of isigma).

Cf. A056550 (corresponding with sigma), A064611 (unitary), A307043 (exponential), A307161 (bi-unitary).

f[p_, e_] := p^(2^(-1 + Position[Reverse @ IntegerDigits[e, 2], ?(# == 1 &)])); isigma[1] = 1; isigma[n] := Times @@ (Flatten @ (f @@@ FactorInteger[n]) + 1); seq = {}; s = 0; Do[s = s + isigma [n]; If[Divisible[s, n], AppendTo[seq, n]], {n, 1, 10^6}]; seq

A327574 Decimal expansion of the constant that appears in the asymptotic formula for average order of the infinitary divisors sum function (A049417).

Original entry on oeis.org

7, 3, 0, 7, 1, 8, 2, 4, 2, 1, 2, 7, 3, 8, 4, 2, 2, 5, 8, 3, 8, 9, 7, 5, 4, 6, 8, 1, 7, 3, 5, 3, 0, 1, 6, 1, 9, 5, 7, 2, 5, 6, 4, 3, 3, 8, 6, 1, 7, 2, 7, 8, 6, 9, 7, 0, 7, 3, 3, 6, 7, 6, 2, 3, 0, 1, 0, 7, 9, 8, 8, 3, 3, 2, 8, 0, 0, 5, 3, 4, 6, 3, 7, 0, 2, 9, 9

Offset: 0

Amiram Eldar, Sep 17 2019

The asymptotic mean of the infinitary abundancy index lim_{n->oo} (1/n) * Sum_{k=1..n} A049417(k)/k = 1.461436... is twice this constant. - Amiram Eldar, Jun 13 2020

			0.730718242127384225838975468173530161957256433861727...

Steven R. Finch, Mathematical Constants II, Cambridge University Press, 2018, section 1.7.5, pp. 53-54.

Graeme L. Cohen and Peter Hagis, Jr., Arithmetic functions associated with infinitary divisors of an integer, International Journal of Mathematics and Mathematical Sciences, Vol. 16, No. 2 (1993), pp. 373-383.

Cf. A049417, A050376, A327566.

Cf. A013661 (corresponding constant for all divisors), A275480 (exponential), A306633 (unitary), A307160 (bi-unitary).

$MaxExtraPrecision = 1000; m = 1000; em = 10; f[x_] := Sum[Log[1 + x^(2^e)/(1 + 1/x^(2^e))], {e, 0, em}]; c = Rest[CoefficientList[Series[f[x], {x, 0, m}], x]*Range[0, m]]; RealDigits[(1/2) * Exp[NSum[Indexed[c, k]*PrimeZetaP[k]/k, {k, 2, m}, NSumTerms -> m, WorkingPrecision -> m]], 10, 100][[1]]

Equals Limit_{k->oo} A327566(k)/k^2.

Equals (1/2) * Product_{P} (1 + 1/(P*(P+1))), where P are numbers of the form p^(2^k) where p is prime and k >= 0 (A050376).

cp's OEIS Frontend

A326488 Numbers m such that A327566(m) = Sum_{k=1..m} isigma(k) is divisible by m, where isigma(k) is the sum of infinitary divisors of k (A049417).

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica