A080359 -id:A080359 - cp's OEIS frontend

A164554 Ramanujan primes A104272(n) for which A104272(n) = A080359(n).

2, 71, 101, 181, 239, 241, 269, 349, 373, 409, 419, 433, 439, 491, 593, 599, 601, 607, 647, 653, 659, 823, 827, 857, 947, 1021, 1031, 1061, 1063, 1091, 1103, 1301, 1427, 1429, 1447, 1451, 1489, 1553, 1559, 1567, 1601, 1607, 1609, 1789, 1867, 1871, 1913, 1999, 2003

Offset: 1

Vladimir Shevelev, Aug 15 2009

nonn

For every n>=1, A104272(n) >= A080359(n), and the sequence shows where the inequality becomes an equality.

Let prime(m) < a(n)/2 < prime(m+1); then there exist primes p

For example, a(2) = 71, 31 < a(2)/2 < 37 and intervals (62,71), (71,74) contain the primes p = 67 and q = 73 respectively.

Let us call a prime p compatible with another prime q, if the intervals (p/2,q/2) and (p,q], if q>p, (or intervals (q/2,p/2) and (q,p], if qVladimir Shevelev, Apr 25 2012]

			a(2)=71, such that 31<71/2<37, and we see that p=67 is in interval (62, 71) and q=73 is in interval (71, 74).

V. Shevelev, On critical small intervals containing primes, arXiv:0908.2319 [math.NT] [From _Vladimir Shevelev_, Aug 20 2009]
V. Shevelev, Ramanujan and Labos primes, their generalizations, and classifications of primes, J. Integer Seq. 15 (2012) Article 12.5.4
J. Sondow, J. W. Nicholson, and T. D. Noe, Ramanujan Primes: Bounds, Runs, Twins, and Gaps, J. Integer Seq. 14 (2011) Article 11.6.2

Cf. A104272, A080359, A164368, A164333, A164288, A164294, A212493, A212541.

nn = 200; t = Table[0, {nn+1}]; s = 0;
Do[If[PrimeQ[k], s++]; If[PrimeQ[k/2], s--]; If[s <= nn && t[[s+1]] == 0, t[[s+1]] = k], {k, Prime[3nn]}
];
A080359 = Rest[t];
R = Table[0, {nn}]; s = 0;
Do[If[PrimeQ[k], s++]; If[PrimeQ[k/2], s--]; If[s < nn, R[[s+1]] = k], {k, Prime[3nn]}
];
A104272 = R+1;
Intersection[A104272, A080359] (* Jean-François Alcover, Oct 28 2018, after T. D. Noe in A104272 *)

All solutions of the equation A104272(x)=A080359(x) are x=pi(a(n))-pi(a(n)/2). - Vladimir Shevelev, Apr 25 2012

Terms beyond 659 from R. J. Mathar, Dec 17 2009

A194217 a(n) = A104272(n)-A080359(n).

Original entry on oeis.org

0, 8, 4, 10, 10, 4, 6, 6, 0, 24, 0, 4, 18, 36, 12, 10, 6, 0, 36, 36, 34, 0, 0, 12, 0, 10, 24, 18, 34, 0, 14, 0, 22, 0, 0, 10, 0, 0, 18, 24, 0, 4, 60, 48, 10, 0, 0, 0, 0, 28, 24, 0, 0, 0, 16, 36, 36, 6, 8, 12, 36, 10, 0, 0, 24, 0, 22, 54, 30, 0, 14, 12, 18, 22

Offset: 1

Vladimir Shevelev, Aug 18 2011

nonn

Conjecture: Asymptotic density of nonzero terms is 3/4.

Alois P. Heinz, Table of n, a(n) for n = 1..1000
V. Shevelev, Ramanujan and Labos primes, their generalizations and classifications of primes
J. Sondow, Ramanujan Prime. (MathWorld)

Cf. A104272, A080359, A193507, A194184, A194186.

nn = 100;
R = Table[0, {nn}]; s = 0;
Do[If[PrimeQ[k], s++]; If[PrimeQ[k/2], s--]; If[s < nn, R[[s+1]] = k], {k, Prime[3nn]}
];
A104272 = R = R + 1;
T = Table[0, {nn + 1}]; s = 0;
Do[If[PrimeQ[k], s++]; If[PrimeQ[k/2], s--]; If[s <= nn && T[[s+1]] == 0, T[[s+1]] = k], {k, Prime[3nn]}
];
A080359 = Rest[T];
A104272 - A080359 (* Jean-François Alcover, Aug 19 2018, after T. D. Noe *)

A164372 Greater of twin primes (A006512) which are not Labos primes (A080359).

Original entry on oeis.org

5, 7, 151, 229, 571, 643, 1051, 1153, 1669, 1723, 1951, 2029, 2131, 2239, 2311, 2593, 2659, 3001, 3121, 3169, 3253, 3583, 3769, 4003, 4219, 4231, 4483, 4549, 4723, 4789, 5641, 6451, 6553, 6661, 6763, 6949, 6961, 7129, 7351, 8011, 9043, 9463, 9631, 10009

