login
Primality of Mersenne number 2^prime(n)-1
0

%I #11 Jul 10 2021 17:00:00

%S 1,1,1,1,0,1,1,1,0,0,1,0,0,0,0,0,0,1,0,0,0,0,0,1,0,0,0,1,0,0,1,0,0,0,

%T 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,

%U 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,1

%N Primality of Mersenne number 2^prime(n)-1

%C Characteristic function of A016027 = indices of prime Mersenne numbers (A001348). See these sequences for further references.

%H Richard K. Guy, <a href="http://www.jstor.org/stable/2322249">The Strong Law of Small Numbers</a>, Example 2.

%H <a href="/index/Ch#char_fns">Index entries for characteristic functions</a>

%t Table[If[PrimeQ[2^n-1],1,0],{n,Prime[Range[120]]}] (* or *) Module[ {mpe = MersennePrimeExponent[Range[15]]},Table[If[MemberQ[mpe,p],1,0],{p,Prime[ Range[120]]}]] (* Requires Mathematica version 10 or later *) (* _Harvey P. Dale_, Jul 10 2021 *)

%o (PARI) for(n=1,199,print1(ispseudoprime(2^prime(n)-1)","))

%K nonn

%O 1

%A _M. F. Hasler_, Mar 23 2014