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-3 of 3 results.

A063974 Number of terms in inverse set of usigma = sum of unitary divisors = A034448.

Original entry on oeis.org

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

Views

Author

Labos Elemer, Sep 05 2001

Keywords

Examples

			usigma(x) = 288, invusigma(288) = {138,154,165,168,213,235,248,253}, so a(288) = 8, the number of all terms in the inverse set.
		

Crossrefs

Formula

Size of set {x; usigma(x) = n}.
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = A308041. - Amiram Eldar, Dec 23 2024

A379513 Numerators of the partial sums of the reciprocals of the sum of unitary divisors function (A034448).

Original entry on oeis.org

1, 4, 19, 107, 39, 61, 259, 817, 853, 97, 301, 307, 2209, 187, 2279, 39583, 121129, 122557, 124699, 126127, 509863, 171541, 173921, 526523, 6930479, 6983519, 7063079, 7118771, 7193027, 802663, 405199, 13495327, 1131701, 30726097, 123670153, 622026437, 11910394103
Offset: 1

Views

Author

Amiram Eldar, Dec 23 2024

Keywords

Examples

			Fractions begin with 1, 4/3, 19/12, 107/60, 39/20, 61/30, 259/120, 817/360, 853/360, 97/40, 301/120, 307/120, ...
		

Crossrefs

Cf. A034448, A064609, A370898, A379514 (denominators), A379515.

Programs

  • Mathematica
    usigma[n_] := Times @@ (1 + Power @@@ FactorInteger[n]); usigma[1] = 1; Numerator[Accumulate[Table[1/usigma[n], {n, 1, 50}]]]
  • PARI
    usigma(n) = {my(f = factor(n)); prod(i = 1, #f~, 1 + f[i, 1]^f[i, 2]);}
    list(nmax) = {my(s = 0); for(k = 1, nmax, s += 1 / usigma(k); print1(numerator(s), ", "))};

Formula

a(n) = numerator(Sum_{k=1..n} 1/A034448(k)).
a(n)/A379514(n) = B * log(n) + D + O(log(n)^(5/3) * log(log(n))^(4/3) / n), where B = A308041, D = B * (gamma + A1 - A2), gamma = A001620, A1 = Sum_{p prime} ((p*C(p)*log(p)/(p-1)) * Sum_{k>=1} (k/(p^k*(p^(k+1)+1)))), A2 = Sum_{p prime} ((C(p)*log(p)/p^2) * Sum_{k>=0} (1/(p^k*(p^(k+1)+1)))), and C(p) = 1 - (p/(p-1)) * Sum_{k>=1} (1/(p^k*(p^(k+1)+1))) (Sita Ramaiah and Suryanarayana, 1980).

A379515 Numerators of the partial alternating sums of the reciprocals of the sum of unitary divisors function (A034448).

Original entry on oeis.org

1, 2, 11, 43, 53, 4, 37, 293, 329, 103, 113, 107, 809, 129, 809, 12913, 41119, 39691, 41833, 8081, 33395, 32443, 33871, 10973, 148361, 48275, 7149, 34861, 108119, 319937, 164941, 1761311, 112361, 662011, 5405483, 26502319, 516671461, 508357441, 3620857237, 3556192637
Offset: 1

Views

Author

Amiram Eldar, Dec 23 2024

Keywords

Examples

			Fractions begin with 1, 2/3, 11/12, 43/60, 53/60, 4/5, 37/40, 293/360, 329/360, 103/120, 113/120, 107/120, ...
		

Crossrefs

Programs

  • Mathematica
    usigma[n_] := Times @@ (1 + Power @@@ FactorInteger[n]); usigma[1] = 1; Numerator[Accumulate[Table[(-1)^(n+1)/usigma[n], {n, 1, 50}]]]
  • PARI
    usigma(n) = {my(f = factor(n)); prod(i = 1, #f~, 1 + f[i, 1]^f[i, 2]);}
    list(nmax) = {my(s = 0); for(k = 1, nmax, s += (-1)^(k+1) / usigma(k); print1(numerator(s), ", "))};

Formula

a(n) = numerator(Sum_{k=1..n} (-1)^(k+1)/A034448(k)).
a(n)/A379516(n) = E * log(n) + F + O(log(n)^(5/3) * log(log(n))^(4/3) / n^u), where u > 0, E = A308041 * (2/(A323482 + 1/2) - 1) = 0.10259754363391420806..., and F is a constant.
Showing 1-3 of 3 results.