Achilles numbers: Difference between revisions
→{{header|RPL}}: Faster version
m (→{{header|RPL}}: add a few comments) |
(→{{header|RPL}}: Faster version) |
||
Line 2,779:
Based on Wikipedia definition: n = p<sub>1</sub><sup>a<sub>1</sub></sup>p<sub>2</sub><sup>a<sub>2</sub></sup>…p<sub>k</sub><sup>a<sub>k</sub></sup> is an Achilles number if min(a<sub>1</sub>, a<sub>2</sub>, …, a<sub>k</sub>) ≥ 2 and gcd(a<sub>1</sub>, a<sub>2</sub>, …, a<sub>k</sub>) = 1.
{{works with|HP|49g}}
≪
1 s SIZE '''FOR''' j▼
SWAP s j GET + '''NEXT''' ▼
▲ '''IF''' DUP SIZE 4 < '''THEN''' DROP2 0
'''ELSE'''
a ≪ MIN ≫ STREAM 2 ≥▼
→ a
≫
'''END'''
Line 2,805 ⟶ 2,799:
≪ { } 1 → n
≪ '''WHILE''' DUP SIZE 20 < '''REPEAT'''
'''IF''' 'n' INCR <span style="color:blue">ACH?</span>
'''THEN''' IF n EULER <span style="color:blue">ACH?</span> '''THEN''' n + '''END END'''
'''END'''
≫ ≫ EVAL
Line 2,813 ⟶ 2,808:
1: {500 864 1944 2000 2592 3456 5000 10125 10368 12348 12500 16875 19652 19773 30375 31104 32000 33275 37044 40500}
</pre>
First 50 Achilles numbers found in
=={{header|Rust}}==
|