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

A271718 Numbers n such that n*(n+1)^n - 1 is prime.

Original entry on oeis.org

2, 3, 7, 14, 43, 81, 943, 1621
Offset: 1

Views

Author

Jeppe Stig Nielsen, Apr 12 2016

Keywords

Comments

The corresponding primes are a subset of the generalized Woodall primes (A210340).

Examples

			14 is a member because 14*15^14 - 1 = 408700964355468749 is a prime number.
		

Crossrefs

Programs

  • Mathematica
    Select[Range[10^3], PrimeQ[# (# + 1)^# - 1] &] (* Michael De Vlieger, Apr 12 2016 *)
  • PARI
    for(n=1,10^10,ispseudoprime(n*(n+1)^n-1)&&print1(n,", "))

A353122 Numbers k such that k^k*(k+1) + 1 is prime.

Original entry on oeis.org

0, 1, 2, 3, 6, 9, 186, 198, 8390
Offset: 1

Views

Author

Juri-Stepan Gerasimov, Apr 24 2022

Keywords

Comments

Corresponding primes start 2, 3, 13, 109, 326593, 3874204891, ...
a(9) > 6000. - Jon E. Schoenfield, Jun 05 2022
a(10) > 18000. - Michael S. Branicky, Aug 08 2024

Examples

			9 is in the sequence because 9^9*(9+1) + 1 = 3874204891, which is prime.
		

Crossrefs

Programs

  • Magma
    [n: n in [0..200] | IsPrime(n^n*(n+1) + 1)];
    
  • Mathematica
    Join[{0}, Select[Range[200], PrimeQ[#^#*(# + 1) + 1] &]] (* Amiram Eldar, Apr 25 2022 *)
  • PARI
    isok(k) = ispseudoprime(k^k*(k+1) + 1); \\ Michel Marcus, May 16 2022

Extensions

a(9) from Michael S. Branicky, Dec 22 2023
Showing 1-2 of 2 results.