HP 39g-Grafenberechner - Handbuch (Seite 159)

Exakte Berechnungen und Mathematik mit HP40G
Arithmetische Programme 159
Man sucht, wenn 2 N teilen kann, dann testet man den ungeraden Teiler D
zwischen 3 und E(
ÖN).
In der Reihe FACT folgt nach jedem Teiler sein Exponent (Potenzexponent):
Decomp(12)=
{2,2,3.,1}.
Funktion facprem(N)
lokal K D FACT
{}->FACT
0 -> K
solange N mod 2 = 0 durchführen
K+1 -> K
N/2 -> N
/solange/
wenn K
¹ 0 so
FACT +
{2 K} -> FACT
/wenn/
3-> D
solange D*D
£ N durchführen
0 -> K
solange N mod D = 0 durchführen
K+1 -> K
N/D -> N
/solange/
wenn K
¹ 0 so
FACT +
{D K} -> FACT
/wenn/
D+2 -> D
/solange/
wenn N
¹ 1 so