R=range def rec_str(a):return"".join(["[",", ".join(rec_str(x)for x in a),"]"])if isinstance(a,list)else str(a) class BIT: def __init__(self,x): if isinstance(x,int): self.N=x self.F=[0]*(x+1) else: self.N=len(x) self.F=[0]*(self.N+1) for j in R(1,self.N+1): i=j-1 self.F[j]=x[i] k=j-(j&-j) while i>k: self.F[j]+=self.F[i] i-=i&-i def copy(self): a=__class__([]) a.N=__self__.N a.F=__self__.F[:] return a def Add(self,i,u): if i<0:return i+=1 while i<=self.N: self.F[i]+=u i+=i&-i def Set(self,i,u):self.Add(i,u-self.Get(i)) def Get(self,i): assert 0<=i=u or j==N j=s=n=0 p=1<<17 #131072 while p: k=j|p p>>=1 if k<=self.N: n+=self.F[k] if nr:return self.F.Add(l,-(l-1)*u) self.F.Add(r+1,r*u) self.G.Add(l,u) self.G.Add(r+1,-u) def Add(self,i,u):self.IntervalAdd(i,i,u) def Set(self,i,u):self.Add(i,u-self.Get(i)) def InitialSegmentSum(self,r): assert -2>(L*L//(N-L+1)): C=[[I]*(L+1)for i in R(L+1)] def G(i,j):return C[i][j]if 0<=i<=L and 0<=j<=L else I C[0][0]=0 for(v,b)in sorted([(v,0)for v in A]+[(v,1)for v in B]): D=[c[:]for c in C] if b: for i in R(L+1): for j in R(L+1):D[i][j]=G(i,j-1)-v if ij else min(G(i,j),G(i,j-1)+v)) else: for i in R(L+1): for j in R(L+1):D[i][j]=G(i-1,j)+v if ij else min(G(i,j),G(i-1,j)+v)) C=D print(G(L,L)) else: C=[IntervalAddBIT([I]*(M-L+1))for i in R(N-L+1)] def G(i,j):return C[i].Get(j)if 0<=i<=N-L and 0<=j<=M-L else I def S(i,j,v): if 0<=i<=N-L and 0<=j<=M-L:C[i].Set(j,v) S(0,0,0) c=d=0 for(v,b)in sorted([(v,0)for v in A]+[(v,1)for v in B]): if b:d+=1 else:c+=1 for i in R(N-L,-1,-1): D=i-c+d if b: l=G(i,D-1) C[i].IntervalAdd(0,D-1,-v) C[i].IntervalAdd(D,M-L,v) S(i,D,min(l,G(i,D))) else: l=G(i-1,D) C[i].IntervalAdd(0,D,v) C[i].IntervalAdd(D+1,M-L,-v) S(i,D,min(l,G(i,D))) print(G(N-L,M-L))