数学中国

 找回密码
 注册
搜索
热搜: 活动 交友 discuz
楼主: ysr

[求助]素性测试

[复制链接]
 楼主| 发表于 2024-1-15 22:33 | 显示全部楼层
本帖最后由 ysr 于 2024-1-15 15:38 编辑

Private Function fenjiecx(sa As String, sb As String) As String '分解程序

Dim a, b, c
m = sa
p = sb
q = zhengchuqy(MCC1(Trim(m), Trim(p)))
m1 = Right(m, 2)
n1 = zhengchuqy(MCC(Trim(p), 2))
sd = 22
d = zhengchuqy(MCC(Trim(m), 4))
n1 = zhengliys2(Trim(n1), Val(sd))
d = zhengliys2(Trim(d), Val(sd))
If m1 Mod 4 = 1 Then
b = zhengchuqy(mcc2(MPC(zhengchuqy(mcc2(MPC(Trim(q), 1), Trim(p), Val(sd))), 1 & String(Val(sd), "0")), 2, Val(sd)))
Else
b = zhengchuqy(mcc2(MPC(zhengchuqy(mcc2(MPC(Trim(q), 3), Trim(p), Val(sd))), 1 & String(Val(sd), "0")), 2, Val(sd)))
End If
x = zhengliys2(1, Len(p) - 1)
bxn1 = mbc2(mbc2(Trim(b), Trim(x), Val(sd)), Trim(n1), Val(sd))
a = 4
b2 = mpc3(mpc2(mpc2(mbc2(8 & String(Val(sd), "0"), Trim(n1), Val(sd)), mbc2(mbc2(4 & String(Val(sd), "0"), Trim(b), Val(sd)), Trim(x), Val(sd))), mbc2(2 & String(Val(sd), "0"), Trim(x), Val(sd))), 2 & String(Val(sd), "0"))
If m1 Mod 4 = 1 Then
c = mpc2(mpc2(mpc3(mpc2(mpc2(mbc2(mbc2(4 & String(Val(sd), "0"), Trim(n1), Val(sd)), Trim(n1), Val(sd)), mbc2(8 & String(Val(sd), "0"), Trim(bxn1), Val(sd))), mbc2(mbc2(4 & String(Val(sd), "0"), Trim(n1), Val(sd)), Trim(x), Val(sd))), mbc2(2 & String(Val(sd), "0"), Trim(n1), Val(sd))), mbc2(mbc2(2 & String(Val(sd), "0"), Trim(b), Val(sd)), Trim(x), Val(sd))), mbc2(2 & String(Val(sd), "0"), Trim(x), Val(sd)))
d = mpc2(mpc2(mpc2(mpc2(mbc2(mbc2(2 & String(Val(sd), "0"), Trim(x), Val(sd)), Trim(d), Val(sd)), mbc2(mbc2(4 & String(Val(sd), "0"), Trim(bxn1), Val(sd)), Trim(n1), Val(sd))), mbc2(mbc2(2 & String(Val(sd), "0"), Trim(n1), Val(sd)), Trim(x), Val(sd))), mbc2(2 & String(Val(sd), "0"), Trim(bxn1), Val(sd))), mbc2(mbc2(2 & String(Val(sd), "0"), Trim(n1), Val(sd)), Trim(x), Val(sd)))
Else

c = mpc2(mpc2(mpc3(mpc2(mpc2(mbc2(mbc2(4 & String(Val(sd), "0"), Trim(n1), Val(sd)), Trim(n1), Val(sd)), mbc2(8 & String(Val(sd), "0"), Trim(bxn1), Val(sd))), mbc2(mbc2(4 & String(Val(sd), "0"), Trim(n1), Val(sd)), Trim(x), Val(sd))), mbc2(2 & String(Val(sd), "0"), Trim(n1), Val(sd))), mbc2(mbc2(2 & String(Val(sd), "0"), Trim(b), Val(sd)), Trim(x), Val(sd))), mbc2(4 & String(Val(sd), "0"), Trim(x), Val(sd)))
d = mpc2(mpc2(mpc2(mpc2(mbc2(mbc2(2 & String(Val(sd), "0"), Trim(x), Val(sd)), Trim(d), Val(sd)), mbc2(mbc2(4 & String(Val(sd), "0"), Trim(bxn1), Val(sd)), Trim(n1), Val(sd))), mbc2(mbc2(2 & String(Val(sd), "0"), Trim(n1), Val(sd)), Trim(x), Val(sd))), mbc2(2 & String(Val(sd), "0"), Trim(bxn1), Val(sd))), mbc2(mbc2(4 & String(Val(sd), "0"), Trim(n1), Val(sd)), Trim(x), Val(sd)))
End If
b2 = shuchujg(Trim(b2), Val(sd))
c = shuchujg(Trim(c), Val(sd))
d = shuchujg(Trim(d), Val(sd))


