add GetSquarefree(poly f) in nforder.lib
use GetSquarefree in is_primitve