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

A284411 Least prime p such that more than half of all integers are divisible by n distinct primes not greater than p.

Original entry on oeis.org

3, 37, 42719, 5737850066077
Offset: 1

Views

Author

Peter Munn, Mar 26 2017

Keywords

Comments

The proportion of all integers that satisfy the divisibility criterion for p=prime(m) is determined using the proportion that satisfy it over any interval of primorial(m)=A002110(m) integers.
a(4) is from De Koninck, 2009; calculation credited to David Grégoire.
a(5) is about 7.887*10^34 assuming the Riemann Hypothesis, and about 7*10^34 unconditionally (De Koninck and Tenenbaum, 2002). - Amiram Eldar, Dec 05 2024

Examples

			Exactly half of the integers are divisible by 2, so a(1)>2. Two-thirds of all integers are divisible by 2 or 3, so a(1) = 3.
		

References

  • Jean-Marie De Koninck, Those Fascinating Numbers, American Mathematical Society, 2009, pp. 13, 216 and 368.

Crossrefs

Formula

a(n) is least p=prime(m) such that 2*Sum_{k=0..n-1} A096294(m,k) < A002110(m).
log(log(a(n))) = n - b + O(1/sqrt(n)), where b = 1/3 + A077761 (De Koninck and Tenenbaum, 2002). - Amiram Eldar, Dec 05 2024

Extensions

Definition edited by N. J. A. Sloane, Apr 01 2017

A342480 a(n) is the denominator of the asymptotic density of numbers whose second smallest prime divisor (A119288) is prime(n).

Original entry on oeis.org

1, 6, 10, 15, 1155, 1365, 12155, 1616615, 37182145, 11849255, 33426748355, 247357937827, 10141675450907, 25652473199353, 2928046583754721, 155186468939000213, 223317113839049087, 558516101711461766587, 796182527971658263007, 241532826894674874877669, 430046252763689411367557
Offset: 1

Views

Author

Amiram Eldar, Mar 13 2021

Keywords

Comments

See A342479 for details.

Crossrefs

Cf. A038110, A038111, A119288, A342479 (numerators).

Programs

  • Mathematica
    f[n_] := Module[{p = Prime[n], q}, q = Select[Range[p - 1], PrimeQ]; Plus @@ (1/(q - 1))*Times @@ ((q - 1)/q)/p]; Denominator @ Array[f, 30]

A378720 a(n) is the numerator of the asymptotic density of numbers whose third smallest prime divisor is prime(n).

Original entry on oeis.org

0, 0, 1, 1, 4, 326, 628, 992, 98304, 125568, 733440, 281163264, 386427322368, 3178249003008, 12454223855616, 6450728943845376, 342348724735967232, 20218431581110665216, 39814891891080560640, 82739188294287768944640, 15336676441718784000, 61298453882755419734016000
Offset: 1

Views

Author

Robert G. Wilson v and Amiram Eldar, Dec 05 2024

Keywords

Comments

The third smallest prime divisor of a number k is the third member in the ordered list of the distinct prime divisors of k. Only numbers in A000977 have a third smallest prime divisor.
The partial sums of the fractions first exceed 1/2 after summing 4467 terms. Therefore, the median value of the distribution of the third prime divisor is prime(4467) = 42719 = A284411(3).

Examples

			The fractions begin with 0/1, 0/1, 1/30, 1/30, 4/165, 326/15015, 628/36465, 992/62985, 98304/7436429, 125568/11849255, ..., .
a(1) = a(2) = 0 since there are no numbers whose third prime divisor is 2 or 3.
a(3)/A378721(3) = 1/30 since the numbers whose third prime divisor is 5 are the numbers that are divisible by 2, 3 and 5, and their density if (1/2)*(1/3)*(1/5) = 1/30.
a(4)/A378721(4) = 1/30 since the numbers whose third prime divisor is 7 are the union of the numbers that are divisible by 2, 3 and 7 and not by 5 whose density is (1/2)*(1/3)*(1-1/5)*(1/7) = 2/105, the numbers that are divisible by 2, 5 and 7 and not by 3 whose density is (1/2)*(1-1/3)*(1/5)*(1/7) = 1/105, and the numbers that are divisible by 3, 5 and 7 and not by 2 whose density is (1-1/2)*(1/3)*(1/5)*(1/7) = 1/210, and 2/105 + 1/105 + 1/210 = 1/30.
		

References

  • József Sándor and Borislav Crstici, Handbook of Number theory II, Kluwer Academic Publishers, 2004, Chapter 4, pp. 337-341.

Crossrefs