Offset: 1

Vladimir Shevelev, Aug 14 2009

nonn

The terms greater than 7 are in A164294.

V. Shevelev, Ramanujan and Labos primes, their generalizations, and classifications of primes, J. Integer Seq. 15 (2012) Article 12.5.4

Cf. A006512, A080359, A164294, A164371, A104272, A001359.

lista(nn)= {my(vlp = readvec("/gp/bfiles/b080359.txt")); forprime (p=3, nn, if (isprime(p-2) && !vecsearch(vlp, p), print1(p, ", ")););} \\ Michel Marcus, Jan 15 2014

More terms from Michel Marcus, Jan 15 2014

A182391 Numbers n for which A104272(n) = A080359(n).

Original entry on oeis.org

1, 9, 11, 18, 22, 23, 25, 30, 32, 34, 35, 37, 38, 41, 46, 47, 48, 49, 52, 53, 54, 63, 64, 66, 70, 75, 76, 79, 80, 82, 84, 94, 98, 99, 101, 102, 105, 108, 109, 110, 113, 114, 115, 124, 127, 128, 131, 135, 136, 139, 140, 148, 149, 150, 151, 154, 156, 158, 160

Offset: 1

Vladimir Shevelev, Apr 27 2012

nonn

Number m is in the sequence iff 1) there exists only composite number k such that 2*k-1 is prime and A060715(k)=m; 2) there is no prime p such that 2*p-1 is prime and A060715(p)=m-1.

Cf. A104272, A080359, A164554, A193507, A194184, A194186, A194217, A182366.

A194217(n)=0.

A182392 Numbers n for which there exists only composite number k such that A060715(k) = n and 2*k-1 is prime, but A104272(n) differs from A080359(n).

Original entry on oeis.org

3, 8, 36, 55, 58, 83, 129, 134, 143, 155, 186, 197, 207, 218, 269, 295, 309, 310, 361, 362, 380, 396, 412, 454, 466, 473, 505, 511, 514, 544, 549, 556, 563, 616, 631, 660, 666, 677, 683, 697, 771, 781, 788, 797, 812, 873, 874, 881, 883, 894, 906, 953

Offset: 1

Vladimir Shevelev and Peter J. C. Moses, Apr 27 2012

nonn

There exists a prime p=p(n) such that 2*p-1 is prime and A060715(p)=a(n)-1 (cf. comment in A182391).

Cf. A104272, A080359, A164554, A193507, A194184, A194186, A194217, A182366, A182391.

A050376 "Fermi-Dirac primes": numbers of the form p^(2^k) where p is prime and k >= 0.

Original entry on oeis.org

2, 3, 4, 5, 7, 9, 11, 13, 16, 17, 19, 23, 25, 29, 31, 37, 41, 43, 47, 49, 53, 59, 61, 67, 71, 73, 79, 81, 83, 89, 97, 101, 103, 107, 109, 113, 121, 127, 131, 137, 139, 149, 151, 157, 163, 167, 169, 173, 179, 181, 191, 193, 197, 199, 211, 223, 227, 229, 233, 239, 241

Offset: 1

Christian G. Bower, Nov 15 1999

