Jacobsthal numbers: Difference between revisions
Content added Content deleted
imported>Maxima enthusiast No edit summary |
|||
Line 2,191: | Line 2,191: | ||
313 5562466239377370006237035693149875298444543026970449921737087520370363869220418099018130434731 |
313 5562466239377370006237035693149875298444543026970449921737087520370363869220418099018130434731 |
||
347 95562442332919646317117537304253622533190207882011713489066201641121786503686867002917439712921903606443</pre> |
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}}== |
=={{header|Nim}}== |