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.

A381642 Numbers k such that F(k) > F(m) for all m < k, where F(k) = A381641(k)/A381639(k).

Original entry on oeis.org

1, 6, 10, 14, 22, 26, 34, 38, 42, 46, 58, 62, 66, 78, 102, 114, 130, 170, 190, 230, 290, 310, 370, 406, 410, 430, 434, 470, 518, 574, 602, 658, 742, 826, 854, 938, 994, 1022, 1106, 1162, 1218, 1302, 1554, 1722, 1806, 1974, 2226, 2478, 2562, 2706, 2814, 2838, 2982
Offset: 1

Views

Author

Amiram Eldar, Mar 03 2025

Keywords

Comments

Called "F-champion numbers" by Erdős and Nicolas (1989).
All the terms are squarefree numbers.

Crossrefs

Programs

  • Mathematica
    f[n_] := Module[{p = FactorInteger[n][[;; , 1]]}, Total[1 - Most[p]/Rest[p]]]; seq[lim_] := Module[{fi, fmax = -1, s = {}}, Do[fi = f[i]; If[fi > fmax, fmax = fi; AppendTo[s, i]], {i, 1, lim}]; s]; seq[3000]
  • PARI
    f(n) = {my(p = factor(n)[,1]); sum(i = 1, #p-1, 1 - p[i]/p[i+1]);}
    list(lim) = {my(fm = -1, f1); for(k = 1, lim, f1 = f(k); if(f1 > fm, print1(k, ", "); fm = f1));}

A381639 Denominators of Sum_{i=1..omega(n)-1} p_{i}/p_{i+1}, where omega(n) = A001221(n) and p_1 < p_2 < ... p_omega(n) are the distinct prime factors of n; a(1) = 1.

Original entry on oeis.org

1, 1, 1, 1, 1, 3, 1, 1, 1, 5, 1, 3, 1, 7, 5, 1, 1, 3, 1, 5, 7, 11, 1, 3, 1, 13, 1, 7, 1, 15, 1, 1, 11, 17, 7, 3, 1, 19, 13, 5, 1, 21, 1, 11, 5, 23, 1, 3, 1, 5, 17, 13, 1, 3, 11, 7, 19, 29, 1, 15, 1, 31, 7, 1, 13, 33, 1, 17, 23, 35, 1, 3, 1, 37, 5, 19, 11, 39, 1
Offset: 1

Views

Author

Amiram Eldar, Mar 03 2025

Keywords

Comments

First differs from A119288 at n = 30.
First differs from {A226040(n-1)} at n = 35.
Also denominators of the fractions whose numerators are A381641.

Crossrefs

Cf. A001221, A119288, A226040, A381638 (numerators), A381640, A381641.

Programs

  • Mathematica
    a[n_] := Module[{p = FactorInteger[n][[;; , 1]]}, Denominator[Total[Most[p]/Rest[p]]]]; Array[a, 100]
  • PARI
    a(n) = {my(p = factor(n)[,1]); denominator(sum(i = 1, #p-1, p[i]/p[i+1]));}

A381638 Numerators of Sum_{i=1..omega(n)-1} p_{i}/p_{i+1}, where omega(n) = A001221(n) and p_1 < p_2 < ... p_omega(n) are the distinct prime factors of n; a(1) = 0.

Original entry on oeis.org

0, 0, 0, 0, 0, 2, 0, 0, 0, 2, 0, 2, 0, 2, 3, 0, 0, 2, 0, 2, 3, 2, 0, 2, 0, 2, 0, 2, 0, 19, 0, 0, 3, 2, 5, 2, 0, 2, 3, 2, 0, 23, 0, 2, 3, 2, 0, 2, 0, 2, 3, 2, 0, 2, 5, 2, 3, 2, 0, 19, 0, 2, 3, 0, 5, 31, 0, 2, 3, 39, 0, 2, 0, 2, 3, 2, 7, 35, 0, 2, 0, 2, 0, 23, 5
Offset: 1

Views

Author

Amiram Eldar, Mar 03 2025

Keywords

Examples

			Fractions begin with 0, 0, 0, 0, 0, 2/3, 0, 0, 0, 2/5, 0, 2/3, ...
		

Crossrefs

Cf. A000961, A001221, A381639 (denominators), A381640, A381641.

Programs

  • Mathematica
    a[n_] := Module[{p = FactorInteger[n][[;; , 1]]}, Numerator[Total[Most[p]/Rest[p]]]]; Array[a, 100]
  • PARI
    a(n) = {my(p = factor(n)[,1]); numerator(sum(i = 1, #p-1, p[i]/p[i+1]));}

Formula

a(n) = 0 if and only if n is a power of a prime (A000961).
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Sum_{q prime} (1/q^2) * Sum_{primes p < q} Product_{primes r, p < r < q} (1-1/r). This sum converges slowly: for primes q that are not exceeding 10^9, 10^10, 10^11, and 10^12, the sums are 0.5399..., 0.5447..., 0.5487..., and 0.5520..., respectively.
Showing 1-3 of 3 results.