|

楼主 |
发表于 2008-9-20 13:37
|
显示全部楼层
偶数分成两个素数的分法数量可以看作一个概率问题而推导出猜想成立的表达式
[这个贴子最后由愚工688在 2008/09/27 00:31pm 第 2 次编辑]
本着实事求是的原则,我把所用的偶数分成两个素数的QBasic 程序的文本发表在此,以供同好者参考,验证。
偶数分成两个素数的QBasic 程序[编程:愚工688]
": r$ = ""
s1 = 0: s2 = 0: R1 = 2: K = 1
IF INT(a / 2) = a / 2 THEN B = 1: ELSE B = 0
FOR x = B TO C STEP 2
GOSUB pp2
NEXT x
GOSUB pp3
s = s1 + s2: E = (d - s1) / s1
d = INT(100 * d + .5) / 100: E = INT(100 * E + .5) / 100
K = INT(100 * K + .5) / 100
PRINT TAB(0); "M="; m; TAB(10); "S(m)="; s; TAB(21); "S1(m)="; s1; TAB(32); "Sp(m)="; d; TAB(47); "E(m)="; E; TAB(58); "K(m)="; K; TAB(70); "r="; R1
PRINT "* Sp("; m$; ")="; d$ + r$; "="; d
PRINT #1, TAB(1); "M="; m; TAB(12); "S(m)="; s; TAB(24); "S1(m)="; s1; TAB(36); "Sp(m)="; d; TAB(51); "E(m)="; E; TAB(62); "K(m)="; K; TAB(74); "r="; R1
PRINT #1, " * Sp("; m$; ")="; d$ + r$; "="; d
PRINT #1, ""
PRINT
NEXT m
po2: END
pp2: READ r
IF SQR(m - 3) < r THEN
GOTO po3
ELSEIF INT((a + x) / r) = (a + x) / r THEN GOTO po4
ELSEIF INT((a - x) / r) = (a - x) / r AND a - x > r THEN GOTO po4
ELSE GOTO pp2
END IF
po3: IF a - x >= r THEN
s1 = s1 + 1
PRINT #1, a - x; "+"; a + x;
ELSE s2 = s2 + 1
PRINT #1, "("; a - x; "+"; a + x; ")";
END IF
s = s1 + s2
PRINT a - x; "+"; a + x;
po4: RESTORE
RETURN
pp3: READ r
IF SQR(m - 3) < r THEN
GOTO po6
END IF
R1 = r
r0$ = STR$(r): R1$ = STR$(r - 1): R2$ = STR$(r - 2)
IF INT(a / r) = a / r THEN
d = d * (r - 1) / r
ELSE
d = d * (r - 2) / r
END IF
IF INT(a / r) = a / r THEN
r$ = r$ + "*(" + R1$ + "/" + r0$ + ")"
ELSE
r$ = r$ + "*(" + R2$ + "/" + r0$ + ")"
END IF
IF INT(a / r) = a / r THEN K = K * (r - 1) / (r - 2)
GOTO pp3
po6: RESTORE
RETURN
DATA 3 , 5 , 7 , 11 , 13 , 17 , 19 , 23 , 29 , 31 , 37 , 41 , 43 , 47 , 53 , 59 , 61 , 67 , 71 , 73 , 79 , 83 , 89 , 97 , 101 , 103 , 107 , 109 , 113 , 127 , 131 , 137 , 139 , 149 , 151 , 157 , 163 , 167 , 173 , 179 , 181 , 191 , 193 , 197 , 199 , 211 , 223 , 227 , 229 , 233 , 239 , 241 , 251 , 257 , 263 , 269 , 271 , 277 , 281 , 283 , 293 , 307 , 311 , 313 , 317 , 331 , 337 , 347 , 349 , 353 , 359 , 367 , 373 , 379 , 383 , 389 , 397 , 401 , 409 , 419 , 421 , 431 , 433 , 439 , 443 , 449 , 457 , 461 , 463 , 467 , 479 , 487 , 491 , 499
DATA 503 , 509 , 521 , 523 , 541 , 547 , 557 , 563 , 569 , 571 , 577 , 587 , 593 , 599 , 601 , 607 , 613 , 617 , 619 , 631 , 641 , 643 , 647 , 653 , 659 , 661 , 673 , 677 , 683 , 691 , 701 , 709 , 719 , 727 , 733 , 739 , 743 , 751 , 757 , 761 , 769 , 773 , 787 , 797 , 809 , 811 , 821 , 823 , 827 , 829 , 839 , 853 , 857 , 859 , 863 , 877 , 881 , 883 , 887 , 907 , 911 , 919 , 929 , 937 , 941 , 947 , 953 , 967 , 971 , 977 , 983 , 991 , 997
DATA 1009 , 1013 , 1019 , 1021 , 1031 , 1033 , 1039 , 1049 , 1051 , 1061 , 1063 , 1069 , 1087 , 1091 , 1093 , 1097 , 1103 , 1109 , 1117 , 1123 , 1129 , 1151 , 1153 , 1163 , 1171 , 1181 , 1187 , 1193 , 1201 , 1213 , 1217 , 1223 , 1229 , 1231 , 1237 , 1249 , 1259 , 1277 , 1279 , 1283 , 1289 , 1291 , 1297 , 1301 , 1303 , 1307 , 1319 , 1321 , 1327 , 1361 , 1367 , 1373 , 1381 , 1399 , 1409 , 1423 , 1427 , 1429 , 1433 , 1439 , 1447 , 1451 , 1453 , 1459 , 1471 , 1481 , 1483 , 1487 , 1489 , 1493 , 1499
注:DATA 素数数据库决定了可以计算的最大偶数为1499*1499+1,若需超出,需依大小次序扩展素数库。
|
|