# hdu 4052 Adding New Machine 扫描线求矩形面积并

690人阅读 评论(0)

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <set>
#define ll long long
#define midt (tr[t].l+tr[t].r>>1)
#define ls t<<1
#define rs t<<1|1
using namespace std;
const int maxn=7e4+9;
int w,h,n,m;
int y[maxn<<2],lon;
struct
{
int xl,xu,yl,yu;
}rec[maxn],now[maxn<<1];

struct D
{
int x,y1,y2,count;
bool operator <(const D &xx) const
{
if(x!=xx.x) return(x<xx.x);
else return count>xx.count;
}
}d[maxn<<1];

struct
{
int l,r,y1,y2,count,len;
}tr[maxn<<5];

void maketree(int t,int l,int r)
{
tr[t].l=l;
tr[t].r=r;
tr[t].y1=y[l];
tr[t].y2=y[r];
tr[t].len=0;
tr[t].count=0;
if(l+1==r) return ;

int mid=midt;
maketree(ls,l,mid);
maketree(rs,mid,r);
}

void insert(int t,int y1,int y2,D &d)
{
if(tr[t].y1==y1&&tr[t].y2==y2)
{
tr[t].count+=d.count;
}
else
{
int mid=midt;
if(y2<=y[mid])
insert(ls,y1,y2,d);
else if(y[mid]<=y1)
insert(rs,y1,y2,d);
else
{
insert(ls,y1,y[mid],d);
insert(rs,y[mid],y2,d);
}
}
if(tr[t].count>0) tr[t].len=tr[t].y2-tr[t].y1;
else if(tr[t].l+1==tr[t].r) tr[t].len=0;
else tr[t].len=tr[ls].len+tr[rs].len;
}

ll solve()
{
set <int> s;
s.clear();
for(int i=1;i<=n;i++)
{
d[i*2-1].x=now[i].xl;
d[i*2-1].y1=now[i].yl;
d[i*2-1].y2=now[i].yu;
d[i*2-1].count=1;

d[i*2].x=now[i].xu;
d[i*2].y1=now[i].yl;
d[i*2].y2=now[i].yu;
d[i*2].count=-1;

s.insert(now[i].yl);
s.insert(now[i].yu);
}
lon=0;
for(set <int> ::iterator k=s.begin();k!=s.end();k++)
{
y[++lon]=*k;
}

if(lon==0) return 0;
sort(d+1,d+1+n+n);

maketree(1,1,lon);
ll ans=0;
d[0].x=-1;
for(int i=1;i<=n*2;i++)
{
if(d[i].y1<=d[i].y2)
{
ans+=(ll)(d[i].x-d[i-1].x)*tr[1].len;
insert(1,d[i].y1,d[i].y2,d[i]);
}
}
return ans;
}

int main()
{
//    freopen("in.txt","r",stdin);
while(scanf("%d%d%d%d",&w,&h,&n,&m)!=EOF)
{
for(int i=1;i<=n;i++)
{
scanf("%d%d",&rec[i].xl,&rec[i].yl);
scanf("%d%d",&rec[i].xu,&rec[i].yu);
rec[i].xl--;
rec[i].yl--;
}
for(int i=1;i<=n;i++)
{
now[i].xl=max(0,rec[i].xl-m+1);
now[i].xu=min(rec[i].xu,w-m+1);
now[i].yl=rec[i].yl;
now[i].yu=rec[i].yu;
}
ll ret=0;
if(w-m+1>0)
ret+=solve();

for(int i=1;i<=n;i++)
{
now[i].xl=rec[i].xl;
now[i].xu=rec[i].xu;
now[i].yl=max(0,rec[i].yl-m+1);
now[i].yu=min(rec[i].yu,h-m+1);
}
if(h-m+1>0)
ret+=solve();

ret=max((ll)(w-m+1)*h,(ll)0)+max((ll)w*(h-m+1),(ll)0)-ret;
if(m==1)
ret/=2;
cout<<ret<<endl;
}
return 0;
}


0
0

【直播】机器学习&数据挖掘7周实训--韦玮
【套餐】系统集成项目管理工程师顺利通关--徐朋
【直播】3小时掌握Docker最佳实战-徐西宁
【套餐】机器学习系列套餐（算法+实战）--唐宇迪
【直播】计算机视觉原理及实战--屈教授
【套餐】微信订阅号+服务号Java版 v2.0--翟东平
【直播】机器学习之矩阵--黄博士
【套餐】微信订阅号+服务号Java版 v2.0--翟东平
【直播】机器学习之凸优化--马博士
【套餐】Javascript 设计模式实战--曾亮

* 以上用户言论只代表其个人观点，不代表CSDN网站的观点或立场
个人资料
• 访问：168622次
• 积分：4321
• 等级：
• 排名：第6942名
• 原创：274篇
• 转载：2篇
• 译文：0篇
• 评论：11条
最新评论