Special neighbor primes: Difference between revisions
add fermat
m (blah) |
(add fermat) |
||
Line 139:
{ 67 71 }
{ 73 79 }
</pre>
=={{header|Fermat}}==
<lang fermat>Func Nextprime(p) =
q:=1;
while not Isprime(p+q)=1 do
q:=q + 1;
od;
p+q.;
for p1 = 3 to 99 by 2 do
p2:=Nextprime(p1);
if p2<100 and Isprime(p1)=1 and Isprime(p1+p2-1) then
!!(p1,' +',p2,' - 1 =',p1+p2-1);
fi;
od;</lang>
{{out}}<pre>
3 + 5 - 1 = 7
5 + 7 - 1 = 11
7 + 11 - 1 = 17
11 + 13 - 1 = 23
13 + 17 - 1 = 29
19 + 23 - 1 = 41
29 + 31 - 1 = 59
31 + 37 - 1 = 67
41 + 43 - 1 = 83
43 + 47 - 1 = 89
61 + 67 - 1 = 127
67 + 71 - 1 = 137
73 + 79 - 1 = 151
</pre>
| |||