# why (4*i+7) in Z[i] not a prime?norm(4*i+7) =65, 65mod4=1

```
K.<a> = NumberField(x^2 +1);
(4*a+7).norm()
65
mod (65,4)
1
factor(4*a+7)
(-3*a - 2) * (a - 2)
(4*a+7).multiplicative_order();
+Infinity
K.elements_of_norm(65);factor(-8*a + 1);K.prime_factors(4*a -7);
-8*a - 1, 4*a - 7, 4*a + 7, -8*a + 1]
(a - 2) * (3*a - 2)
[Fractional ideal (3*a - 2), Fractional ideal (-a - 2)]
[Fractional ideal (a - 2), Fractional ideal (3*a - 2)]
```

add a comment