clc,clear,close all
for n = 3:200
    s = 4;
    num = 2^n-1;
    for i = 1:n-2
        s = mod(s^2-2,num);
        if s == 0
             disp(n)
             break
        end
    end
end

---

相关阅读资料
*-卢卡斯-莱默检验法

---

2021年11月5日20:55:12