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.

Showing 1-4 of 4 results.

A327837 Decimal expansion of the asymptotic mean of the number of exponential divisors function (A049419).

Original entry on oeis.org

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

Views

Author

Amiram Eldar, Sep 27 2019

Keywords

Examples

			1.602317102305418052349626315621161003776939495785572...
		

Crossrefs

Cf. A059956 (constant for unitary divisors), A306071 (bi-unitary), A327576 (infinitary).

Programs

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

Formula

Equals lim_{k->oo} A145353(k)/k.
Equals Product_{p prime} (1 + Sum_{e >= 2} p^(-e) * (d(e) - d(e-1))), where d(e) is the number of divisors of e (A000005).
Equals Product_{p prime} (1 - 1/p) * (2 - (log(p-1) + QPolyGamma(0, 1, 1/p)) / log(p)). - Vaclav Kotesovec, Feb 27 2023
From Amiram Eldar, Dec 24 2024: (Start)
Equals lim_{m->oo} (1/m) * Sum_{k=1..m} k/uphi(k) = lim_{m->oo} (1/m) * Sum_{k=1..m} A319677(k)/A319676(k), where uphi(k) is the unitary totient function (A047994).
Equals lim_{m->oo} (1/log(m)) * Sum_{k=1..m} 1/uphi(k) = lim_{m->oo} (1/log(m)) * A379517(m)/A379518(m).
Equals lim_{m->oo} (1/m) * Sum_{k=1..m} A361967(k).
Equals Product_{p prime} ((1-1/p) * (1 + Sum_{k>=1} 1/(p^k-1))).
Equals Product_{p prime} (1 + (1-1/p) * Sum_{k>=1} 1/(p^k*(p^k-1))). (End)

Extensions

More digits from Vaclav Kotesovec, Jun 13 2021

A348341 a(n) is the number of noninfinitary divisors of n.

Original entry on oeis.org

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

Views

Author

Amiram Eldar, Oct 13 2021

Keywords

Examples

			a(4) = 1 since 4 has one noninfinitary divisor, 2.
		

Crossrefs

Programs

  • Mathematica
    a[1] = 0; a[n_] := DivisorSigma[0, n] - Times @@ Flatten[2^DigitCount[#, 2, 1] & /@ FactorInteger[n][[;; , 2]]]; Array[a, 100]
  • PARI
    A348341(n) = (numdiv(n)-factorback(apply(a -> 2^hammingweight(a), factorint(n)[, 2]))); \\ Antti Karttunen, Oct 13 2021

Formula

a(n) = A000005(n) - A037445(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).
Sum_{k=1..n} a(k) ~ c * n * log(n), where c = (1 - 2 * A327576) = 0.266749... . - Amiram Eldar, Dec 09 2022

A327573 Partial sums of the number of infinitary divisors function: a(n) = Sum_{k=1..n} id(k), where id is A037445.

Original entry on oeis.org

1, 3, 5, 7, 9, 13, 15, 19, 21, 25, 27, 31, 33, 37, 41, 43, 45, 49, 51, 55, 59, 63, 65, 73, 75, 79, 83, 87, 89, 97, 99, 103, 107, 111, 115, 119, 121, 125, 129, 137, 139, 147, 149, 153, 157, 161, 163, 167, 169, 173, 177, 181, 183, 191, 195, 203, 207, 211, 213, 221
Offset: 1

Views

Author

Amiram Eldar, Sep 17 2019

Keywords

Comments

Differs from A306069 at n >= 16.

References

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

Crossrefs

Cf. A006218 (all divisors), A064608 (unitary), A306069 (bi-unitary), A145353 (exponential).

Programs

  • Mathematica
    f[p_, e_] := 2^DigitCount[e, 2, 1]; id[1] = 1; id[n_] := Times @@ (f @@@ FactorInteger[n]); Accumulate[Array[id, 100]]

Formula

a(n) ~ 2 * c * n * log(n), where c = 0.366625... (A327576). [Corrected by Amiram Eldar, May 07 2021]

A344274 a(n) is the least k such that the average number of infinitary divisors of {1..k} is >= n.

Original entry on oeis.org

1, 6, 24, 105, 385, 1554, 6063, 23688, 92610, 362112, 1416360, 5539296, 21663378, 84725487, 331362185, 1295952084, 5068450464, 19822658688
Offset: 1

Views

Author

Amiram Eldar, May 13 2021

Keywords

Examples

			a(2) = 6 since the average number of infinitary divisors of {1..6} is A327573(6)/6 = 13/6 > 2.
		

Crossrefs

The infinitary version of A085829.
Similar sequences: A328331, A336304, A338891, A338943, A344272, A344273.

Programs

  • Mathematica
    f[p_, e_] := 2^DigitCount[e, 2, 1]; idivnum[1] = 1; idivnum[n_] := Times @@ (f @@@ FactorInteger[n]); seq={}; s = 0; k = 1; Do[While[s = s + idivnum[k]; s < k*n, k++]; AppendTo[seq, k]; k++, {n, 1, 10}]; seq

Formula

Lim_{n->oo} a(n+1)/a(n) = exp(1/(2*c)) = 3.9109891037..., where c is A327576.
Showing 1-4 of 4 results.