|
|
A333368
|
|
Primes of the form k*m^(k*m) - 1 with m > 1.
|
|
1
|
|
|
3, 31, 191, 5119, 131071, 524287, 3758096383, 4353564671, 1356446145697, 1618481116086271, 2058911320946489, 1046695266054721074427023041, 847823165504324070285888664019, 5359447279004780799548150067050349330431, 2817103802133904744169307240538184064530443801964688726052818649087
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
LINKS
|
|
|
EXAMPLE
|
31 appears in this sequence because 31=2*2^(2*2)-1 and 31 is prime.
|
|
PROG
|
(PARI) lista(nn) = {my(k, n=2, v=List([]), x=4, y); while(x<nn, k=1; y=x-1; while(y<nn, if(ispseudoprime(y), listput(v, y)); k++; y=k*x^k-1); n++; x=n^n); Set(v); } \\ Jinyuan Wang, Mar 18 2020
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|