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

A196668 The Chebyshev primes of index 2.

Original entry on oeis.org

17, 29, 41, 53, 61, 71, 83, 89, 101, 103, 113, 127, 137, 149, 151, 157, 193, 211, 239, 241, 251, 257, 269, 293, 307, 311, 313, 317, 331, 353, 359, 373, 379, 389, 397, 401, 433, 439, 443, 457, 461, 463, 479, 499, 503, 509, 521, 523, 569, 571, 577, 587, 593, 599
Offset: 1

Views

Author

Michel Planat, Oct 05 2011

Keywords

Comments

The sequence consists of such odd prime numbers p that satisfy li[psi(p^2)]-li[psi(p^2-1)]<1/2, where li(x) is the logarithmic integral and psi(x) is the Chebyshev's psi function.

Links

Crossrefs

Cf. A196667, A196669, A196670.

Programs

  • Magma
    Mangoldt:=function(n);
if #Factorization(n) eq 1 then return Log(Factorization(n)[1][1]); else return 0; end if;
end function;
tcheb:=function(n);
x:=0;
for i in [1..n] do
x:=x+Mangoldt(i);
end for;
return(x);
end function;
jump2:=function(n);
x:=LogIntegral(tcheb(NthPrime(n)^2))-LogIntegral(tcheb(NthPrime(n)^2-1));
return x;
end function;
Set2:=[];
for i in [2..1000] do
if jump2(i)-1/2 lt 0 then Set2:=Append(Set2,NthPrime(i)); NthPrime(i); end if;
end for;
Set2;
  • Maple
    # The function PlanatSole(n,r) is in A196667.
A196668 := n -> PlanatSole(n,2); # Peter Luschny, Oct 23 2011
  • Mathematica
    ChebyshevPsi[n_] := Log[LCM @@ Range[n]];
Reap[Do[If[LogIntegral[ChebyshevPsi[p^2]] - LogIntegral[ChebyshevPsi[p^2 - 1]] < 1/2, Print[p]; Sow[p]], {p, Prime[Range[2, 200]]}]][[2, 1]] (* Jean-François Alcover, Jul 14 2018, updated Dec 06 2018 *)
  • Perl
    use ntheory ":all"; forprimes { say if 2 * (LogarithmicIntegral(chebyshev_psi($**2)) - LogarithmicIntegral(chebyshev_psi($**2-1))) < 1 } 3, 1000; # Dana Jacobsen, Dec 29 2015
  • Sage
    def A196668(n) : return PlanatSole(n,2)
# The function PlanatSole(n,r) is in A196667.
# Peter Luschny, Oct 23 2011

A196670 The Chebyshev primes of index 4.

Original entry on oeis.org

5, 7, 17, 19, 31, 37, 41, 43, 53, 59, 67, 73, 79, 83, 101, 103, 107, 127, 131, 149, 157, 163, 179, 181, 197, 199, 211, 223, 227, 257, 269, 277, 281, 317, 331, 337, 347, 353, 379, 389, 419, 421, 439, 461, 463, 467, 479, 491, 499, 509, 541, 563, 569, 577, 617
Offset: 1

Views

Author

Michel Planat, Oct 05 2011

Keywords

Comments

The sequence consists of such odd prime numbers p that satisfy li(psi(p^4)) - li(psi(p^4-1)) < 1/4, where li(x) is the logarithmic integral and psi(x) is the Chebyshev psi function.

Links

Crossrefs

Cf. A196667, A196668, A196669.

Programs

  • Magma
    Mangoldt:=function(n);
if #Factorization(n) eq 1 then return Log(Factorization(n)[1][1]); else return 0; end if;
end function;
tcheb:=function(n);
x:=0;
for i in [1..n] do
x:=x+Mangoldt(i);
end for;
return(x);
end function;
jump4:=function(n);
x:=LogIntegral(tcheb(NthPrime(n)^4))-LogIntegral(tcheb(NthPrime(n)^4-1));
return x;
end function;
Set4:=[];
for i in [2..1000] do
if jump4(i)-1/4 lt 0 then Set4:=Append(Set4,NthPrime(i)); NthPrime(i); end if;
end for;
Set4;
  • Maple
    # The function PlanatSole(n,r) is in A196667.
