gaussiant 0.8.0

Gaussian integers
Documentation 
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

//! See <https://en.wikipedia.org/wiki/Modular_arithmetic#Properties>.

use gaussiant::{gaussint, GaussianInt};

fn main() {
    let a = gaussint!(2, -1);
    let b = gaussint!(2, 3);
    let n = gaussint!(1, 1);

    let a1 = a;
    let b1 = b;
    let a2 = gaussint!(4, 1);
    let b2 = gaussint!(4, -3);

    // a ≡ b (mod n)
    assert!(a.congruent(b, n));
    // a₂ ≡ b₂ (mod n)
    assert!(a2.congruent(b2, n));

    // The congruence relation satisfies all the conditions of an equivalence relation:

    // a ≡ a (mod n) (reflexivity)
    assert!(a.congruent(a, n));

    // a ≡ b (mod n) if b ≡ a (mod n) for all a, b, and n. (symmetry)
    assert_eq!(a.congruent(b, n), b.congruent(a, n));

    // a ≡ b (mod n) and b ≡ c (mod n) ⇒ a ≡ c (mod n) (transitivity)
    let c = gaussint!(5, 4);
    assert!(b.congruent(c, n));
    assert!(a.congruent(c, n));

    // If a₁ ≡ b₁ (mod n) and a₂ ≡ b₂ (mod n), or if a ≡ b (mod n), then:

    // a₁ + a₂ ≡ b₁ + b₂ (mod n) (compatibility with addition)
    assert!((a1 + a2).congruent(b1 + b2, n));

    // a₁ – a₂ ≡ b₁ – b₂ (mod n) (compatibility with subtraction)
    assert!((a1 - a2).congruent(b1 - b2, n));

    // a + k ≡ b + k (mod n) for any integer k (compatibility with translation)
    let k = gaussint!(10);

    assert!((a + k).congruent(b + k, n));
    // k a ≡ k b (mod n) for any integer k (compatibility with scaling)
    assert!((a * k).congruent(b * k, n));

    // a₁ + a₂ ≡ b₁ + b₂ (mod n) (compatibility with addition)
    assert!((a1 + a2).congruent(b1 + b2, n));

    // a₁ – a₂ ≡ b₁ – b₂ (mod n) (compatibility with subtraction)
    assert!((a1 + a2).congruent(b1 + b2, n));

    // a₁ a₂ ≡ b₁ b₂ (mod n) (compatibility with multiplication)
    assert!((a1 * a2).congruent(b1 * b2, n));

    // a^k ≡ b^k (mod n) for any non-negative integer k (compatibility with exponentiation)
    // TODO

    // p(a) ≡ p(b) (mod n), for any polynomial p(x) with integer coefficients
    // (compatibility with polynomial evaluation)
    // TODO?
}