login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A243472 Primes p such that p^6 - p^5 - 1 is prime. 2
2, 31, 101, 151, 181, 199, 229, 277, 307, 317, 379, 439, 479, 491, 647, 691, 797, 911, 997, 1039, 1051, 1181, 1291, 1367, 1381, 1471, 1511, 1549, 1657, 1709, 1847, 1867, 1987, 2081, 2099, 2111, 2207, 2467, 2621, 2707, 3041, 3221, 3259, 3541, 3571, 3581, 3769 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,1

LINKS

K. D. Bajpai, Table of n, a(n) for n = 1..4699

EXAMPLE

31 appears in the sequence because it is prime and 31^6 - 31^5 - 1 = 858874529 is also prime.

101 appears in the sequence because it is prime and 101^6 - 101^5  - 1 = 1051010050099 is also prime.

MAPLE

A243472 := proc() local a, b; a:=ithprime(n); b:= a^6-a^5-1; if isprime (b) then RETURN (a); fi; end: seq(A243472 (), n=1..2000);

MATHEMATICA

c = 0;  Do[k=Prime[n]; If[PrimeQ[k^6-k^5-1], c++; Print[c, " ", k]], {n, 1, 200000}];

Select[Prime[Range[600]], PrimeQ[#^6-#^5-1]&] (* Harvey P. Dale, Jan 21 2015 *)

PROG

(PARI) s=[]; forprime(p=2, 4000, if(isprime(p^6-p^5-1), s=concat(s, p))); s \\ Colin Barker, Jun 06 2014

CROSSREFS

Cf. A000040, A091567, A053182, A053183, A053184, A238136.

Sequence in context: A141978 A212569 A101254 * A107170 A273403 A291052

Adjacent sequences:  A243469 A243470 A243471 * A243473 A243474 A243475

KEYWORD

nonn

AUTHOR

K. D. Bajpai, Jun 05 2014

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 26 20:47 EDT 2019. Contains 321535 sequences. (Running on oeis4.)