Input: 
n=integer  to be factored
B=smoothness bound

Output:
The product of the factors p of n, with p-1 B-smooth


pminus(n,B)=
d=gcd(lift(Mod(random,n)^(B!)-1),n);\
print(d,"*",n/d)