/* Predicate function that checks k-almost primality for given integer n and parameter k */ k_almost_primep(n,k):=if integerp((n)^(1/k)) and primep((n)^(1/k)) then true else lambda([x],(length(ifactors(x))=k and unique(map(second,ifactors(x)))=[1]) or (length(ifactors(x))