hdu 4052 线段树求面积并

解题思路:

我们可以将这个问题转换成看一个坐标向上或者想右（分别表示了竖着和横着）可以放一个新的机器，那么原来旧的机器所占的位置就可以用两个计算面积方式来求面积并，最后用总面积减去机器所占面积就行了。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mx = 1e5+10;
int W,H,N,M;
ll mark[mx<<2];//记录某个区间的下底边个数
ll sum[mx<<2];//记录某个区间的下底边总长度
ll has[mx];//对x进行离散化,否则x为浮点数且很大无法进行线段树

//以横坐标作为线段(区间),对横坐标线段进行扫描
//扫描的作用是每次更新下底边总长度和下底边个数,增加新面积
struct seg{//线段
    int l,r,h;
    int d;
    seg(){}
    seg(int x1,int x2,int H,int c):l(x1),r(x2),h(H),d(c){}
    bool operator<(const seg &a)const{
        return h<a.h;
    }
}s[mx],num[mx];

void Upfather(int n,int left,int right){
    if(mark[n]) sum[n]=has[right+1]-has[left];//表示该区间整个线段长度可以作为底边
    else if(left == right)sum[n]=0;//叶子结点则底边长度为0(区间内线段长度为0)
    else sum[n]=sum[n<<1]+sum[n<<1|1];
}

void Update(int L,int R,int d,int n,int left,int right){
    if(L<=left && right<=R){//该区间是当前扫描线段的一部分,则该区间下底边总长以及上下底边个数差更新
        mark[n]+=d;//更新底边相差差个数
        Upfather(n,left,right);//更新底边长
        return;
    }
    int mid=left+right>>1;
    if(L<=mid)Update(L,R,d,n<<1,left,mid);
    if(R>mid)Update(L,R,d,n<<1|1,mid+1,right);
    Upfather(n,left,right);
}
ll get_ans(int k)
{
    sort(has,has+k);
    sort(s,s+k);
    int m = 1;
    for(int i=1;i<k;++i)//去重复端点
    if(has[i] != has[i-1])has[m++]=has[i];
    ll ans = 0;
    for(int i=0;i<k;++i){//扫描线段
        int L=lower_bound(has,has+m,s[i].l)-has;
        int R=lower_bound(has,has+m,s[i].r)-has-1;
        Update(L,R,s[i].d,1,0,m-1);//扫描线段时更新底边长度和底边相差个数
        ans+=sum[1]*(s[i+1].h-s[i].h);//新增加面积
    }
    return ans;
}
int main()
{
    while(~scanf("%d%d%d%d",&W,&H,&N,&M))
    {
        int l1,l2 = W,r1 = H,r2;
        l1 = W - M + 1;
        r2 = H - M + 1;
        int x1,y1,x2,y2,k = 0;
        for(int i=0;i<N;i++){
            scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
            x1--,y1--;
            num[i] = seg(x1,y1,x2,y2);
        }
        for(int i=0;i<N;i++){
            int t1 = max(0,num[i].l-M+1),t2 = min(num[i].h,l1);
            has[k] = t1;
            s[k++] = seg(t1,t2,num[i].r,1);
            has[k] = t2;
            s[k++] = seg(t1,t2,num[i].d,-1);
        }
        ll ans = 0;
        if(l1>0) ans += 1ll*l1*r1-get_ans(k);
        if(M==1) { printf("%lld\n",ans); continue;  }
        k = 0;
        for(int i=0;i<N;i++){
            int t1 = max(0,num[i].r-M+1),t2 = min(num[i].d,r2);
            has[k] = num[i].l;
            s[k++] = seg(num[i].l,num[i].h,t1,1);
            has[k] = num[i].h;
            s[k++] = seg(num[i].l,num[i].h,t2,-1);
        }
        if(r2>0) ans += 1ll*l2*r2-get_ans(k);
        printf("%lld\n",ans);
    }
    return 0;
}