[求助]素性测试

ysr

2^1277=有385位，用时0秒2601983048666099770481310081841021384653815561816676201329778087600902014918340074503059860433081046210605403488570251947845891562080866227034976651419330190731032377347305086443295837415395887618239855136922452802923419286887119716740625346109565072933087221327790207134604146257063901166556207972729700461767055550785130256674608872183239507219512717434046725178680177638925792182272

ysr

2601983048666099770481310081841021384653815561816676201329778087600902014918340074503059860433081046210605403488570251947845891562080866227034976651419330190731032377347305086443295837415395887618239855136922452802923419286887119716740625346109565072933087221327790207134604146257063901166556207972729700461767055550785130256674608872183239507219512717434046725178680177638925792182271~2601983048666099770481310081841021384653815561816676201329778087600902014918340074503059860433081046210605403488570251947845891562080866227034976651419330190731032377347305086443295837415395887618239855136922452802923419286887119716740625346109565072933087221327790207134604146257063901166556207972729700461767055550785130256674608872183239507219512717434046725178680177638925792182271之间的素数有个：（用时3178.844秒）
这个数是合数！但目前无法分解

ysr

输入m=14585113 p=975输出a=4  b=2800.04989035974  c=-120338.165417003 d=17207854.860784(这是一元三次方程的系数)

ysr

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

ysr 发表于 2024-1-5 11:18
输入m=14585113 p=975输出a=4  b=2800.04989035974  c=-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，是巧合吧？

ysr

输入m=6.958000001675E+21 p=68683288434输出a=4  b=245233816433.822  c=2.69128261764411E+21 d=2.35869705497168E+31

ysr

输入m=14585111 p=975输出a=4  b=830.102564102564  c=-2040316.05128205 d=48088500.0000001

实际14585111=173*84307

ysr

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

ysr

输入m=6958000001674999998647 p=68683288435输出a=4  b=265233816439.41597945753600000000  c=4064948386561471129439.27193266738900000000 d=-11203029455730576652764507051908.18176563222500000000

ysr

输入m=6958000001674999998647 p=68683288434输出a=4  b=265233816438.40409812681800000000  c=4064948386491971792140.08185368641800000000 d=-11203029456923937410196226156380.83528051903000000000