2^1277=有385位，用时0秒2601983048666099770481310081841021384653815561816676201329778087600902014918340074503059860433081046210605403488570251947845891562080866227034976651419330190731032377347305086443295837415395887618239855136922452802923419286887119716740625346109565072933087221327790207134604146257063901166556207972729700461767055550785130256674608872183239507219512717434046725178680177638925792182272
2601983048666099770481310081841021384653815561816676201329778087600902014918340074503059860433081046210605403488570251947845891562080866227034976651419330190731032377347305086443295837415395887618239855136922452802923419286887119716740625346109565072933087221327790207134604146257063901166556207972729700461767055550785130256674608872183239507219512717434046725178680177638925792182271~2601983048666099770481310081841021384653815561816676201329778087600902014918340074503059860433081046210605403488570251947845891562080866227034976651419330190731032377347305086443295837415395887618239855136922452802923419286887119716740625346109565072933087221327790207134604146257063901166556207972729700461767055550785130256674608872183239507219512717434046725178680177638925792182271之间的素数有个：（用时3178.844秒）
这个数是合数！但目前无法分解
本帖最后由 ysr 于 2024-1-5 11:27 编辑
输入m=14585113 p=975输出a=4b=2800.04989035974c=-120338.165417003 d=17207854.860784(这是一元三次方程的系数)
本帖最后由 ysr 于 2024-1-5 11:47 编辑
Private Sub Command1_Click()
Dim a, b, c
m = Val(Text1)
p = Val(Text2)
q = Val(m \ p)
If m Mod 4 = 1 Then
b = ((q - 1) / p - 1) / 2
Else
b = ((q - 3) / p - 1) / 2
End If
b1 = b - 1
d = m \ 4
n1 = p \ 2
n2 = (d / (2 * b1 + 1)) ^ (1 / 2)
x = n2 - n1
If m Mod 4 = 1 Then
a = 4
b2 = 8 * n1 - 4 * b * x - 2 * x + 2
c = 4 * n1 ^ 2 - 8 * b * x * n1 - 4 * n1 * x + 2 * n1 - 2 * b * x - 2 * x
d = 2 * x * d - 4 * b * x * n1 ^ 2 - 2 * n1 * x - 2 * b * x * n1 - 2 * n1 * x
Else
a = 4
b2 = 8 * n1 - 4 * b * x - 2 * x + 2
c = 4 * n1 ^ 2 - 8 * b * x * n1 - 4 * n1 * x + 2 * n1 - 2 * b * x - 4 * x
d = 2 * x * d - 4 * b * x * n1 ^ 2 - 2 * n1 * x - 2 * b * x * n1 - 4 * n1 * x
End If
Text3 = "输入m=" & m & " p=" & p & "输出a=" & a & "b=" & b2 & "c=" & c & " d=" & d
End Sub
Private Sub Command2_Click()
Text1 = ""
Text2 = ""
Text3 = ""
End Sub
ysr 发表于 2024-1-5 11:18
输入m=14585113 p=975输出a=4b=2800.04989035974c=-120338.165417003 d=17207854.860784(这是一元三次 ...
输入1:a=4,b=2800.04989035974,c=-120338.165417003,d=17207854.860784;输出结果1: x1=-747.9267789072,x2=23.9571531586+71.9576698187i,x3=23.9571531586-71.9576698187i m=-3966902711.9967554016 n=1797674006.1047658716
而975+23=998，实际因子是997，是巧合吧？
输入m=6.958000001675E+21 p=68683288434输出a=4b=245233816433.822c=2.69128261764411E+21 d=2.35869705497168E+31
输入m=14585111 p=975输出a=4b=830.102564102564c=-2040316.05128205 d=48088500.0000001
实际14585111=173*84307
Private Sub Command1_Click()
Dim a, b, c
m = Val(Text1)
p = Val(Text2)
q = Val(m / p)
If InStr(q, ".") > 0 Then
q = Mid(q, 1, InStr(q, ".") - 1)
Else
q = q
End If
m1 = Right(m, 2)
If m1 Mod 4 = 1 Then
b = ((q - 1) / p - 1) / 2
d = (m - 1) / 4
Else
b = ((q - 3) / p - 1) / 2
d = (m - 3) / 4
End If
n1 = (p - 1) / 2
X = zhengliys2(1, Len(p) - 1)
If m1 Mod 4 = 1 Then
a = 4
b2 = 8 * n1 - 4 * b * X - 2 * X + 2
c = 4 * n1 ^ 2 - 8 * b * X * n1 - 4 * n1 * X + 2 * n1 - 2 * b * X - 2 * X
d = 2 * X * d - 4 * b * X * n1 ^ 2 - 2 * n1 * X - 2 * b * X * n1 - 2 * n1 * X
Else
a = 4
b2 = 8 * n1 - 4 * b * X - 2 * X + 2
c = 4 * n1 ^ 2 - 8 * b * X * n1 - 4 * n1 * X + 2 * n1 - 2 * b * X - 4 * X
d = 2 * X * d - 4 * b * X * n1 ^ 2 - 2 * n1 * X - 2 * b * X * n1 - 4 * n1 * X
End If
Text3 = "输入m=" & m & " p=" & p & "输出a=" & a & "b=" & b2 & "c=" & c & " d=" & d
End Sub
Private Sub Command2_Click()
Text1 = ""
Text2 = ""
Text3 = ""
End Sub
输入m=6958000001674999998647 p=68683288435输出a=4b=265233816439.41597945753600000000c=4064948386561471129439.27193266738900000000 d=-11203029455730576652764507051908.18176563222500000000
输入m=6958000001674999998647 p=68683288434输出a=4b=265233816438.40409812681800000000c=4064948386491971792140.08185368641800000000 d=-11203029456923937410196226156380.83528051903000000000