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

A205617 Number of decompositions of 2n into an unordered sum of two non-Ramanujan primes (A174635).

Original entry on oeis.org

0, 0, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 2, 2, 1, 2, 2, 2, 1, 3, 0, 1, 3, 1, 2, 2, 1, 3, 2, 0, 1, 3, 2, 1, 3, 1, 3, 4, 1, 2, 4, 1, 4, 2, 0, 3, 2, 3, 2, 3, 2, 3, 5, 1, 3, 4, 0, 5, 1, 0, 4, 3, 3, 1, 4, 3, 5, 4, 0, 4, 3, 1, 4, 2, 2, 6, 2, 3, 4, 4, 1, 3
Offset: 1

Views

Author

Keywords

Comments

There are 15 zeros in the first 10^8 terms. a(n) > 0 for n from 315 to 10^8.

Examples

			a(25) = 3. 2*25 = 50 = 7+43 = 13+37 = 19+31 (7, 13, 19, 31, 37, and 43 are all non-Ramanujan primes (A174635)). 50 is the unordered sum of two non-Ramanujan primes in three ways.
		

Crossrefs

A205616 Even numbers that are not the sum of two non-Ramanujan primes (A174635).

Original entry on oeis.org

2, 4, 52, 70, 100, 124, 130, 148, 208, 232, 238, 292, 352, 418, 628
Offset: 1

Views

Author

Keywords

Comments

No other terms < 2*10^8. Conjectured to be complete.

Examples

			70 is a term because no two non-Ramanujan primes (A174635) sum to 70. 72 is not a term because 19 + 53 = 72. 19 and 53 are both non-Ramanujan primes.
		

Crossrefs

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

Views

Author

Jonathan Sondow, Feb 27 2005

Keywords

Comments

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

Examples

			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]
		

References

  • 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.

Crossrefs

Cf. A006992 (Bertrand primes), A056171 (pi(n) - pi(n/2)).
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. A185004 - A185007 ("modular" Ramanujan primes).
Not to be confused with the Ramanujan numbers or Ramanujan tau function, A000594.

Programs

  • Maple
    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
  • Mathematica
    (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 *)
  • PARI
    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
  • Perl
    use ntheory ":all"; my $r = ramanujan_primes(1000); say "[@$r]"; # Dana Jacobsen, Sep 06 2015
    

Formula

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

A174641 Smallest prime that begins a run of n consecutive primes that are not Ramanujan primes.

Original entry on oeis.org

3, 3, 3, 73, 191, 191, 509, 2539, 2539, 5279, 9901, 9901, 9901, 11593, 11593, 55343, 55343, 55343, 55343, 55343, 174929, 174929, 174929, 225977, 225977, 225977, 225977, 225977, 534889, 534889, 534889, 534889, 534889, 534889, 534889, 534889, 2492317, 2492317
Offset: 1

Views

Author

T. D. Noe, Nov 29 2010

Keywords

Comments

The run of 10 consecutive non-Ramanujan primes was mentioned by Sondow.

Crossrefs

Cf. A104272 (Ramanujan primes), A174635 (non-Ramanujan primes).
Cf. A174602 (runs of Ramanujan primes).

