#include #include #include using namespace std; long op(long a,long b){return min(a,b);} long e(){return(long)1e18;} long mp(long f,long x){return f+x;} long id(){return 0L;} int N,K; long L,U,D; long A[2<<17]; struct dat{ long t; long b,sgn; long l,r; bool operator<(const dat&rhs)const{return t>N>>K>>L>>U; D=U-L; vectorOP; OP.reserve(2*N); for(int i=0;i>A[i]; A[i]+=D; if(A[i]>D) {//pos long a=A[i]-D; long tl=0,tr=(a+K-1)/K; long b=a/K,sgn=-1; long al=0,ar=a%K; OP.push_back((dat){tl,b,sgn,al,ar}); OP.push_back((dat){tr,b+sgn*(tr-tl),sgn*2,al,ar}); } {//neg long tl=A[i]/K,tr=2e12; long b=0,sgn=1; long al=A[i]%K+1,ar=K; OP.push_back((dat){tl,b,sgn,al,ar}); OP.push_back((dat){tr,b+sgn*(tr-tl),sgn*2,al,ar}); } } atcoder::lazy_segtreeseg(vector(K,0L)); sort(OP.begin(),OP.end()); long prvt=0; long B=0,SGN=0; long ans=1e18; for(dat d:OP) { if(prvt=0)ans=min(ans,B+seg.all_prod()); else ans=min(ans,B+SGN*(d.t-prvt-1)+seg.all_prod()); B+=SGN*(d.t-prvt); prvt=d.t; } if(abs(d.sgn)==1) { seg.apply(d.l,d.r,1L); B+=d.b; SGN+=d.sgn; } else { seg.apply(d.l,d.r,-1L); B-=d.b; SGN-=d.sgn/2; } } cout<