module Data.Numbers
(
primeFactors,
numOfFactors,
factorSum,
factors) where
import Data.List
import Data.Numbers.Primes
primeFactors :: Integer -> [Integer]
primeFactors n = primeFactors' n n [] primes where
primeFactors' n m l (p:ps)
| p>n = l
| mo == 0 = primeFactors' n d (p:l) (p:ps)
| otherwise = primeFactors' n n l ps where
(d,mo) = divMod m p
numOfFactors :: Integer -> Int
numOfFactors n = product (map ((+1).length) (group (primeFactors n)))
factorSum :: Integer -> Integer
factorSum = product.(map (\p@(x:xs) -> (product p * x 1) `div` (x1))).(group.primeFactors)
factors :: Integer -> [Integer]
factors n = 1:factors' n [] [2..] where
factors' n l (m:ms)
| m*m==n= m:l
| m*m>n = l
| r == 0 = factors' n (d:m:l) ms
| otherwise = factors' n l ms where
(d,r) = divMod n m