Programs

  • Mathematica
    nn=10000; t=Table[0, {nn}]; len=Prime[3*nn]; s=0; Do[If[PrimeQ[k], s++]; If[PrimeQ[k/2], s--]; If[s
    				
  • Perl
    use ntheory ":all";
    my($k, $max, $start, $end, $inc, $p, $q, $r, $pi)
       = (0, 0, 0, 10, 1e9, 0, 2, [], prime_iterator(3));
    while (1) {
      if (!@$r) {
        ($start, $end) = ($end+1, $end+$inc);
        $r = ramanujan_primes($start, $end);
      }
      ($p, $q, $k) = ($q, shift(@$r), 0);
      # $k = prime_count($p+1,$q-1);
      $k++ while $pi->() < $q;
      say ++$max," ",next_prime($p)   while $k > $max;
    }
    # Dana Jacobsen, Jul 14 2016

A205618 Last occurrence of n partitions in A205617.

Original entry on oeis.org

314, 629, 959, 1154, 1424, 4619, 1922, 4094, 2549, 3884, 3989, 5774, 4724, 5669, 6404, 5879, 7664, 5594, 8609, 9239, 9029, 8714, 10562, 10394, 9869, 11549, 9764, 12239, 11444, 11969, 11654, 14279, 14489, 12209, 13229, 15014, 13859, 14804, 15584, 16979, 19634
Offset: 0

Views

Author

Keywords

Examples

			a(0) = 314 because the last occurrence of a zero in A205617 is at a(314).
		

Crossrefs

A233866 Ramanujan primes A104272 that are primes p(k) such that (k+1)*p(k)>k*p(k+1).

Original entry on oeis.org

2, 11, 17, 29, 41, 59, 71, 101, 107, 127, 149, 179, 227, 229, 281, 307, 311, 347, 349, 419, 431, 439, 461, 487, 569, 599, 641, 643, 659, 739, 769, 809, 821, 823, 853, 857, 937, 967, 983, 1009, 1019, 1031, 1049, 1061, 1087, 1091, 1151, 1187, 1217, 1229
Offset: 1

Views

Author

Irina Gerasimova, Dec 17 2013

Keywords

Comments

Non-Ramanujan primes (A174635) that are primes p(k) such that (k+1)*p(k)>k*p(k+1): 137, 163, 191, 197, 223, 277, 379, 397, 457, 499, 521, 613, 617, 673, 757, 859, 877, 907, 1093, 1181, 1213, 1223, 1231,...
Primes p(k) such that (k+1)*p(k)> k*p(k+1): 2, 11, 17, 29, 41, 59, 71,...
primes p(k) such that (k+1)*p(k) - k*p(k+1)=1: 2, 11,...
Primes p(m) such that (m+1)*p(m) < m*p(m+1): 3, 5, 7, 13, 19, 23, 31,...

Examples

			Ramanujan prime 2 is in this sequence because 2 = p(1) such that (1 + 1)*p(1) = 2*2 = 4 > 1*p(1 + 1) = 3;
Ramanujan prime 11 is in this sequence because 11 = p(5) such that (5 + 1)*p(5) = 6*11 = 66 > 5*p(5 + 1) = 65.
		

A267549 Primes prime(k) such that floor( (prime(k)/k)^2 ) <= prime(k+1) - prime(k).

Original entry on oeis.org

3, 5, 7, 13, 23, 113
Offset: 1

Views

Author

John W. Nicholson, Jan 16 2016

Keywords

Comments

Prime index A000720 is: 2, 3, 4, 6, 9, and 30.
floor( (prime(k)/k)^2 ) is: 2, 2, 3, 4, 6, and 14.
Similarly, ceiling( (prime(k)/k)^2 ) > prime(k+1) - prime(k) holds for all prime(k) < 10^8 with the exception of prime(k) = 7. For prime(k) = 7, 4 = ceiling((prime(k) / k)^2) = prime(k+1) - prime(k).
Stronger than Firoozbakht's conjecture which states that the sequence prime(k)^(1/k) is strictly decreasing.
Conjecture: list is complete. If so, subsequence of A124147 and A174635.
Andrew Granville conjectures that lim sup (prime(n+1)-prime(n))/log(prime(n))^2 >= 2/e^gamma = 1.1229189.... If so (or at least if the lim sup is greater than 1) then this sequence is infinite. - Charles R Greathouse IV, Feb 18 2016

Examples

			For a(3) = 7, floor((7 / 4)^2) = 3 < 4 = 11 - 7. Note that all other a(n) use = instead of <.
		

Crossrefs

Programs

  • Mathematica
    Select[Prime@ Range[10^5], Floor[(#/PrimePi@ #)^2] <= NextPrime@ # - # &] (* Michael De Vlieger, Jan 21 2016 *)
  • PARI
    L=10^11;p=2;forprime(q=3,L,a=floor((p/primepi(p))^2.);if(a<=q-p, print1(p, ", "));p=q)

Formula

Floor((A000040(k) / k)^2) <= A000040(k+1)-A000040(k), where k = A000720.
Floor(A001248(k) / A000290(k)) <= A001223(k), where k = A000720.
Showing 1-7 of 7 results.