aoc2019_22_part2_try2.py

N=119315717514047
P=101741582076661
def G(a,b):
 if a<1:return(b,0,1)
 g,x,y=G(b%a,a);return(g,y-(b/a)*x,x)
I=lambda x:G(x,N)[1]%N
m,a,s,c=1,0,1,0
for l in open("input.txt"):
 w=l.split()
 if'c'==l[0]:a+=m*int(w[-1])
 elif'w'==l[5]:m*=I(int(w[-1]))
 else:m=-m;a+=m
while P:
 if P&1:s,c=m*s,m*c+a
 P>>=1;m,a=m*m%N,m*a+a
print(2020*s+c)%N