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

A063983 Least k such that k*2^n +/- 1 are twin primes.

Original entry on oeis.org

4, 2, 1, 9, 12, 6, 3, 9, 57, 30, 15, 99, 165, 90, 45, 24, 12, 6, 3, 69, 132, 66, 33, 486, 243, 324, 162, 81, 90, 45, 345, 681, 585, 375, 267, 426, 213, 429, 288, 144, 72, 36, 18, 9, 147, 810, 405, 354, 177, 1854, 927, 1125, 1197, 666, 333, 519, 1032, 516, 258, 129, 72
Offset: 0

Views

Author

Robert G. Wilson v, Sep 06 2001

Keywords

Comments

Excluding the first three terms, all remaining terms have digital root 3, 6, or 9. - J. W. Helkenberg, Jul 24 2013

Examples

			a(3) = 9 because 9*2^3 = 72 and 71 and 73 are twin primes.
a(6) = 3 because 3*2^6 = 192 and {191, 193} are twin primes.
a(71) = 630 because 630*2^71 = 1487545442103938242314240 and {1487545442103938242314239, 1487545442103938242314241} are twin primes.
		

References

  • Richard Crandall and Carl Pomerance, 'Prime Numbers: A Computational Perspective,' Springer-Verlag, NY, 2001, page 12.

Crossrefs

Cf. A071256, A060210, A060256. For records see A125848, A125019.
Cf. A076806 (requires odd k).

Programs

  • Mathematica
    Table[Do[s=(2^j)*k; If[PrimeQ[s-1]&&PrimeQ[s+1],Print[{j,k}]], {k,1,2*j^2}],{j,0,100}]; (* outprint of a[j]=k *)
    Do[ k = 1; While[ ! PrimeQ[ k*2^n + 1 ] || ! PrimeQ[ k*2^n - 1 ], k++ ]; Print[ k ], {n, 0, 50} ]
    f[n_] := Block[{k = 1},While[Nand @@ PrimeQ[{-1, 1} + 2^n*k], k++ ];k];Table[f[n], {n, 0, 60}] (* Ray Chandler, Jan 09 2009 *)

Extensions

More terms from Labos Elemer, May 24 2002
Edited by N. J. A. Sloane, Jul 03 2008 at the suggestion of R. J. Mathar

A045753 Numbers n such that 4n-1 and 4n+1 are both primes.

Original entry on oeis.org

1, 3, 15, 18, 27, 45, 48, 57, 60, 78, 87, 105, 108, 150, 165, 207, 255, 258, 273, 288, 330, 357, 363, 372, 402, 405, 417, 447, 468, 483, 507, 522, 528, 567, 585, 648, 672, 678, 750, 780, 792, 813, 825, 840, 843, 867, 882, 885, 918, 942, 963, 1005, 1023
Offset: 1

Views

Author

Keywords

Examples

			3 belongs to the sequence because 4*3+1 and 4*3-1 are both primes.
		

Crossrefs

