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

A075714 1+n+n^s is a prime, s=18.

Original entry on oeis.org

1, 2, 9, 24, 27, 44, 80, 251, 263, 311, 332, 356, 366, 371, 458, 515, 546, 548, 561, 566, 597, 599, 608, 650, 674, 713, 717, 722, 746, 762, 855, 867, 909, 969, 989, 993, 1010, 1011, 1022, 1052, 1064, 1191, 1245, 1269, 1275, 1284, 1355, 1376, 1431, 1473
Offset: 1

Views

Author

Zak Seidov, Oct 03 2002

Keywords

Comments

For s = 5,8,11,14,17,20,..., n_s=1+n+n^s is always composite for any n>1. Also at n=1, n_s=3 is a prime for any s. So it is interesting to consider only the cases of s =/= 5,8,11,14,17,20,... and n>1. Here i consider the case s=18 and find several first n's making n_s a prime (or a probable prime).

Examples

			2 is OK because at s=18, n=2, n_s=1+n+n^s=262147 is a prime.
		

Crossrefs

Programs

  • Magma
    [n: n in [0..1600] | IsPrime(s) where s is 1+n+n^18]; // Vincenzo Librandi, Jul 28 2014
  • Mathematica
    Select[Range[1800], PrimeQ[1 + # + #^18] &]  (* Harvey P. Dale, Mar 24 2011 *)
  • PARI
    for(n=1,1000,if(isprime(1+n+n^18),print1(n",")))
    

Extensions

More terms from Ralf Stephan, Apr 05 2003

A075716 1+n+n^s is a prime, s=15.

Original entry on oeis.org

1, 2, 30, 32, 35, 54, 62, 77, 101, 120, 138, 161, 171, 186, 210, 234, 269, 285, 311, 341, 362, 368, 374, 467, 476, 486, 531, 567, 578, 720, 737, 740, 780, 806, 824, 932, 990, 1035, 1037, 1041, 1049, 1067, 1089, 1136, 1137, 1146, 1167, 1202, 1251, 1269
Offset: 1

Views

Author

Zak Seidov, Oct 03 2002

Keywords

Comments

For s = 5,8,11,14,17,20,..., n_s=1+n+n^s is always composite for any n>1. Also at n=1, n_s=3 is a prime for any s. So it is interesting to consider only the cases of s =/= 5,8,11,14,17,20,... and n>1. Here i consider the case s=15 and find several first n's making n_s a prime (or a probable prime).

Examples

			2 is OK because at s=15, n=2, n_s=1+n+n^s=32771 is a prime.
		

Crossrefs

Programs

  • Magma
    [n: n in [0..2000] | IsPrime(s) where s is 1+n+n^15]; // Vincenzo Librandi, Jul 28 2014
  • Mathematica
    Select[Range[1500], PrimeQ[1 + # + #^15] &] (* Harvey P. Dale, Dec 13 2010 *)
    Select[Range[2000], PrimeQ[Total[#^Range[1, 15, 14]] + 1] &] (* Vincenzo Librandi, Jul 28 2014 *)
  • PARI
    for(n=1,1000,if(isprime(1+n+n^15),print1(n",")))
    

Extensions

More terms from Ralf Stephan, Apr 05 2003

A245476 Least number k > 1 such that k^n + k + 1 is prime, or 0 if no such number exists.

Original entry on oeis.org

2, 2, 2, 2, 0, 2, 2, 0, 3, 3, 0, 2, 5, 0, 2, 2, 0, 2, 8, 0, 6, 3, 0, 6, 15, 0, 6, 2, 0, 2, 23, 0, 23, 56, 0, 15, 114, 0, 14, 11, 0, 3, 14, 0, 29, 110, 0, 21, 9, 0, 53, 59, 0, 6, 2, 0, 3, 29, 0, 71, 21, 0, 146, 17, 0, 35, 2, 0, 9, 6, 0, 77, 41, 0, 27, 176, 0, 153, 21, 0, 39, 32, 0, 2, 314, 0, 3, 5, 0, 66, 44, 0, 234
Offset: 1

Views

Author

Derek Orr, Jul 23 2014

Keywords

Comments

Except for a(2), a(n) = 0 if n == 2 mod 3 (A016789).
It appears that this is an "if and only if".
a(n) = 2 if and only if n is in A057732.
Many terms in the linked table correspond to probable primes. If n == 2 mod 3 then k^2+k+1 divides k^n+k+1. This is why a(n) = 0 if n > 2 and n == 2 mod 3. - Jens Kruse Andersen, Jul 28 2014

Examples

			2^9 + 2 + 1 = 515 is not prime. 3^9 + 3 + 1 = 19687 is prime. Thus a(9) = 3.
		

Crossrefs

Cf. Numbers n such that n^s + n + 1 is prime: A005097 (s = 1), A002384 (s = 2), A049407 (s = 3), A049408 (s = 4), A075723 (s = 6), A075722 (s = 7), A075720 (s = 9), A075719 (s = 10), A075718 (s = 12), A075717 (s = 13), A075716 (s = 15), A075715 (s = 16), A075714 (s = 18), A075713 (s = 19).

Programs

  • Maple
    f:= proc(n) local k;
       if n mod 3 = 2 and n > 2 then return 0 fi;
       for k from 2 to 10^6 do
          if isprime(k^n+k+1) then return k fi
       od:
      error("no solution found for n = %1",n);
    end proc:
    seq(f(n),n=1..100); # Robert Israel, Jul 27 2014
  • PARI
    a(n) = if(n>2&&n==Mod(2, 3), return(0)); k=2; while(!ispseudoprime(k^n+k+1), k++); k
    vector(150, n, a(n)) \\ Derek Orr with corrections and improvements from Colin Barker, Jul 23 2014
Showing 1-3 of 3 results.