WC2018集训 吉老师的军训练

WC2018集训 吉老师的军训练

#include<bits/stdc++.h>
#define RG register
#define IL inline
#define _ 200005
#define X 100000000
#define ll unsigned long long
using namespace std;

IL int gi(){
    RG int data = 0 , m = 1; RG char ch = 0;
    while(ch != '-' && (ch<'0' || ch > '9')) ch = getchar();
    if(ch == '-'){m = 0; ch = getchar();}
    while(ch>='0' && ch<='9'){data = (data<<1) + (data<<3) + ch - '0' ;  ch = getchar();}
    return (m) ? data : -data ; 
}

struct HJT{int ls,rs; ll sumK,sumB,tagK,tagB;}t[40*_] ; 
struct YCB{
    int l,r,ps; ll dk,db;
    bool operator < (const YCB &B) const{
        return ps < B.ps ; 
    }
}q[_<<1];
int tot,n,m,Q,X1,X2,Y1,Y2,d,xx,yy,Y[_],rt[_],oo,yoy; ll S,ans ; 

void Update(int &o,ll l,ll r,int ql,int qr,ll dk,ll db){
    t[++oo] = t[o]; o = oo ; 
    if(ql <= l && r <= qr){
        t[o].tagK += dk ; t[o].tagB += db ;
        t[o].sumK += 1ll * (r - l + 1) * dk ;
        t[o].sumB += 1ll * (r - l + 1) * db ; 
        return ; 
    }RG int mid = (l + r) >> 1;
    if(ql <= mid) Update(t[o].ls , l , mid , ql , qr , dk , db) ;
    if(qr  > mid) Update(t[o].rs , mid + 1 , r , ql , qr , dk , db) ;
    t[o].sumK = t[t[o].ls].sumK + t[t[o].rs].sumK + (r-l+1) * t[o].tagK ;
    t[o].sumB = t[t[o].ls].sumB + t[t[o].rs].sumB + (r-l+1) * t[o].tagB ; 
}
ll Query(int &o,int l,int r,int ql,int qr,ll x){
    if(!o) return 0;
    if(ql == l && r == qr) return 1ll * t[o].sumK * x + t[o].sumB ;
    RG int mid = (l + r) >> 1;
    RG ll Data = (qr-ql+1) * ( t[o].tagK * x  + t[o].tagB );
    if(qr <= mid) return Data + Query(t[o].ls,l,mid,ql,qr,x) ;
    else if(ql > mid) return Data + Query(t[o].rs,mid+1,r,ql,qr,x) ;
    else return
             Data +
             Query(t[o].ls,l,mid,ql,mid,x) + Query(t[o].rs,mid+1,r,mid+1,qr,x) ;
    return 0;
}

int main(){
    freopen("c.in","r",stdin) ;
    freopen("c.out","w",stdout) ; 
    n = gi(); m = gi(); d = gi(); Q = gi();
    for(RG int i = 1; i <= d; i ++){
        X1 = gi(); X2 = gi(); Y1 = gi(); Y2 = gi(); S = gi();
        q[++tot] = (YCB){X1 , X2 , Y1 , S , 1ll*S*(1-Y1)} ; 
        q[++tot] = (YCB){X1 , X2 , Y2+1 , -S , 1ll*S*Y2 } ;
        Y[++yoy] = Y1 ; Y[++yoy] = Y2 + 1;
    }
    sort(q + 1 , q + tot + 1) ;
    sort(Y + 1 , Y + yoy + 1) ;
    rt[0] = ++ oo ;
    for(RG int i = 1; i <= tot; i ++)
        rt[i] = rt[i-1] , Update(rt[i] , 1 , X , q[i].l , q[i].r , q[i].dk , q[i].db) ;
    ans = 0;
    while(Q --){
        xx = gi(); yy = gi();
        X1 = ans % n + 1; X2 = (ans + xx) % n + 1 ;
        Y1 = ans % m + 1; Y2 = (ans + yy) % m + 1 ;
        if(X1 > X2) swap(X1 , X2) ;
        if(Y1 > Y2) swap(Y1 , Y2) ;
        xx = upper_bound(Y + 1 , Y + yoy + 1 , Y1 - 1) - Y - 1 ;
        yy = upper_bound(Y + 1 , Y + yoy + 1 , Y2) - Y - 1 ;
        ans = 0;
        ans = ans + Query(rt[yy] , 1 , X , X1 , X2 , Y2) ;
        ans = ans - Query(rt[xx] , 1 , X , X1 , X2 , Y1-1) ;
        printf("%llu",ans) ; puts("");
    }return 0;
}