Programs

  • Mathematica
    a[n_] := Block[{p, q = Prime@ Range@ n}, p = Fold[Times, 1, q]; q = Most@ q; Plus @@ Times @@@ Subsets[q -1, {n -3}]/p]; a[1] = 0; Numerator@ Array[a, 22]
  • PARI
    a(n) = {my(v = primes(n), q = vecextract(apply(x -> x-1, v),"^-1"), p = vecprod(v), prd = vecprod(q)/p, sm = 0, sb); forsubset([#q, 2], s, sb = vecextract(q, s); sm += 1/vecprod(sb)); numerator(prd * sm);}

Formula

a(n)/A378721(n) = (1/prime(n)#) * (Product_{k=1..n-1} (prime(k) - 1)) * Sum_{j=1..n-1, i=1..j-1} 1/((prime(i)-1)*(prime(j)-1)), where prime(n)# = A002110(n) is the n-th primorial number.
Sum_{n>=1} a(n)/A378721(n) = 1.
Sum_{n=1..m} a(n)/A378721(n) > 1/2 for m >= 4467 = primepi(A284411(3)).

A378721 a(n) is the denominator of the asymptotic density of numbers whose third smallest prime divisor is prime(n).

Original entry on oeis.org

1, 1, 30, 30, 165, 15015, 36465, 62985, 7436429, 11849255, 73465381, 33426748355, 50708377254535, 436092044389001, 1863302371480277, 1086305282573001491, 64092011671807087969, 3909612711980232366109, 8449808119441147371913, 18598027670889965365580513, 3543193335582015099413
Offset: 1

Views

Author

Robert G. Wilson v and Amiram Eldar, Dec 05 2024

Keywords

Comments

See A378720 for more details.

Crossrefs

Cf. A000040, A038110, A038111, A342479, A342480, A378720 (numerators).

Programs

  • Mathematica
    a[n_] := Block[{p, q = Prime@ Range@ n}, p = Fold[Times, 1, q]; q = Most@ q; Plus @@ Times @@@ Subsets[q -1, {n -3}]/p]; a[1] = 0; Denominator@ Array[a, 21]
  • PARI
    a(n) = {my(v = primes(n), q = vecextract(apply(x -> x-1, v),"^-1"), p = vecprod(v), prd = vecprod(q)/p, sm = 0, sb); forsubset([#q, 2], s, sb = vecextract(q, s); sm += 1/vecprod(sb)); denominator(prd * sm);}

A356093 a(n) = numerator((prime(n)-1)/prime(n)#), where prime(n)# = A002110(n) is the n-th primorial.

Original entry on oeis.org

1, 1, 2, 1, 1, 2, 8, 3, 1, 2, 1, 6, 4, 1, 1, 2, 1, 2, 1, 1, 12, 1, 1, 4, 16, 10, 1, 1, 18, 8, 3, 1, 4, 1, 2, 5, 2, 27, 1, 2, 1, 6, 1, 32, 14, 3, 1, 1, 1, 2, 4, 1, 8, 25, 128, 1, 2, 9, 2, 4, 1, 2, 3, 1, 4, 2, 1, 8, 1, 2, 16, 1, 1, 2, 9, 1, 2, 6, 40, 4, 1, 2, 1
Offset: 1

Views

Author

Amiram Eldar, Jul 26 2022

Keywords

Comments

f(n) = a(n)/A356094(n) is the asymptotic density of numbers k such that prime(n) = A053669(k) is the smallest prime not dividing k.
The asymptotic mean of A053669 is 2.9200509773... (A249270) which is the weighted mean of the primes with f(n) as weights. The corresponding asymptotic standard deviation, which can be evaluated from the second raw moment Sum_{n>=1} f(n) * prime(n)^2, is 2.8013781465... .

Examples

			Fractions begin with 1/2, 1/3, 2/15, 1/35, 1/231, 2/5005, 8/255255, 3/1616615, 1/10140585, 2/462120945, ...
		

Crossrefs

Cf. A002110, A006093, A039787, A053669, A249270, A356094 (denominators).
Similar sequences: A038110, A338559, A340818, A341431, A342450, A342479.

Programs

  • Mathematica
    primorial[n_] := Product[Prime[i], {i, 1, n}]; Numerator[Table[(Prime[i] - 1)/primorial[i], {i, 1, 100}]]
  • PARI
    a(n) = numerator((prime(n)-1)/factorback(primes(n))); \\ Michel Marcus, Jul 26 2022
    
  • Python
    from math import gcd
    from sympy import prime, primorial
    def A356093(n): return (p:=prime(n)-1)//gcd(p,primorial(n)) # Chai Wah Wu, Jul 26 2022

Formula

a(n) = 1 iff prime(n) is in A039787.
Let f(n) = a(n)/A356094(n):
f(n) = A006093(n)/A002110(n).
Sum_{n>=1} f(n) = 1.
Sum_{n>=1} f(n) * prime(n) = A249270.
Showing 1-5 of 5 results.