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

A103199 Primes p such that p-1 has more divisors than any smaller prime-1.

Original entry on oeis.org

2, 3, 5, 7, 13, 31, 37, 61, 181, 241, 421, 1009, 1321, 1801, 2161, 2521, 6301, 7561, 12601, 15121, 20161, 30241, 45361, 55441, 100801, 110881, 196561, 332641, 498961, 786241, 982801, 1108801, 1580041, 1940401, 1995841, 2402401, 3880801, 4324321, 11476081, 11531521
Offset: 1

Views

Author

Don Reble, Mar 19 2005

Keywords

Comments

There are infinitely many primes p such that d(p-1) > exp(c*log(p)/log(log(p))), where d(k) is the number of divisors of k, and c > 0 is a constant (Prachar, 1955). Therefore, this sequence is infinite. - Amiram Eldar, Apr 16 2024

Crossrefs

Programs

  • Mathematica
    seq[pmax_] := Module[{d, dm = 0, s = {}, p = 1}, While[p < pmax, p = NextPrime[p]; d = DivisorSigma[0, p-1]; If[d > dm, dm = d; AppendTo[s, p]]]; s]; seq[10^6] (* Amiram Eldar, Apr 16 2024 *)
  • PARI
    lista(pmax) = {my(dm = 0, d); forprime(p = 1, pmax, d = numdiv(p-1); if(d > dm, dm = d; print1(p, ", ")));} \\ Amiram Eldar, Apr 16 2024

Extensions

a(38)-a(40) added and name clarified by Amiram Eldar, Apr 16 2024

A293059 Numbers k such that sigma(phi(k))/k > sigma(phi(m))/m for all m < k, where sigma is the sum of divisors function (A000203) and phi is Euler's totient function (A000010).

Original entry on oeis.org

1, 5, 7, 13, 31, 37, 61, 181, 241, 421, 899, 1321, 1333, 1763, 2161, 2521, 5183, 7561, 12601, 15121, 28187, 30241, 55441, 110881, 167137, 278263, 332641, 555911, 666917, 722473, 1443853, 2165407, 3607403, 4324321, 7212581, 8654539, 10817761, 21631147, 36768847
Offset: 1

Views

Author

Amiram Eldar, Oct 15 2017

Keywords

Comments

Alaoglu and Erdős proved that lim sup sigma(phi(n))/n = oo, thus this sequence is infinite.

Crossrefs

Programs

  • Mathematica
    a={}; rm=0; Do[r = DivisorSigma[1, EulerPhi[n]]/n; If[r>rm, rm=r; AppendTo[a,n]],{n,1,100000}]; a
  • PARI
    lista(nn) = {my(rmax = 0); for (n=1, nn, if ((r=sigma(eulerphi(n))/n) > rmax, rmax = r; print1(n, ", ")););} \\ Michel Marcus, Oct 18 2017
Showing 1-2 of 2 results.