isprime.asv

来自「用MATLAB实现的RSA。界面友好」· ASV 代码 · 共 21 行

%判断num是不是素数，如果是，返回1，否则返回0
function prime=IsPrime(num)
k=0;n=num-1;
while mod(n,2)==0
    k=k+1;
    n=n/2;
end
m=n;
a=unidrnd(n);
b=mod(a^m, n);
if b==1
    prime=1;
end
for i=0:1:k-1
    if mod(b,n)==n-1
        prime=1;
    else
        b=b*b;
    end
end
prime=0;