Anonymous user
Jacobsthal numbers: Difference between revisions
no edit summary
imported>Maxima enthusiast No edit summary |
|||
Line 2,191:
313 5562466239377370006237035693149875298444543026970449921737087520370363869220418099018130434731
347 95562442332919646317117537304253622533190207882011713489066201641121786503686867002917439712921903606443</pre>
=={{header|Maxima}}==
<syntaxhighlight lang="maxima">
jacobstahl(n):=(2^n-(-1)^n)/3$
jacobstahl_lucas(n):=2^n+(-1)^n$
jacobstahl_oblong(n):=jacobstahl(n)*jacobstahl(n+1)$
/* Function that returns a list of the first len jacobstahl primes */
jacobstahl_primes_count(len):=block(
[i:0,count:0,result:[]],
while count<len do (if primep(jacobstahl(i)) then (result:endcons(jacobstahl(i),result),count:count+1),i:i+1),
result)$
/* Test cases */
makelist(jacobstahl(i),i,0,29);
makelist(jacobstahl_lucas(i),i,0,29);
makelist(jacobstahl_oblong(i),i,0,19);
jacobstahl_primes_count(10);
</syntaxhighlight>
{{out}}
<pre>
[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]
[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]
[0,1,3,15,55,231,903,3655,14535,58311,232903,932295,3727815,14913991,59650503,238612935,954429895,3817763271,15270965703,61084037575]
[3,5,11,43,683,2731,43691,174763,2796203,715827883]
</pre>
=={{header|Nim}}==
| |||