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.

A370387 a(n) is the least prime p such that p + 6*k*(k+1) is prime for 0 <= k <= n-1 but not for k=n.

Original entry on oeis.org

2, 19, 5, 67, 7, 281, 1051, 6791, 11, 115599457, 365705201, 79352440891, 286351937491, 5810592517241, 17, 1942721697854617
Offset: 1

Views

Author

J.W.L. (Jan) Eerland, Mar 12 2024

Keywords

Comments

a(10), ..., a(14) > 10^7, a(15) = 17, a(16), ..., a(20) > 10^7.
a(29) = 31. - Chai Wah Wu, Apr 10 2024

Crossrefs

Programs

  • Maple
    f:= proc(p) local k;
      for k from 1 while isprime(p+k*(k+1)*6) do od:
      k
    end proc:
    A:= Vector(12): count:= 0:
    for i from 1 while count < 12 do
      v:= f(ithprime(i));
      if A[v] = 0 then count:= count+1; A[v]:= ithprime(i) fi
    od:
    convert(A,list);
  • Mathematica
    Table[p=1;m=6;Monitor[Parallelize[While[True,If[And[MemberQ[PrimeQ[Table[p+m*k*(k+1),{k,0,n-1}]],False]==False,PrimeQ[p+m*n*(n+1)]==False],Break[]];p++];p],p],{n,1,10}]
  • PARI
    isok(p, n) = for (k=0, n-1, if (! isprime(p + 6*k*(k+1)), return(0))); return (!isprime(p + 6*n*(n+1)));
    a(n) = my(p=2); while (!isok(p, n), p=nextprime(p+1)); p;
    
  • Perl
    use ntheory qw(:all); sub a { my $n = $[0]; my $lo = 2; my $hi = 2*$lo; while (1) { my @terms = grep { !is_prime($ + 6*$n*($n+1)) } sieve_prime_cluster($lo, $hi, map { 6*$*($+1) } 1 .. $n-1); return $terms[0] if @terms; $lo = $hi+1; $hi = 2*$lo; } }; $| = 1; for my $n (1..100) { print a($n), ", " } # Daniel Suteu, Dec 30 2024

Extensions

a(10)-a(11) from Chai Wah Wu, Apr 10 2024
a(12) from Chai Wah Wu, Apr 11 2024
a(13)-a(14) from David A. Corneth, Apr 11 2024
a(15) from J.W.L. (Jan) Eerland, Mar 12 2024
a(16) from Daniel Suteu, Dec 30 2024

A376675 a(n) is the least prime p such that p + 7*k*(k+1) is prime for 0 <= k <= n-1 but not for k=n.

Original entry on oeis.org

2, 3, 59, 5, 89, 599, 3329, 617, 269, 21107, 9833477, 19497833669, 215830859597, 111338387, 251704297005767, 17
Offset: 1

Views

Author

J.W.L. (Jan) Eerland, Oct 01 2024

Keywords

Crossrefs

Programs

  • Maple
    f:= proc(p) local k;
      for k from 1 while isprime(p+k*(k+1)*7) do od:
      k
    end proc:
    A:= Vector(12): count:= 0:
    for i from 1 while count < 12 do
      v:= f(ithprime(i));
      if A[v] = 0 then count:= count+1; A[v]:= ithprime(i) fi
    od:
    convert(A,list);
  • Mathematica
    Table[p=1;m=7;Monitor[Parallelize[While[True,If[And[MemberQ[PrimeQ[Table[p+m*k*(k+1),{k,0,n-1}]],False]==False,PrimeQ[p+m*n*(n+1)]==False],Break[]];p++];p],p],{n,1,10}]
  • PARI
    isok(p, n) = for (k=0, n-1, if (! isprime(p + 7*k*(k+1)), return(0))); return (!isprime(p + 7*n*(n+1)));
    a(n) = my(p=2); while (!isok(p, n), p=nextprime(p+1)); p;
    
  • Perl
    use ntheory qw(:all); sub a { my $n = $[0]; my $lo = 2; my $hi = 2*$lo; while (1) { my @terms = grep { !is_prime($ + 7*$n*($n+1)) } sieve_prime_cluster($lo, $hi, map { 7*$*($+1) } 1 .. $n-1); return $terms[0] if @terms; $lo = $hi+1; $hi = 2*$lo; } }; $| = 1; for my $n (1..100) { print a($n), ", " }; # Daniel Suteu, Oct 04 2024

