diff --git a/python/pdtd.py b/python/pdtd.py
new file mode 100644
index 0000000..9de6964
--- /dev/null
+++ b/python/pdtd.py
@@ -0,0 +1,11 @@
+for num in range(1, 5001):
+  total = 0
+  rem = num
+
+  while rem > 0:
+    digit = rem % 10
+    total += pow(digit, digit)
+    rem //= 10
+
+  if total == num:
+    print(total)