Hi, for a given $N,$ I want to compute a product of the type in SageMath
$$ \prod_{p\nmid N} \left(1-\frac{1}{p}\right)\left(1+\frac{1}{p}^{-1}\right) $$
For this I need to generate a list of these primes which does not divide $N.$ I do not know any method or command which does that. How can I compute this in sage? Thanks in advance.
