In [3]:
# Step 1
y = 2001

# Step 2
a = y % 19

# Step 3
b = y // 100
c = y % 10

# Step 4
d = b // 4
e = b % 4

# Step 5
g = (8 * b + 13) // 25

# Step 6
h = (19 * a + b - d - g + 15) % 30

# Step 7
j = c // 4
k = c % 4

# Step 8 
m = (a + 11 * h) // 319

# Step 9
r = (2 * e + 2 * j - k - h + m + 32) % 7

# Step 10
n = (h - m + r + 90) // 25

# Step 11
p = (h - m + r + n + 19) % 32

print(p)
 
15
In [ ]: