[求助]素性测试

ysr

48159457367+24078843941=72238301308.

ysr

39097039560272928010164175021985880602580272792922567487344554734374127201273854611544066056097115357697238578058106247273461503061095065912451866010325996513800613924935742199093690953691140303632422462719709060227310529936492064349560839847109726836021275336252370239606709890960860239645355109048787
经过程序判断，结果是：*这是合数，有302位，用时3401.309秒

ysr

2571548369866401154762124488623114006103344329106851940612135327634004653071260863682696945178204681400011946477794407710938848804017401062757414714439486114198518710341540080693877842223985143550250435065198963213728593736084588180107843186915782390310479016413374955446436206039程序判断结果：
*这是合数，有280位，用时3952.102秒

ysr

'修正程序，细化公式的一种程序代码：
Private Sub Command1_Click()
Dim m, n
m = Trim(Text1)
n = Trim(Text2)
b16 = "8325594683879298"
q = zhengchuqy(MCC1(Trim(m), Trim(n)))
q = q & String(10, "0")
B1 = zhengchuqy(MCC1(MPC(Trim(q), 1 & String(10, "0")), MPC(Trim(n), 1)))
b2 = zhengchuqy(MCC1(MPC(Trim(B1), 1 & String(10, "0")), 2))
b3 = MbC(Trim(b2), Trim(b16))
b4 = Mid(b3, 1, Len(b3) - 16)
b5 = MPC1(1 & String(10, "0"), MbC(Trim(b4), 2))
m1 = m & String(10, "0")
b10 = zhengchuqy(MCC1(Trim(m1), Trim(b5)))
n2 = zhengchuqy(MBBC(Trim(b10)))

Text3 = n2




End Sub


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

ysr

输入：6958000001674999998647，
72238301309，
输出：1.33336599364964879897，
288960283210.64579957836343063979，
15655500003533973748207.12575346059581765915，
1739500000376610400733.60772629647099185986，
-4538679201004402756564633339781693092072378.77668261584207350467

输入：6958000001674999998647，
62238301309，
输出：1.79625902391500316968，
335388331074.70409419257292914772，
15655500003544866736242.15216383711423240783，
1739500000375241396745.07758291697614104653，
-3369068949487167329425824452753062874092234.72899176895173032750.

ysr

输入:  a=1.33336599364964879897,  b=288960283210.64579957836343063979,  c=15655500003533973748207.12575346059581765915,  d=1739500000376610400733.60772629647099185986，e=-4538679201004402756564633339781693092072378.77668261584207350467;  输出结果:    x1,2=-105258918697.7571014210+ -18059575327.6119410465,  x3,4=-3098533264.7956239397+ -18059575327.0906165853


输入:  a=1.79625902391500316968,  b=335388331074.70409419257292914772,  c=15655500003544866736242.15216383711423240783,  d=1739500000375241396745.07758291697614104653，e=-3369068949487167329425824452753062874092234.72899176895173032750;  输出结果:    x1,2=-90687850885.7647971919+ -15559575327.8433433063,  x3,4=-2669601076.7353675467+ -15559575327.1126397862

ysr

-3098533264.7956239397   +   18059575327.0906165853 =14,961,042,062.2949926456，
14961042062.2949926456   -   72238301309 =-57,277,259,246.7050073544

57277259246.7050073544   +   72238301309 =129,515,560,555.7050073544，
129515560555.7050073544   ÷   2 =64,757,780,277.8525036772，
而这个值比较接近实际。

或者：
-105258918697.757101421   -   18059575327.6119410465 =-123,318,494,025.3690424675，
-123318494025.3690424675   +   72238301309 =-51,080,192,716.3690424675


又或者：

-2669601076.7353675467   +   15559575327.1126397862 =12,889,974,250.3772722395,
12889974250.3772722395   +   62238301309 =75,128,275,559.3772722395.

75128275559.3772722395   +   62238301309 =137,366,576,868.3772722395，
137366576868.3772722395   ÷   2 =68,683,288,434.18863611975，
而这个值更接近实际。

ysr

Private Sub Command1_Click()
Dim a, b, c, d, f, g
Dim ja, jb, jc
k = DeleteSpace(Text1.Text)
ja = DeleteSpace(Text2.Text)
jb = DeleteSpace(Text3.Text)
jc = DeleteSpace(Text4.Text)
jd = DeleteSpace(Text5.Text)
sd1 = DeleteSpace(Text14.Text)
If Abs(Val(sd1)) <= 10 Then
sd = 10 + 2
Else
If Len(sd1) > 3 And Abs(Val(sd1)) > 300 Then
sd = 300 + 2
Else
sd = Int(Abs(Val(sd1))) + 2
End If
End If

a2 = zhengliys2(Trim(ja), Val(sd))
b2 = zhengliys2(Trim(jb), Val(sd))

c2 = zhengliys2(Trim(jc), Val(sd))
D2 = zhengliys2(Trim(jd), Val(sd))
k2 = zhengliys2(Trim(k), Val(sd))
If MBJC(Trim(k2), 0) = 0 Then
Text4.Text = "a 不能为 0"
   Text5.Text = "a 为 0可能已不是1元4次方程"
  Else
   
a3 = mcc2(Trim(a2), Trim(k2), Val(sd))
b3 = mcc2(Trim(b2), Trim(k2), Val(sd))
c3 = mcc2(Trim(c2), Trim(k2), Val(sd))
d3 = mcc2(Trim(D2), Trim(k2), Val(sd))
jk = 8 & String(sd, "0")
ja1 = mbc2("-4" & String(sd, "0"), Trim(b3), Val(sd))
jb1 = mpc2(mbc2(2 & String(sd, "0"), mbc2(Trim(a3), Trim(c3), Val(sd)), Val(sd)), mbc2(8 & String(sd, "0"), Trim(d3), Val(sd)))
jc1 = mpc2(mbc2(Trim(d3), mpc2(mbc2(4 & String(sd, "0"), Trim(b3), Val(sd)), mbc2(Trim(a3), Trim(a3), Val(sd))), Val(sd)), mbc2(Trim(c3), Trim(c3), Val(sd)))
y = jie3cifc(Trim(ja1), Trim(jb1), Trim(jc1), Trim(jk), Val(sd))
Y1 = zhengliys2(Trim(y), Val(sd))
End If
Text6 = zhengliys2(Trim(y), Val(sd))
Text8 = y
za = 1 & String(sd, "0")
z2 = mpc2(mpc3(mbc2(8 & String(sd, "0"), Trim(Y1), Val(sd)), mbc2(Trim(a3), Trim(a3), Val(sd))), mbc2(4 & String(sd, "0"), Trim(b3), Val(sd)))
z2 = mbbc2(qdfh(Trim(z2)), Val(sd))

zb1 = mcc2(mpc3(Trim(a3), Trim(z2)), 2 & String(sd, "0"), Val(sd))
zb2 = mcc2(mpc2(Trim(a3), Trim(z2)), 2 & String(sd, "0"), Val(sd))
zc1 = mpc3(Trim(Y1), mcc2(mpc2(mbc2(Trim(a3), Trim(Y1), Val(sd)), Trim(c3)), Trim(z2), Val(sd)))
zc2 = mpc2(Trim(Y1), mcc2(mpc2(mbc2(Trim(a3), Trim(Y1), Val(sd)), Trim(c3)), Trim(z2), Val(sd)))
y3 = jie2cifc(Trim(za), Trim(zb1), Trim(zc1), Val(sd))
y4 = jie2cifc(Trim(za), Trim(zb2), Trim(zc2), Val(sd))
Text9 = y3
Text10 = y4

Text15.Text = Text15.Text & "  输入" & ":  " & "a=" & k & ",  b=" & ja & ",  c=" & jb & ",  d=" & jc _
  & "，e=" & jd & ";  输出结果" & ":    " & "x1,2=" & Text9.Text & ",  x3,4=" & Text10.Text
End Sub


Private Sub Command2_Click()
Text1 = ""
Text2 = ""
Text3 = ""
Text4 = ""
Text5 = ""

Text6 = ""
Text7 = ""
Text8 = ""
Text9 = ""

Text10 = ""
Text11 = ""
Text12 = ""
Text13 = ""
Text14 = ""
Text15 = ""

End Sub

ysr

输入：6958000001674999998647，
68683288434，
输出：68749242188
68937973266
69247974510
69574961778
70036451573
70148794420
70190094151
70413727759
70910113965
71049692356
71108853442

ysr

输入：6958000001674999998647，
70910113965，
输出：71000008221
非常接近于实际

进一步改进程序得到：
输入：6958000001674999998647，
70910113965，
输出：70999998350
非常接近于实际