素数筛选及伪素数及素数通项公式资料.docVIP

素数筛选分析及伪素数通项公式 （1） 素数筛选程序 Option Explicit Sub 素数筛选() Dim MyArray() As Variant Dim Arr() As String Dim Temp() As String Dim i As Variant Dim r As Variant Dim j As Variant Dim n As Variant n = Sheet1.Range(B1) 声明动态数组。 ReDim MyArray(n) For i = 1 To n MyArray(i) = pu(i) Next i On Error Resume Next For j = 1 To n Temp = Filter(Arr, MyArray(j)) If UBound(Temp) 0 Then r = r + 1 ReDim Preserve Arr(1 To r) Arr(r) = MyArray(j) End If Next Sheet1.Range(a1).Resize(r, 1) = Application.Transpose(Arr) Sheet1.Range(C1) = UBound(Arr) Sheet1.Range(a1).Resize(1, r) = Arr Sheet1.Range(f1) = Join(Arr, ,) End Sub Function SZU(n As Variant, k As Variant) As Variant Dim MyArray() As Variant Dim Arr() As String Dim Temp() As String Dim i As Variant Dim r As Variant Dim j As Variant Dim m As Variant Dim d As Variant n = n + Int(0.25 * n) ReDim MyArray(n) For i = 1 To n MyArray(i) = pu(i) Next i On Error Resume Next For j = 1 To n Temp = Filter(Arr, MyArray(j)) If UBound(Temp) 0 Then r = r + 1 ReDim Preserve Arr(1 To r) Arr(r) = MyArray(j) End If Next Sheet1.Range(a1).Resize(r, 1) = Application.Transpose(Arr) Sheet1.Range(C1) = UBound(Arr) Sheet1.Range(a1).Resize(1, r) = Arr SZU = Arr(k) End Function Function pu(Y As Variant) As Variant Dim d As Variant d = 3.14159265358979 pu = 3 * Y + 3 * (Sin(0.5 * d * Y)) ^ 2 + (-1) ^ Y If qp(pu) = 0 Then pu = 5 Else pu = 3 * Y + 3 * (Sin(0.5 * d * Y)) ^ 2 + (-1) ^ Y End If End Function Function p(m As Variant, W As Variant) As Variant If W = m Then p = 1 End If If m / W = Int(m / W) And W m Then p = 0 Else p = 1 End If End Function Function qp(m As Variant) As Variant Dim s As Long, s1 As L