Jacobsthal numbers: Difference between revisions

=={{header|Quackery}}==

<code>isprime</code> is defined at [[Primality by trial division#Quackery]].

<syntaxhighlight lang="Quackery"> [ 2 over ** -1 rot ** - 3 / ] is j ( n --> n )

[ 2 over ** -1 rot ** + ] is jl ( n --> n )

[ dup 1+ j swap j * ] is jo ( n --> n )

say "First 30 Jacobsthal numbers:"

cr

30 times [ i^ j echo sp ]

cr cr

say "First 30 Jacobsthal-Lucas numbers:"

cr

30 times [ i^ jl echo sp ]

cr cr

say "First 20 Jacobsthal oblong numbers:"

cr

20 times [ i^ jo echo sp ]

cr cr

say "First 10 Jacobsthal primes:"

cr

[] 0

[ dup j dup isprime iff

[ swap dip join ]

else drop

1+

over size 10 = until ]

drop

witheach [ echo sp ]</syntaxhighlight>

{{out}}

<pre>First 30 Jacobsthal numbers:

0 1 1 3 5 11 21 43 85 171 341 683 1365 2731 5461 10923 21845 43691 87381 174763 349525 699051 1398101 2796203 5592405 11184811 22369621 44739243 89478485 178956971

First 30 Jacobsthal-Lucas numbers:

2 1 5 7 17 31 65 127 257 511 1025 2047 4097 8191 16385 32767 65537 131071 262145 524287 1048577 2097151 4194305 8388607 16777217 33554431 67108865 134217727 268435457 536870911

First 20 Jacobsthal oblong numbers:

0 1 3 15 55 231 903 3655 14535 58311 232903 932295 3727815 14913991 59650503 238612935 954429895 3817763271 15270965703 61084037575

First 10 Jacobsthal primes:

3 5 11 43 683 2731 43691 174763 2796203 715827883

</pre>