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

A088768 a(n) = (A087682(n)-1)/2.

Original entry on oeis.org

5, 7, 10, 19, 22, 25, 37, 40, 52, 79, 82, 85, 94, 109, 115, 124, 142, 169, 172, 187, 190, 220, 235, 247, 274, 277, 289, 292, 304, 319, 325, 334, 367, 382, 409, 415, 472, 487, 502, 520, 547, 550, 589, 604, 610, 649, 655, 715, 739, 745, 775, 787, 802, 814, 850
Offset: 1

Views

Author

Ray Chandler, Oct 26 2003

Keywords

Crossrefs

Programs

  • Magma
    [(n-1)/2: n in [6..2000] |IsPrime(n+8) and IsPrime(n-8)]; // Vincenzo Librandi, May 20 2017
  • Mathematica
    Select[Range[6, 2000], PrimeQ[2 # + 7] && PrimeQ[2 # - 9] &] - 1 (* Vincenzo Librandi, May 20 2017 *)

A082467 Least k>0 such that n-k and n+k are both primes.

Original entry on oeis.org

1, 2, 1, 4, 3, 2, 3, 6, 1, 6, 3, 2, 3, 6, 1, 12, 3, 2, 9, 6, 5, 6, 3, 4, 9, 12, 1, 12, 9, 4, 3, 6, 5, 6, 9, 2, 3, 12, 1, 24, 3, 2, 15, 6, 5, 12, 3, 8, 9, 6, 7, 12, 3, 4, 15, 12, 1, 18, 9, 4, 3, 6, 5, 6, 15, 2, 3, 12, 1, 6, 15, 4, 3, 6, 5, 18, 9, 2, 15, 24, 5, 12, 3, 14, 9, 18, 7, 12, 9, 4, 15, 6, 7, 30, 9
Offset: 4

Views

Author

Benoit Cloitre, Apr 27 2003

Keywords

Comments

The existence of k>0 for all n >= 4 is equivalent to the strong Goldbach Conjecture that every even number >= 8 is the sum of two distinct primes.
n and k are coprime, because otherwise n + k would be composite. So the rational sequence r(n) = a(n)/n = k/n is injective. - Jason Kimberley, Sep 21 2011
Because there are arbitrarily many composites from m!+2 to m!+m, there are also arbitrarily large a(n) but they increase very slowly. The twin prime conjecture implies that infinitely many a(n) are 1. - Juhani Heino, Apr 09 2020

Examples

			n=10: k=3 because 10-3 and 10+3 are both prime and 3 is the smallest k such that n +/- k are both prime.
		

Crossrefs

Cf. A129301 (records), A129302 (where records occur).
Cf. A047160 (allows k=0).
Cf. A078611 (subset for prime n).

Programs

  • Magma
    A082467 := func; [A082467(n):n in [4..98]]; // Jason Kimberley, Sep 03 2011
  • Maple
    A082467 := proc(n) local k; k := 1+irem(n,2);
    while n > k do if isprime(n-k) then if isprime(n+k)
    then RETURN(k) fi fi; k := k+2 od; print("Goldbach erred!") end:
    seq(A082467(i),i=4..90); # Peter Luschny, Sep 21 2011
  • Mathematica
    f[n_] := Block[{k}, If[OddQ[n], k = 2, k = 1]; While[ !PrimeQ[n - k] || !PrimeQ[n + k], k += 2]; k]; Table[ f[n], {n, 4, 98}] (* Robert G. Wilson v, Mar 28 2005 *)
  • PARI
    a(n)=if(n<0,0,k=1; while(isprime(n-k)*isprime(n+k) == 0,k++); k)
    

Formula

A078496(n)-a(n) = A078587(n)+a(n) = n.

Extensions

Entries checked by Klaus Brockhaus, Apr 08 2007

A087697 Numbers k such that k + 7 and k - 7 are both prime.

Original entry on oeis.org

10, 12, 24, 30, 36, 54, 60, 66, 90, 96, 120, 144, 156, 174, 186, 204, 234, 264, 270, 276, 300, 324, 360, 366, 390, 426, 450, 456, 516, 564, 570, 594, 600, 606, 624, 654, 666, 684, 726, 750, 780, 804, 816, 846, 870, 960, 984, 990
Offset: 1

Views

Author

Zak Seidov, Sep 27 2003

Keywords

Crossrefs

Programs

  • Magma
    [n: n in [5..1000] | IsPrime(n-7) and IsPrime(n+7)]; // Vincenzo Librandi, Jul 23 2018
  • Maple
    select(t -> isprime(t+7) and isprime(t-7), [seq(i,i=8..1000,2)]); # Robert Israel, Jul 22 2018
  • Mathematica
    Rest[Select[Range[1000], PrimeQ[# - 7] && PrimeQ[# + 7] &]] (* Vincenzo Librandi, Jul 23 2018 *)
  • PARI
    isok(n) = isprime(n-7) && isprime(n+7); \\ Michel Marcus, Jul 23 2018
    

A087683 Numbers n such that n + 10 and n - 10 are both prime.

Original entry on oeis.org

13, 21, 27, 33, 51, 57, 63, 69, 93, 99, 117, 141, 147, 183, 189, 201, 261, 267, 273, 303, 321, 327, 357, 363, 369, 399, 411, 429, 453, 477, 489, 513, 531, 567, 597, 603, 609, 651, 663, 729, 819, 849, 867, 873, 897, 957, 981, 987
Offset: 1

Views

Author

Zak Seidov, Sep 27 2003

Keywords

Comments

3 divides every term except the first. - T. D. Noe, May 14 2008

Crossrefs

Programs

  • Mathematica
    f[n_]:=PrimeQ[n-10]&&PrimeQ[n+10]; lst={};Do[If[f[n],AppendTo[lst,n]],{n,9,8!,2}];lst (* Vladimir Joseph Stephan Orlovsky, Oct 09 2009 *)
    Select[Range[10,1000],AllTrue[#+{10,-10},PrimeQ]&] (* Harvey P. Dale, Jul 29 2024 *)

A087711 a(n) = smallest number k such that both k-n and k+n are primes.

Original entry on oeis.org

2, 4, 5, 8, 7, 8, 11, 10, 11, 14, 13, 18, 17, 16, 17, 22, 21, 20, 23, 22, 23, 26, 25, 30, 29, 28, 33, 32, 31, 32, 37, 36, 35, 38, 37, 38, 43, 42, 41, 44, 43, 48, 47, 46, 57, 52, 51, 50, 53, 52, 53, 56, 55, 56, 59, 58, 75, 70, 69, 72, 67, 66, 65, 68, 67, 72, 71, 70, 71, 80, 81, 78
Offset: 0

Views

Author

Zak Seidov, Sep 28 2003

Keywords

Comments

Let b(n), c(n) and d(n) be respectively, smallest number m such that phi(m-n) + sigma(m+n) = 2n, smallest number m such that phi(m+n) + sigma(m-n) = 2n and smallest number m such that phi(m-n) + sigma(m+n) = phi(m+n) + sigma(m-n), we conjecture that for each positive integer n, a(n)=b(n)=c(n)=d(n). Namely we conjecture that for each positive integer n, a(n) < A244446(n), a(n) < A244447(n) and a(n) < A244448(n). - Jahangeer Kholdi and Farideh Firoozbakht, Sep 05 2014

Examples

			n=10: k=13 because 13-10 and 13+10 are both prime and 13 is the smallest k such that k +/- 10 are both prime
4-1=3, prime, 4+1=5, prime; 5-2=3, 5+2=7; 8-3=5, 8+3=11; 9-4=5, 9+4=13, ...
		

Crossrefs

Programs

  • Magma
    distance:=function(n); k:=n+2; while not IsPrime(k-n) or not IsPrime(k+n) do k:=k+1; end while; return k; end function; [ distance(n): n in [1..71] ]; /* Klaus Brockhaus, Apr 08 2007 */
    
  • Maple
    Primes:= select(isprime,{seq(2*i+1,i=1..10^3)}):
    a[0]:= 2:
    for n from 1 do
      Q:= Primes intersect map(t -> t-2*n,Primes);
      if nops(Q) = 0 then break fi;
      a[n]:= min(Q) + n;
    od:
    seq(a[i],i=0..n-1); # Robert Israel, Sep 08 2014
  • Mathematica
    s = ""; k = 0; For[i = 3, i < 22^2, If[PrimeQ[i - k] && PrimeQ[i + k], s = s <> ToString[i] <> ","; k++ ]; i++ ]; Print[s] (* Vladimir Joseph Stephan Orlovsky, Apr 03 2008 *)
    snk[n_]:=Module[{k=n+1},While[!PrimeQ[k+n]||!PrimeQ[k-n],k++];k]; Array[ snk,80,0] (* Harvey P. Dale, Dec 13 2020 *)
  • PARI
    a(n)=my(k);while(!isprime(k-n) || !isprime(k+n),k++);return(k) \\ Edward Jiang, Sep 05 2014

Formula

a(n) = A020483(n)+n for n >= 1. - Robert Israel, Sep 08 2014

Extensions

Entries checked by Klaus Brockhaus, Apr 08 2007
Showing 1-5 of 5 results.