Добавлена операция возведения в степень по модулю

1 files changed, 21 insertions, 0 deletions

diff --git a/sem2/src/mpn.rs b/sem2/src/mpn.rs
index 97ba55d..991e813 100644
--- a/sem2/src/mpn.rs
+++ b/sem2/src/mpn.rs
@@ -343,6 +343,27 @@ impl Rem for Number {
     }
 }
 
+impl Number {
+    pub fn pow_mod(&self, power: &Self, modulo: &Self) -> Result<Self, String> {
+        if modulo.is_zero() {
+            return Err("Модуль должен быть отличен от нуля".to_string());
+        }
+        let one: Number = Number::from_digits(&vec![1], self.radix).unwrap();
+        let two: Number = Number::from_digits(&vec![2], self.radix).unwrap();
+        let mut acc = one.clone();
+        let mut p = power.clone();
+        let mut base = self.clone();
+        while !p.is_zero() {
+            if p.clone() % two.clone() == one.clone() {
+                acc = (acc * base.clone()) % modulo.clone();
+            }
+            p = p / two.clone();
+            base = (base.clone() * base) % modulo.clone();
+        }
+        return Ok(acc);
+    }
+}
+
 impl PartialEq for Number {
     fn eq(&self, other: &Self) -> bool {
         if self.len() != other.len() {