AKS test for primes: Difference between revisions
add RPL: improved display of (x-1)^p
(add RPL) |
(add RPL: improved display of (x-1)^p) |
||
Line 4,770:
≫ ≫ ‘'''XPM1'''’ STO
≪ DUP SIZE → coeffs size
≪ 0 1 size '''FOR''' j
coeffs j GET 'X' size j - ^ * + '''NEXT'''
≫ ≫ ‘'''SHOWA'''’ STO
≪ DUP '''XPM1''' → p coeffs
≪ 1
Line 4,783 ⟶ 4,788:
end loop
return coeffs
'''SHOWA''' ''( { coeffs } -- 'polynom' )''
loop for each coeff
append jth polynomial term
return polynom
'''PRIM?''' ''( p -- boolean ) ''
Line 4,792 ⟶ 4,802:
{{in}}
<pre>
≪ { } 1 7 FOR n n XPM1 SHOWA + NEXT ≫
≪ { } 2 35 FOR n IF n PRIM? THEN n + END NEXT ≫
</pre>
{{out}}
<pre>
8:
7:
6:
5:
4:
3:
2:
1: { 2 3 5 7 11 13 17 19 23 29 31 }
</pre>
| |||