Smarandache-Wellin primes: Difference between revisions

Line 303:

7th: index 435 digits 1543 last prime 3037

8th: index 1429 digits 5719 last prime 11927

</pre>

=={{header|Phix}}==

Rather slow...

with javascript_semantics

include mpfr.e

string sw = ""

mpz n = mpz_init(), dn = mpz_init()

sequence swp = {}, swdp = {}, dc = repeat(0,10)

integer np = 0

atom t0 = time(), t1 = time()+1

while length(swp)<iff(platform()=JS?5:8) do

np += 1

integer p = get_prime(np)

string sp = sprintf("%d",p)

sw &= sp

mpz_set_str(n,sw,10)

if mpz_prime(n) then

swp = append(swp,{np,p,shorten(sw)})

end if

for c in sp do

dc[c-'0'+1] += 1

end for

integer dc9 = dc[10]

if odd(dc9) and remainder(dc9,5)!=0 then

string swdps = join(dc,"",fmt:="%d")

mpz_set_str(dn,swdps,10)

if mpz_prime(dn) then

swdp = append(swdp,{np,swdps})

end if

if platform()!=JS and time()>t1 then

progress("processing prime %d...\r",{np})

t1 = time()+1

end if

end while

printf(1,"Smarandache-Whellen primes:\n")

for i,s in swp do

printf(1," %d%s:",{i,ord(i)})

printf(1," Index: %4d, Last prime %5d, %s\n",s)

end for

printf(1,"Smarandache-Whellen derived primes:\n")

for i,s in swdp do

printf(1," %d%s:",{i,ord(i)})

printf(1," Index: %4d, %s\n",s)

end for

?elapsed(time()-t0)

<pre>

Smarandache-Whellen primes:

1st: Index: 1, Last prime 2, 2

2nd: Index: 2, Last prime 3, 23

3rd: Index: 4, Last prime 7, 2357

4th: Index: 128, Last prime 719, 23571113171923293137...73677683691701709719 (355 digits)

5th: Index: 174, Last prime 1033, 23571113171923293137...10131019102110311033 (499 digits)

6th: Index: 342, Last prime 2297, 23571113171923293137...22732281228722932297 (1,171 digits)

7th: Index: 435, Last prime 3037, 23571113171923293137...30013011301930233037 (1,543 digits)

8th: Index: 1429, Last prime 11927, 23571113171923293137...11903119091192311927 (5,719 digits)

Smarandache-Whellen derived primes:

1st: Index: 32, 4194123321127

2nd: Index: 72, 547233879626521

3rd: Index: 73, 547233979727521

4th: Index: 134, 13672766322929571043

5th: Index: 225, 3916856106393739943689

6th: Index: 303, 462696313560586013558131

7th: Index: 309, 532727113760586013758133

8th: Index: 363, 6430314317473636515467149

9th: Index: 462, 8734722823685889120488197

10th: Index: 490, 9035923128899919621189209

11th: Index: 495, 9036023329699969621389211

12th: Index: 522, 9337023533410210710923191219

13th: Index: 538, 94374237357103109113243102223

14th: Index: 624, 117416265406198131121272110263

15th: Index: 721, 141459282456260193137317129313

16th: Index: 738, 144466284461264224139325131317

17th: Index: 790, 156483290479273277162351153339

18th: Index: 852, 164518312512286294233375158359

19th: Index: 1087, 208614364610327343341589284471

20th: Index: 1188, 229667386663354357356628334581

"13 minutes and 13s"

</pre>

@@ Line 303: / Line 303: @@
 th: index  435  digits 1543  last prime  3037
 th: index 1429  digits 5719  last prime 11927
