fermattest.vbs

可以进行任意多位的大整数的乘方、开方、四则运算等

'----------------------------------------------------------
' Probabilistic Primality Test Script
'------------------------------------
' This VBScript code checks the primality of strHex for 4 
' times using Fermat Primality Test.
' 
' strHex = "10 9f e4 57 14 86 6e 56 " + _
'          "fd d4 ad 9b 6b 68 6d f2 " + _
'          "72 24 af b7 86 8c f4 f0 " + _
'          "cb b7 94 52 69 32 85 3c " + _
'          "bf 0b ee a6 15 94 16 66 " + _
'          "54 d1 3c d9 fe 0d 9d a5 " + _
'          "94 a9 7e e2 02 30 f1 2f " + _
'          "b5 43 4d e7 3f b4 f8 10 " + _
'          "27 25 a0 16 22 b3 1b 1e " + _
'          "a4 2e 3a 26 50 19 03 9a " + _
'          "c1 df 31 86 9b d9 79 30 " + _
'          "d7 92 fb 72 cd aa 97 1d " + _
'          "8a 80 15 af"
' 
' is a known 800-bit prime.
'
'----------------------------------------------------------

Dim strHex, nPrime

strHex = "10 9f e4 57 14 86 6e 56 " + _
         "fd d4 ad 9b 6b 68 6d f2 " + _
         "72 24 af b7 86 8c f4 f0 " + _
         "cb b7 94 52 69 32 85 3c " + _
         "bf 0b ee a6 15 94 16 66 " + _
         "54 d1 3c d9 fe 0d 9d a5 " + _
         "94 a9 7e e2 02 30 f1 2f " + _
         "b5 43 4d e7 3f b4 f8 10 " + _
         "27 25 a0 16 22 b3 1b 1e " + _
         "a4 2e 3a 26 50 19 03 9a " + _
         "c1 df 31 86 9b d9 79 30 " + _
         "d7 92 fb 72 cd aa 97 1d " + _
         "8a 80 15 af"

nPrime = Hpmbmath.FermatTest(strHex, 4)
If nPrime = 1 Then
    Hpmbmath.Output("strHex is a prime.")
Else 
  If nPrime = 0 Then
    Hpmbmath.Output("strHex is not a prime.")
  Else 
    Hpmbmath.Output("internal error.")
  End If
End If