ff = jie3cifc(Trim(b), Trim(c), Trim(d), Trim(a), 22)
If InStr(ff, ".") = 0 Then
ff = ff
Else
ff = Mid(ff, 1, InStr(ff, ".") - 1)
End If
p = MPC3(Trim(p), Trim(ff))
fenjiecx = Trim(p)
End Function

Public Function jie3cifc(a2 As String, b2 As String, c2 As String, k2 As String, sd As String) As String '3次方程
a3 = mcc2(Trim(a2), Trim(k2), Val(sd))
b3 = mcc2(Trim(b2), Trim(k2), Val(sd))
c3 = mcc2(Trim(c2), Trim(k2), Val(sd))


   
   
   ' m = 36 * Val(ja) * Val(jb) - 8 * Val(ja) ^ 3 - 108 * Val(jc)
   If mbjc2(Trim(a3), 0) = 0 And mbjc2(Trim(b3), 0) = 0 Then
   m = mbc2(-108 & String(sd, "0"), Trim(c3), Val(sd))
   
   n4 = qdfh(mbc2(Trim(m), Trim(m), Val(sd)))
   Else
   
   m1 = mbc2(mbc2(Trim(a3), Trim(b3), Trim(sd)), Val(36) & String(sd, "0"), Val(sd))
   m2 = mbc2(mbc2(mbc2(Trim(a3), Trim(a3), Trim(sd)), Trim(a3), Val(sd)), Val(8) & String(sd, "0"), Val(sd))
   m3 = mbc2(Trim(c3), Val(108) & String(sd, "0"), Val(sd))
   m = mpc2(mpc2(Trim(m1), Trim(m2)), Trim(m3))
  Print "m="; m
  
'n1 = Val(m) ^ 2 + (12 * Val(B) - 4 * Val(A) ^ 2) ^ 3

  n1 = mbc2(Trim(m), Trim(m), Val(sd))
    n2 = mpc2(mbc2(Trim(b3), Val(12) & String(sd, "0"), Val(sd)), mbc2(mbc2(Trim(a3), Trim(a3), Val(sd)), Val(4) & String(sd, "0"), Val(sd)))
     n3 = mbc2(mbc2(Trim(n2), Trim(n2), Val(sd)), Trim(n2), Val(sd))
     n4 = mpc3(Trim(n1), Trim(n3))
     End If
     fn4 = fhys(Trim(n4))
    n5 = mbbc2(qqdl(qdfh(Trim(n4))), Val(sd))
   Print "n5="; n5
   
    'If n1 < 0 And m < 0 Then
    If Val(fn4) < 0 And Val(fhys(Trim(m))) < 0 Then
    'n=(93312R^6+311040R^5+285120R^4+221824R^3+27020R^2+3712R+320)^(1/2)
    'a=-(2+81k)
    'm=216r^3+360r^2+128r-8
    'k=(n^2/m^2)/27
    ''a1=(M/8*(9t-1))^(1/3),,b1=a1*(t)^(1/2)
    '令m=19683N^6+1215N^4M^2+17N^2M^4-M^6,
