mersenneforum.org

mersenneforum.org (https://www.mersenneforum.org/index.php)
-   And now for something completely different (https://www.mersenneforum.org/forumdisplay.php?f=119)
-   -   Smarandache prime(s) (https://www.mersenneforum.org/showthread.php?t=20527)

Batalov 2015-10-05 21:10

Smarandache prime(s)
 
[QUOTE="Neil J. A. Sloane, OEIS Foundation"]To Number Theory List,

Consider the sequence with nth term equal to the
[URL="http://mathworld.wolfram.com/ConsecutiveNumberSequences.html"]concatenation of the decimal numbers 1234...n[/URL] ([URL]https://oeis.org/A007908[/URL]).
[B]When is the first prime?[/B] The comments in A007908 say that there should be infinitely many primes, and that there
are no primes among the first 64000 terms.
If you would like to help with this search, you could leave a comment in A007908 saying that there are no primes among terms X through Y, or, of course, that n = Z gives a (probable) prime, which would be pretty exciting.

Best regards, Neil[/QUOTE]
There is a PRPNet server now, working on candidates n>200,000. Everyone is welcome to attach their prpclients.

[B]Instructions:[/B]
1a. Get PRPclient 5.4.0 for your platform from [URL="http://prpnet.primegrid.com/"]here[/URL]. It contains all programs, bundled in.
1b. (optionally) Get PFGW from [URL="http://sourceforge.net/projects/openpfgw/files/"]here[/URL] if you already have prpclient but not PFGW
2. Edit prpclient.ini : add email, username, machineid, server and pfgwexe :
[CODE]//============= in prpclient.ini =============
email=youremail@somewhere.com
userid=UserId
machineid=PhenomII
instanceid=1
//...
server=SmarPRP:100:1:smarandache.ddns.net:1200
//...
pfgwexe=pfgw64.exe
[/CODE][COLOR=Blue]The current status is here: [url]http://smarandache.ddns.net:1200/pending_tests.html[/url] ( or [url]http://99.121.249.54:1200/[/url] )
The current typical test duration is [B]16[/B] hours (on a Broadwell), [B]10-24[/B] hours on any modern CPU.[/COLOR]

The credit for the find will go to "A,B,PFGW", where A is the finder, B is me and PFGW is Mark+George and other authors of PFGW.

[B]Happy hunting! [/B]And read the whole thread, too (optional but important).
[COLOR=LemonChiffon].[/COLOR]

paulunderwood 2015-10-05 21:20

Are you going to take a punt at it, Serge :question:

science_man_88 2015-10-05 21:39

well as far as I can tell only prime indexes need to be checked because for example 123456 = 1001 *123 +333 and both 123 and 333 divide 3 and so on for the other primes but this is likely well known already, in fact it sounds familiar ! not much else I can think of that's useful to help though. edit: doh I forgot that it expands the number of digits faster than doubling for doubling n for the most part.

wombatman 2015-10-05 21:54

[QUOTE]There are no primes in the first 77000 terms. - Max Alekseyev, Oct 03 2015[/QUOTE]

Just wanted to this statement out as it's buried at the end of the comments on the OEIS page.

science_man_88 2015-10-05 22:08

I realize now that 1,12,123 .. means each number adds it's value mod 3 to the sum of digits and so the sum of digits is only not divisible by 3 when you are at an index that is 3x+1 and this is only odd if x is even so we get 6y+1 and then when you follow this up you'll see the sequence 1,7,13,19,25 (25 can't work as it ends in 5), mod 30 and then if any of the endings can't work with an even number in front of them I'd say try that to eliminate tests but once again I'm probably randomly ranting again.

danaj 2015-10-05 23:15

sm88, a half-decent primality test will find these cases (divisible by small values) quite easily, even at hundreds of thousands of digits. Filtering out values divisible by the first 10k or more small primes is also quite fast -- a trivial amount of time compared to running a PRP test. I'd be leery of excluding values to test unless one was quite sure of the math.

science_man_88 2015-10-05 23:31

[QUOTE=danaj;412045]sm88, a half-decent primality test will find these cases (divisible by small values) quite easily, even at hundreds of thousands of digits. Filtering out values divisible by the first 10k or more small primes is also quite fast -- a trivial amount of time compared to running a PRP test. I'd be leery of excluding values to test unless one was quite sure of the math.[/QUOTE]

well of my second post I can be certain because 12 is 0 mod 3 and 123 is 0 mod 3 and the numbers appended have sum digits that are 0 1 and 2 mod 3 0+1+2 =3 if you have a number that is like 123 appending a number that is 1 mod 3 lands you at 1 mod 3 and appending a number that is 2 mod 3 on that will land you at 0 mod 3 and likewise when you append a 0 mod 3 number after that so both 2 and 0 mod 3 are divisible by 3 that leaves only 1 mod 3 or 3x+1 which can only be odd ( the only indexes that can append an odd last digit when x is even leading to 6y+1 and then any of those that have y that are 4 mod 5 lead to a number that will end in 5 so only if y is 0,1,2,3 mod 5 will 6y+1 be a number that doesn't end in 5.

Batalov 2015-10-06 01:45

[QUOTE=wombatman;412043]Just wanted to this statement out as it's buried at the end of the comments on the OEIS page.[/QUOTE]
Yep, that's why I additionally linked to Mathworld (except for the link, the rest is verbatim from N.J.A.S.' posting to NUBTHRY).

I've previously taken a bite on up-and-down-concatentated primes. And just like for them, I'd like to note that for every 10[SUP]k[/SUP] <= n < 10[SUP]k+1[/SUP], there is an explicit summation formula. Additionally, Capt.Obvious reports that only n [FONT=&quot]≡[/FONT] 1,3,7,9 (mod 10) need to be considered. And n>=77000. :rolleyes:

LaurV 2015-10-06 02:12

Sum of digits of any 3 consecutive numbers is divisible by 3, so only 2/3 of the 1,3,5,7 (mod 10) have to be considered :razz:
This is what sm88 tries to tell you hehe, and you don't get it...

edit: pari goes quite fast to few thousands and stays there. Interesting that the "isprime" version is faster than other "probable prime" (modular) tests, because probably the builtin "isprime" filters the numbers which are obviously not primes (like even, etc). It can be seen how the "counter" progresses "in steps", jumping over big chunks of integers.

[code]
gp> i=1; n=1; while(!isprime(n),print(i++); n=eval(concat(Str(n),Str(i))))
[/code]

edit2, this is also faster to generate the sequence than the version on oeis, because it doesn't do any calculus (one should print n instead of i).

Batalov 2015-10-06 02:19

1 Attachment(s)
[QUOTE=LaurV;412052]This is what sm88 tries to tell you hehe, and you don't get it...[/QUOTE]
By G*d, you are right!! I never [I]get[/I] what sm88 tries to whatever... because...

LaurV 2015-10-06 02:29

Bad Batalov! Bad! The [STRIKE]mods[/STRIKE] Gods should not put [STRIKE]people[/STRIKE] their sheep on ignore list, they should read everything, that is why they are [STRIKE]mods[/STRIKE] Gods, to protect the other sheep, etc... with the function comes the responsibility :wink: :leaving:

[YOUTUBE]T-K8qi_AoXI[/YOUTUBE]
[SPOILER]Just re-watched Jackie Brown, in which particular case Michael Keaton's name was the [STRIKE]Birdman[/STRIKE] Lord, as everyone most probably knows. It is a so much more mellow film :rolleyes: [/SPOILER]


All times are UTC. The time now is 17:18.

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