Every number n is a product of a unique subset of these numbers. This is sometimes called the Fermi-Dirac factorization of n (see A182979). Proof: In the prime factorization n = Product_{j>=1} p(j)^e(j) expand every exponent e(j) as binary number and pick the terms of this sequence corresponding to the positions of the ones in binary (it is clear that both n and n^2 have the same number of factors in this sequence, and that each factor appears with exponent 1 or 0).
Or, a(1) = 2; for n>1, a(n) = smallest number which cannot be obtained as the product of previous terms. This is evident from the unique factorization theorem and the fact that every number can be expressed as the sum of powers of 2. - Amarnath Murthy, Jan 09 2002
Except for the first term, same as A084400. - David Wasserman, Dec 22 2004
The least number having 2^n divisors (=A037992(n)) is the product of the first n terms of this sequence according to Ramanujan.
According to the Bose-Einstein distribution of particles, an unlimited number of particles may occupy the same state. On the other hand, according to the Fermi-Dirac distribution, no two particles can occupy the same state (by the Pauli exclusion principle). Unique factorizations of the positive integers by primes (A000040) and over terms of A050376 one can compare with two these distributions in physics of particles. In the correspondence with this, the factorizations over primes one can call "Bose-Einstein factorizations", while the factorizations over distinct terms of A050376 one can call "Fermi-Dirac factorizations". - Vladimir Shevelev, Apr 16 2010
The numbers of the form p^(2^k), where p is prime and k >= 0, might thus be called the "Fermi-Dirac primes", while the classic primes might be called the "Bose-Einstein primes". - Daniel Forgues, Feb 11 2011
In the theory of infinitary divisors, the most natural name of the terms is "infinitary primes" or "i-primes". Indeed, n is in the sequence, if and only if it has only two infinitary divisors. Since 1 and n are always infinitary divisors of n>1, an i-prime has no other infinitary divisors. - Vladimir Shevelev, Feb 28 2011
{a(n)} is the minimal set including all primes and closed with respect to squaring. In connection with this, note that n and n^2 have the same number of factors in their Fermi-Dirac representations. - Vladimir Shevelev, Mar 16 2012
In connection with this sequence, call an integer compact if the factors in its Fermi-Dirac factorization are pairwise coprime. The density of such integers equals (6/Pi^2)*Product_{prime p} (1+(Sum_{i>=1} p^(-(2^i-1))/(p+1))) = 0.872497... It is interesting that there exist only 7 compact factorials listed in A169661. - Vladimir Shevelev, Mar 17 2012
The first k terms of the sequence solve the following optimization problem:
Let x_1, x_2,..., x_k be integers with the restrictions: 2<=x_1A064547(Product{i=1..k} x_i) >= k.  Let the goal function be Product_{i=1..k} x_i. Then the minimal value of the goal function is Product_{i=1..k} a(i). - Vladimir Shevelev, Apr 01 2012
From Joerg Arndt, Mar 11 2013: (Start)
Similarly to the first comment, for the sequence "Numbers of the form p^(3^k) or p^(2*3^k) where p is prime and k >= 0" one obtains a factorization into distinct factors by using the ternary expansion of the exponents (here n and n^3 have the same number of such factors).
The generalization to base r would use "Numbers of the form p^(r^k), p^(2*r^k), p^(3*r^k), ..., p^((r-1)*r^k) where p is prime and k >= 0" (here n and n^r have the same number of (distinct) factors).  (End)
The first appearance of this sequence as a multiplicative basis in number theory with some new notions, formulas and theorems may have been in my 1981 paper (see the Abramovich reference). - Vladimir Shevelev, Apr 27 2014
Numbers n for which A064547(n) = 1. - Antti Karttunen, Feb 10 2016
Lexicographically earliest sequence of distinct nonnegative integers such that no term is a product of 2 or more distinct terms. Removing the distinctness requirement, the sequence becomes A000040 (the prime numbers); and the equivalent sequence where the product is of 2 distinct terms is A026416 (without its initial term, 1). - Peter Munn, Mar 05 2019
The sequence was independently developed as a multiplicative number system in 1985-1986 (and first published in 1995, see the Uhlmann reference) using a proof method involving representations of positive integers as sums of powers of 2. This approach offers an arguably simpler and more flexible means for analyzing the sequence. - Jeffrey K. Uhlmann, Nov 09 2022

			Prime powers which are not terms of this sequence:
  8 = 2^3 = 2^(1+2), 27 = 3^3 = 3^(1+2), 32 = 2^5 = 2^(1+4),
  64 = 2^6 = 2^(2+4), 125 = 5^3 = 5^(1+2), 128 = 2^7 = 2^(1+2+4)
"Fermi-Dirac factorizations":
  6 = 2*3, 8 = 2*4, 24 = 2*3*4, 27 = 3*9, 32 = 2*16, 64 = 4*16,
  108 = 3*4*9, 120 = 2*3*4*5, 121 = 121, 125 = 5*25, 128 = 2*4*16.

V. S. Abramovich, On an analog of the Euler function, Proceeding of the North-Caucasus Center of the Academy of Sciences of the USSR (Rostov na Donu) (1981) No. 2, 13-17 (Russian; MR0632989(83a:10003)).
S. Ramanujan, Highly Composite Numbers, Collected Papers of Srinivasa Ramanujan, p. 125, Ed. G. H. Hardy et al., AMS Chelsea 2000.
V. S. Shevelev, Multiplicative functions in the Fermi-Dirac arithmetic, Izvestia Vuzov of the North-Caucasus region, Nature sciences 4 (1996), 28-43 (in Russian; MR 2000f: 11097, pp. 3912-3913).
J. K. Uhlmann, Dynamic map building and localization: new theoretical foundations, Doctoral Dissertation, University of Oxford, Appendix 16, 1995.

T. D. Noe and Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
Steven R. Finch, Unitarism and Infinitarism, February 25, 2004. [Cached copy, with permission of the author]
Simon Litsyn and Vladimir Shevelev, On factorization of integers with restrictions on the exponent, INTEGERS: Electronic Journal of Combinatorial Number Theory, 7 (2007), #A33, 1-36.
OEIS Wiki, "Fermi-Dirac representation" of n.
Vladimir Shevelev, Compact integers and factorials, Acta Arith. 126 (2007), no.3, 195-236.
J. K. Uhlmann, Appendix 16, Doctoral Dissertation, University of Oxford, page 243, 1995.

