mersenneforum.org  

Go Back   mersenneforum.org > Extra Stuff > Blogorrhea > sweety439

Reply
 
Thread Tools
Old 2016-12-28, 20:00   #45
Batalov
 
Batalov's Avatar
 
"Serge"
Mar 2008
Phi(4,2^7658614+1)/2

22×2,281 Posts
Default

Quote:
Originally Posted by sweety439 View Post
(269*10^n+1)/9 seems to be tested to n=10000 (I ran the program of 25 hours and 38 minutes!!!), still no (probable) prime found, base released.
Mr Winston Wolf once said: "That's 30 minutes away. I'll be there in 10."

In comparison, it takes you 25 hours and 38 minutes to do a job that takes 10 minutes.
That doesn't bother you? Not at all?
Batalov is offline   Reply With Quote
Old 2016-12-28, 20:26   #46
gd_barnes
 
gd_barnes's Avatar
 
May 2007
Kansas; USA

23·3·52·17 Posts
Default

Quote:
Originally Posted by Batalov View Post
Mr Winston Wolf once said: "That's 30 minutes away. I'll be there in 10."

In comparison, it takes you 25 hours and 38 minutes to do a job that takes 10 minutes.
That doesn't bother you? Not at all?
Hey now! How ruuuuuuude! All of my machines are more than 4 years old. I ran a 5+ year old version of PFGW (version 3.3.6) with no sieving (trial-factoring only) on my newest 4-year-old machine with 5 other cores running at the same time and it took 18 minutes.

So now I'm sad that you make fun of my old machines and software that take almost twice as long as you say it should. (lol)

Anyway...a modern machine running a modern version of PFGW with some sieving done should likely take ~5 minutes.

I don't think anything bothers Sweety. We can insult him all day long and he just keeps posting trivial results.

Last fiddled with by gd_barnes on 2016-12-28 at 20:28
gd_barnes is offline   Reply With Quote
Old 2016-12-28, 22:43   #47
Batalov
 
Batalov's Avatar
 
"Serge"
Mar 2008
Phi(4,2^7658614+1)/2

22·2,281 Posts
Cool

Quote:
Originally Posted by gd_barnes View Post
I ran a 5+ year old version of PFGW (version 3.3.6) with no sieving (trial-factoring only)
Well for these PRPs with a denominator one should run only LLR or P95. They are a few times faster. And there are <300 tests for n<10,000 and each is much faster than a second, ...so 10 minutes was said only to match the Wolf quote.

Quote:
Originally Posted by Mr.Wolf
I'm Winston Wolf. I solve problems.
Batalov is offline   Reply With Quote
Old 2016-12-29, 00:24   #48
gd_barnes
 
gd_barnes's Avatar
 
May 2007
Kansas; USA

23·3·52·17 Posts
Default

Quote:
Originally Posted by Batalov View Post
Well for these PRPs with a denominator one should run only LLR or P95. They are a few times faster. And there are <300 tests for n<10,000 and each is much faster than a second, ...so 10 minutes was said only to match the Wolf quote.
Hum....How might one run LLR with a denominator? I've never tried it before and just now tried this:

ABC ($a*10^$b$c)/9 // Sieved to 1000000000 with srsieve
269 47 +1
269 51 +1
269 84 +1
269 96 +1
269 105 +1
(etc.)

LLR would not run. It runs fine without the denominator. I'm running LLR 3.8.13...not the latest version but it should be recent enough.
gd_barnes is offline   Reply With Quote
Old 2016-12-29, 02:29   #49
Batalov
 
Batalov's Avatar
 
"Serge"
Mar 2008
Phi(4,2^7658614+1)/2

22·2,281 Posts
Default

Code:
ABC ($a*$b^$c$d)/$e
269 10 47 +1 9
...
You will be very pleasantly surprised. Essentially the speed is the same as if without $e.
Batalov is offline   Reply With Quote
Old 2016-12-29, 07:29   #50
gd_barnes
 
gd_barnes's Avatar
 
May 2007
Kansas; USA

1020010 Posts
Default

