Product of primes not dividing $N$ in SAGE.

asked 2017-07-07 00:25:32 +0200

Manoj
5 ●1 ●1 ●3

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.

edit retag flag offensive close merge delete

add a comment

1 Answer

Sort by » oldest newest most voted

answered 2017-07-07 03:10:31 +0200

kcrisman
12222 ●41 ●135 ●254

This is a bit ill-defined - the infinitely many such primes? But here is something you can probably modify.

sage: N = 123
sage: prod(1/p for p in set(primes(N)).difference(set(N.prime_factors())))
1/256992314149736648684107368223932215392477130

edit flag offensive delete link

Comments

Thank you ver much.

Manoj ( 2017-07-07 22:17:48 +0200 )edit

add a comment

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Question Tools

Stats

Asked: 2017-07-07 00:25:32 +0200

Seen: 440 times

Last updated: Jul 07 '17

Product of primes not dividing $N$ in SAGE. edit