mersenne.m

来自「关于质数的一些程序」· M 代码 · 共 7 行

function y=mersenne(n)
m=0;
for k=1:n
    if(isprime(2^k-1)) m=m+1;
        y(m,1)=k;y(m,2)=2^k-1;
    end        
end