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.

A249566 Numbers n such that A182134(n) = 4, i.e., there exist exactly four primes p with prime(n) < p < prime(n)^(1+1/n).

Original entry on oeis.org

17, 19, 24, 26, 32, 33, 35, 36, 37, 38, 40, 42, 43, 47, 50, 51, 52, 58, 62, 63, 64, 76, 77, 78, 79, 90, 91, 93, 95, 121, 123, 124, 125, 126, 134, 135, 137, 150, 153, 185, 186, 187, 188, 189, 201, 203, 213, 218, 219, 238, 239, 259, 263, 278, 279, 289, 293
Offset: 1

Views

Author

Robert Price, Nov 01 2014

Keywords

Comments

See A246782 for a more complete description of this sequence.
a(1136) > 10^12.
It is interesting that three consecutive integers n = 20004097201301075, n + 1 and n + 2 are in the sequence. Conjecture: The sequence is infinite. - Farideh Firoozbakht, Nov 01 2014

Crossrefs

Programs

  • Haskell
    a249566 n = a249566_list !! (n-1)
    a249566_list = filter ((== 4) . a182134) [1..]
    -- Reinhard Zumkeller, Nov 17 2014
  • Mathematica
    np[n_]:=(a = Prime[n]; b = a^(1 + 1/n); Length[Select[Range[a+1,b], PrimeQ]]); Do[If[np[n] == 4,Print[n]], {n, 293}]
    np[n_]:=(a = Prime[n]; b = a^(1 + 1/n); Length[Select[Range[a+1,b], PrimeQ]]); Select[Range[293], np[#]==4&] (* Farideh Firoozbakht, Nov 01 2014 *)
  • PARI
    for(n=1,9e9,primepi(prime(n)^(1+1/n))-n==4&&print1(n",")) \\ M. F. Hasler, Nov 03 2014