zoj 3540

/*
题意：给你n*m的棋盘，有些棋盘被一些矩形覆盖，求放进一个长为M，宽为1的矩形的方案数

离散，线段树 +扫描线 
分析：最直接得想法，对于每一层统计空白连续的长度然后统计放的方案数，不过编程复杂度有点高
要分别记录下每个区间左边连续最长空白长度和右边最长连续空白长度，
然后count[rt]=count[r<<1]+count[rt<<1|1]+{跨两段的个数} 

一个好想法：当以y为分层标准时，统计竖着放的方案数，那么把每个矩形都向下延伸M-1个单位，
那么剩下的每一个空位都对应一种方法，直接矩形并求面积就好，问题得到转化；

一些细节：扫描线的细节处非常有意思;
首先题目给我们的数据是对应一段区间的，比如矩形[2,2,2,2]表示长宽都为1的矩形,为了方便，我们把它转化成
[1,1,2,2]这样在统计的时候直接相减就可以了， 
*/
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<vector>
#include<set>
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define Find(i) lower_bound(LX.begin(),LX.begin()+n1,i)-LX.begin()
using namespace std;
const int NN=100000+10;
typedef long long LL;
struct Edge{
    int hi,ls,rs,s;
    Edge(){}
    Edge(int a,int b,int c,int e):hi(a),ls(b),rs(c),s(e){}
    bool operator < (const Edge &p)const{
        return hi<p.hi;
    }
}; 
vector<Edge> E;
vector<int> LX;
int len[NN<<2],cov[NN<<2];
int W,H,N,M,n1;
int x[NN/2],y[NN/2],x2[NN/2],y2[NN/2];
LL ans;
void pushup(int l,int r,int rt){
    if (cov[rt]>=1){
        len[rt]=(LX[r+1]-LX[l]);
    }else if (l==r) len[rt]=0;
    else len[rt]=len[rt<<1]+len[rt<<1|1]; 
}
void update(int L,int R,int v,int l,int r,int rt){
    if (L<=l && r<=R){
        cov[rt]+=v;
        if (cov[rt]>=1) len[rt]=(LX[r+1]-LX[l]);
        else if (l==r) len[rt]=0;
        else len[rt]=len[rt<<1]+len[rt<<1|1];
        return ;
    }
    int m=(l+r)>>1;
    if (L<=m) update(L,R,v,lson);
    if (m< R) update(L,R,v,rson);
    pushup(l,r,rt);
}

void work(){
    E.clear();LX.clear();
    E.push_back(Edge(max(0,H-M+1),0,W,1));
    E.push_back(Edge(H,0,W,-1));
    LX.push_back(0);LX.push_back(W);//最后面要加入一个大矩形    
    for (int i=0;i<N;i++){
        E.push_back(Edge(max(0,y[i]-M),x[i]-1,x2[i],1));//矩形往下延伸 
        E.push_back(Edge(y2[i],x[i]-1,x2[i],-1));
        LX.push_back(x[i]-1);LX.push_back(x2[i]);

    }
    sort(E.begin(),E.end());
    sort(LX.begin(),LX.end());
    n1=unique(LX.begin(),LX.end())-LX.begin();
    
    LL ret=0;
    memset(len,0,sizeof(len));
    memset(cov,0,sizeof(cov));
//    cout<<"**** "<<endl;
    for (int i=0;i<E.size();i++){
    //    cout<<E[i].hi<<" "<<E[i].ls<<" "<<E[i].rs<<" "<<E[i].s<<endl;
        int l=Find(E[i].ls),r=Find(E[i].rs)-1;
        if (l<=r) update(l,r,E[i].s,0,n1-1,1);
        if (E[i].hi!=E[i+1].hi){
        //    cout<<len[1]<<endl;
            ret+=(LL)len[1]*(E[i+1].hi-E[i].hi);
        }
    }
    ans+=(LL)W*H-ret;
}
int main(){
    while (~scanf("%d%d%d%d",&W,&H,&N,&M)){
        for (int i=0;i<N;i++){
            scanf("%d%d%d%d",x+i,y+i,x2+i,y2+i);
        }
        ans=0;
        work();
        if (M!=1) {
            swap(H,W);
            for (int i=0;i<N;i++){
                swap(x[i],y[i]);swap(x2[i],y2[i]);
            }
            work();
        }
        printf("%lld\n",ans);
    }
}