A196670 := n -> PlanatSole(n,4); # Peter Luschny, Oct 23 2011
  • Mathematica
    ChebyshevPsi[n_] := Log[LCM @@ Range[n]];
Reap[Do[If[LogIntegral[ChebyshevPsi[p^4]] - LogIntegral[ChebyshevPsi[p^4 - 1]] < 1/4, Print[p]; Sow[p]], {p, Prime[Range[2, 120]]}]][[2, 1]] (* Jean-François Alcover, Jul 14 2018, updated Dec 06 2018 *)
  • Perl
    use ntheory ":all"; forprimes { say if 4 *(LogarithmicIntegral(chebyshev_psi($**4)) - LogarithmicIntegral(chebyshev_psi($**4-1))) < 1 } 3,100; # Dana Jacobsen, Dec 29 2015
  • Sage
    def A196670(n) : return PlanatSole(n,4)
# The function PlanatSole(n,r) is in A196667.
# Peter Luschny, Oct 23 2011

Extensions

More terms from Dana Jacobsen, Dec 29 2015

A197187 The Riemann primes of the psi type and index 3.

Original entry on oeis.org

2, 3, 5, 7, 11, 13, 17, 29, 59, 67, 97, 103, 149, 151, 233, 251, 277, 311, 313, 479, 643, 719, 919, 967, 1039, 1373, 1489, 1553, 1847, 1973, 1979, 2663, 2953, 3323, 3677, 3691, 4651, 4663, 4789
Offset: 1

Views

Author

Michel Planat, Oct 11 2011

Keywords

Comments

The sequence consists of the prime numbers p that are champions (left to right maxima) of the function |psi(p^3)-p^3|, where psi(p) is the Chebyshev psi function.

Links

Crossrefs

Cf. A196669, A197185, A197186, A197188.

Programs

  • Mathematica
    ChebyshevPsi[n_] := Range[n] // MangoldtLambda // Total;
Reap[For[max=0; p=2, p < 1000, p = NextPrime[p], f = Abs[ChebyshevPsi[p^3] - p^3]; If[f > max, max = f; Print[p]; Sow[p]]]][[2, 1]] (* Jean-François Alcover, Dec 03 2018 *)
  • Perl
    use ntheory ":all"; my($max,$f)=(0); forprimes { $f=abs(chebyshev_psi($**3)-$**3); if ($f > $max) { say; $max=$f; } } 1000; # Dana Jacobsen, Dec 28 2015

Extensions

More terms from Dana Jacobsen, Dec 28 2015

A197299 The Riemann primes of the theta type and index 3.

Original entry on oeis.org

2, 3, 5, 7, 13, 17, 23, 31, 37, 41, 43, 47, 53, 59, 67, 73, 83, 89, 101, 103, 137, 163, 167, 179, 197, 211, 223, 239, 251, 277, 331, 379, 397, 431, 463, 467, 521, 577, 593, 601, 613, 617, 719, 809, 881, 919, 967, 1091, 1123, 1129, 1237, 1249, 1289
Offset: 1

Views

Author

Michel Planat, Oct 13 2011

Keywords

References

  • The sequence consists of the prime numbers p that are champions (left to right maxima) of the function |theta(p^3)-p^3|, where theta(p) is the Chebyshev theta function.

Links

Crossrefs

Cf. A196669, A197187, A197297, A197298, A197300.

Programs

  • Perl
    use ntheory ":all"; my($max,$f)=(0); forprimes { $f=abs(chebyshev_theta($**3)-$**3); if ($f > $max) { say; $max=$f; } } 1000; # Dana Jacobsen, Dec 28 2015

Extensions

More terms from Dana Jacobsen, Dec 28 2015
Showing 1-4 of 4 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.