hdu 4052 Adding New Machine(线段树+扫描线)

最新推荐文章于 2026-06-15 09:03:01 发布

原创最新推荐文章于 2026-06-15 09:03:01 发布 · 849 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

标签

#datastruct

数据结构同时被 2 个专栏收录

91 篇文章

订阅专栏

线段树

50 篇文章

订阅专栏

本文探讨了一个关于在W×H的矩形区域内摆放1×m大小设备的问题，旨在计算所有可能的摆放方式数量，同时确保任意两个设备不会相互重叠。文章通过横纵两向分别枚举并排除原有设备占用的空间，最终利用矩形面积并集的方法计算有效摆放区域。

题意：给出一个W×H的矩形，有些地方有一些设备，现在要添加一个1×m的设备，可以横着或竖着放，问有多少种放法，并且任意两个矩形相互覆盖。

思路：可以想到分为横着和竖着放都枚举一下，对于原有的设备，将矩形扩大一下，标记出来不能放置的位置，最后用总的面积减去不能放的面积就行了，这就转化为矩形面积并，本来以为挺好写的，由于矩形面积并掌握的不是很好，被华丽丽的坑了，果然自己YY的不可信啊。。。

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<map>
#include<queue>
#include<stack>
#include<cmath>
#include<vector>
#define inf 0x3f3f3f3f
#define Inf 0x3FFFFFFFFFFFFFFFLL
#define eps 1e-9
#define pi acos(-1.0)
using namespace std;
typedef long long ll;
const int maxn=100000+10;
struct Node
{
    int x,L,R,c;
    Node(){}
    Node(int x,int L,int R,int c):x(x),L(L),R(R),c(c){}
    bool operator <(const Node &a) const
    {
        return x<a.x;
    }
}node[maxn];
struct Mat
{
    int x1,y1,x2,y2;
    void read()
    {
        scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
        if(x1>x2) swap(x1,x2);
        if(y1>y2) swap(y1,y2);
    }
}mat[maxn];
map<int,int>mp;
int num[maxn];
int sum[maxn<<2],cover[maxn<<2];
inline void PushUp(int l,int r,int rt)
{
    if(cover[rt]) sum[rt]=num[r]-num[l-1];
    else if(l==r) sum[rt]=0;
    else sum[rt]=sum[rt<<1]+sum[rt<<1|1];
}
void build(int l,int r,int rt)
{
    sum[rt]=cover[rt]=0;
    if(l==r) return ;
    int m=(l+r)>>1;
    build(l,m,rt<<1);
    build(m+1,r,rt<<1|1);
}
void Update(int L,int R,int l,int r,int rt,int v)
{
    if(l>=L&&r<=R)
    {
        cover[rt]+=v;
        PushUp(l,r,rt);
        return ;
    }
    int m=(l+r)>>1;
    if(m>=L) Update(L,R,l,m,rt<<1,v);
    if(m<R) Update(L,R,m+1,r,rt<<1|1,v);
    PushUp(l,r,rt);
}
ll solve(int W,int H,int n,int m)
{
    if(m>H) return 0;
    mp.clear();
    int tot=0,L,R,N;
    for(int i=0;i<n;++i)
    {
        L=mat[i].y1;
        R=min(H,mat[i].y2+m-1)+1;
        num[tot++]=L;
        num[tot++]=R;
        node[i]=Node(mat[i].x1,L,R,1);
        node[i+n]=Node(mat[i].x2+1,L,R,-1);
    }
    L=1;
    R=min(H,m-1)+1;
    if(L<=R)
    {
        node[tot]=Node(1,L,R,1);
        num[tot++]=L;
        node[tot]=Node(W+1,L,R,-1);
        num[tot++]=R;
    }
    sort(num,num+tot);
    sort(node,node+tot);
    ll ans=0;
    N=unique(num,num+tot)-num;
    for(int i=0;i<N;++i)
        mp[num[i]]=i+1;
    if(N) build(1,N,1);
    for(int i=0;i<tot-1;++i)
    {
        L=mp[node[i].L];
        R=mp[node[i].R]-1;
        if(L<=R) Update(L,R,1,N,1,node[i].c);
        ans+=(ll)sum[1]*(node[i+1].x-node[i].x);
    }
    ans=(ll)H*W-ans;
    return ans;
}
int main()
{
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    int W,H,n,m;
    while(~scanf("%d%d%d%d",&W,&H,&n,&m))
    {
        for(int i=0;i<n;++i)
            mat[i].read();
        ll ans=solve(W,H,n,m);
        for(int i=0;i<n;++i)
        {
            swap(mat[i].x1,mat[i].y1);
            swap(mat[i].x2,mat[i].y2);
        }
        swap(W,H);
        ans+=solve(W,H,n,m);
        if(m==1) ans/=2;
        printf("%I64d\n",ans);
    }
    return 0;
}