几个vb小程序

ysr

1与800之间的素数打头有4组差为2和4096和2的4生素数对: (用时0.015625秒）
/29/31/4127/4129
/59/61/4157/4159
/239/241/4337/4339
/419/421/4517/4519

ysr

1与800之间的素数打头有10组差为2和10和2的4生素数对: (用时1.953125E-02秒）
/5/7/17/19
/17/19/29/31
/29/31/41/43
/59/61/71/73
/137/139/149/151
/179/181/191/193
/227/229/239/241
/269/271/281/283
/419/421/431/433
/809/811/821/823

ysr

1与800之间的素数打头有3组差为2和22和2的4生素数对: (用时0.015625秒）
/5/7/29/31
/17/19/41/43
/617/619/641/643
1与800之间的素数打头有11组差为2和28和2的4生素数对: (用时0.0234375秒）
/11/13/41/43
/29/31/59/61
/41/43/71/73
/71/73/101/103
/107/109/137/139
/149/151/179/181
/197/199/227/229
/239/241/269/271
/281/283/311/313
/431/433/461/463
/569/571/599/601
1与800之间的素数打头有5组差为2和34和2的4生素数对: (用时0.015625秒）
/5/7/41/43
/71/73/107/109
/101/103/137/139
/191/193/227/229
/311/313/347/349

ysr

1与800之间的素数打头有13组差为2和40和2的4生素数对: (用时0.0234375秒）
/17/19/59/61
/29/31/71/73
/59/61/101/103
/107/109/149/151
/137/139/179/181
/149/151/191/193
/197/199/239/241
/227/229/269/271
/239/241/281/283
/269/271/311/313
/419/421/461/463
/599/601/641/643
/617/619/659/661

ysr

'超大整数的加减乘除运算及x进制转换

Option Explicit

Public Function IIf(ByVal blnExp, vtTrue, vtFalse)
    If blnExp Then
        IIf = vtTrue
    Else
        IIf = vtFalse
    End If
End Function

Public Function ChangeType(vtData, vtType)
    Dim ret
    Select Case vtType
    Case vbEmpty
    Case vbNull
        ret = Null
    Case vbInteger
        ret = ChangeType(vtData, vbDouble)
        If ret >= -32768 And ret <= 32767 Then
            ret = CInt(ret)
        Else
            ret = 0
        End If
    Case vbLong
        ret = ChangeType(vtData, vbDouble)
        If ret >= -2147483648 And ret <= 2147483647 Then
            ret = CLng(ret)
        Else
            ret = CLng(0)
        End If
    Case vbSingle
        If IsNumeric(vtData) Then
            ret = CSng(vtData)
        ElseIf VarType(vtData) = vbDecimal Then
            ret = CSng(vtData)
        Else
            ret = CSng(0)
        End If
    Case vbDouble
        If IsNumeric(vtData) Then
            ret = CDbl(vtData)
        ElseIf VarType(vtData) = vbDecimal Then
            ret = CDbl(vtData)
        Else
            ret = CDbl(0)
        End If
    Case vbCurrency
        ret = ChangeType(vtData, vbDouble)
        If ret >= -922337203685477.5808 And ret <= 922337203685477.5807 Then
            ret = CCur(ret)
        Else
            ret = CCur(0)
        End If
    Case vbDate
        If IsDate(vtData) Then
            ret = CDate(vtData)
        End If
    Case vbString
        If Not IsNull(vtData) Then
            ret = CStr(vtData)
        Else
            ret = Empty
        End If
    Case vbBoolean
        ret = ChangeType(vtData, vbDouble)
        ret = CBool(Not ret = 0)
    Case vbByte
        ret = ChangeType(vtData, vbDouble)
        If ret >= 0 And ret <= 255 Then
            ret = CByte(ret)
        Else
            ret = CByte(0)
        End If
    Case Else
        If VarType(vtData) = vbObject Then
            Set ret = vtData
        Else
            ret = vtData
        End If
    End Select
    ChangeType = ret
End Function

Public Function atos(vtData)
    atos = ChangeType(vtData, vbString)
End Function

Public Function atoi(vtData)
    atoi = ChangeType(vtData, vbInteger)
End Function

Public Function atol(vtData)
    atol = ChangeType(vtData, vbLong)
End Function

Public Function atof(vtData)
    atof = ChangeType(vtData, vbDouble)
End Function

Class ImplNumber
Public Function MyFormat(ByVal n)
    Dim p, r, l, i
    p = atos(n)
    l = Len(p)
    ReDim r(l - 1)
    For i = 1 To l
        r(i - 1) = (Asc(Mid(p, i, 1)) And &HFF) - 48
        If r(i - 1) < 0 Or r(i - 1) > 10 Then
            Err.Raise vbObjectError + 1, "ImplNumber.Format", "第" & (i - 1) & "位字符(" & Chr(r(i - 1) + 48) & ")非数字"
        End if
    Next
    MyFormat = r
End Function

Public Function MyFix(ByVal n)
    Dim p
    p = atos(n)
    Do While Left(p, 1) = "0"
        p = Mid(p, 2)
    Loop
    If p = "" Then p = "0"
    MyFix = p
End Function

Private Function Compare(ByVal n1, ByVal n2)
    Dim p1, p2
    Dim l1, l2
    Dim i, i1, i2
    p1 = atos(n1)
    p2 = atos(n2)
    l1 = Len(p1)
    l2 = Len(p2)
    If l1 > l2 Then
        Compare = 1
    ElseIf l1 < l2 Then
        Compare = -1
    ElseIf p1 = p2 Then
        Compare = 0
    Else
        For i = 1 To l1 Step 8
            i1 = CLng(Mid(p1, i, 8))
            i2 = CLng(Mid(p2, i, 8))
            If i1 > i2 Then
                Compare = 1
                Exit For
            ElseIf i1 < i2 Then
                Compare = -1
                Exit For
            End If
        Next
    End If
End Function

Private Function MyNumber(ByVal l)
    Dim r, i
    ReDim r(l)
    For i = 0 To l
        r(i) = 0
    Next
    MyNumber = r
End Function

Private Function MyAdd(ByVal n1, ByVal n2)
    Dim p1, p2, p3
    Dim l1, l2, l3
    Dim i, t
    p1 = MyFormat(n1)
    p2 = MyFormat(n2)
    l1 = UBound(p1)
    l2 = UBound(p2)
    l3 = IIf(l1 > l2, l1, l2) + 1
    p3 = MyNumber(l3)
    t = 0
    For i = 0 To l3
        If l1 - i >= 0 Then t = t + p1(l1 - i)
        If l2 - i >= 0 Then t = t + p2(l2 - i)
        p3(l3 - i) = IIf(t > 9, t - 10, t)
        t = IIf(t > 9, 1, 0)
    Next
    MyAdd = MyFix(Join(p3, ""))
    Erase p1
    Erase p2
    Erase p3
End Function

Private Function MySubtract(ByVal n1, ByVal n2)
    Dim p1, p2, p3, sign
    Dim i, t, l1, l2, l3
    Select Case Compare(n1, n2)
    Case -1
        p1 = MyFormat(n2)
        p2 = MyFormat(n1)
        sign = "-"
    Case 0
        MySubtract = "0"
        Exit Function
    Case 1
        p1 = MyFormat(n1)
        p2 = MyFormat(n2)
    End Select
    l1 = UBound(p1)
    l2 = UBound(p2)
    l3 = l1
    p3 = MyNumber(l3)
    t = 0
    For i = 0 To l3
       If l1 - i >= 0 Then t = p1(l1 - i) - t
       If l2 - i >= 0 Then t = t - p2(l2 - i)
       p3(l3 - i) = IIf(t < 0, t + 10, t)
       t = IIf(t < 0, 1, 0)
    Next
    MySubtract = sign & MyFix(Join(p3, ""))
    Erase p1
    Erase p2
    Erase p3
End Function

'加法
Public Function Add(ByVal n1, ByVal n2)
    Dim s1, s2
    Dim p1, p2
    p1 = MyFix(n1)
    p2 = MyFix(n2)
    s1 = Left(p1, 1)
    s2 = Left(p2, 1)
    If s1 = "-" Then p1 = MyFix(Mid(p1, 2))
    If s2 = "-" Then p2 = MyFix(Mid(p2, 2))
    If s1 = "-" Then
        If s2 = "-" Then
            Add = "-" & MyAdd(p1, p2)
        Else
            Add = MySubstract(p2, p1)
        End If
    Else
        If s2 = "-" Then
            Add = MySubtract(p1, p2)
        Else
            Add = MyAdd(p1, p2)
        End If
    End If
End Function

'减法
Public Function Subtract(ByVal n1, ByVal n2)
    Dim s1, s2
    Dim p1, p2
    p1 = MyFix(n1)
    p2 = MyFix(n2)
    s1 = Left(p1, 1)
    s2 = Left(p2, 1)
    If s1 = "-" Then p1 = MyFix(Mid(p1, 2))
    If s2 = "-" Then p2 = MyFix(Mid(p2, 2))
    If s1 = "-" Then
        If s2 = "-" Then
            Subtract = MySubstract(p2, p1)
        Else
            Subtract = "-" & MyAdd(p1, p2)
        End If
    Else
        If s2 = "-" Then
            Subtract = MyAdd(p1, p2)
        Else
            Subtract = MySubtract(p1, p2)
        End If
    End If
End Function

Private Function MyMultiply(ByVal n1, ByVal n2)
    Dim p1, p2, p3, p4
    Dim l1, l2, l3
    Dim i, k, t
    If Compare(n1, n2) = 1 Then
        p1 = MyFormat(n2)
        p2 = MyFormat(n1)
    Else
        p1 = MyFormat(n1)
        p2 = MyFormat(n2)
    End If
    l1 = UBound(p1)
    l2 = UBound(p2)
    p4 = "0"
    For i = 0 To l1
        l3 = l2 + i + 1
        p3 = MyNumber(l3)
        t = 0
        For k = 0 To l2
            t = t + p1(l1 - i) * p2(l2 - k)
            p3(l3 - i - k) = IIf(t > 9, (t Mod 10), t)
            t = IIf(t > 9, t / 10, 0)
        Next
        If t > 0 Then
            p3(l3 - i - k) = t
        End If
        p4 = MyAdd(p4, MyFix(Join(p3, "")))
        Erase p3
    Next
    MyMultiply = p4
    Erase p1
    Erase p2
End Function

'乘法
Public Function Multiply(ByVal n1, ByVal n2)
    Dim s1, s2
    Dim p1, p2
    p1 = MyFix(n1)
    p2 = MyFix(n2)
    s1 = Left(p1, 1)
    s2 = Left(p2, 1)
    If s1 = "-" Then p1 = MyFix(Mid(p1, 2))
    If s2 = "-" Then p2 = MyFix(Mid(p2, 2))
    If p1 = "0" Or p2 = "0" Then
        Multiply = "0"
    ElseIf s1 = "-" Then
        If s2 = "-" Then
            Multiply = MyMultiply(p1, p2)
        Else
            Multiply = "-" & MyMultiply(p1, p2)
        End If
    Else
        If s2 = "-" Then
            Multiply = "-" & MyMultiply(p1, p2)
        Else
            Multiply = MyMultiply(p1, p2)
        End If
    End If
End Function

Private Function MyDiv(ByVal n1, ByVal n2)
    Dim p(1), n3, i
    n3 = MySubtract(n1, n2)
    i = 1
    Do While Compare(n3, n2) <> -1
        n3 = MySubtract(n3, n2)
        i = i + 1
    Loop
    p(0) = i
    p(1) = n3
    MyDiv = p
End Function

Private Function MyDivision(ByVal n1, ByVal n2, ByVal decimal, ByVal sign)
    Dim p1, p2, p3(1), p4, p5
    Dim i, cmp, l1, l2, lx
    If decimal > 0 Then
        p1 = n1 & String(decimal, "0")
    Else
        p1 = n1
    End If
    p2 = n2
    cmp = Compare(p1, p2)
    If cmp = -1 Then
        p3(0) = 0
        p3(1) = n1
        MyDivision = p3
        Exit Function
    End If
    If cmp = 0 Then
        If decimal > 0 Then
            p3(0) = sign & "0." & String(decimal - 1, "0") & 1
            p3(1) = n1
        Else
            p3(0) = sign & "1"
            p3(1) = 0
        End If
        MyDivision = p3
        Exit Function
    End If
    l1 = Len(p1)
    l2 = Len(p2)
    lx = Len(n1)
    p4 = Mid(p1, 1, l2)
    i = l2
    p3(0) = sign
    If decimal > 0 And i > lx Then
        p3(0) = p3(0) & "." & String(i - lx - 1, "0")
    End If
    Do While i <= l1
        If Compare(p4, p2) <> -1 Then
            p5 = MyDiv(p4, p2)
            p4 = p5(1)
            p3(0) = p3(0) & p5(0)
        ElseIf i = l1 Then
            If i = lx Then p3(1) = MyFix(p4)
            Exit Do
        Else
            If i = lx Then
                p3(1) = MyFix(p4)
                If decimal > 0 Then p3(0) = p3(0) & "."
            End If
            i = i + 1
            p4 = MyFix(p4 & Mid(p1, i, 1))
            If Compare(p4, p2) = -1 Then p3(0) = p3(0) & "0"
        End If
    Loop
    MyDivision = p3
End Function

'除法
'decimal = 小数点后的位数
'函数返回拥有两个元素的数组
'元素0 = 商
'元素1 = 余数
Public Function Division(ByVal n1, ByVal n2, ByVal decimal)
    Dim s1, s2
    Dim p1, p2
    p1 = MyFix(n1)
    p2 = MyFix(n2)
    s1 = Left(p1, 1)
    s2 = Left(p2, 1)
    If s1 = "-" Then p1 = MyFix(Mid(p1, 2))
    If s2 = "-" Then p2 = MyFix(Mid(p2, 2))
    If p1 = "0" Then
        Division = Array(0, 0)
    ElseIf p2 = "0" Then
        Err.Raise vbObjectError + 1, "ImplNumber.Division", "被零除"
    ElseIf s1 = "-" Then
        If s2 = "-" Then
            Division = MyDivision(p1, p2, decimal, "")
        Else
            Division = MyDivision(p1, p2, decimal, "-")
        End If
    Else
        If s2 = "-" Then
            Division = MyDivision(p1, p2, decimal, "-")
        Else
            Division = MyDivision(p1, p2, decimal, "")
        End If
    End If
End Function

'将一个10进制整数进行(2 - 36)进制的转换
Public Function BaseX(ByVal n, ByVal x)
    Dim s, i, p
    If x < 2 Then
        Err.Raise vbObjectError + 1, "ImplNumber.BaseX", "错误的进制"
    End If
    If Compare(n, "0") = 1 Then
        p = Division(n, x, 0)
        s = s & BaseX(p(0), x)
        i = CInt(p(1))
        If i < 10 Then
            s = s & i
        Else
            s = s & Chr(i + 55)
        End If
    End If
    BaseX = s
End Function

'将一个(2 - 36)进制的字符转换成10进制的整数
Public Function ConvertX(ByVal s, ByVal x)
    Dim i
    Dim n, p, t
    If x < 2 Then
        Err.Raise vbObjectError + 1, "ImplNumber.ConvertX", "错误的进制"
    End If
    n = 0
    p = 1
    For i = Len(s) To 1 Step -1
        t = Asc(Mid(s, i, 1)) And &HFF
        If t >= 48 And t <= 57 Then
            t = t - 48
        ElseIf t >= 65 And t < 55 + x And t <= 90 Then
            t = t - 55
        ElseIf t >= 97 And t < 87 + x And t <= 122 Then
            t = t - 87
        Else
            Err.Raise vbObjectError + 1, "ImplNumber.ConvertX", "错误的进制字符串"
        End If
        n = Add(n, Multiply(t, p))
        p = Multiply(p, x)
    Next
    ConvertX = n
End Function
End Class

'范例
Dim num, i, x
Set num = New ImplNumber
WScript.Echo num.Add(784921795923989, 5215632421426)
WScript.Echo num.Subtract(784921795923989, 5215632421426)
WScript.Echo num.Multiply(784921795923989, 5215632421426)
WScript.Echo Join(num.Division(784921795923989, 5215632421426, 12), " - ")
For i = 2 To 36
    x = num.BaseX(784921795923989, i)
    WScript.Echo "Base" & i & "(" & num.ConvertX(x, i) & ") = " & x
Next
Set num = Nothing

ysr

连乘积公式结果： 偶数9866 其方根为99.3277403347121  其方根内最大素数97 方根内的素数个数m=25  每m-1个中的平均值3.75350895581372  总个数为94.4596518857898  方根内能产生的素数对个数：0.950989638619403

我们知道99内有25个素数，其中有24个奇素数，所以，9866的哥德巴赫猜想解的个数的下限为94-24=70，实际为86个，如下为实际值（哟与偶数的哥德巴赫猜想解的个数是波动式上升的，其波动性让我们不容易验证，而9866的哥德巴赫猜想的个数是波谷，所以，用这个典型可以验证该个数公式就是下限），所以这个确实是下限：
9866的方根为99.3277403347121,方根内有4个总数有86个:
9866=7+ 9859
37+ 9829
79+ 9787
97+ 9769
127+ 9739
223+ 9643
433+ 9433
463+ 9403
523+ 9343
547+ 9319
709+ 9157
733+ 9133
739+ 9127
757+ 9109
823+ 9043
853+ 9013
859+ 9007
937+ 8929
1063+ 8803
1087+ 8779
1129+ 8737
1153+ 8713
1237+ 8629
1303+ 8563
1327+ 8539
1399+ 8467
1423+ 8443
1447+ 8419
1489+ 8377
1549+ 8317
1579+ 8287
1597+ 8269
1657+ 8209
1699+ 8167
1777+ 8089
1873+ 7993
1933+ 7933
1987+ 7879
1993+ 7873
1999+ 7867
2113+ 7753
2143+ 7723
2179+ 7687
2293+ 7573
2377+ 7489
2389+ 7477
2473+ 7393
2557+ 7309
2647+ 7219
2659+ 7207
2689+ 7177
2707+ 7159
2797+ 7069
2917+ 6949
3037+ 6829
3163+ 6703
3187+ 6679
3229+ 6637
3259+ 6607
3313+ 6553
3319+ 6547
3469+ 6397
3499+ 6367
3529+ 6337
3637+ 6229
3733+ 6133
3793+ 6073
3823+ 6043
3943+ 5923
4027+ 5839
4129+ 5737
4177+ 5689
4219+ 5647
4243+ 5623
4297+ 5569
4339+ 5527
4363+ 5503
4423+ 5443
4447+ 5419
4519+ 5347
4639+ 5227
4657+ 5209
4759+ 5107
4789+ 5077
4909+ 4957
4933+ 4933

ysr

1与9866之间有204对孪生素数对:
3和 5  孪中4
5和 7  孪中6
11和 13  孪中12
17和 19  孪中18
29和 31  孪中30
41和 43  孪中42
59和 61  孪中60
71和 73  孪中72
101和 103  孪中102
107和 109  孪中108
137和 139  孪中138
149和 151  孪中150
179和 181  孪中180
191和 193  孪中192
197和 199  孪中198
227和 229  孪中228
239和 241  孪中240
269和 271  孪中270
281和 283  孪中282
311和 313  孪中312
347和 349  孪中348
419和 421  孪中420
431和 433  孪中432
461和 463  孪中462
521和 523  孪中522
569和 571  孪中570
599和 601  孪中600
617和 619  孪中618
641和 643  孪中642
659和 661  孪中660
809和 811  孪中810
821和 823  孪中822
827和 829  孪中828
857和 859  孪中858
881和 883  孪中882
1019和 1021  孪中1020
1031和 1033  孪中1032
1049和 1051  孪中1050
1061和 1063  孪中1062
1091和 1093  孪中1092
1151和 1153  孪中1152
1229和 1231  孪中1230
1277和 1279  孪中1278
1289和 1291  孪中1290
1301和 1303  孪中1302
1319和 1321  孪中1320
1427和 1429  孪中1428
1451和 1453  孪中1452
1481和 1483  孪中1482
1487和 1489  孪中1488
1607和 1609  孪中1608
1619和 1621  孪中1620
1667和 1669  孪中1668
1697和 1699  孪中1698
1721和 1723  孪中1722
1787和 1789  孪中1788
1871和 1873  孪中1872
1877和 1879  孪中1878
1931和 1933  孪中1932
1949和 1951  孪中1950
1997和 1999  孪中1998
2027和 2029  孪中2028
2081和 2083  孪中2082
2087和 2089  孪中2088
2111和 2113  孪中2112
2129和 2131  孪中2130
2141和 2143  孪中2142
2237和 2239  孪中2238
2267和 2269  孪中2268
2309和 2311  孪中2310
2339和 2341  孪中2340
2381和 2383  孪中2382
2549和 2551  孪中2550
2591和 2593  孪中2592
2657和 2659  孪中2658
2687和 2689  孪中2688
2711和 2713  孪中2712
2729和 2731  孪中2730
2789和 2791  孪中2790
2801和 2803  孪中2802
2969和 2971  孪中2970
2999和 3001  孪中3000
3119和 3121  孪中3120
3167和 3169  孪中3168
3251和 3253  孪中3252
3257和 3259  孪中3258
3299和 3301  孪中3300
3329和 3331  孪中3330
3359和 3361  孪中3360
3371和 3373  孪中3372
3389和 3391  孪中3390
3461和 3463  孪中3462
3467和 3469  孪中3468
3527和 3529  孪中3528
3539和 3541  孪中3540
3557和 3559  孪中3558
3581和 3583  孪中3582
3671和 3673  孪中3672
3767和 3769  孪中3768
3821和 3823  孪中3822
3851和 3853  孪中3852
3917和 3919  孪中3918
3929和 3931  孪中3930
4001和 4003  孪中4002
4019和 4021  孪中4020
4049和 4051  孪中4050
4091和 4093  孪中4092
4127和 4129  孪中4128
4157和 4159  孪中4158
4217和 4219  孪中4218
4229和 4231  孪中4230
4241和 4243  孪中4242
4259和 4261  孪中4260
4271和 4273  孪中4272
4337和 4339  孪中4338
4421和 4423  孪中4422
4481和 4483  孪中4482
4517和 4519  孪中4518
4547和 4549  孪中4548
4637和 4639  孪中4638
4649和 4651  孪中4650
4721和 4723  孪中4722
4787和 4789  孪中4788
4799和 4801  孪中4800
4931和 4933  孪中4932
4967和 4969  孪中4968
5009和 5011  孪中5010
5021和 5023  孪中5022
5099和 5101  孪中5100
5231和 5233  孪中5232
5279和 5281  孪中5280
5417和 5419  孪中5418
5441和 5443  孪中5442
5477和 5479  孪中5478
5501和 5503  孪中5502
5519和 5521  孪中5520
5639和 5641  孪中5640
5651和 5653  孪中5652
5657和 5659  孪中5658
5741和 5743  孪中5742
5849和 5851  孪中5850
5867和 5869  孪中5868
5879和 5881  孪中5880
6089和 6091  孪中6090
6131和 6133  孪中6132
6197和 6199  孪中6198
6269和 6271  孪中6270
6299和 6301  孪中6300
6359和 6361  孪中6360
6449和 6451  孪中6450
6551和 6553  孪中6552
6569和 6571  孪中6570
6659和 6661  孪中6660
6689和 6691  孪中6690
6701和 6703  孪中6702
6761和 6763  孪中6762
6779和 6781  孪中6780
6791和 6793  孪中6792
6827和 6829  孪中6828
6869和 6871  孪中6870
6947和 6949  孪中6948
6959和 6961  孪中6960
7127和 7129  孪中7128
7211和 7213  孪中7212
7307和 7309  孪中7308
7331和 7333  孪中7332
7349和 7351  孪中7350
7457和 7459  孪中7458
7487和 7489  孪中7488
7547和 7549  孪中7548
7559和 7561  孪中7560
7589和 7591  孪中7590
7757和 7759  孪中7758
7877和 7879  孪中7878
7949和 7951  孪中7950
8009和 8011  孪中8010
8087和 8089  孪中8088
8219和 8221  孪中8220
8231和 8233  孪中8232
8291和 8293  孪中8292
8387和 8389  孪中8388
8429和 8431  孪中8430
8537和 8539  孪中8538
8597和 8599  孪中8598
8627和 8629  孪中8628
8819和 8821  孪中8820
8837和 8839  孪中8838
8861和 8863  孪中8862
8969和 8971  孪中8970
8999和 9001  孪中9000
9011和 9013  孪中9012
9041和 9043  孪中9042
9239和 9241  孪中9240
9281和 9283  孪中9282
9341和 9343  孪中9342
9419和 9421  孪中9420
9431和 9433  孪中9432
9437和 9439  孪中9438
9461和 9463  孪中9462
9629和 9631  孪中9630
9677和 9679  孪中9678
9719和 9721  孪中9720
9767和 9769  孪中9768
9857和 9859  孪中9858
用时1.59200000000419秒

ysr

1与4933之间有125对孪生素数对:
3和 5  孪中4
5和 7  孪中6
11和 13  孪中12
17和 19  孪中18
29和 31  孪中30
41和 43  孪中42
59和 61  孪中60
71和 73  孪中72
101和 103  孪中102
107和 109  孪中108
137和 139  孪中138
149和 151  孪中150
179和 181  孪中180
191和 193  孪中192
197和 199  孪中198
227和 229  孪中228
239和 241  孪中240
269和 271  孪中270
281和 283  孪中282
311和 313  孪中312
347和 349  孪中348
419和 421  孪中420
431和 433  孪中432
461和 463  孪中462
521和 523  孪中522
569和 571  孪中570
599和 601  孪中600
617和 619  孪中618
641和 643  孪中642
659和 661  孪中660
809和 811  孪中810
821和 823  孪中822
827和 829  孪中828
857和 859  孪中858
881和 883  孪中882
1019和 1021  孪中1020
1031和 1033  孪中1032
1049和 1051  孪中1050
1061和 1063  孪中1062
1091和 1093  孪中1092
1151和 1153  孪中1152
1229和 1231  孪中1230
1277和 1279  孪中1278
1289和 1291  孪中1290
1301和 1303  孪中1302
1319和 1321  孪中1320
1427和 1429  孪中1428
1451和 1453  孪中1452
1481和 1483  孪中1482
1487和 1489  孪中1488
1607和 1609  孪中1608
1619和 1621  孪中1620
1667和 1669  孪中1668
1697和 1699  孪中1698
1721和 1723  孪中1722
1787和 1789  孪中1788
1871和 1873  孪中1872
1877和 1879  孪中1878
1931和 1933  孪中1932
1949和 1951  孪中1950
1997和 1999  孪中1998
2027和 2029  孪中2028
2081和 2083  孪中2082
2087和 2089  孪中2088
2111和 2113  孪中2112
2129和 2131  孪中2130
2141和 2143  孪中2142
2237和 2239  孪中2238
2267和 2269  孪中2268
2309和 2311  孪中2310
2339和 2341  孪中2340
2381和 2383  孪中2382
2549和 2551  孪中2550
2591和 2593  孪中2592
2657和 2659  孪中2658
2687和 2689  孪中2688
2711和 2713  孪中2712
2729和 2731  孪中2730
2789和 2791  孪中2790
2801和 2803  孪中2802
2969和 2971  孪中2970
2999和 3001  孪中3000
3119和 3121  孪中3120
3167和 3169  孪中3168
3251和 3253  孪中3252
3257和 3259  孪中3258
3299和 3301  孪中3300
3329和 3331  孪中3330
3359和 3361  孪中3360
3371和 3373  孪中3372
3389和 3391  孪中3390
3461和 3463  孪中3462
3467和 3469  孪中3468
3527和 3529  孪中3528
3539和 3541  孪中3540
3557和 3559  孪中3558
3581和 3583  孪中3582
3671和 3673  孪中3672
3767和 3769  孪中3768
3821和 3823  孪中3822
3851和 3853  孪中3852
3917和 3919  孪中3918
3929和 3931  孪中3930
4001和 4003  孪中4002
4019和 4021  孪中4020
4049和 4051  孪中4050
4091和 4093  孪中4092
4127和 4129  孪中4128
4157和 4159  孪中4158
4217和 4219  孪中4218
4229和 4231  孪中4230
4241和 4243  孪中4242
4259和 4261  孪中4260
4271和 4273  孪中4272
4337和 4339  孪中4338
4421和 4423  孪中4422
4481和 4483  孪中4482
4517和 4519  孪中4518
4547和 4549  孪中4548
4637和 4639  孪中4638
4649和 4651  孪中4650
4721和 4723  孪中4722
4787和 4789  孪中4788
4799和 4801  孪中4800
4931和 4933  孪中4932
用时0.716000000000349秒

ysr

哥德巴赫猜想解的个数和孪生素数对的个数可以有个统一的下限公式（下限公式有很多，网上看到的大多是二者各不相同的，这个无所谓），哥德巴赫猜想解的个数与孪生素数对个数比较的话注意一点，二者有个等效区间，偶数2A的哥德巴赫猜想解与孪生素数对个数的等效区间，是A内的孪生素数对个数而不是2A内的孪生素数对个数。

如果不是等效区间，那就无法统一公式，无法比较的。孪生素数对个数的增长曲线是不波动的，没有低谷值可言，而哥德巴赫猜想解的个数是波动式上升的，比较的只能是哥德巴赫猜想解的个数的低谷值的增长曲线，就是下限公式，二者的下限公式可以统一。

ysr

比如偶数9866的实际哥德巴赫猜想解的个数是86，而其等效的孪生素数对的区间是4933内的孪生素数对个数，实际是125对，二者数量级一致，都略大于其一个理论下限值，该值可以统一。