Programs

  • Magma
    [n: n in [1..2000] | IsPrime(4*n+1) and IsPrime(4*n-1)] // Vincenzo Librandi, Nov 18 2010
    
  • Mathematica
    Select[Range[1023], And @@ PrimeQ[{-1, 1} + 4# ] &] (* Ray Chandler, Dec 06 2006 *)
  • PARI
    list(lim)=my(v=List(),p=2); forprime(q=3,4*lim+1, if(q-p==2 && p%4==3, listput(v,q\4)); p=q); Vec(v) \\ Charles R Greathouse IV, Dec 03 2016

Extensions

More terms from Erich Friedman

A124519 Numbers k such that 12*k - 1 and 12*k + 1 are twin primes.

Original entry on oeis.org

1, 5, 6, 9, 15, 16, 19, 20, 26, 29, 35, 36, 50, 55, 69, 85, 86, 91, 96, 110, 119, 121, 124, 134, 135, 139, 149, 156, 161, 169, 174, 176, 189, 195, 216, 224, 226, 250, 260, 264, 271, 275, 280, 281, 289, 294, 295, 306, 314, 321, 335, 341, 344, 355, 356, 379, 399
Offset: 1

Views

Author

Artur Jasinski, Nov 04 2006

Keywords

Examples

			1 is in the sequence since 12*1 - 1 = 11 and 12*1 + 1 = 13 are twin primes.
		

Crossrefs

Programs

  • Mathematica
    Select[Range[400], And @@ PrimeQ[{-1, 1} + 12# ] &] (* Ray Chandler, Nov 16 2006 *)

Extensions

Extended by Ray Chandler, Nov 16 2006

A124522 a(n) = smallest k such that 2nk-1 and 2nk+1 are primes.

Original entry on oeis.org

2, 1, 1, 9, 3, 1, 3, 12, 1, 3, 9, 3, 12, 15, 1, 6, 3, 2, 6, 6, 1, 15, 3, 4, 3, 6, 2, 48, 6, 1, 21, 3, 3, 15, 6, 1, 27, 3, 4, 3, 15, 5, 12, 15, 2, 9, 3, 2, 9, 6, 1, 3, 60, 1, 6, 24, 2, 3, 9, 2, 129, 12, 7, 9, 15, 5, 12, 27, 1, 3, 9, 3, 42, 45, 1, 90, 3, 2, 66, 21, 5, 63, 27, 16, 6, 6, 2, 12, 24, 1, 6
Offset: 1

Views

Author

Artur Jasinski, Nov 04 2006

Keywords

Crossrefs

Programs

  • Maple
    isA001359 := proc(n) RETURN( isprime(n) and isprime(n+2)) ; end: A124522 := proc(n) local k; k :=1 ; while true do if isA001359(2*n*k-1) then RETURN(k) ; fi ; k := k+1 ; od ; end: for n from 1 to 60 do printf("%d,",A124522(n)) ; od ; # R. J. Mathar, Nov 06 2006
  • Mathematica
    f[n_] := Block[{k = 1},While[Nand @@ PrimeQ[{-1, 1} + 2n*k], k++ ];k];Table[f[n], {n, 91}] (* Ray Chandler, Nov 16 2006 *)
    skp[n_]:=Module[{k=1},While[AnyTrue[2n k+{1,-1},CompositeQ],k++];k]; Join[{2},Array[skp,100,2]] (* Harvey P. Dale, Mar 30 2024 *)
  • PARI
    {for(n=1,91,k=1;while(!isprime(2*n*k-1)||!isprime(2*n*k+1),k++);print1(k, ","))}

Extensions

Edited and extended by Klaus Brockhaus and R. J. Mathar, Nov 06 2006

A124518 Numbers k such that 10k-1 and 10k+1 are twin primes.

Original entry on oeis.org

3, 6, 15, 18, 24, 27, 42, 57, 60, 66, 81, 102, 105, 123, 129, 132, 162, 195, 213, 231, 234, 255, 273, 279, 297, 300, 312, 330, 333, 336, 339, 354, 393, 402, 405, 423, 426, 465, 480, 501, 510, 528, 552, 564, 585, 588, 609, 627, 630, 636, 645, 657, 666, 669, 678
Offset: 1

Views

Author

Artur Jasinski, Nov 04 2006

Keywords

Comments

All terms are divisible by 3. - Robert Israel, Apr 07 2019

Crossrefs

Programs

  • Maple
    select(t -> isprime(10*t+1) and isprime(10*t-1), [seq(i,i=3..1000,3)]); # Robert Israel, Apr 07 2019
  • Mathematica
    Select[Range[678], And @@ PrimeQ[{-1, 1} + 10# ] &] (* Ray Chandler, Nov 16 2006 *)

A124520 Numbers k such that 14*k - 1 and 14*k + 1 are twin primes.

Original entry on oeis.org

3, 30, 33, 63, 75, 78, 93, 102, 123, 138, 153, 162, 165, 192, 195, 240, 252, 273, 297, 303, 342, 387, 393, 420, 435, 438, 450, 468, 483, 522, 525, 540, 588, 630, 633, 660, 663, 717, 738, 747, 750, 765, 798, 825, 837, 855, 957, 978, 993, 1023, 1032, 1062
Offset: 1

Views

Author

Artur Jasinski, Nov 04 2006

Keywords

Examples

			3 is in the sequence since 14*3 - 1 = 41 and 14*3 + 1 = 43 are twin primes.
		

Crossrefs

Programs

  • Mathematica
    Select[Range[1100], And @@ PrimeQ[{-1, 1} + 14# ] &] (* Ray Chandler, Nov 16 2006 *)

Extensions

Extended by Ray Chandler, Nov 16 2006

A124521 Numbers k such that 16*k - 1 and 16*k + 1 are twin primes.

Original entry on oeis.org

12, 15, 27, 72, 93, 117, 132, 162, 168, 195, 198, 210, 258, 267, 300, 327, 330, 345, 435, 468, 642, 765, 813, 855, 903, 912, 960, 978, 993, 1128, 1143, 1182, 1290, 1350, 1353, 1365, 1392, 1398, 1440, 1632, 1680, 1713, 1737, 1797, 1848, 1860, 1947, 1953, 1962
Offset: 1

Views

Author

Artur Jasinski, Nov 04 2006

Keywords

Examples

			12 is in the sequence since 16*12 - 1 = 191 and 16*12 + 1 = 193 are twin primes.
		

Crossrefs

Programs

  • Maple
    A124521:=n->`if`(isprime(16*n-1) and isprime(16*n+1), n, NULL): seq(A124521(n), n=1..2000); # Wesley Ivan Hurt, Oct 10 2014
  • Mathematica
    Select[Range[2000], And @@ PrimeQ[{-1, 1} + 16# ] &] (* Ray Chandler, Nov 16 2006 *)

Extensions

Extended by Ray Chandler, Nov 16 2006

A123200 Numbers k such that 1000000*k-1 and 1000000*k+1 are twin primes.

Original entry on oeis.org

24, 30, 198, 345, 348, 432, 438, 471, 492, 609, 669, 774, 777, 858, 864, 1032, 1083, 1125, 1218, 1395, 1536, 1824, 1914, 1929, 2088, 2139, 2301, 2334, 2376, 2418, 2448, 2460, 2544, 2763, 2832, 2970, 3021, 3297, 3369, 3384, 3495, 3528, 3540, 3633, 3777
Offset: 1

Views

Author

Jonathan Vos Post, Nov 05 2006

Keywords

Examples

			a(6) = 432 because 431999999 and 432000001 are primes.
		

Crossrefs

Programs

  • Maple
    a:=proc(n) if isprime(10^6*n-1)=true and isprime(10^6*n+1)=true then n else fi end: seq(a(n),n=1..4500); # Emeric Deutsch, Nov 16 2006
  • Mathematica
    Select[Range[3800],AllTrue[#*10^6+{1,-1},PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, Nov 13 2017 *)

Extensions

More terms from Emeric Deutsch, Nov 16 2006
Showing 1-8 of 8 results.