'm1=243N^4+3N^2M^2+M^4,
'n=(m^2-m1^3)^(1/2),
'则:t=(((m+n)^(1/3)+(m-n)^(1/3))+(2*M^2+3*N^2))/(3*M^2),
''其中M≠m,N≠n,
   jq = mbbc3(mbbc2(qqdl(MPC1(mbc2(qdfh(Trim(m)), qdfh(Trim(m)), Val(sd)), mbc2(qdfh(Trim(n5)), qdfh(Trim(n5)), Val(sd)))), Val(sd)), Val(sd))
   jq3 = mcc2(Trim(n5), Trim(m), Val(sd))
   jq1 = jsfanzq(mcc2(Trim(n5), Trim(m), Val(sd)), Val(sd))
   Print "jq3="; jq3
   Print "jq1="; jq1
   q1 = mbc2(Trim(jq), jsyuxian(mcc2(Trim(jq1), 3 & String(sd, "0"), Val(sd)), Val(sd)), Val(sd))
   
   Print "q1="; q1
   
   If Val(fhys(Trim(q1))) > 0 Then
  q = tjfh(Trim(q1), Val(-1))
  Else
  q = qdfh(Trim(q1))
  End If
   jq3 = mbbc2(MPC(1 & String(sd, "0"), mbc2(jsyuxian(mcc2(Trim(jq1), 3 & String(sd, "0"), Val(sd)), Val(sd)), jsyuxian(mcc2(Trim(jq1), 3 & String(sd, "0"), Val(sd)), Val(sd)), _
   Val(sd))), Val(sd))
   jq2 = tjfh(Trim(jq3), Val(-1))
   Print "jq2="; jq2
   
   qa = mbc2(Trim(jq), Trim(jq2), Val(sd))
   p = q
   pa = qa
  Else
  'If n1 < 0 And m > 0 Then
If Val(fn4) < 0 And MBJC(qdfh(Trim(m)), 0) > 0 Then

jq = mbbc3(mbbc2(qqdl(MPC1(mbc2(qdfh(Trim(m)), qdfh(Trim(m)), Val(sd)), mbc2(qdfh(Trim(n5)), qdfh(Trim(n5)), Val(sd)))), Val(sd)), Val(sd))
   jq1 = jsfanzq(mcc2(Trim(n5), Trim(m), Val(sd)), Val(sd))
   jq3 = mcc2(Trim(n5), Trim(m), Val(sd))
   q = mbc2(Trim(jq), jsyuxian(mcc2(Trim(jq1), 3 & String(sd, "0"), Val(sd)), Val(sd)), Val(sd))
   jq2 = mbbc2(MPC(1 & String(sd, "0"), mbc2(jsyuxian(mcc2(Trim(jq1), 3 & String(sd, "0"), Val(sd)), Val(sd)), jsyuxian(mcc2(Trim(jq1), 3 & String(sd, "0"), Val(sd)), Val(sd)), _
   Val(sd))), Val(sd))
   Print "jq3="; jq3
   Print "jq1="; jq1
   Print "q"; q
   Print "jq2="; jq2
   qa = mbc2(Trim(jq), Trim(jq2), Val(sd))
   p = q
   pa = qa
Else
'If Val(m) + Val(n2) < 0 Then
If fhys(mpc3(Trim(m), Trim(n5))) = -1 Then
p = "-" & mbbc3(qdfh(mpc3(Trim(m), Trim(n5))), Val(sd))
Else
p = mbbc3(mpc3(Trim(m), Trim(n5)), Val(sd))

End If

'If Val(m) < Val(n2) Then
If mbjc2(Trim(m), Trim(n5)) = -1 Then
'q = -(Val(n2) - Val(m)) ^ (1 / 3)
q = "-" & mbbc3(qdfh(mpc2(Trim(n5), Trim(m))), Val(sd))
Else
q = mbbc3(mpc2(Trim(m), Trim(n5)), Val(sd))

End If



  End If
  
  
  End If
  
  'If a = 0 And b = 0 Then
  If Val(fn4) > 0 And mbjc2(Trim(n4), 0) <> 0 Then  'Val(n4) > 0 Or mbjc2(Trim(n4), 0) = 0
  tx1 = mcc2(mpc2(mpc3(Trim(p), Trim(q)), mbc2(Trim(a3), 2 & String(sd, "0"), Val(sd))), 6 & String(sd, "0"), Val(sd))
  ftx1 = fhys(Trim(tx1))
  tx4 = mcc2(mpc3(mpc3(Trim(p), Trim(q)), mbc2(Trim(a3), 4 & String(sd, "0"), Val(sd))), 12 & String(sd, "0"), Val(sd))
  If Val(fhys(Trim(tx4))) > 0 Then
  tx2 = tjfh(Trim(tx4), Val(-1))
  Else
  tx2 = qdfh(Trim(tx4))
  End If
  ' -(Val(p) + Val(q) + 4 * Val(A)) / 12 & "+" & Abs((3 ^ (1 / 2) * (Val(p) - Val(q)))) / 12 & "i数据太大有误差"
