《数论探秘》中复制黏贴过来的筛选大素数的代码

ysr

如下是从《数论探秘》中复制黏贴过来的筛选大素数的代码及其结果：
170141183460469231731687303715884105757~170141183460469231731687303715884105773之间的素数有2个(用时13.60156秒 )：
170141183460469231731687303715884105757   170141183460469231731687303715884105773   

（代码仅仅发主程序）
Private Sub Command1_Click()
Dim a, n
n = Trim(Text1)
n3 = n
m = Trim(Text2)
t = Timer
Do While MBJC(Trim(n), 3) < 0
If MBJC(Trim(n), 2) = 0 Then
m1 = m1 & n & "  "
s = s + 1
Else
m1 = m1
End If
n = n + 1
Loop
If Right(n, 1) Mod 2 = 0 Then
n = MPC1(Trim(n), 1)
Else
n = n
End If
Do While MBJC(Trim(n), Trim(m)) <= 0
If Len(n) < 11 Then
a1 = fenjieyinzi(Trim(n))
If InStr(a1, "*") = 0 Then
s = s + 1
If s Mod 10 = 0 Then
m1 = m1 & n & vbCrLf
Else
m1 = m1 & n & Space(3)
End If
Else
m1 = m1
End If
Else
n1 = MPC(Trim(n), 1)
a = 123
'a为明文
a1 = zzxc(Trim(n), Trim(a))
If Val(a1) > 1 Then
m1 = m1
Else
c = 999
'c为公约
Do While zzxc(Trim(n1), Trim(c)) > 1
c = Val(c - 1)
Loop
d = qniyuan(Trim(c), Trim(n1))
'd为逆元为私钥
a2 = qksmimo(Trim(a), Trim(c), Trim(n))
'a2为密文
a3 = qksmimo(Trim(a2), Trim(d), Trim(n))
If MBJC(Trim(a3), Trim(a)) = 0 Then
s = s + 1
If s Mod 5 = 0 Then
m1 = m1 & n & vbCrLf
Else
m1 = m1 & n & Space(3)
End If
Else
m1 = m1
End If
End If
End If
n = MPC1(Trim(n), 2)
Loop
Text3 = n3 & "~" & m & "之间的素数有" & s & "个(用时" & Timer - t & "秒 )：" & vbCrLf & m1
End Sub
Private Sub Command2_Click()
Text1 = ""
Text2 = ""
Text3 = ""
End Sub

ysr

73425633894622870812912264279427304662415407498217399761627450589750429468816261632146658967160832384997092708273473~73425633894622870812912264279427304662415407498217399761627450589750429468816261632146658967160832384997092708273673之间的素数有1个(用时1782.273秒 )：
73425633894622870812912264279427304662415407498217399761627450589750429468816261632146658967160832384997092708273663   
终于找到了一个素数（有116位）（筛选核判断利用的是欧拉原理）

ysr

我的方法是串行程序，不能象分布程序一样由多台电脑连网并行运算而减少时间，只能是单台计算，循环程序语句不多，但每一次循环结果就是下一次循环的输入，所以，不能并行运算，而循环次数多时间就长，能缩短时间的是程序采用了快速算法，大整数加减乘除的快速算法，以及快速幂摸算法等。判断原理是利用的欧拉原理。