1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
use crate::{Domain, EuclideanDomain, Field, IntegralDomain, UnitaryRing};
#[derive(Clone, Debug, Default)]
pub struct QuotientField<R: EuclideanDomain> {
base: R,
modulo: R::Elem,
}
impl<R: EuclideanDomain> QuotientField<R> {
pub fn new(base: R, modulo: R::Elem) -> Self {
assert!(base.contains(&modulo));
let one = base.one();
assert!(base.is_zero(&base.rem(&one, &one)));
QuotientField { base, modulo }
}
pub fn base(&self) -> &R {
&self.base
}
pub fn modulo(&self) -> &R::Elem {
&self.modulo
}
}
impl<R: EuclideanDomain> Domain for QuotientField<R> {
type Elem = R::Elem;
fn contains(&self, elem: &Self::Elem) -> bool {
self.base.is_reduced(elem, &self.modulo)
}
}
impl<R: EuclideanDomain> UnitaryRing for QuotientField<R> {
fn zero(&self) -> Self::Elem {
self.base.zero()
}
fn neg(&self, elem: &Self::Elem) -> Self::Elem {
self.base.rem(&self.base.neg(elem), &self.modulo)
}
fn add(&self, elem1: &Self::Elem, elem2: &Self::Elem) -> Self::Elem {
self.base.rem(&self.base.add(elem1, elem2), &self.modulo)
}
fn one(&self) -> Self::Elem {
self.base.one()
}
fn mul(&self, elem1: &Self::Elem, elem2: &Self::Elem) -> Self::Elem {
self.base.rem(&self.base.mul(elem1, elem2), &self.modulo)
}
}
impl<R: EuclideanDomain> IntegralDomain for QuotientField<R> {}
impl<R: EuclideanDomain> EuclideanDomain for QuotientField<R> {
fn quo_rem(&self, elem1: &Self::Elem, elem2: &Self::Elem) -> (Self::Elem, Self::Elem) {
if self.is_zero(elem2) {
(self.zero(), elem1.clone())
} else {
(self.div(elem1, elem2), self.zero())
}
}
fn associate_repr(&self, elem: &Self::Elem) -> (Self::Elem, Self::Elem) {
if self.is_zero(elem) {
(self.zero(), self.one())
} else {
(self.one(), self.div(&self.one(), elem))
}
}
}
impl<R: EuclideanDomain> Field for QuotientField<R> {
fn inv(&self, elem: &Self::Elem) -> Self::Elem {
let (g, _, mut r) = self.base.extended_gcd(&self.modulo, elem);
if !self.base.is_one(&g) {
let (a, b) = self.base.quo_rem(&self.base.one(), &g);
assert!(self.base.is_zero(&b), "modulo was not irreducible");
r = self.base.mul(&r, &a);
}
self.base.rem(&r, &self.modulo)
}
}
#[cfg(test)]
mod tests {
use super::*;
use crate::CheckedInts;
#[test]
fn field_1721() {
let ring: CheckedInts<i32> = Default::default();
let field = QuotientField::new(ring, 1721);
for a in -860..860 {
assert!(field.contains(&a));
if a != 0 {
let b = field.inv(&a);
assert!(field.contains(&b));
println!("{} {}", a, b);
assert!(field.is_one(&field.mul(&a, &b)));
}
}
}
}