![]() |
1 Attachment(s)
Karsten:
I finished 751 to 1000. Here is the file: [ATTACH]5640[/ATTACH] Data for 576 coming up: 197*576^576+1 351*576^576+1 438*576^576+1 1088*576^576+1 1400*576^576+1 2657*576^576+1 3658*576^576+1 3833*576^576+1 4238*576^576+1 4373*576^576+1 5196*576^576+1 5897*576^576+1 5923*576^576+1 6540*576^576+1 7101*576^576+1 7896*576^576+1 8276*576^576+1 569: No hitter |
[QUOTE=3.14159;228071]For each base value, I sieved to 10[sup]9[/sup] and tested from there.[/QUOTE]
[QUOTE=3.14159;228075]I increased the trial-division bound to 131072.[/QUOTE] Which? |
[QUOTE=Charles]Which?
[/QUOTE] 60 to 750, sieved to 10[sup]9[/sup] + 1 751 to 930, trial factoring to 10000 931 to 1000, trial factoring to 131072. |
thats not increasing the bound,thats decreasing the bound
|
[QUOTE=3.14159;228095]60 to 750, sieved to 10[sup]9[/sup] + 1
751 to 930, trial factoring to 10000 931 to 1000, trial factoring to 131072.[/QUOTE] Why TD and not sieving? |
I'm just testing another script (slightly changed one for SG-hunting found [url=http://www.mersenneforum.org/showpost.php?p=222587&postcount=13]here[/url]) which uses cnewpgen and cLLR.
Calling like "kbbp 250 500 1 10000" to search for k*b^b+1 in the range 250<=b<=500 and 1<=k<=10000. Primes found will be written in separate file, sieved-files and LLR-results will be saved (for later checking) and a log file gives timestamps and found primes per base. Here first timings: (cnewpgen sieves to p=1e9) [code] 03:27 ---------- 501.RES: 5 03:29 ---------- 502.RES: 7 03:30 ---------- 503.RES: 6 03:32 ---------- 504.RES: 12 03:33 ---------- 505.RES: 5 03:34 ---------- 506.RES: 6 03:35 ---------- 507.RES: 6 03:36 ---------- 508.RES: 4 [/code] So about 1 to 2 minutes per base here! Will run above range over night and give results later. |
[QUOTE=Charles]Why TD and not sieving?
[/QUOTE] Because it requires manual work, because I have no script for it. |
The whole range for k*b^b+1, k=1-10000, b=500-750 took 6h 35m with my new script!
But: I'll test with the old pfgw-script to check this: I run pfgw with the -a1 option to reduce "round off" errors in some cases. The comparison of the cnewpgen/cllr results with the pfgw results showed: pfgw -q"9238*619^619+1" 9238*619^619+1 is composite: RES64: [AF5E2DAE777932BE] (0.1462s+0.0003s) pfgw -a1 -q"9238*619^619+1" 9238*619^619+1 is 3-PRP! (0.4438s+0.0507s) cllr says: 9238*619^619+1 is not prime. RES64: 724FCB9240CF6A8B. OLD64: 56EF62B6C26E3F9E Time : 113.720 ms. I've posted this in the software-thread. Also I have to complete the pfgw-script for checking against cLLR! |
Using -q9238*619^619+1
Special modular reduction using all-complex FFT length 768 on 9238*619^619+1 9238*619^619+1 is composite: RES64: [AF5E2DAE777932BE] (0.0741s+0.0002s) Using -a1 flag in addition gets me the same problem. Program bug! |
Karsten: Apply an N-1 test to those candidates!
To settle that one: The N-1 says it's prime, too! [code] Primality testing 9238*619^619+1 [N-1, Brillhart-Lehmer-Selfridge] Running N-1 test using base 2 Special modular reduction using zero-padded FFT length 1024 on 9238*619^619+1 Calling Brillhart-Lehmer-Selfridge with factored part 99.77% 9238*619^619+1 is prime! (0.0949s+0.0007s)[/code] Proth says it is prime! [code]9238*619^619 + 1 may be prime. (a = 2) 9238*619^619 + 1 is prime! (a = 5) [1733 digits][/code] Alpertron's says it is prime! (It passed 10 bases of the Miller-Rabin) PARI's BPSW says it is prime! [code](11:06) gp > 9238 * 619^619 + 1 %1 = 10521791737525042862401806764632398042923657646716674823891306969544323994101205573714689176655991243919696037400605876253090814150192720905005250600278527578603715897533064247720049370458189693194742739471331386062060226231152068132449381714203037477244507010150743857086315973307931445506463765180636402221325082260574514705043119532103371205796084602880167838795160083882742990382044393993937877942109159969176893137731505511443594195632129250974266102509197167398347241644619582544739326773184631176946557614228911246124279968467124335096832877276494708071836502597194675835760025207070058819623791589873538017514699104839516369724057476973317779549907820921238467968708101157773715775229137147445624253490812363202905314850203123750279184140101599810266564876794130176650942133667692045560517578711657310396773969859130897909422717222474450785887484430721340768732580606878583663529307704013058922465888983847570246684125064774635387489495048705818313022904661092221971580003670900196543072813582391612419312990433035693932318434028018475159448232131979475843627419327094411710193137723609876562170798883462462890788358773191400558060227927082377847828746744722363109408973769759022853500831936828573713477526188280062405873595063014082515959578625328121022025303038079835479877190289631640250889417325517360526822004037166559595833373545081201533456667326717244754058156866313102720668337785735523213374418933780267085847731441131646783734282382026888566190859947369071339203752307037823982765686606116810496498400500831792697383650409814354309432990915988249213952145419643054072557205187975497633901417779944914352553689998866949563500653804695584455804439178530705791578972065873683664923178708868256986592481613851621055203 (11:06) gp > ispseudoprime(%) %2 = 1[/code] So, I must conclude: It is prime! |
Pari certifies this as prime (4.9 seconds) with Pocklington-Lehmer:
[code]> isprime(9238 * 619^619 + 1, 1) [2 2 1] [31 2 1] [149 2 1] [619 2 1][/code] |
| All times are UTC. The time now is 06:20. |
Powered by vBulletin® Version 3.8.11
Copyright ©2000 - 2021, Jelsoft Enterprises Ltd.