《数论探秘》电子版

ysr

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可选大数据高精度的"
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))
  If fhys(Trim(tx1)) >= 0 Then
   jie3cifc = tx1
   ElseIf fhys(Trim(tx12)) >= 0 Then
   jie3cifc = tx12
   Else
   jie3cifc = tx13
   End If
  
  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))
  If fhys(Trim(tx1)) >= 0 Then
   jie3cifc = tx1
   ElseIf fhys(Trim(tx12)) >= 0 Then
   jie3cifc = tx12
   Else
   jie3cifc = tx13
   End If

  End If
  
End If


  


  

End Function

ysr

146376311461999081269634564951与146376311461999081269634584951之间有14对差大于200的素数对:
146376311461999081269634566007/210
146376311461999081269634566239/232
146376311461999081269634567517/210
146376311461999081269634569919/218
146376311461999081269634573531/470
146376311461999081269634574569/230
146376311461999081269634575593/292
146376311461999081269634576591/222
146376311461999081269634577809/320
146376311461999081269634579277/264
146376311461999081269634580471/282
146376311461999081269634583419/438
146376311461999081269634583713/206
146376311461999081269634584751/270
用时7181.339秒

ysr

146376311461999081269634623673与146376311461999081269634643673之间有11对差大于200的素数对:
146376311461999081269634625267/208
146376311461999081269634630259/250
146376311461999081269634631119/336
146376311461999081269634632367/348
146376311461999081269634634141/228
146376311461999081269634638781/240
146376311461999081269634639127/300
146376311461999081269634639673/214
146376311461999081269634640651/208
146376311461999081269634641139/360
146376311461999081269634641787/224
用时-58000.787秒

ysr

146376311461999081269634641787与146376311461999081269634661787之间有14对差大于200的素数对:
146376311461999081269634644853/284
146376311461999081269634645067/214
146376311461999081269634646699/250
146376311461999081269634647143/262
146376311461999081269634648523/214
146376311461999081269634648879/270
146376311461999081269634649291/412
146376311461999081269634653209/396
146376311461999081269634654133/204
146376311461999081269634654667/216
146376311461999081269634658231/270
146376311461999081269634658669/376
146376311461999081269634659343/374
146376311461999081269634659817/254
用时4645.508秒

ysr

146376311461999081269634659817与146376311461999081269634699817之间有28对差大于200的素数对:
146376311461999081269634664033/300
146376311461999081269634664363/216
146376311461999081269634665043/254
146376311461999081269634666009/240
146376311461999081269634666589/406
146376311461999081269634667759/228
146376311461999081269634670039/240
146376311461999081269634672003/204
146376311461999081269634672547/264
146376311461999081269634672807/260
146376311461999081269634673639/228
146376311461999081269634673887/248
146376311461999081269634674097/204
146376311461999081269634675681/224
146376311461999081269634677061/260
146376311461999081269634678471/338
146376311461999081269634682537/246
146376311461999081269634683731/270
146376311461999081269634685021/288
146376311461999081269634687579/276
146376311461999081269634690627/210
146376311461999081269634691077/278
146376311461999081269634694351/214
146376311461999081269634694651/214
146376311461999081269634695089/402
146376311461999081269634697101/222
146376311461999081269634698301/240
146376311461999081269634699763/246
用时22384.119秒

ysr

146376311461999081269634699763与146376311461999081269634709763之间有7对差大于200的素数对:
146376311461999081269634700447/208
146376311461999081269634701073/276
146376311461999081269634701389/226
146376311461999081269634703273/414
146376311461999081269634704419/238
146376311461999081269634707587/298
146376311461999081269634709129/342
用时3872.872秒

大海捞针一样找不到那个大间距

ysr

146376311462001376818808030489与146376311462001376818808070489之间有36对差大于200的素数对:
146376311462001376818808033201/204
146376311462001376818808033559/210
146376311462001376818808034101/212
146376311462001376818808034999/222
146376311462001376818808035677/208
146376311462001376818808035937/216
146376311462001376818808036591/282
146376311462001376818808037263/204
146376311462001376818808039249/278
146376311462001376818808039961/274
146376311462001376818808040497/204
146376311462001376818808041767/208
146376311462001376818808044349/308
146376311462001376818808044983/276
146376311462001376818808045303/314
146376311462001376818808047761/262
146376311462001376818808048213/212
146376311462001376818808049209/218
146376311462001376818808050181/218
146376311462001376818808051193/306
146376311462001376818808053403/272
146376311462001376818808054553/204
146376311462001376818808056811/204
146376311462001376818808059107/454
146376311462001376818808059529/206
146376311462001376818808059763/230
146376311462001376818808060411/234
146376311462001376818808061171/202
146376311462001376818808061429/258
146376311462001376818808064507/274
146376311462001376818808064899/248
146376311462001376818808065289/236
146376311462001376818808065799/266
146376311462001376818808067627/226
146376311462001376818808067843/204
146376311462001376818808069601/534
用时10681.743秒

ysr

146376311462001376818808069601与146376311462001376818808109601之间有20对差大于200的素数对:
146376311462001376818808072759/308
146376311462001376818808073777/286
146376311462001376818808075411/204
146376311462001376818808076143/272
146376311462001376818808076623/366
146376311462001376818808079861/298
146376311462001376818808080397/318
146376311462001376818808082137/260
146376311462001376818808086949/236
146376311462001376818808087343/252
146376311462001376818808088839/218
146376311462001376818808089727/338
146376311462001376818808089973/210
146376311462001376818808091831/220
146376311462001376818808100581/234
146376311462001376818808100827/216
146376311462001376818808103857/374
146376311462001376818808104643/212
146376311462001376818808107701/210
146376311462001376818808109711/202
用时-49784.057秒

ysr

146376311462001376818808109711与146376311462001376818808119711之间有11对差大于200的素数对:
146376311462001376818808112203/422
146376311462001376818808113757/374
146376311462001376818808114409/276
146376311462001376818808115791/272
146376311462001376818808116781/230
146376311462001376818808117091/204
146376311462001376818808117363/272
146376311462001376818808117747/218
146376311462001376818808118279/208
146376311462001376818808118507/214
146376311462001376818808119433/326
用时5326.065秒

ysr

146376311462001376818808126949与146376311462001376818808136949之间有5对差大于200的素数对:
146376311462001376818808130417/310
146376311462001376818808130843/280
146376311462001376818808134271/222
146376311462001376818808134923/364
146376311462001376818808135829/216
用时4395.957秒