' -(Val(p) + Val(q) + 4 * Val(A)) / 12 & "-" & Abs((3 ^ (1 / 2) * (Val(p) - Val(q)))) / 12 & "i裳〈笫莞呔鹊&#65533;"
tx3 = zhengchuqy(MCC1(qdfh(mbc2(mpc2(Trim(p), Trim(q)), zhengchuqy(MBBC(3 & String(2 * Val(sd), "0"))), Val(sd))), 12))
jie3cifc = shuchujg(Trim(tx1), Val(sd))

  Else
  If Mid(fn4, 1, 1) = "-" And Trim(m) = 0 Or Trim(n5) = 0 Then
tx1 = mcc2(mpc2(mpc3(Trim(p), Trim(q)), mbc2(Trim(a3), 2 & String(sd, "0"), Val(sd))), 6 & String(sd, "0"), Val(sd))
  ftx1 = fhys(Trim(tx1))
  tx4 = mcc2(mpc3(mpc3(Trim(p), Trim(q)), mbc2(Trim(a3), 4 & String(sd, "0"), Val(sd))), 12 & String(sd, "0"), Val(sd))
  ' -(Val(p) + Val(q) + 4 * Val(A)) / 12 & "+" & Abs((3 ^ (1 / 2) * (Val(p) - Val(q)))) / 12 & "i数据太大有误差"
' -(Val(p) + Val(q) + 4 * Val(A)) / 12 & "-" & Abs((3 ^ (1 / 2) * (Val(p) - Val(q)))) / 12 & "i可选大数据高精度的"
tx3 = mcc2(mbc2(mpc2(Trim(p), Trim(q)), zhengchuqy(MBBC(3 & String(2 * Val(sd), "0"))), Val(sd)), 12 & String(sd, "0"), Val(sd))
  If Val(fhys(Trim(tx4))) > 0 Then
  tx2 = tjfh(Trim(tx4), Val(-1))
  Else
  tx2 = qdfh(Trim(tx4))
  End If
  tx1 = shuchujg(Trim(tx1), Val(sd))
tx12 = shuchujg(mpc3(Trim(tx2), Trim(tx3)), Val(sd))
  tx13 = shuchujg(mpc2(Trim(tx2), Trim(tx3)), Val(sd))
  tx14 = tx1 & "/" & tx12 & "/" & tx13
  jie3cifc = tx1
  
  Else
  
  'd = mcc2(mpc2(mpc3(Trim(p), Trim(q)), mbc2(Trim(a3), 2 & String(sd, "0"), Val(sd))), 6 & String(sd, "0"), Val(sd))
  'd = mcc2(MPC(qdfh(mbc2(Trim(a3), 2 & String(sd, "0"), Val(sd))), mpc3(Trim(p), Trim(q))), 6 & String(sd, "0"), Val(sd))
  'mpc2减法器有问题?
  'fd = fhys(Trim(d))
  tx1 = mcc2(mpc2(mpc3(Trim(p), Trim(q)), mbc2(Trim(a3), 2 & String(sd, "0"), Val(sd))), 6 & String(sd, "0"), Val(sd))
  ftx1 = fhys(Trim(tx1))
  tx4 = mcc2(mpc3(mpc3(Trim(p), Trim(q)), mbc2(Trim(a3), 4 & String(sd, "0"), Val(sd))), 12 & String(sd, "0"), Val(sd))
  ' -(Val(p) + Val(q) + 4 * Val(A)) / 12 & "+" & Abs((3 ^ (1 / 2) * (Val(p) - Val(q)))) / 12 & "i数据太大有误差"
' -(Val(p) + Val(q) + 4 * Val(A)) / 12 & "-" & Abs((3 ^ (1 / 2) * (Val(p) - Val(q)))) / 12 & "i可选大数据高精度的"
tx3 = mcc2(mbc2(mpc3(Trim(pa), Trim(qa)), zhengchuqy(MBBC(3 & String(2 * Val(sd), "0"))), Val(sd)), 12 & String(sd, "0"), Val(sd))
If Val(fhys(Trim(tx4))) > 0 Then
  tx2 = tjfh(Trim(tx4), Val(-1))
  Else
  tx2 = qdfh(Trim(tx4))
  End If
  tx1 = shuchujg(Trim(tx1), Val(sd))