Quote:
Originally Posted by Batalov View Post
Code:
ABC ($a*$b^$c$d)/$e
269 10 47 +1 9
...
You will be very pleasantly surprised. Essentially the speed is the same as if without $e.
Thanks! Not bad. About 8 mins on that same old machine with 5 other cores running, which includes about 1 minute to over-sieve to 1G. If done efficiently on a modern machine...likely 3 mins.
gd_barnes is offline   Reply With Quote
Old 2016-12-29, 08:52   #51
LaurV
Romulan Interpreter
 
LaurV's Avatar
 
Jun 2011
Thailand

873410 Posts
Default

Quote:
Originally Posted by gd_barnes View Post
ABC ($a*10^$b$c)/9 // Sieved to 1000000000 with srsieve
269 47 +1
269 51 +1
269 84 +1
269 96 +1
269 105 +1
(etc.)
What version of srsieve can do that sieving?
LaurV is offline   Reply With Quote
Old 2016-12-29, 09:22   #52
axn
 
axn's Avatar
 
Jun 2003

469310 Posts
Default

Quote:
Originally Posted by LaurV View Post
What version of srsieve can do that sieving?
Are you worried about the denominator 9? Just skip 3, and srsieve will take care of correctly sieving the numerator k*b^n+/-c
axn is offline   Reply With Quote
Old 2016-12-29, 10:19   #53
gd_barnes
 
gd_barnes's Avatar
 
May 2007
Kansas; USA

23·3·52·17 Posts
Default

Quote:
Originally Posted by LaurV View Post
What version of srsieve can do that sieving?
Quote:
Originally Posted by axn View Post
Are you worried about the denominator 9? Just skip 3, and srsieve will take care of correctly sieving the numerator k*b^n+/-c
Just to clarify:

I did not actually sieve ($a*10^$b$c)/9. I sieved $a*10^$b$c accomplishing exactly what axn said with -p 4 -P 1e9. After the sieving was done I then added the parens and denominator. Btw, this was a technique that I learned from Serge less than a week ago...one that I was not aware of because I've never searched repeat-digit type forms.

A note about this: In this case the sieving still leaves terms with factors of 3^q after dividing by 9. When running PFGW I get around this by trial factoring to 1% with the -f1 switch so it doesn't do a full primality test. It generally takes less than 1/10th of a second to trial factor to 1%.

This begs a question: Does LLR look for teeny factors before doing a primality test? I could not tell with the short test that I ran on it here. If not I may stick with PFGW because it can quickly get rid of very small factors that are left when sieving a form with a denominator like this.

Last fiddled with by gd_barnes on 2016-12-29 at 10:21
gd_barnes is offline   Reply With Quote
Old 2016-12-29, 10:43   #54
LaurV
Romulan Interpreter
 
LaurV's Avatar
 
Jun 2011
Thailand

2×11×397 Posts
Default

Yes, I did exactly the same in the other thread of sweety, I just created a folder (called bp and copied his bases there, in a vector), filled it with sr_bxxx.pfgw files using pari:

Code:
gp> \r bp\bases
gp > (16:13:53) gp > v
%1 = [184, 185, 200, 210, 269, 281, 306, 311, 326, 331, 371, 380, 384, 385, 394, 396, 452, 465, 485, 511, 522, 570, 574, 598, 601, 629, 631, 632, 636, 640, 649, 670, 684, 691, 693, 711, 713, 731, 752, 759, 771, 795, 820, 861, 866, 872, 881, 932, 938, 948, 951, 956, 963, 996, 1005, 1015]
gp > for(i=1,#v,write("BP\\nrep_b"v[i]".pfgw","ABC $a*"v[i]"^$b$c // Sieved to "nextprime(v[i])" with srsieve");forprime(p=3,10^4,write("BP\\nrep_b"v[i]".pfgw","1 "p" -1")))
time = 1min, 4,259 ms.
This just took a minute and the result was: (random sample)
Code:
ABC $a*636^$b$c // Sieved to 641 with srsieve
1 3 -1
1 5 -1
1 7 -1
1 11 -1
1 13 -1
1 17 -1
1 19 -1
1 23 -1
...etc...
This can be sieved with primes larger than 641 (inclusive) and a simple dos for command solves the problem:

Code:
for %p in (nrep_b*.pfgw) do srsieve -S 2 -m 1e8 -w "%p"
In minutes, all the files with even base are reduced to less then half candidates, and files like "sr_bxxx.pfgw" are created from "nrep_bxxx.fgw". For the odd bases, srsieve gives an error that all terms are even, which I don't know how to avoid.

Then a perl command will replace all headers:
Code:
for %p in (sr_*.pfgw) do perl -i.bak  -p -e "s/ABC \$a\*(\d*)\^\$b\$c/ABC (\$a*\1^\$b\$c)\/(\1-1)/g;" %p
(this is heavy artillery, hehe)

And of course, you will want pfgw to stop when a prime is found, therefore
Code:
for %p in (*.pfgw) do perl -i.bak  -p -e "s/Sieved to \d* with srsieve/{number_primes,\$a,1}/g;" %p
Result:
Code:
ABC ($a*210^$b$c)/(210-1) // {number_primes,$a,1}
1 7 -1
1 11 -1
1 17 -1
1 31 -1
1 103 -1
1 107 -1
1 157 -1
1 173 -1
1 227 -1
1 257 -1
1 269 -1
1 331 -1
1 347 -1
1 353 -1
1 383 -1
1 397 -1
...etc...
This feed to pfgw with another dos for command for all files. I know my tools, trust me, I was just wondering if there is a newer version of srsieve. Especially because the version which is marketed like 1.0.6 on the web, shows 1.0.5 when is executed (I downloaded from many places, and yes, all are the same). So I was more like wondering if there is another, real version 1.0.6

Thanks anyhow...
LaurV is offline   Reply With Quote
Old 2016-12-29, 11:19   #55
gd_barnes
 
gd_barnes's Avatar
 
May 2007
Kansas; USA

237308 Posts
Default

Quote:
Originally Posted by LaurV View Post
For the odd bases, srsieve gives an error that all terms are even, which I don't know how to avoid.
There is no way around it except for forms that are c=-1 or +1 and even the solution for those is not ideal. The program would have to be changed to avoid checking for the error. Both srsieve and sr2sieve have the same problem.

For c=-1 or +1 forms, you can use sr1sieve but it will only sieve where the factors are greater than the base so you'd have to figure a way to remove smaller factors or just do like I did and run PFGW with factoring set to 1% with -f1. Also if you are sieving multiple k's at once, which I am frequently doing on some of these conjectures, you have to sieve them one k at a time...not a great solution.

Would someone want to get ahold of Mark (rogue) and see if he can remove the error check from srsieve and sr2sieve where it says all terms are even? That becomes a major issue when sieving near repeat digit forms for odd bases.

Last fiddled with by gd_barnes on 2016-12-29 at 11:28
gd_barnes is offline   Reply With Quote
Reply

Thread Tools


Similar Threads
Thread Thread Starter Forum Replies Last Post
The reverse Sierpinski/Riesel problem sweety439 sweety439 20 2020-07-03 17:22
Semiprime and n-almost prime candidate for the k's with algebra for the Sierpinski/Riesel problem sweety439 sweety439 10 2018-12-14 21:59
The dual Sierpinski/Riesel problem sweety439 sweety439 12 2017-12-01 21:56
Sierpinski/ Riesel bases 6 to 18 robert44444uk Conjectures 'R Us 139 2007-12-17 05:17
Sierpinski/Riesel Base 10 rogue Conjectures 'R Us 11 2007-12-17 05:08

All times are UTC. The time now is 20:42.

Fri Sep 18 20:42:20 UTC 2020 up 8 days, 17:53, 1 user, load averages: 2.15, 1.79, 1.79

Powered by vBulletin® Version 3.8.11
Copyright ©2000 - 2020, Jelsoft Enterprises Ltd.

This forum has received and complied with 0 (zero) government requests for information.

Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.2 or any later version published by the Free Software Foundation.
A copy of the license is included in the FAQ.