Euclid-Mullin sequence: Difference between revisions
add fermat
(add parigp) |
(add fermat) |
||
Line 44:
23003
</pre>
=={{header|Fermat}}==
<lang fermat>Func Firstfac(n) =
j := 3;
up := Sqrt(n);
while j <= up do
if Divides(j,n) then Return(j) fi;
j:=j+2;
od;
Return(n).;
Array eu[16];
eu[1]:=2;
!(eu[1],' ');
for i=2 to 16 do
eu[i]:=Firstfac(1+Prod<k=1,i-1>[eu[k]]);
!(eu[i],' ');
od;</lang>
{{out}}<pre> 2 3 7 43 13 53 5 6221671 38709183810571 139 2801 11 17 5471 52662739 23003</pre>
=={{header|Julia}}==
| |||