diff --git a/ChaCha20_Poly1305_64/sim/do_poly_1305.py b/ChaCha20_Poly1305_64/sim/do_poly_1305.py
index 7a63286..dc0377e 100644
--- a/ChaCha20_Poly1305_64/sim/do_poly_1305.py
+++ b/ChaCha20_Poly1305_64/sim/do_poly_1305.py
@@ -4,13 +4,13 @@ def mask_r(r: int) -> int:
     r_bytes = r.to_bytes(16, "little")
 
     r_masked = bytearray(r_bytes)
-    r_masked[3] &= 15;
-    r_masked[7] &= 15;
-    r_masked[11] &= 15;
-    r_masked[15] &= 15;
-    r_masked[4] &= 252;
-    r_masked[8] &= 252;
-    r_masked[12] &= 252;
+    r_masked[3] &= 15
+    r_masked[7] &= 15
+    r_masked[11] &= 15
+    r_masked[15] &= 15
+    r_masked[4] &= 252
+    r_masked[8] &= 252
+    r_masked[12] &= 252
 
 
     r_masked = int.from_bytes(r_masked, "little")
@@ -39,6 +39,8 @@ def poly1305(message: bytes, r: int, s: int):
 def parallel_poly1305(message: bytes, r: int, s: int, lanes: int):
     r = mask_r(r)
     p = 2**130-5
+
+    r_powers = [r**i % p for i in range(lanes+1)]
     
     acc = [0]*lanes
 
@@ -54,7 +56,7 @@ def parallel_poly1305(message: bytes, r: int, s: int, lanes: int):
             byte_length = (lane.bit_length() + 7) // 8
             lane += 1 << (8*byte_length)
 
-            acc[j] = ((acc[j] + lane)*(r**power)) % p
+            acc[j] = ((acc[j] + lane)*(r_powers[power])) % p
 
     combined_acc = sum(acc) % p
     combined_acc += s