Quadrat special primes: Difference between revisions

Content added Content deleted
(→‎{{header|Ring}}: cleaned up columns, added additional degrees)
m (→‎{{header|Phix}}: added limit array, oct..dec)
Line 227:
=={{header|Phix}}==
<!--<lang Phix>-->
<span style="color: #004080008080;">atomconstant</span> <span style="color: #000000;">Ndesc</span> <span style="color: #0000FF;">=</span> <span style="color: #0000007060A8;">split</span><span style="color: #0000FF;">(</span><span style="color: #008000;">"linear quadratic cubic quartic quintic sextic septic octic nonic decic"</span><span style="color: #0000FF;">16000),</span>
<span style="color: #004080000000;">sequencelimits</span> <span style="color: #0000FF;">=</span> <span style="color: #0000FF;">{</span><span style="color: #000000;">desc1</span><span style="color: #0000FF;">,</span> <span style="color: #000000;">16000</span><span style="color: #0000FF;">,</span> <span style="color: #000000;">15000</span><span style="color: #0000FF;">,</span> <span style="color: #7060A8000000;">split14e9</span><span style="color: #0000FF;">(,</span> <span style="color: #008000000000;">8025e5</span><span style="linearcolor: quadratic#0000FF;">,</span> cubic<span quarticstyle="color: quintic#000000;">25e12</span><span sexticstyle="color: septic#0000FF;">,</span> <span style="color: #000000;">195e12</span><span style="color: #0000FF;">,</span><span style="color: #000000;">75e11</span><span style="color: #0000FF;">,</span> <span style="color: #000000;">3e9</span><span style="color: #0000FF;">,</span> <span style="color: #000000;">11e8</span><span style="color: #0000FF;">)}</span>
<span style="color: #008080;">for</span> <span style="color: #000000;">p</span><span style="color: #0000FF;">=</span><span style="color: #000000;">2</span> <span style="color: #008080;">to</span> <span style="color: #7060A8;">length</span><span style="color: #0000FF;">(</span><span style="color: #000000;">desc</span><span style="color: #0000FF;">)</span> <span style="color: #008080;">do</span>
<span style="color: #004080;">integeratom</span> <span style="color: #000000;">lastnN</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">limits</span><span style="color: #0000FF;">[</span><span style="color: #000000;">p</span><span style="color: #0000FF;">],</span> <span style="color: #000000;">lastn</span> <span style="color: #0000FF;">=</span> <span style="color: #7060A8;">ceil</span><span style="color: #0000FF;">(</span><span style="color: #7060A8;">power</span><span style="color: #0000FF;">(</span><span style="color: #000000;">N</span><span style="color: #0000FF;">,</span><span style="color: #000000;">1</span><span style="color: #0000FF;">/</span><span style="color: #000000;">p</span><span style="color: #0000FF;">))</span>
<span style="color: #004080;">sequence</span> <span style="color: #000000;">res</span> <span style="color: #0000FF;">=</span> <span style="color: #0000FF;">{</span><span style="color: #000000;">2</span><span style="color: #0000FF;">}</span>
<span style="color: #004080;">bool</span> <span style="color: #000000;">done</span> <span style="color: #0000FF;">=</span> <span style="color: #004600;">false</span>
Line 246:
<span style="color: #008080;">end</span> <span style="color: #008080;">while</span>
<span style="color: #004080;">string</span> <span style="color: #000000;">r</span> <span style="color: #0000FF;">=</span> <span style="color: #7060A8;">join_by</span><span style="color: #0000FF;">(</span><span style="color: #7060A8;">apply</span><span style="color: #0000FF;">(</span><span style="color: #004600;">true</span><span style="color: #0000FF;">,</span><span style="color: #7060A8;">sprintf</span><span style="color: #0000FF;">,{{</span><span style="color: #008000;">"%,6d"</span><span style="color: #0000FF;">},</span><span style="color: #000000;">res</span><span style="color: #0000FF;">}),</span><span style="color: #000000;">1</span><span style="color: #0000FF;">,</span><span style="color: #000000;">p</span><span style="color: #0000FF;">+</span><span style="color: #000000;">5</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;">"Found %d %s special primes &lt; %,dg:\n%s\n"</span><span style="color: #0000FF;">,{</span><span style="color: #7060A8;">length</span><span style="color: #0000FF;">(</span><span style="color: #000000;">res</span><span style="color: #0000FF;">),</span><span style="color: #000000;">desc</span><span style="color: #0000FF;">[</span><span style="color: #000000;">p</span><span style="color: #0000FF;">],</span><span style="color: #000000;">N</span><span style="color: #0000FF;">,</span><span style="color: #000000;">r</span><span style="color: #0000FF;">})</span>
<span style="color: #008080;">end</span> <span style="color: #008080;">for</span>
<!--</lang>-->
{{out}}
<pre>
Found 49 quadratic special primes < 16,00016000:
2 3 7 11 47 83 227
263 587 911 947 983 1,019 1,163
Line 260:
12,791 13,367 13,691 14,591 14,627 14,771 15,671
 
Found 23 cubic special primes < 16,00015000:
2 3 11 19 83 1,811 2,027 2,243
2,251 2,467 2,531 2,539 3,539 3,547 4,547 5,059
10,891 12,619 13,619 13,627 13,691 13,907 14,419
 
Found 39 quartic special primes < 16,0001.4e+10:
2 3 19 160,019 1,049,920,019 1,050,730,019 1,051,540,019 12,910,750,019 13,960,510,019
 
Found 29 quintic special primes < 16,0008.025e+8:
2 3 32,771 32,803 33,827 41,603 579,427 778,179,427 802,479,427
2 3
 
Found 46 sextic special primes < 16,0002.5e+13:
2 3 67 131 2,176,782,467 22,485,250,805,891
 
Found 37 septic special primes < 16,0001.95e+14:
2 3 131 194,871,710,000,131 194,893,580,000,131 194,893,580,280,067 194,971,944,444,163
 
Found 4 octic special primes < 7.5e+12:
2 3 65,539 6,553,600,065,539
 
Found 6 nonic special primes < 3e+9:
2 3 262,147 10,339,843 20,417,539 1,020,417,539
 
Found 4 decic special primes < 1.1e+9:
2 3 1,073,741,827 1,073,742,851
</pre>
 
Retrieved from "https://rosettacode.org/wiki/Quadrat_special_primes"