Maxima Function
ifactors (n)
For a positive integer n returns the factorization of n. If
n=p1^e1..pk^nk
is the decomposition of n into prime
factors, ifactors returns [[p1, e1], ... , [pk, ek]]
.
Factorization methods used are trial divisions by primes up to 9973, Pollard's rho method and elliptic curve method.
(%i1) ifactors(51575319651600); (%o1) [[2, 4], [3, 2], [5, 2], [1583, 1], [9050207, 1]] (%i2) apply("*", map(lambda([u], u[1]^u[2]), %)); (%o2) 51575319651600