Cf. A000040 (primes, is a subsequence), A026416, A026477, A037992 (partial products), A050377-A050380, A052330, A064547, A066724, A084400, A176699, A182979.
Cf. A268388 (complement without 1).
Cf. A124010, subsequence of A000028, A000961, A213925, A223490.
Cf. A228520, A186945 (Fermi-Dirac analog of Ramanujan primes, A104272, and Labos primes, A080359).
Cf. also A268385, A268391, A268392.

a050376 n = a050376_list !! (n-1)
a050376_list = filter ((== 1) . a209229 . a100995) [1..]
-- Reinhard Zumkeller, Mar 19 2013

isA050376 := proc(n)
    local f,e;
    f := ifactors(n)[2] ;
    if nops(f) = 1 then
        e := op(2,op(1,f)) ;
        if isA000079(e) then
            true;
        else
            false;
        end if;
    else
        false;
    end if;
end proc:
A050376 := proc(n)
    option remember ;
    local a;
    if n = 1 then
        2 ;
    else
        for a from procname(n-1)+1 do
            if isA050376(a) then
                return a;
            end if;
        end do:
    end if;
end proc: # R. J. Mathar, May 26 2017

nn = 300; t = {}; k = 1; While[lim = nn^(1/k); lim > 2,  t = Join[t, Prime[Range[PrimePi[lim]]]^k]; k = 2 k]; t = Union[t] (* T. D. Noe, Apr 05 2012 *)