+</pre>
+=={{header|Phix}}==
+Rather slow...
+<!--<syntaxhighlight lang="phix">(phixonline)-->
+ <span style="color: #008080;">with</span> <span style="color: #008080;">javascript_semantics</span>
+ <span style="color: #008080;">include</span> <span style="color: #004080;">mpfr</span><span style="color: #0000FF;">.</span><span style="color: #000000;">e</span>
+ <span style="color: #004080;">string</span> <span style="color: #000000;">sw</span> <span style="color: #0000FF;">=</span> <span style="color: #008000;">""</span>
+ <span style="color: #004080;">mpz</span> <span style="color: #000000;">n</span> <span style="color: #0000FF;">=</span> <span style="color: #7060A8;">mpz_init</span><span style="color: #0000FF;">(),</span> <span style="color: #000000;">dn</span> <span style="color: #0000FF;">=</span> <span style="color: #7060A8;">mpz_init</span><span style="color: #0000FF;">()</span>
+ <span style="color: #004080;">sequence</span> <span style="color: #000000;">swp</span> <span style="color: #0000FF;">=</span> <span style="color: #0000FF;">{},</span> <span style="color: #000000;">swdp</span> <span style="color: #0000FF;">=</span> <span style="color: #0000FF;">{},</span> <span style="color: #000000;">dc</span> <span style="color: #0000FF;">=</span> <span style="color: #7060A8;">repeat</span><span style="color: #0000FF;">(</span><span style="color: #000000;">0</span><span style="color: #0000FF;">,</span><span style="color: #000000;">10</span><span style="color: #0000FF;">)</span>
+ <span style="color: #004080;">integer</span> <span style="color: #000000;">np</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">0</span>
+ <span style="color: #004080;">atom</span> <span style="color: #000000;">t0</span> <span style="color: #0000FF;">=</span> <span style="color: #7060A8;">time</span><span style="color: #0000FF;">(),</span> <span style="color: #000000;">t1</span> <span style="color: #0000FF;">=</span> <span style="color: #7060A8;">time</span><span style="color: #0000FF;">()+</span><span style="color: #000000;">1</span>
+ <span style="color: #008080;">while</span> <span style="color: #7060A8;">length</span><span style="color: #0000FF;">(</span><span style="color: #000000;">swp</span><span style="color: #0000FF;">)<</span><span style="color: #7060A8;">iff</span><span style="color: #0000FF;">(</span><span style="color: #7060A8;">platform</span><span style="color: #0000FF;">()=</span><span style="color: #004600;">JS</span><span style="color: #0000FF;">?</span><span style="color: #000000;">5</span><span style="color: #0000FF;">:</span><span style="color: #000000;">8</span><span style="color: #0000FF;">)</span> <span style="color: #008080;">do</span>
+     <span style="color: #000000;">np</span> <span style="color: #0000FF;">+=</span> <span style="color: #000000;">1</span>
+     <span style="color: #004080;">integer</span> <span style="color: #000000;">p</span> <span style="color: #0000FF;">=</span> <span style="color: #7060A8;">get_prime</span><span style="color: #0000FF;">(</span><span style="color: #000000;">np</span><span style="color: #0000FF;">)</span>
+     <span style="color: #004080;">string</span> <span style="color: #000000;">sp</span> <span style="color: #0000FF;">=</span> <span style="color: #7060A8;">sprintf</span><span style="color: #0000FF;">(</span><span style="color: #008000;">"%d"</span><span style="color: #0000FF;">,</span><span style="color: #000000;">p</span><span style="color: #0000FF;">)</span>
+     <span style="color: #000000;">sw</span> <span style="color: #0000FF;">&=</span> <span style="color: #000000;">sp</span>
+     <span style="color: #7060A8;">mpz_set_str</span><span style="color: #0000FF;">(</span><span style="color: #000000;">n</span><span style="color: #0000FF;">,</span><span style="color: #000000;">sw</span><span style="color: #0000FF;">,</span><span style="color: #000000;">10</span><span style="color: #0000FF;">)</span>
+     <span style="color: #008080;">if</span> <span style="color: #7060A8;">mpz_prime</span><span style="color: #0000FF;">(</span><span style="color: #000000;">n</span><span style="color: #0000FF;">)</span> <span style="color: #008080;">then</span>
+         <span style="color: #000000;">swp</span> <span style="color: #0000FF;">=</span> <span style="color: #7060A8;">append</span><span style="color: #0000FF;">(</span><span style="color: #000000;">swp</span><span style="color: #0000FF;">,{</span><span style="color: #000000;">np</span><span style="color: #0000FF;">,</span><span style="color: #000000;">p</span><span style="color: #0000FF;">,</span><span style="color: #7060A8;">shorten</span><span style="color: #0000FF;">(</span><span style="color: #000000;">sw</span><span style="color: #0000FF;">)})</span>
+     <span style="color: #008080;">end</span> <span style="color: #008080;">if</span>
+     <span style="color: #008080;">for</span> <span style="color: #000000;">c</span> <span style="color: #008080;">in</span> <span style="color: #000000;">sp</span> <span style="color: #008080;">do</span>
+         <span style="color: #000000;">dc</span><span style="color: #0000FF;">[</span><span style="color: #000000;">c</span><span style="color: #0000FF;">-</span><span style="color: #008000;">'0'</span><span style="color: #0000FF;">+</span><span style="color: #000000;">1</span><span style="color: #0000FF;">]</span> <span style="color: #0000FF;">+=</span> <span style="color: #000000;">1</span>
+     <span style="color: #008080;">end</span> <span style="color: #008080;">for</span>
+     <span style="color: #004080;">integer</span> <span style="color: #000000;">dc9</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">dc</span><span style="color: #0000FF;">[</span><span style="color: #000000;">10</span><span style="color: #0000FF;">]</span>
+     <span style="color: #008080;">if</span> <span style="color: #7060A8;">odd</span><span style="color: #0000FF;">(</span><span style="color: #000000;">dc9</span><span style="color: #0000FF;">)</span> <span style="color: #008080;">and</span> <span style="color: #7060A8;">remainder</span><span style="color: #0000FF;">(</span><span style="color: #000000;">dc9</span><span style="color: #0000FF;">,</span><span style="color: #000000;">5</span><span style="color: #0000FF;">)!=</span><span style="color: #000000;">0</span> <span style="color: #008080;">then</span>
+         <span style="color: #004080;">string</span> <span style="color: #000000;">swdps</span> <span style="color: #0000FF;">=</span> <span style="color: #7060A8;">join</span><span style="color: #0000FF;">(</span><span style="color: #000000;">dc</span><span style="color: #0000FF;">,</span><span style="color: #008000;">""</span><span style="color: #0000FF;">,</span><span style="color: #000000;">fmt</span><span style="color: #0000FF;">:=</span><span style="color: #008000;">"%d"</span><span style="color: #0000FF;">)</span>
+         <span style="color: #7060A8;">mpz_set_str</span><span style="color: #0000FF;">(</span><span style="color: #000000;">dn</span><span style="color: #0000FF;">,</span><span style="color: #000000;">swdps</span><span style="color: #0000FF;">,</span><span style="color: #000000;">10</span><span style="color: #0000FF;">)</span>
+         <span style="color: #008080;">if</span> <span style="color: #7060A8;">mpz_prime</span><span style="color: #0000FF;">(</span><span style="color: #000000;">dn</span><span style="color: #0000FF;">)</span> <span style="color: #008080;">then</span>
+             <span style="color: #000000;">swdp</span> <span style="color: #0000FF;">=</span> <span style="color: #7060A8;">append</span><span style="color: #0000FF;">(</span><span style="color: #000000;">swdp</span><span style="color: #0000FF;">,{</span><span style="color: #000000;">np</span><span style="color: #0000FF;">,</span><span style="color: #000000;">swdps</span><span style="color: #0000FF;">})</span>
+         <span style="color: #008080;">end</span> <span style="color: #008080;">if</span>
+     <span style="color: #008080;">end</span> <span style="color: #008080;">if</span>
+     <span style="color: #008080;">if</span> <span style="color: #7060A8;">platform</span><span style="color: #0000FF;">()!=</span><span style="color: #004600;">JS</span> <span style="color: #008080;">and</span> <span style="color: #7060A8;">time</span><span style="color: #0000FF;">()></span><span style="color: #000000;">t1</span> <span style="color: #008080;">then</span>
+         <span style="color: #7060A8;">progress</span><span style="color: #0000FF;">(</span><span style="color: #008000;">"processing prime %d...\r"</span><span style="color: #0000FF;">,{</span><span style="color: #000000;">np</span><span style="color: #0000FF;">})</span>
+         <span style="color: #000000;">t1</span> <span style="color: #0000FF;">=</span> <span style="color: #7060A8;">time</span><span style="color: #0000FF;">()+</span><span style="color: #000000;">1</span>
+     <span style="color: #008080;">end</span> <span style="color: #008080;">if</span>
+ <span style="color: #008080;">end</span> <span style="color: #008080;">while</span>
+ <span style="color: #7060A8;">printf</span><span style="color: #0000FF;">(</span><span style="color: #000000;">1</span><span style="color: #0000FF;">,</span><span style="color: #008000;">"Smarandache-Whellen primes:\n"</span><span style="color: #0000FF;">)</span>
+ <span style="color: #008080;">for</span> <span style="color: #000000;">i</span><span style="color: #0000FF;">,</span><span style="color: #000000;">s</span> <span style="color: #008080;">in</span> <span style="color: #000000;">swp</span> <span style="color: #008080;">do</span>
+     <span style="color: #7060A8;">printf</span><span style="color: #0000FF;">(</span><span style="color: #000000;">1</span><span style="color: #0000FF;">,</span><span style="color: #008000;">" %d%s:"</span><span style="color: #0000FF;">,{</span><span style="color: #000000;">i</span><span style="color: #0000FF;">,</span><span style="color: #7060A8;">ord</span><span style="color: #0000FF;">(</span><span style="color: #000000;">i</span><span style="color: #0000FF;">)})</span>
+     <span style="color: #7060A8;">printf</span><span style="color: #0000FF;">(</span><span style="color: #000000;">1</span><span style="color: #0000FF;">,</span><span style="color: #008000;">" Index: %4d, Last prime %5d, %s\n"</span><span style="color: #0000FF;">,</span><span style="color: #000000;">s</span><span style="color: #0000FF;">)</span>
+ <span style="color: #008080;">end</span> <span style="color: #008080;">for</span>
+ <span style="color: #7060A8;">printf</span><span style="color: #0000FF;">(</span><span style="color: #000000;">1</span><span style="color: #0000FF;">,</span><span style="color: #008000;">"Smarandache-Whellen derived primes:\n"</span><span style="color: #0000FF;">)</span>
+ <span style="color: #008080;">for</span> <span style="color: #000000;">i</span><span style="color: #0000FF;">,</span><span style="color: #000000;">s</span> <span style="color: #008080;">in</span> <span style="color: #000000;">swdp</span> <span style="color: #008080;">do</span>
+     <span style="color: #7060A8;">printf</span><span style="color: #0000FF;">(</span><span style="color: #000000;">1</span><span style="color: #0000FF;">,</span><span style="color: #008000;">" %d%s:"</span><span style="color: #0000FF;">,{</span><span style="color: #000000;">i</span><span style="color: #0000FF;">,</span><span style="color: #7060A8;">ord</span><span style="color: #0000FF;">(</span><span style="color: #000000;">i</span><span style="color: #0000FF;">)})</span>
+     <span style="color: #7060A8;">printf</span><span style="color: #0000FF;">(</span><span style="color: #000000;">1</span><span style="color: #0000FF;">,</span><span style="color: #008000;">" Index: %4d, %s\n"</span><span style="color: #0000FF;">,</span><span style="color: #000000;">s</span><span style="color: #0000FF;">)</span>
+ <span style="color: #008080;">end</span> <span style="color: #008080;">for</span>
+ <span style="color: #0000FF;">?</span><span style="color: #7060A8;">elapsed</span><span style="color: #0000FF;">(</span><span style="color: #7060A8;">time</span><span style="color: #0000FF;">()-</span><span style="color: #000000;">t0</span><span style="color: #0000FF;">)</span>
+<!--</syntaxhighlight>-->
+{{out}}
+<pre>
+Smarandache-Whellen primes:
+st: Index:    1, Last prime     2, 2
+nd: Index:    2, Last prime     3, 23
+rd: Index:    4, Last prime     7, 2357
+th: Index:  128, Last prime   719, 23571113171923293137...73677683691701709719 (355 digits)
+th: Index:  174, Last prime  1033, 23571113171923293137...10131019102110311033 (499 digits)
+th: Index:  342, Last prime  2297, 23571113171923293137...22732281228722932297 (1,171 digits)
+th: Index:  435, Last prime  3037, 23571113171923293137...30013011301930233037 (1,543 digits)
+th: Index: 1429, Last prime 11927, 23571113171923293137...11903119091192311927 (5,719 digits)
+Smarandache-Whellen derived primes:
+st: Index:   32, 4194123321127
+nd: Index:   72, 547233879626521
+rd: Index:   73, 547233979727521
+th: Index:  134, 13672766322929571043
+th: Index:  225, 3916856106393739943689
+th: Index:  303, 462696313560586013558131
+th: Index:  309, 532727113760586013758133
+th: Index:  363, 6430314317473636515467149
+th: Index:  462, 8734722823685889120488197
+th: Index:  490, 9035923128899919621189209
+th: Index:  495, 9036023329699969621389211
+th: Index:  522, 9337023533410210710923191219
+th: Index:  538, 94374237357103109113243102223
+th: Index:  624, 117416265406198131121272110263
+th: Index:  721, 141459282456260193137317129313
+th: Index:  738, 144466284461264224139325131317
+th: Index:  790, 156483290479273277162351153339
+th: Index:  852, 164518312512286294233375158359
+th: Index: 1087, 208614364610327343341589284471
+th: Index: 1188, 229667386663354357356628334581
+"13 minutes and 13s"
 </pre>