Motzkin numbers: Difference between revisions

m
Line 85:
40 66,368,199,913,921,497 ❌
41 192,137,918,101,841,817 ❌
</pre>
 
=={{header|Julia}}==
<lang julia>using Primes
 
function motzkin(N)
m = zeros(Int, N)
m[1] = m[2] = 1
for i in 3:N
m[i] = (m[i - 1] * (2i - 1) + m[i - 2] * (3i - 6)) ÷ (i + 1)
end
return m
end
 
println(" n M[n] Prime?\n-----------------------------------")
for (i, m) in enumerate(motzkin(42))
println(lpad(i, 2), lpad(m, 20), lpad(isprime(m), 8))
end
</lang>{{out}}
<pre>
n M[n] Prime?
-----------------------------------
1 1 false
2 1 false
3 2 true
4 4 false
5 9 false
6 21 false
7 51 false
8 127 true
9 323 false
10 835 false
11 2188 false
12 5798 false
13 15511 true
14 41835 false
15 113634 false
16 310572 false
17 853467 false
18 2356779 false
19 6536382 false
20 18199284 false
21 50852019 false
22 142547559 false
23 400763223 false
24 1129760415 false
25 3192727797 false
26 9043402501 false
27 25669818476 false
28 73007772802 false
29 208023278209 false
30 593742784829 false
31 1697385471211 false
32 4859761676391 false
33 13933569346707 false
34 40002464776083 false
35 114988706524270 false
36 330931069469828 false
37 953467954114363 true
38 2750016719520991 false
39 7939655757745265 false
40 22944749046030949 false
41 66368199913921497 false
42 192137918101841817 false
</pre>
 
4,108

edits