In [13]:
def main():
    y =int(input("Enter year==> "))
     #1. Let y be the year (such as 1800 or 2001). (This would be your input)   
        
    a = y % 19
    #2. Divide y by 19 and call the remainder a . Ignore the quotient.
    
    b = y // 100
    c = y % 100
      #3. Divide y by 100 to get a quotient b and a remainder c .
        
    d = b // 4 
    e = b % 4
    #4. Divide b by 4 to get a quotient d and a remainder e .
    
    g = (8 * b + 13) // 25
    #5. Divide 8 * b + 13 by 25 to get a quotient g . Ignore the remainder.
    
    h = (19 * a + b - g + 15) % 30
    #6. Divide 19 * a + b - d - g + 15 by 30 to get a remainder h . Ignore the quotient.
    
    j = c // 4
    k = c % 4
    #7. Divide c by 4 to get a quotient j and a remainder k .

    m = (a + 11 * h) // 319
    #8. Divide a + 11 * h by 319 to get a quotient m . Ignore the remainder.

    r = (2 * e + 2 * j - k - h + m + 32) % 7
    #9. Divide 2 * e + 2 * j - k - h + m + 32 by 7 to get a remainder r . Ignore the quotient.

    n = (h - m + r + 90) // 25
    #10. Divide h - m + r + 90 by 25 to get a quotient n . Ignore the remainder.

    p = (h - m + r + n + 19) % 32
    #11. Divide h - m + r + n + 19 by 32 to get a remainder p . Ignore the quotient
     
    print(n,p)
    print(a,b,c,d,e,g,h,j,m,r)
    
main()
Enter year==> 2021
4 4
7 20 21 5 0 6 12 5 0 1
In [ ]: