数学中国

 找回密码
 注册
搜索
热搜: 活动 交友 discuz
楼主: 蔡家雄

判定梅森质数的卢卡斯序列

  [复制链接]
发表于 2020-8-23 13:13 | 显示全部楼层
本帖最后由 ysr 于 2020-8-23 05:27 编辑

修改后的结果:(按照你的示例,p算到1~30足够了,下面就是1~30的)
48/50/52/54/56/58/60/64/66/68/70/72/74/76/78/80/82/84/86/88/
90/94/96/100/102/104/106/108/110/112/114/116/118/120/122/124/126/130/132/134/
136/138/140/142/144/146/148/150/152/154/156/160/162/164/166/168/170/174/178/180/
184/186/188/190/192/194/196/198/200/202/204/206/208/210/214/216/218/220/222/224/
226/228/230/234/236/238/240/242/244/246/248/250/252/254/256/258/264/268/270/274/
276/278/280/282/284/286/288/292/294/298/300/304/306/310/314/316/318/320/322/324/
326/328/330/334/336/340/342/344/348/350/354/358/360/364/366/370/374/376/378/380/
384/386/390/392/394/396/400/402/404/406/408/410/412/414/416/418/420/422/426/430/
432/436/438/442/444/446/448/450/452/454/456/460/462/466/468/470/472/474/476/478/
480/484/486/490/492/494/496/498/500/502/504/508/510/514/516/520/522/526/528/532/
534/538/544/550/578/582/584/588/590/594/598/600/604/606/608/610/612/614/616/618/
620/622/624/628/630/634/636/638/640/642/644/646/648/650/654/656/658/660/664/666/
668/670/672/674/676/678/680/684/686/688/690/694/696/698/700/702/704/706/708/712/
714/720/746/750/752/756/762/764/768/770/774/776/780/782/784/786/790/792/794/796/
798/800/802/804/810/816/834/838/840/844/846/850/852/856/860/862/864/866/868/870/
872/874/876/880/882/886/888/890/894/896/898/900/904/906/910/912/914/916/918/920/
922/924/928/930/934/936/940/944/948/950/952/954/958/960/964/966/970/974/978/980/
982/984/988/990/994/996/1000/1004/1006/1008/1010/1012/1014/1020/1024/1026/1028/1030/1032/1034/
1036/1038/1042/1044/1046/1050/1052/1056/1058/1060/1062/1064/1066/1068/1070/1072/1074/1076/1078/1080/
1082/1084/1086/1090/1092/1098/1100/1102/1104/1106/1108/1110/1114/1116/1120/1122/1124/1128/1130/1134/
1136/1140/1146/1152/1160/1164/1166/1170/1176/1182/1188/1192/1194/1198/1200/1204/1206/1208/1210/1212/
1214/1216/1218/1222/1224/1228/1230/1234/1236/1238/1240/1242/1244/1246/1248/1252/1254/1260/1266/1270/
1272/1276/1282/1286/1288/1290/1292/1296/1300/1302/1306/1308/1312/1314/1318/1320/1324/1326/1328/1330/
1332/1334/1336/1338/1340/1342/1344/1348/1350/1356/1436/1440/1442/1446/1450/1452/1456/1458/1460/1462/
1464/1466/1468/1470/1472/1476/1482/1488/1490/1492/1494/1496/1498/1500/1502/1504/1506/1510/1512/1518/
1522/1524/1528/1530/1534/1536/1540/1546/1552/1560/1564/1566/1570/1576/1582/1586/1590/1592/1594/1596/
1600/1602/1604/1606/1608/1610/1612/1614/1616/1618/1620/1622/1624/1626/1630/1632/1634/1636/1638/1640/
1642/1644/1648/1650/1654/1656/1660/1664/1666/1668/1670/1674/1678/1680/1684/1686/1690/1696/1700/1704/
1706/1708/1710/1712/1714/1716/1720/1722/1726/1728/1730/1734/1736/1740/1746/1752/1760/1764/1766/1770/
1776/1782/1784/1788/1790/1794/1796/1800/1802/1806/1812/1818/1838/1842/1844/1848/1854/1860/1868/1872/
1874/1878/1884/1888/1890/1894/1900/1906/1908/1912/1914/1918/1920/1924/1930/1936/1938/1942/1944/1948/
1954/1960/1986/1990/1992/1996/2002/2008/2010/2014/2016/2020/2024/2026/2028/2030/2032/2034/2036/2038/
2040/2042/2044/2046/2050/2052/2056/2058/2076/2080/2082/2086/2090/2092/2094/2096/2098/2100/2106/2110/
2112/2116/2118/2120/2122/2124/2126/2128/2130/2134/2136/2140/2142/2146/2148/2150/2152/2154/2156/2158/
2160/2164/2166/2168/2170/2172/2174/2178/2184/2190/2244/2248/2250/2254/2256/2258/2260/2262/2264/2266/
2268/2272/2274/2276/2278/2280/2282/2284/2286/2288/2290/2292/2294/2296/2298/2302/2304/2308/2310/2314/
2318/2320/2322/2324/2326/2328/2334/2340/2346/2348/2350/2352/2354/2356/2358/2362/2364/2368/2370/2378/
2382/2384/2388/2390/2392/2394/2398/2400/2404/2406/2410/2418/2422/2424/2428/2430/2434/2436/2440/2446/
2448/2452/2454/2458/2460/2464/2470/2474/2476/2478/2480/2484/2488/2490/2494/2496/2500/2506/2510/2514/
2516/2520/2526/2532/2558/2562/2564/
回复 支持 反对

使用道具 举报

发表于 2020-8-23 13:29 | 显示全部楼层
有了160,下面再算1~6000的,这样可能就全覆盖大于等于64的偶数了。时间长等会吧。
回复 支持 反对

使用道具 举报

发表于 2020-8-23 16:37 | 显示全部楼层
程序结果还没有出来明天算吧,上班了!
回复 支持 反对

使用道具 举报

发表于 2020-8-23 16:39 | 显示全部楼层
1与6000之间有由间距为-2060的2对差为30的素数组成的296215组4生素数对:
/44/7/37/7/37
/48/7/37/11/41
/50/7/37/13/43
/54/7/37/17/47
/60/7/37/23/53
/66/7/37/29/59
/68/7/37/31/61
/74/7/37/37/67
/78/7/37/41/71
/80/7/37/43/73
/90/7/37/53/83
/96/7/37/59/89
/104/7/37/67/97
/108/7/37/71/101
/110/7/37/73/103
/116/7/37/79/109
/120/7/37/83/113
/134/7/37/97/127
/138/7/37/101/131
/144/7/37/107/137
/146/7/37/109/139
/164/7/37/127/157
/174/7/37/137/167
/186/7/37/149/179
/188/7/37/151/181
/200/7/37/163/193
/204/7/37/167/197
/218/7/37/181/211
/230/7/37/193/223
/234/7/37/197/227
/236/7/37/199/229
/248/7/37/211/241
/264/7/37/227/257
/270/7/37/233/263
/276/7/37/239/269
/278/7/37/241/271
/288/7/37/251/281
/300/7/37/263/293
/314/7/37/277/307
/318/7/37/281/311
/320/7/37/283/313
/344/7/37/307/337
/354/7/37/317/347
/374/7/37/337/367
/386/7/37/349/379
/390/7/37/353/383
/396/7/37/359/389
/404/7/37/367/397
/416/7/37/379/409
/426/7/37/389/419
/438/7/37/401/431
/446/7/37/409/439
/456/7/37/419/449
/468/7/37/431/461
/470/7/37/433/463
/486/7/37/449/479
/494/7/37/457/487
/498/7/37/461/491
/516/7/37/479/509
/528/7/37/491/521
/578/7/37/541/571
/584/7/37/547/577
/594/7/37/557/587
/600/7/37/563/593
/606/7/37/569/599
/608/7/37/571/601
/614/7/37/577/607
/624/7/37/587/617
/638/7/37/601/631
/650/7/37/613/643
/654/7/37/617/647
/668/7/37/631/661
/680/7/37/643/673
/684/7/37/647/677
/690/7/37/653/683
/698/7/37/661/691
/746/7/37/709/739
/764/7/37/727/757
/776/7/37/739/769
/780/7/37/743/773
/794/7/37/757/787
/834/7/37/797/827
/846/7/37/809/839
/860/7/37/823/853
/864/7/37/827/857
/866/7/37/829/859
/890/7/37/853/883
/894/7/37/857/887
/914/7/37/877/907
/918/7/37/881/911
/944/7/37/907/937
/948/7/37/911/941
/974/7/37/937/967
/978/7/37/941/971
/984/7/37/947/977
/990/7/37/953/983
/1004/7/37/967/997
/1020/7/37/983/1013(程序结果出来了,后面还有不发了)
回复 支持 反对

使用道具 举报

发表于 2020-8-23 16:39 | 显示全部楼层
48/50/52/54/56/58/60/64/66/68/70/72/74/76/78/80/82/84/86/88/
90/92/94/96/98/100/102/104/106/108/110/112/114/116/118/120/122/124/126/128/
130/132/134/136/138/140/142/144/146/148/150/152/154/156/158/160/162/164/166/168/
170/172/174/176/178/180/182/184/186/188/190/192/194/196/198/200/202/204/206/208/
210/212/214/216/218/220/222/224/226/228/230/232/234/236/238/240/242/244/246/248/
250/252/254/256/258/260/262/264/266/268/270/272/274/276/278/280/282/284/286/288/
290/292/294/296/298/300/302/304/306/308/310/312/314/316/318/320/322/324/326/328/
330/332/334/336/338/340/342/344/346/348/350/352/354/356/358/360/362/364/366/368/
370/372/374/376/378/380/382/384/386/388/390/392/394/396/398/400/402/404/406/408/
410/412/414/416/418/420/422/424/426/428/430/432/434/436/438/440/442/444/446/448/
450/452/454/456/458/460/462/464/466/468/470/472/474/476/478/480/482/484/486/488/
490/492/494/496/498/500/502/504/506/508/510/512/514/516/518/520/522/524/526/528/
530/532/534/536/538/540/542/544/546/548/550/552/554/556/558/560/562/564/566/568/
570/572/574/576/578/580/582/584/586/588/590/592/594/596/598/600/602/604/606/608/
610/612/614/616/618/620/622/624/626/628/630/632/634/636/638/640/642/644/646/648/
650/652/654/656/658/660/662/664/666/668/670/672/674/676/678/680/682/684/686/688/
690/692/694/696/698/700/702/704/706/708/710/712/714/716/718/720/722/724/726/728/
730/732/734/736/738/740/742/744/746/748/750/752/754/756/758/760/762/764/766/768/
770/772/774/776/778/780/782/784/786/788/790/792/794/796/798/800/802/804/806/808/
810/812/814/816/818/820/822/824/826/828/830/832/834/836/838/840/842/844/846/848/
850/852/854/856/858/860/862/864/866/868/870/872/874/876/878/880/882/884/886/888/
890/892/894/896/898/900/902/904/906/908/910/912/914/916/918/920/922/924/926/928/
930/932/934/936/938/940/942/944/946/948/950/952/954/956/958/960/962/964/966/968/
970/972/974/976/978/980/982/984/986/988/990/992/994/996/998/1000/1002/1004/1006/1008/
1010/1012/1014/1016/1018/1020/1022/1024/1026/1028/1030/1032/1034/1036/1038/1040/1042/1044/1046/1048/
1050/1052/1054/1056/1058/1060/1062/1064/1066/1068/1070/1072/1074/1076/1078/1080/1082/1084/1086/1088/
1090/1092/1094/1096/1098/1100/1102/1104/1106/1108/1110/1112/1114/1116/1118/1120/1122/1124/1126/1128/
1130/1132/1134/1136/1138/1140/1142/1144/1146/1148/1150/1152/1154/1156/1158/1160/1162/1164/1166/1168/
1170/1172/1174/1176/1178/1180/1182/1184/1186/1188/1190/1192/1194/1196/1198/1200/1202/1204/1206/1208/
1210/1212/1214/1216/1218/1220/1222/1224/1226/1228/1230/1232/1234/1236/1238/1240/1242/1244/1246/1248/
1250/1252/1254/1256/1258/1260/1262/1264/1266/1268/1270/1272/1274/1276/1278/1280/1282/1284/1286/1288/
1290/1292/1294/1296/1298/1300/1302/1304/1306/1308/1310/1312/1314/1316/1318/1320/1322/1324/1326/1328/
1330/1332/1334/1336/1338/1340/1342/1344/1346/1348/1350/1352/1354/1356/1358/1360/1362/1364/1366/1368/
1370/1372/1374/1376/1378/1380/1382/1384/1386/1388/1390/1392/1394/1396/1398/1400/1402/1404/1406/1408/
1410/1412/1414/1416/1418/1420/1422/1424/1426/1428/1430/1432/1434/1436/1438/1440/1442/1444/1446/1448/
1450/1452/1454/1456/1458/1460/1462/1464/1466/1468/1470/1472/1474/1476/1478/1480/1482/1484/1486/1488/
1490/1492/1494/1496/1498/1500/1502/1504/1506/1508/1510/1512/1514/1516/1518/1520/1522/1524/1526/1528/
1530/1532/1534/1536/1538/1540/1542/1544/1546/1548/1550/1552/1554/1556/1558/1560/1562/1564/1566/1568/
1570/1572/1574/1576/1578/1580/1582/1584/1586/1588/1590/1592/1594/1596/1598/1600/1602/1604/1606/1608/
1610/1612/1614/1616/1618/1620/1622/1624/1626/1628/1630/1632/1634/1636/1638/1640/1642/1644/1646/1648/
1650/1652/1654/1656/1658/1660/1662/1664/1666/1668/1670/1672/1674/1676/1678/1680/1682/1684/1686/1688/
1690/1692/1694/1696/1698/1700/1702/1704/1706/1708/1710/1712/1714/1716/1718/1720/1722/1724/1726/1728/
1730/1732/1734/1736/1738/1740/1742/1744/1746/1748/1750/1752/1754/1756/1758/1760/1762/1764/1766/1768/
1770/1772/1774/1776/1778/1780/1782/1784/1786/1788/1790/1792/1794/1796/1798/1800/1802/1804/1806/1808/
1810/1812/1814/1816/1818/1820/1822/1824/1826/1828/1830/1832/1834/1836/1838/1840/1842/1844/1846/1848/
1850/1852/1854/1856/1858/1860/1862/1864/1866/1868/1870/1872/1874/1876/1878/1880/1882/1884/1886/1888/
1890/1892/1894/1896/1898/1900/1902/1904/1906/1908/1910/1912/1914/1916/1918/1920/1922/1924/1926/1928/
1930/1932/1934/1936/1938/1940/1942/1944/1946/1948/1950/1952/1954/1956/1958/1960/1962/1964/1966/1968/
1970/1972/1974/1976/1978/1980/1982/1984/1986/1988/1990/1992/1994/1996/1998/2000/2002/2004/2006/2008/
2010/2012/2014/2016/2018/2020/2022/2024/2026/2028/2030/2032/2034/2036/2038/2040/2042/2044/2046/2048/
2050/2052/2054/2056/2058/2060/2062/2064/2066/2068/2070/2072/2074/2076/2078/2080/2082/2084/2086/2088/
2090/2092/2094/2096/2098/2100/2102/2104/2106/2108/2110/2112/2114/2116/2118/2120/2122/2124/2126/2128/
2130/2132/2134/2136/2138/2140/2142/2144/2146/2148/2150/2152/2154/2156/2158/2160/2162/2164/2166/2168/
2170/2172/2174/2176/2178/2180/2182/2184/2186/2188/2190/2192/2194/2196/2198/2200/2202/2204/2206/2208/
2210/2212/2214/2216/2218/2220/2222/2224/2226/2228/2230/2232/2234/2236/2238/2240/2242/2244/2246/2248/
2250/2252/2254/2256/2258/2260/2262/2264/2266/2268/2270/2272/2274/2276/2278/2280/2282/2284/2286/2288/
2290/2292/2294/2296/2298/2300/2302/2304/2306/2308/2310/2312/2314/2316/2318/2320/2322/2324/2326/2328/
2330/2332/2334/2336/2338/2340/2342/2344/2346/2348/2350/2352/2354/2356/2358/2360/2362/2364/2366/2368/
2370/2372/2374/2376/2378/2380/2382/2384/2386/2388/2390/2392/2394/2396/2398/2400/2402/2404/2406/2408/
2410/2412/2414/2416/2418/2420/2422/2424/2426/2428/2430/2432/2434/2436/2438/2440/2442/2444/2446/2448/
2450/2452/2454/2456/2458/2460/2462/2464/2466/2468/2470/2472/2474/2476/2478/2480/2482/2484/2486/2488/
2490/2492/2494/2496/2498/2500/2502/2504/2506/2508/2510/2512/2514/2516/2518/2520/2522/2524/2526/2528/
2530/2532/2534/2536/2538/2540/2542/2544/2546/2548/2550/2552/2554/2556/2558/2560/2562/2564/
看看是否有缺项?
回复 支持 反对

使用道具 举报

发表于 2020-8-23 16:40 | 显示全部楼层
大于等于64的偶数已经无缺项了。
回复 支持 反对

使用道具 举报

发表于 2020-8-23 16:43 | 显示全部楼层
代码如下:
Private Sub Command1_Click()
Dim a, b
a = Val(Text1)
a1 = a
q = Val(Text2)
m = Sqr(q)
If Right(a, 1) Mod 2 = 0 Then
a = a + 1
Else
a = a
End If
s = 0
a2 = a

Do While a2 < q
n = a2 + 1 + 30
Do While n < 100
B1 = a2

b2 = B1 + 30
B3 = n - B1 - 30
b4 = B3 + 30

c1 = fenjieyinzi0(Val(B1))
C2 = fenjieyinzi0(Val(b2))
c3 = fenjieyinzi0(Val(B3))
D1 = fenjieyinzi0(Val(b4))

If InStr(c1, "*") = 0 And InStr(D1, "*") = 0 And InStr(C2, "*") = 0 And InStr(c3, "*") = 0 And n < 2506 + 60 Then
s = s + 1
Print B1, b2, B3, b4
Text3 = Text3 & "/" & n & "/" & B1 & "/" & b2 & "/" & B3 & "/" & b4 & vbCrLf

s187 = s187 & "/" & n
Else
s = s
End If
n = n + 2
Loop
a2 = a2
s1 = s
Do While n < 10000
B1 = a2

b2 = B1 + 30
B3 = n - B1 - 30
b4 = B3 + 30
Print B3, b4, B1, b2
c1 = fenjieyinzi0(Val(B1))
C2 = fenjieyinzi0(Val(b2))
c3 = fenjieyinzi0(Val(B3))
D1 = fenjieyinzi0(Val(b4))

If InStr(c1, "*") = 0 And InStr(D1, "*") = 0 And InStr(C2, "*") = 0 And InStr(c3, "*") = 0 And n < 2506 + 60 Then
s1 = s1 + 1
Print B1, b2, B3, b4
Text3 = Text3 & "/" & n & "/" & B1 & "/" & b2 & "/" & B3 & "/" & b4 & vbCrLf

s187 = s187 & "/" & n
Else
s1 = s1
End If
n = n + 2
Loop
a2 = a2 + 2
s183 = s183 + s1
Loop
Combo1 = a1 & "与" & q & "之间有由间距为" & B3 - b2 & "的2对差为30的素数组成的" & s183 & "组4生素数对:" & vbCrLf & Text3

s103 = s187
s103 = Mid(s103, 1)
Dim i As Integer
Dim ak(), s105, cr(), f
Set f = CreateObject("Scripting.Dictionary")
s105 = Split(s103, "/")
   j1 = UBound(s105)
   Print j1
   For k = 1 To j1
      n1 = n1 + 1
       ReDim Preserve ak(1 To n1)
      ak(n1) = s105(n1)
    Next
    Print ak(1)
     n = 0
           For i = 1 To j1 - 1
             n = n + 1
             ReDim Preserve cr(1 To n)
            m = Val(ak(i + 1))
            f(m) = ""
           
      Next
      n = 0
      m = f.Keys
      For i = 0 To f.Count - 1
          ReDim Preserve cr(1 To i + 1)
          cr(i + 1) = m(i)
      Next
     For i = 1 To UBound(cr) - 1
        For j = i + 1 To UBound(cr)
            If cr(i) > cr(j) Then
                temp = cr(j)
                cr(j) = cr(i)
                cr(i) = temp  'c数组是排序好的
            End If
        Next j
        
       ' If i Mod 20 = 0 Then
       ' s104 = s104 & temp & "/" & vbCrLf
       ' Else
       ' s104 = s104 & temp & "/"
       ' End If
    Next i
   
      For i = 1 To UBound(cr)
        If i Mod 20 = 0 Then
          s104 = s104 & cr(i) & "/" & vbCrLf
        Else
          s104 = s104 & cr(i) & "/"
        End If
     Next
         Print temp
         MsgBox "ok"
     MsgBox s104  '显示数组
