Giuga numbers: Difference between revisions
Content added Content deleted
(add RPL) |
imported>Maxima enthusiast No edit summary |
||
Line 881: | Line 881: | ||
24423128562 (elapsed: 432.06500005722046) |
24423128562 (elapsed: 432.06500005722046) |
||
432.066249 seconds (235 allocations: 12.523 KiB) |
432.066249 seconds (235 allocations: 12.523 KiB) |
||
</pre> |
|||
=={{header|Maxima}}== |
|||
<syntaxhighlight lang="maxima"> |
|||
/* Predicate function that checks wether an integer is a Giuga number or not */ |
|||
giugap(n):=if not primep(n) then block(ifactors(n),map(lambda([x],mod((n/x)-1,x)=0),map(first,%%)), |
|||
if length(unique(%%))=1 and apply(lhs,unique(%%))=0 then true)$ |
|||
/* Function that returns a list of the first len Giuga integers */ |
|||
giuga_count(len):=block( |
|||
[i:1,count:0,result:[]], |
|||
while count<len do (if giugap(i) then (result:endcons(i,result),count:count+1),i:i+1), |
|||
result)$ |
|||
/* Test case */ |
|||
giuga_count(4); |
|||
</syntaxhighlight> |
|||
{{out}} |
|||
<pre> |
|||
[30,858,1722,66198] |
|||
</pre> |
</pre> |
||