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.

A220462 Chebyshev numbers C_v(n) for v=3/2: a(n) is the smallest number such that if x>=a(n), then theta(x)-theta(2*x/3)>=n*log(x), where theta(x) = sum_{prime p<=x} log p.

Original entry on oeis.org

13, 37, 41, 67, 73, 97, 127, 137, 173, 179, 181, 211, 229, 239, 263, 307, 311, 347, 367, 379, 431, 433, 443, 449, 479, 487, 541, 563, 587, 599, 607, 641, 643, 673, 739, 757, 787, 797, 809, 823, 827, 859, 937, 967, 997, 1019, 1031, 1039, 1049, 1061, 1087
Offset: 1

Views

Author

Keywords

Comments

All terms are primes.
Up to a(97)=2333, only four terms of the sequence (a(33)=643, a(34)=673, a(76)=1721 and a(77)=1741) are not (3/2)-Ramanujan numbers as in Shevelev's link; up to 2333, the only (3/2)-Ramanujan numbers missing from the sequence are 2, 617, 653, 709, 1709, 1733, and 1747.

Crossrefs

Cf. A220293.

Programs

  • Mathematica
    (* Assuming range of x is from a(n) to 2*a(n) *) theta[x_] := Sum[Log[p], {p, Table[Prime[k], {k, 1, PrimePi[x]}]}]; Clear[a]; a[0] = 2; a[n_] := a[n] = (t = Table[{an, x >= an && theta[x] - theta[2*(x/3)] >= n*Log[x]}, {an, a[n - 1], Prime[4*(n + 1)]}, {x, an, 2*an}]; sp = t // Flatten[#, 1] & // Sort // Split[#, #1[[1]] == #2[[1]] &] &; Select[sp, And @@ (#[[All, 2]]) &] // First // First // First); Table[Print[a[n]]; a[n], {n, 1, 51}] (* Jean-François Alcover, Jan 24 2013 *)

Formula

a(n)<=prime(4*(n+1)).