Jump to content

Wieferich primes: Difference between revisions

m
m (→‎{{header|Julia}}: simpler version)
m (→‎{{header|Phix}}: simplified)
Line 42:
<!--<lang Phix>(notonline)-->
<span style="color: #008080;">include</span> <span style="color: #7060A8;">mpfr</span><span style="color: #0000FF;">.</span><span style="color: #000000;">e</span>
<span style="color: #7060A8;">mpz_sub_uimpz</span><span style="color: #0000FF;">(</span><span style="color: #000000;">p2pm1m1z</span> <span style="color: #0000FF;">,</span><span style="color: #000000;">p2pm1m1</span><span style="color: #0000FF;">,</span><span style="color: #0000007060A8;">1mpz_init</span><span style="color: #0000FF;">()</span>
<span style="color: #008080;">function</span> <span style="color: #000000;">weiferich</span><span style="color: #0000FF;">(</span><span style="color: #004080;">integer</span> <span style="color: #000000;">p</span><span style="color: #0000FF;">)</span>
<span style="color: #7060A8;">mpzmpz_set_str</span><span style="color: #0000FF;">(</span><span style="color: #000000;">p2pm1m1z</span><span style="color: #0000FF;">,</span><span style="color: #7060A8;">repeat</span><span style="color: #0000FF;">(</span><span style="color: #008000;">'1'</span><span style="color: #0000FF;">,</span><span style="color: #7060A8000000;">mpz_initp</span><span style="color: #0000FF;">-</span><span style="color: #000000;">1</span><span style="color: #0000FF;">),</span><span style="color: #000000;">2</span><span style="color: #0000FF;">()</span>
<span style="color: #008080;">return</span> <span style="color: #7060A8;">mpz_ui_pow_uimpz_fdiv_q_ui</span><span style="color: #0000FF;">(</span><span style="color: #000000;">p2pm1m1z</span><span style="color: #0000FF;">,</span><span style="color: #000000;">2z</span><span style="color: #0000FF;">,</span><span style="color: #000000;">p</span><span style="color: #0000FF;">-*</span><span style="color: #000000;">1p</span><span style="color: #0000FF;">)=</span><span style="color: #000000;">0</span>
<span style="color: #7060A8;">mpz_sub_ui</span><span style="color: #0000FF;">(</span><span style="color: #000000;">p2pm1m1</span><span style="color: #0000FF;">,</span><span style="color: #000000;">p2pm1m1</span><span style="color: #0000FF;">,</span><span style="color: #000000;">1</span><span style="color: #0000FF;">)</span>
<span style="color: #008080;">return</span> <span style="color: #7060A8;">mpz_fdiv_q_ui</span><span style="color: #0000FF;">(</span><span style="color: #000000;">p2pm1m1</span><span style="color: #0000FF;">,</span><span style="color: #000000;">p2pm1m1</span><span style="color: #0000FF;">,</span><span style="color: #000000;">p</span><span style="color: #0000FF;">*</span><span style="color: #000000;">p</span><span style="color: #0000FF;">)=</span><span style="color: #000000;">0</span>
<span style="color: #008080;">end</span> <span style="color: #008080;">function</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;">"Weiferich primes less than 5000: %V\n"</span><span style="color: #0000FF;">,{</span><span style="color: #7060A8;">filter</span><span style="color: #0000FF;">(</span><span style="color: #7060A8;">get_primes_le</span><span style="color: #0000FF;">(</span><span style="color: #000000;">5000</span><span style="color: #0000FF;">),</span><span style="color: #000000;">weiferich</span><span style="color: #0000FF;">)})</span>
7,820

edits

Cookies help us deliver our services. By using our services, you agree to our use of cookies.