Text4 = s104

End Sub

Private Sub Command2_Click()
Text1 = ""
Text2 = ""
Text3 = ""
Combo1 = ""
Form1.Cls
End Sub
回复 支持 反对

使用道具 举报

 楼主| 发表于 2020-8-24 14:57 | 显示全部楼层
例 2n=9900, p=13, 41, 67, 83, 151, 181, 239, 251, 349, 379, 409, 433, 479, 577, 643, 661, 743, 857, 937, 1033, 1039, 1063, 1163, 1193, 1201, 1481, 1583, 1607, 1627, 1637, 1753, 1831, 2111, 2309, 2311, 2341, 2393, 2411, 2549, 2617, 2657, 2663, 2857, 2887, 2909, 3037, 3137, 3191, 3301, 3511, 3527, 3541, 3583, 3613, 3671, 3697, 3833, 3889, 4019, 4021, 4049, 4211, 4231, 4421, 4451, 4463, 4483, 4519, 4567, 4591, 4673, 4691, 4703,

使 p与p+30 及 2n-p与2n-p-30 均为素数,

则 2n=素数(p)+素数(2n-p)=素数(p+30)+素数(2n-p-30) 均有解。


评分

参与人数 2威望 +20 收起 理由
cz1 + 10 赞一个!
wlc1 + 10 赞一个!

查看全部评分

回复 支持 反对

使用道具 举报

 楼主| 发表于 2020-8-24 15:52 | 显示全部楼层
例 2n=2000, p=211, 277, 307, 541, 547, 571, 877, 907, 937,

使 p与p+30 及 2n-p与2n-p-30 均为素数,

则 2n=素数(p)+素数(2n-p)=素数(p+30)+素数(2n-p-30) 均有解。
回复 支持 反对

使用道具 举报

发表于 2020-8-24 23:58 | 显示全部楼层
定理:4生素数组p,p+2,p+2+6n+4,p+4+6n+4有无穷多组。
此定理很重要,很有趣,可以用来证明李明波大哥的孪中差猜想和孪中和猜想!证明都是成立的!
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|手机版|小黑屋|数学中国 ( 京ICP备05040119号 )

GMT+8, 2024-3-29 20:56 , Processed in 0.089844 second(s), 16 queries .

Powered by Discuz! X3.4

Copyright © 2001-2020, Tencent Cloud.

快速回复 返回顶部 返回列表