Motzkin numbers: Difference between revisions
add fermat
(add parigp) |
(add fermat) |
||
Line 411:
41 192,137,918,101,841,817
</pre>
=={{header|Fermat}}==
<lang fermat>
Array m[42];
m[1]:=1;
m[2]:=2;
!!(1,0); {precompute and print m[0] thru m[2]}
!!(1,0);
!!(2,1);
for n=3 to 41 do
m[n]:=(m[n-1]*(2*n+1) + m[n-2]*(3*n-3))/(n+2);
!!(m[n],Isprime(m[n]));
od;
; {built-in Isprime function returns 0 for 1, 1 for primes, and}
; {the smallest prime factor for composites, so this actually gives}
; {slightly more information than requested}
</lang>
=={{header|FreeBASIC}}==
| |||