{a(n)= local(m, c, k, p); if(n<=1, 2*(n==1), n--; c=0; m=2; while( cMichael Somos, Apr 15 2005; edited by Michel Marcus, Aug 07 2021

lst(lim)=my(v=primes(primepi(lim)),t); forprime(p=2,sqrt(lim),t=p; while((t=t^2)<=lim,v=concat(v,t))); vecsort(v) \\ Charles R Greathouse IV, Apr 10 2012

is_A050376(n)=2^#binary(n=isprimepower(n))==n*2 \\ M. F. Hasler, Apr 08 2015

ispow2(n)=n && n>>valuation(n,2)==1
is(n)=ispow2(isprimepower(n)) \\ Charles R Greathouse IV, Sep 18 2015

isok(n)={my(e=isprimepower(n)); e && !bitand(e,e-1)} \\ Andrew Howroyd, Oct 16 2024

from sympy import isprime, perfect_power
def ok(n):
  if isprime(n): return True
  answer = perfect_power(n)
  if not answer: return False
  b, e = answer
  if not isprime(b): return False
  while e%2 == 0: e //= 2
  return e == 1
def aupto(limit):
  alst, m = [], 1
  for m in range(1, limit+1):
    if ok(m): alst.append(m)
  return alst
print(aupto(241)) # Michael S. Branicky, Feb 03 2021

from sympy import primepi, integer_nthroot
def A050376(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        kmin = kmax >> 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x): return n+x-sum(primepi(integer_nthroot(x,1<Chai Wah Wu, Feb 18-19 2025

(define A050376 (MATCHING-POS 1 1 (lambda (n) (= 1 (A064547 n)))))
;; Requires also my IntSeq-library. - Antti Karttunen, Feb 09 2016

From Vladimir Shevelev, Mar 16 2012: (Start)
Product_{i>=1} a(i)^k_i = n!, where k_i = floor(n/a(i)) - floor(n/a(i)^2) + floor(n/a(i)^3) - floor(n/a(i)^4) + ...
Denote by A(x) the number of terms not exceeding x.
Then A(x) = pi(x) + pi(x^(1/2)) + pi(x^(1/4)) + pi(x^(1/8)) + ...
Conversely, pi(x) = A(x) - A(sqrt(x)). For example, pi(37) = A(37) - A(6) = 16-4 = 12. (End)
A209229(A100995(a(n))) = 1. - Reinhard Zumkeller, Mar 19 2013
From Vladimir Shevelev, Aug 31 2013: (Start)
A Fermi-Dirac analog of Euler product: Zeta(s) = Product_{k>=1} (1+a(k)^(-s)), for s > 1.
In particular, Product_{k>=1} (1+a(k)^(-2)) = Pi^2/6. (End)
a(n) = A268385(A268392(n)). [By their definitions.] - Antti Karttunen, Feb 10 2016
A000040 union A001248 union A030514 union A179645 union A030635 union .... - R. J. Mathar, May 26 2017

Edited by Charles R Greathouse IV, Mar 17 2010

More examples from Daniel Forgues, Feb 09 2011

A104272 Ramanujan primes R_n: a(n) is the smallest number such that if x >= a(n), then pi(x) - pi(x/2) >= n, where pi(x) is the number of primes <= x.

Original entry on oeis.org

2, 11, 17, 29, 41, 47, 59, 67, 71, 97, 101, 107, 127, 149, 151, 167, 179, 181, 227, 229, 233, 239, 241, 263, 269, 281, 307, 311, 347, 349, 367, 373, 401, 409, 419, 431, 433, 439, 461, 487, 491, 503, 569, 571, 587, 593, 599, 601, 607, 641, 643, 647, 653, 659

Offset: 1

Jonathan Sondow, Feb 27 2005

Referring to his proof of Bertrand's postulate, Ramanujan (1919) states a generalization: "From this we easily deduce that pi(x) - pi(x/2) >= 1, 2, 3, 4, 5, ..., if x >= 2, 11, 17, 29, 41, ..., respectively." Since the a(n) are prime (by their minimality), I call them "Ramanujan primes."
See the additional references and links mentioned in A143227.
2n log 2n < a(n) < 4n log 4n for n >= 1, and prime(2n) < a(n) < prime(4n) if n > 1. Also, a(n) ~ prime(2n) as n -> infinity.
Shanta Laishram has proved that a(n) < prime(3n) for all n >= 1.
a(n) - 3n log 3n is sometimes positive, but negative with increasing frequency as n grows since a(n) ~ 2n log 2n. There should be a constant m such that for n >= m we have a(n) < 3n log 3n.
A good approximation to a(n) = R_n for n in [1..1000] is A162996(n) = round(k*n * (log(k*n)+1)), with k = 2.216 determined empirically from the first 1000 Ramanujan primes, which approximates the {k*n}-th prime number which in turn approximates the n-th Ramanujan prime and where abs(A162996(n) - R_n) < 2 * sqrt(A162996(n)) for n in [1..1000]. Since R_n ~ prime(2n) ~ 2n * (log(2n)+1) ~ 2n * log(2n), while A162996(n) ~ prime(k*n) ~ k*n * (log(k*n)+1) ~ k*n * log(k*n), A162996(n) / R_n ~ k/2 = 2.216/2 = 1.108 which implies an asymptotic overestimate of about 10% (a better approximation would need k to depend on n and be asymptotic to 2). - Daniel Forgues, Jul 29 2009
Let p_n be the n-th prime. If p_n >= 3 is in the sequence, then all integers (p_n+1)/2, (p_n+3)/2, ..., (p_(n+1)-1)/2 are composite numbers. - Vladimir Shevelev, Aug 12 2009
Denote by q(n) the prime which is the nearest from the right to a(n)/2. Then there exists a prime between a(n) and 2q(n). Converse, generally speaking, is not true, i.e., there exist primes outside the sequence, but possess such property (e.g., 109). - Vladimir Shevelev, Aug 14 2009
The Mathematica program FasterRamanujanPrimeList uses Laishram's result that a(n) < prime(3n).
See sequence A164952 for a generalization we call a Ramanujan k-prime. - Vladimir Shevelev, Sep 01 2009
From Jonathan Sondow, May 22 2010: (Start)
About 46% of primes < 19000 are Ramanujan primes. About 78% of the lesser of twin primes < 19000 are Ramanujan primes.
About 15% of primes < 19000 are the lesser of twin primes. About 26% of Ramanujan primes < 19000 are the lesser of twin primes.
A reason for the jumps is in Section 7 of "Ramanujan primes and Bertrand's postulate" and in Section 4 of "Ramanujan Primes: Bounds, Runs, Twins, and Gaps". (See the arXiv link for a corrected version of Table 1.)
See Shapiro 2008 for an exposition of Ramanujan's proof of his generalization of Bertrand's postulate. (End)
The (10^n)-th R prime: 2, 97, 1439, 19403, 242057, 2916539, 34072993, 389433437, .... - Robert G. Wilson v, May 07 2011, updated Aug 02 2012
The number of R primes < 10^n: 1, 10, 72, 559, 4459, 36960, 316066, 2760321, .... - Robert G. Wilson v, Aug 02 2012
a(n) = R_n = R_{0.5,n} in "Generalized Ramanujan Primes."
All Ramanujan primes are in A164368. - Vladimir Shevelev, Aug 30 2011
If n tends to infinity, then limsup(a(n) - A080359(n-1)) = infinity; conjecture: also limsup(a(n) - A080359(n)) = infinity (cf. A182366). - Vladimir Shevelev, Apr 27 2012
Or the largest prime x such that the number of primes in (x/2,x] equals n. This equivalent definition underlines an important analogy between Ramanujan and Labos primes (cf. A080359). - Vladimir Shevelev, Apr 29 2012
Research questions on R_n - prime(2n) are at A233739, and on n-Ramanujan primes at A225907. - Jonathan Sondow, Dec 16 2013
The questions on R_n - prime(2n) in A233739 have been answered by Christian Axler in "On generalized Ramanujan primes". - Jonathan Sondow, Feb 13 2014
Srinivasan's Lemma (2014): prime(k-n) < prime(k)/2 if R_n = prime(k) and n > 1. Proof: By the minimality of R_n, the interval (prime(k)/2,prime(k)] contains exactly n primes and so prime(k-n) < prime(k)/2. - Jonathan Sondow, May 10 2014
For some n and k, we see that A168421(k) = a(n) so as to form a chain of primes similar to a Cunningham chain. For example (and the first example), A168421(2) = 7, links a(2) = 11 = A168421(3), links a(3) = 17 = A168421(4), links a(4) = 29 = A168421(6), links a(6) = 47. Note that the links do not have to be of a form like q = 2*p+1 or q = 2*p-1. - John W. Nicholson, Feb 22 2015
Extending Sondow's 2010 comments: About 48% of primes < 10^9 are Ramanujan primes. About 76% of the lesser of twin primes < 10^9 are Ramanujan primes. - Dana Jacobsen, Sep 06 2015
Sondow, Nicholson, and Noe's 2011 conjecture that pi(R_{m*n}) <= m*pi(R_n) for m >= 1 and n >= N_m (see A190413, A190414) was proved for n > 10^300 by Shichun Yang and Alain Togbé in 2015. - Jonathan Sondow, Dec 01 2015
Berliner, Dean, Hook, Marr, Mbirika, and McBee (2016) prove in Theorem 18 that the graph K_{m,n} is prime for n >= R_{m-1}-m; see A291465. - Jonathan Sondow, May 21 2017
Okhotin (2012) uses Ramanujan primes to prove Lemma 8 in "Unambiguous finite automata over a unary alphabet." - Jonathan Sondow, May 30 2017
Sepulcre and Vidal (2016) apply Ramanujan primes in Remark 9 of "On the non-isolation of the real projections of the zeros of exponential polynomials." - Jonathan Sondow, May 30 2017
Axler and Leßmann (2017) compute the first k-Ramanujan prime for k >= 1 + epsilon; see A277718, A277719, A290394. - Jonathan Sondow, Jul 30 2017

			a(1) = 2 is Bertrand's postulate: pi(x) - pi(x/2) >= 1 for all x >= 2.
a(2) = 11 because a(2) < 8 log 8 < 17 and pi(n) - pi(n/2) > 1 for n = 16, 15, ..., 11 but pi(10) - pi(5) = 1.
Consider a(9)=71. Then the nearest prime > 71/2 is 37, and between a(9) and 2*37, that is, between 71 and 74, there exists a prime (73). - _Vladimir Shevelev_, Aug 14 2009 [corrected by _Jonathan Sondow_, Jun 17 2013]

Srinivasa Ramanujan, Collected Papers of Srinivasa Ramanujan (Ed. G. H. Hardy, S. Aiyar, P. Venkatesvara and B. M. Wilson), Amer. Math. Soc., Providence, 2000, pp. 208-209.
Harold N. Shapiro, Ramanujan's idea, Section 9.3B in Introduction to the Theory of Numbers, Dover, 2008.

T. D. Noe, Table of n, a(n) for n = 1..10000
N. Amersi, O. Beckwith, S. J. Miller, R. Ronan, and Jonathan Sondow, Generalized Ramanujan primes, arXiv:1108.0475 [math.NT], 2011.
N. Amersi, O. Beckwith, S. J. Miller, R. Ronan, and Jonathan Sondow, Generalized Ramanujan primes, Combinatorial and Additive Number Theory, Springer Proc. in Math. & Stat., CANT 2011 and 2012, Vol. 101 (2014), 1-13.
Christian Axler, Über die Primzahl-Zählfunktion, die n-te Primzahl und verallgemeinerte Ramanujan-Primzahlen, Ph.D. thesis 2013, in German, English summary.
Christian Axler, On generalized Ramanujan primes, arXiv:1401.7179 [math.NT], 2014.
Christian Axler, On generalized Ramanujan primes, Ramanujan J. 39 (1) (2016) 1-30.
Christian Axler and Thomas Leßmann, An explicit upper bound for the first k-Ramanujan prime, arXiv:1504.05485 [math.NT], 2015.
Christian Axler and Thomas Leßmann, On the first k-Ramanujan prime, Amer. Math. Monthly, 124 (2017), 642-646.
Adam H. Berliner, N. Dean, J. Hook, A. Marr, A. Mbirika, and C. McBee, Coprime and prime labelings of graphs, arXiv preprint arXiv:1604.07698 [math.CO], 2016; Journal of Integer Sequences, Vol. 19 (2016), #16.5.8.
Paul Erdős, A theorem of Sylvester and Schur, J. London Math. Soc., 9 (1934), 282-288.
Peter Hegarty, Why should one expect to find long runs of (non)-Ramanujan primes?, arXiv:1201.3847 [math.NT], 2012.
Ernest G. Hibbs, Component Interactions of the Prime Numbers, Ph. D. Thesis, Capitol Technology Univ. (2022), see p. 33.
Shanta Laishram, On a conjecture on Ramanujan primes, Int. J. Number Theory, 6 (2010), 1869-1873.
Catherine Lee, Minimum coprime graph labelings, arXiv:1907.12670 [math.CO], 2019.
Jaban Meher and M. Ram Murty, Ramanujan's proof of Bertrand's postulate, Amer. Math. Monthly, Vol. 120, No. 7 (2013), pp. 650-653.
Alexander Okhotin, Unambiguous finite automata over a unary alphabet, Inf. Comput., 212 (2012), 15-36.
Murat Baris Paksoy, Derived Ramanujan primes: R'_n, arXiv:1210.6991 [math.NT], 2012.
PlanetMath, Ramanujan prime
Srinivasa Ramanujan, A proof of Bertrand's postulate, J. Indian Math. Soc., 11 (1919), 181-182.
Juan Matias Sepulcre and Tomás Vidal, On the non-isolation of the real projections of the zeros of exponential polynomials, J. Math. Anal. Appl., 437 (2016) No. 1, 513-525.
Vladimir Shevelev, On critical small intervals containing primes, arXiv:0908.2319 [math.NT], 2009.
Vladimir Shevelev, Ramanujan and Labos primes, their generalizations and classifications of primes, arXiv:0909.0715 [math.NT], 2009-2011.
Vladimir Shevelev, Ramanujan and Labos primes, their generalizations, and classifications of primes, J. Integer Seq. 15 (2012) Article 12.5.4
Vladimir Shevelev, Charles R. Greathouse IV, and Peter J. C. Moses, On intervals (kn, (k+1)n) containing a prime for all n>1, Journal of Integer Sequences, Vol. 16 (2013), Article 13.7.3. arXiv:1212.2785
Jonathan Sondow, Ramanujan primes and Bertrand's postulate, arXiv:0907.5232 [math.NT], 2009-2010.
Jonathan Sondow, Ramanujan primes and Bertrand's postulate, Amer. Math. Monthly, 116 (2009), 630-635. Zentralblatt review.
Jonathan Sondow, J. W. Nicholson, and T. D. Noe, Ramanujan Primes: Bounds, Runs, Twins, and Gaps, arXiv:1105.2249 [math.NT] 2011; J. Integer Seq. 14 (2011) Article 11.6.2.
Jonathan Sondow, Ramanujan Prime, Eric Weisstein's MathWorld.
Jonathan Sondow and Eric Weisstein, Bertrand's Postulate, MathWorld.
Anitha Srinivasan, An upper bound for Ramanujan primes, Integers, 14 (2014), #A19.
Anitha Srinivasan and John W. Nicholson, An improved upper bound for Ramanujan primes, Integers, 15 (2015), #A52.
Wikipedia, Bertrand's postulate.
Wikipedia, Ramanujan prime.
Shichun Yang and Alain Togbé, On the estimates of the upper and lower bounds of Ramanujan primes, Ramanujan J., online 14 August 2015, 1-11.

Cf. A006992 (Bertrand primes), A056171 (pi(n) - pi(n/2)).
Cf. A000720, A007053, A014085, A060715, A084139, A084140, A143223, A143224, A143225, A143226, A143227, A080360, A080359, A164368, A164288, A164554, A164333, A164294, A164371, A190303.
Cf. A162996 (Round(kn * (log(kn)+1)), with k = 2.216 as an approximation of R_n = n-th Ramanujan Prime).
Cf. A163160 (Round(kn * (log(kn)+1)) - R_n, where k = 2.216 and R_n = n-th Ramanujan prime).
Cf. A178127 (Lesser of twin Ramanujan primes), A178128 (Lesser of twin primes if it is a Ramanujan prime).
Cf. A181671 (number of Ramanujan primes less than 10^n).
Cf. A174635 (non-Ramanujan primes), A174602, A174641 (runs of Ramanujan and non-Ramanujan primes).
Cf. A189993, A189994 (lengths of longest runs).
Cf. A190124 (constant of summation: 1/a(n)^2).
Cf. A192820 (2- or derived Ramanujan primes R'_n), A192821, A192822, A192823, A192824, A225907.
Cf. A193761 (0.25-Ramanujan primes), A193880 (0.75-Ramanujan primes).
Cf. A190413, A190414, A212493, A212541, A233739, A233822, A277718, A277719, A164952, A290394, A291465.
Cf. A185004 - A185007 ("modular" Ramanujan primes).
Not to be confused with the Ramanujan numbers or Ramanujan tau function, A000594.

A104272 := proc(n::integer)
    local R;
    if n = 1 then
        return 2;
    end if;
    R := ithprime(3*n-1) ; # upper limit Laishram's thrm Thrm 3 arXiv:1105.2249
    while true do
        if A056171(R) = n then # Defn. 1. of Shevelev JIS 14 (2012) 12.1.1
            return R ;
        end if;
        R := prevprime(R) ;
    end do:
end proc:
seq(A104272(n),n=1..200) ; # slow downstream search <= p(3n-1) R. J. Mathar, Sep 21 2017

(RamanujanPrimeList[n_] := With[{T=Table[{k,PrimePi[k]-PrimePi[k/2]}, {k,Ceiling[N[4*n*Log[4*n]]]}]}, Table[1+First[Last[Select[T,Last[ # ]==i-1&]]],{i,1,n}]]; RamanujanPrimeList[54]) (* Jonathan Sondow, Aug 15 2009 *)
(FasterRamanujanPrimeList[n_] := With[{T=Table[{k,PrimePi[k]-PrimePi[k/2]}, {k,Prime[3*n]}]}, Table[1+First[Last[Select[T,Last[ # ]==i-1&]]],{i,1,n}]]; FasterRamanujanPrimeList[54])
nn=1000; R=Table[0,{nn}]; s=0; Do[If[PrimeQ[k], s++]; If[PrimeQ[k/2], s--]; If[sT. D. Noe, Nov 15 2010 *)

ramanujan_prime_list(n) = {my(L=vector(n), s=0, k=1); for(k=1, prime(3*n)-1, if(isprime(k), s++); if(k%2==0 && isprime(k/2), s--); if(sSatish Bysany, Mar 02 2017

use ntheory ":all"; my $r = ramanujan_primes(1000); say "[@$r]"; # Dana Jacobsen, Sep 06 2015

a(n) = 1 + max{k: pi(k) - pi(k/2) = n - 1}.
a(n) = A080360(n-1) + 1 for n > 1.
a(n) >= A080359(n). - Vladimir Shevelev, Aug 20 2009
A193761(n) <= a(n) <= A193880(n).
a(n) = 2*A084140(n) - 1, for n > 1. - Jonathan Sondow, Dec 21 2012
a(n) = prime(2n) + A233739(n) = (A233822(n) + a(n+1))/2. - Jonathan Sondow, Dec 16 2013
a(n) = max{prime p: pi(p) - pi(p/2) = n} (see Shevelev 2012). - Jonathan Sondow, Mar 23 2016
a(n) = A000040(A179196(n)). - R. J. Mathar, Sep 21 2017
Sum_{n>=1} (-1)^(n+1)/a(n) = A190303. - Amiram Eldar, Nov 20 2020

Showing 1-10 of 43 results. Next

The primes of A080359 larger than 3 all have the property that the integers in the interval selected by halving the value of the preceding prime and halving their own value are all composite. This sequence here collects the primes that are not in A080359 but still share this property of the prime-free subinterval.

Conjecture: The sequence is unbounded.
Records are 0, 18, 36, 48, 52, ... with indices 1, 14, 20, 44, 96, ...
The places of nonzero terms correspond to places of those terms of A193507 which are in A164294. Conjecture: The asymptotic density of nonzero terms is 2/(e^2+1). A heuristic proof follows from the comment to A193507 and the first reference there.

Every greater of twin primes (A006512), beginning with 13, is in the sequence.
A very simple sieve for the generation of the terms is the following: Let p_n be the n-th prime. Consider consecutive intervals of the form (2p_n, 2p_{n+1}), n=1,2,... From every interval containing at least one prime we take the first one and remove it from the set of all primes. Then all remaining primes form the sequence. Let us demonstrate this sieve: For primes 2,3,5,7,11,... consider intervals (4,6), (6,10), (10,14), (14,22), (22,26), (26,34), ... . Removing from the set of all primes the first prime of each interval, i.e., 5,7,11,17,23,29,... ,we obtain 2,3,13,19,31, etc.
This sequence and A164368 are the mutually wrapping up sequences:
   a(1) <= A164368(1) <= a(2) <= A164368(2) <= ...
Following the steps to generate T(n,1) in A229608 provides an alternate method of generating this sequence. - Bob Selcoe, Oct 27 2015

cp's OEIS Frontend

A164294 Primes prime(k) such that all integers in [(prime(k-1)+1)/2,(prime(k)-1)/2] are composite, excluding those primes in A080359.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A194186 a(n) = A193507(n) - A080359(n).

Views

Author

Keywords

Comments

Links

Crossrefs

A194598 Union of A080359 and A164294.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Formula

A164554 Ramanujan primes A104272(n) for which A104272(n) = A080359(n).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

Extensions

A194217 a(n) = A104272(n)-A080359(n).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

A164372 Greater of twin primes (A006512) which are not Labos primes (A080359).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

Extensions

A182391 Numbers n for which A104272(n) = A080359(n).

Views

Author

Keywords

Comments

Crossrefs

Formula

A182392 Numbers n for which there exists only composite number k such that A060715(k) = n and 2*k-1 is prime, but A104272(n) differs from A080359(n).

Views

Author

Keywords

Comments

Crossrefs

A050376 "Fermi-Dirac primes": numbers of the form p^(2^k) where p is prime and k >= 0.

Views

Author

Keywords

Comments