Find prime numbers of the form n*n*n+2: Difference between revisions
Content added Content deleted
(Added Algol 68) |
(Added Algol W) |
||
Line 85: | Line 85: | ||
OD |
OD |
||
END</lang> |
END</lang> |
||
{{out}} |
|||
<pre> |
|||
1: 3 |
|||
3: 29 |
|||
5: 127 |
|||
29: 24391 |
|||
45: 91127 |
|||
63: 250049 |
|||
65: 274627 |
|||
69: 328511 |
|||
71: 357913 |
|||
83: 571789 |
|||
105: 1157627 |
|||
113: 1442899 |
|||
123: 1860869 |
|||
129: 2146691 |
|||
143: 2924209 |
|||
153: 3581579 |
|||
171: 5000213 |
|||
173: 5177719 |
|||
189: 6751271 |
|||
</pre> |
|||
=={{header|ALGOL W}}== |
|||
<lang algolw>begin % Find n such that n^3 + 2 is a prime for n < 200 % |
|||
for n := 1 until 199 do begin |
|||
integer candidate; |
|||
logical isPrime; |
|||
candidate := ( n * n * n ) + 2; |
|||
% there will only be 199 candidates, so a primality check by trial % |
|||
% division should be OK % |
|||
isPrime := true; |
|||
for f := 2 until entier( sqrt( candidate ) ) do begin |
|||
isPrime := candidate rem f not = 0; |
|||
if not isPrime then goto endPrimalityCheck |
|||
end for_f ; |
|||
endPrimalityCheck: |
|||
if isPrime then begin |
|||
% n^3 + 2 is prime % |
|||
write( i_w := 4, s_w := 0, n, ": ", i_w := 8, candidate ) |
|||
end if_isPrime |
|||
end for_n |
|||
end.</lang> |
|||
{{out}} |
{{out}} |
||
<pre> |
<pre> |