From c113d0f2d8a2326ec4a6ee29eef5586e8536974c Mon Sep 17 00:00:00 2001
From: VicRen <vic.ren@me.com>
Date: Fri, 20 Nov 2020 10:05:28 +0800
Subject: [PATCH] 48: self powers

---
 48_self_powers/euler48.py | 17 +++++++++++++++++
 1 file changed, 17 insertions(+)
 create mode 100644 48_self_powers/euler48.py

diff --git a/48_self_powers/euler48.py b/48_self_powers/euler48.py
new file mode 100644
index 0000000..7a95cb7
--- /dev/null
+++ b/48_self_powers/euler48.py
@@ -0,0 +1,17 @@
+def modpow(b, e, m):
+  r = 1
+  while e > 0:
+    if e & 1:
+      r = (r * b) % m
+    b = (b * b) % m
+    e = e >> 1
+  return r
+
+def powersum(n, m):
+  s = 0
+  for i in range(1, n + 1):
+    s+= modpow(i, i, m)
+    s%= m
+  return s
+
+print powersum(1000, 10000000000)
\ No newline at end of file