tx12 = shuchujg(mpc3(Trim(tx2), Trim(tx3)), Val(sd))
  tx13 = shuchujg(mpc2(Trim(tx2), Trim(tx3)), Val(sd))
  tx14 = tx1 & "/" & tx12 & "/" & tx13
  jie3cifc = tx1

  End If
  
End If


  


  

End Function

Private Function max(ByVal lists As String) As String
Dim temp As String
Dim a() As String
a = Split(lists, "/")
Dim b As Long
temp = a(0)
For b = 0 To UBound(a)
If Abs(temp) < Abs(a(b)) Then temp = a(b)
Next
max = temp
End Function
回复 支持 反对

使用道具 举报

 楼主| 发表于 2024-1-15 23:04 | 显示全部楼层
Private Function qdqiandao0(sa As String) As String
Do While Left(sa, 1) = "0"
sa = Mid(sa, 2)
Loop
If Len(sa) = 0 Then
qdqiandao0 = 0
Else
qdqiandao0 = sa
End If
End Function

Private Function jsmssushu(sa As String) As String
Dim B
B = 1
Do While c <= sa - 1
B = MbC(Trim(B), 2)
c = c + 1
Loop
jsmssushu = MPC(Trim(B), 1)
End Function

Private Function fenjieyinzi0(A As String) As String
X = 3
Do While InStr(MCC(Trim(A), 2), "/") = 0 And MBJC(Trim(A), 4) >= 0

  
If B = 0 Then
fenjieyinzi0 = fenjieyinzi0 & "2"
B = 1
Else
fenjieyinzi0 = fenjieyinzi0 & "*2"
End If
A = MCC(Trim(A), 2)
k = A
  
Loop

Do While MBJC(Trim(A), "0001") > 0
Do While MBJC(Trim(X), zhengchuqy(MBBC(Trim(A)))) <= 0
Do While InStr(MCC1(Trim(A), Trim(X)), "/") = 0 And MBJC(Trim(A), MbC(Trim(X), Trim(X))) >= 0
  
If B = 0 Then
fenjieyinzi0 = fenjieyinzi0 & X
B = 1
Else
fenjieyinzi0 = fenjieyinzi0 & "*" & X
End If
A = MCC1(Trim(A), Trim(X))
Loop
  
X = MPC1(Trim(X), 2)
Loop
  
k = A
A = 1
Loop
  
If B = 1 Then
fenjieyinzi0 = fenjieyinzi0 & "*" & k
Else
fenjieyinzi0 = "这是一个质数"
End If
End Function



Private Function fenjieyinzi(sa As String) As String
Dim X, A, B, k As String
A = Val(sa)

X = 3
If A <= 1 Or A > Int(A) Then
If A = 1 Then
fenjieyinzi = "它既不是质数,也不是合数"

Else
MsgBox "error"
End If
  
Else

Do While A / 2 = Int(A / 2) And A >= 4
  
If B = 0 Then
fenjieyinzi = fenjieyinzi & "2"
B = 1
Else
fenjieyinzi = fenjieyinzi & "*2"
End If
A = A / 2
k = A
  
Loop

Do While A > 1
Do While X <= Sqr(A)
Do While A / X = Int(A / X) And A >= X * X
  
If B = 0 Then
fenjieyinzi = fenjieyinzi & X
B = 1
Else
fenjieyinzi = fenjieyinzi & "*" & X
End If
A = A / X
Loop
  
X = X + 2
Loop
  
k = A
A = 1
Loop
  
If B = 1 Then
fenjieyinzi = fenjieyinzi & "*" & k
Else
fenjieyinzi = "这是一个质数"
End If
  
  
  
  

End If

End Function

Private Function zhengchuqyushu(sa As String) As String
If InStr(sa, "/") = 0 Then
zhengchuqyushu = 0
Else
zhengchuqyushu = Mid(sa, InStr(sa, "/") + 1)
End If


End Function
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|手机版|小黑屋|数学中国 ( 京ICP备05040119号 )

GMT+8, 2024-4-25 16:05 , Processed in 0.056640 second(s), 14 queries .

Powered by Discuz! X3.4

Copyright © 2001-2020, Tencent Cloud.

快速回复 返回顶部 返回列表