Extensions

a(11)-a(12) from Hugo Pfoertner, Oct 01 2024
a(13)-a(14) from Hugo Pfoertner, Oct 03 2024
a(15)-a(16) from Daniel Suteu, Oct 04 2024

A378839 a(n) is the least prime p such that p + 8*k*(k+1) is prime for 0 <= k <= n-1 but not for k=n.

Original entry on oeis.org

2, 3, 151, 181, 13, 811, 23671, 92221, 45417481, 5078503, 4861, 20379346831, 12180447943, 31, 10347699089473
Offset: 1

Views

Author

J.W.L. (Jan) Eerland, Dec 09 2024

Keywords

Comments

No further terms < 2.5*10^11. - Michael S. Branicky, Dec 16 2024

Crossrefs

Programs

  • Maple
    f:= proc(p) local k;
      for k from 1 while isprime(p+k*(k+1)*8) do od:
      k
    end proc:
    A:= Vector(12): count:= 0:
    for i from 1 while count < 12 do
      v:= f(ithprime(i));
      if A[v] = 0 then count:= count+1; A[v]:= ithprime(i) fi
    od:
    convert(A,list);
  • Mathematica
    Table[p=1;m=8;Monitor[Parallelize[While[True,If[And[MemberQ[PrimeQ[Table[p+m*k*(k+1),{k,0,n-1}]],False]==False,PrimeQ[p+m*n*(n+1)]==False],Break[]];p++];p],p],{n,1,10}]
  • PARI
    isok(p, n) = for (k=0, n-1, if (! isprime(p + 8*k*(k+1)), return(0))); return (!isprime(p + 8*n*(n+1)));
    a(n) = my(p=2); while (!isok(p, n), p=nextprime(p+1)); p;
    
  • Perl
    use ntheory qw(:all); sub a { my $n = $[0]; my $lo = 2; my $hi = 2*$lo; while (1) { my @terms = grep { !is_prime($ + 8*$n*($n+1)) } sieve_prime_cluster($lo, $hi, map { 8*$*($+1) } 1 .. $n-1); return $terms[0] if @terms; $lo = $hi+1; $hi = 2*$lo; } }; $| = 1; for my $n (1..100) { print a($n), ", " }; #

Extensions

a(12)-a(14) from Michael S. Branicky, Dec 15 2024
a(15) from Daniel Suteu, Dec 17 2024

A378841 a(n) is the least prime p such that p + 9*k*(k+1) is prime for 0 <= k <= n-1 but not for k=n.

Original entry on oeis.org

2, 11, 13, 5, 19, 173, 3163, 83, 21013, 878359, 3676219, 239, 43, 5201390418463, 86927887467919
Offset: 1

Views

Author

J.W.L. (Jan) Eerland, Dec 09 2024

Keywords

Crossrefs

Programs

  • Maple
    f:= proc(p) local k;
      for k from 1 while isprime(p+k*(k+1)*9) do od:
      k
    end proc:
    A:= Vector(12): count:= 0:
    for i from 1 while count < 12 do
      v:= f(ithprime(i));
      if A[v] = 0 then count:= count+1; A[v]:= ithprime(i) fi
    od:
    convert(A,list);
  • Mathematica
    Table[p=1;m=9;Monitor[Parallelize[While[True,If[And[MemberQ[PrimeQ[Table[p+m*k*(k+1),{k,0,n-1}]],False]==False,PrimeQ[p+m*n*(n+1)]==False],Break[]];p++];p],p],{n,1,10}]
  • PARI
    isok(p, n) = for (k=0, n-1, if (! isprime(p + 9*k*(k+1)), return(0))); return (!isprime(p + 9*n*(n+1)));
    a(n) = my(p=2); while (!isok(p, n), p=nextprime(p+1)); p;
    
  • Perl
    use ntheory qw(:all); sub a { my $n = $[0]; my $lo = 2; my $hi = 2*$lo; while (1) { my @terms = grep { !is_prime($ + 9*$n*($n+1)) } sieve_prime_cluster($lo, $hi, map { 9*$*($+1) } 1 .. $n-1); return $terms[0] if @terms; $lo = $hi+1; $hi = 2*$lo; } }; $| = 1; for my $n (1..100) { print a($n), ", " }; #

Extensions

a(14) from Daniel Suteu, Dec 17 2024
a(15) from Daniel Suteu, Dec 22 2024
Showing 1-4 of 4 results.