判定梅森质数的卢卡斯序列

蔡家雄

设 30k+7 是素数，

且 (30k+7)^9*4+1 是素数，

则 10 是素数 (30k+7)^9*4+1 的原根。

有 30k+7=397, 457, 2557, 2917, 3547, ......

ysr

蔡家雄 发表于 2023-3-26 01:08
王兄：用你的快速幂运算程序，

a=10^9345586206421825384305800813114

余数是：18691172412843650768611601626228

ysr

蔡家雄 发表于 2023-3-26 01:25
王兄：用你的快速幂运算程序，

a=10^30581434818121773544530677133194

余数是61162869636243547089061354266388

ysr

蔡家雄 发表于 2023-3-26 01:35
王兄：用你的快速幂运算程序，

a=10^ 177738164913441570541606003148534

余数是：355476329826883141083212006297068

ysr

回复您的点评：王守恩：我这个(学习版)只能做十位数，Mod[10^12345678910,123456789]80219890，再多就显示“Overflow”  发表于 2023-3-26 01:54

内置求余数函数超过10位就会溢出，要调用自编的大数运算程序，我有快速幂模程序，幂的模可以先求模再求幂，幂的模等于模的幂（求幂的时候，奇数指数求一次模指数减小1，而指数是偶数的时候求一次模指数减小一半就是除以2，所以降幂速度很快，直到指数为0），这样速度就快了，中间数据小了一点。

蔡家雄

由 a=10^ 487138448432162610794

与 b= 974276896864325221589

则 a/b 的余数==

ysr

蔡家雄 发表于 2023-3-26 20:26
由 a=10^ 487138448432162610794

与 b= 974276896864325221589

由 a=10^ 487138448432162610794

与 b= 974276896864325221589

则 a/b 的余数=974276896864325221588

蔡家雄

由 a=10^ 154671165511282659934024706261946791821975634

与 b= 309342331022565319868049412523893583643951269

则 a/b 的余数==

ysr

蔡家雄 发表于 2023-3-27 21:49
由 a=10^ 154671165511282659934024706261946791821975634

与 b= 30934233102256531986804941252389358 ...

由 a=10^ 154671165511282659934024706261946791821975634

与 b= 309342331022565319868049412523893583643951269

则 a/b 的余数=309342331022565319868049412523893583643951268

蔡家雄

由 a=10^ 531989031935723044148800549586752242668697605054

与 b= 1063978063871446088297601099173504485337395210109

则 a/b 的余数==