|
%I #26 May 26 2020 22:00:55
%S 3,31,191,5119,131071,524287,3758096383,4353564671,1356446145697,
%T 1618481116086271,2058911320946489,1046695266054721074427023041,
%U 847823165504324070285888664019,5359447279004780799548150067050349330431,2817103802133904744169307240538184064530443801964688726052818649087
%N Primes of the form k*m^(k*m) - 1 with m > 1.
%H Sean A. Irvine, <a href="https://github.com/archmageirvine/joeis/blob/master/src/irvine/oeis/a333/A333368.java">Java program</a> (github)
%e 31 appears in this sequence because 31=2*2^(2*2)-1 and 31 is prime.
%o (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
%Y Cf. A000040, A050918.
%K nonn,easy
%O 1,1
%A _Eder Vanzei_, Mar 17 2020
%E Corrected by _Michel Marcus_ and _Jinyuan Wang_, Mar 17 2020
|
