hdu1828(线段树+扫描线)

时间:2022-10-24 08:59:32

又知道了线段树的一种用法,除了单点更新,区间更新,还有这种在一段线段上标号但不往下推。 真是神奇

hdu1828

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <math.h>
#include <algorithm>
#include <string>
#include <queue>
#include <map>
#include <stdlib.h>
using namespace std;
#define N 5050 struct RECT
{
int x1,y1,x2,y2;
}g[N]; struct LINE
{
int k,b,d;
int flag;
}line[*N]; int ans=;
int n; int l[],r[],mark[],yb[],yd[],num[]; int cmp(LINE x,LINE y)
{
return x.k<y.k;
} void build(int tl,int tr,int s)
{
l[s]=tl; r[s]=tr;
mark[s]=;yb[s]=; yd[s]=; num[s]=;
if(tl==tr) return ;
int mid=(tl+tr)/;
build(tl,mid,*s);
build(mid+,tr,*s+);
} void up(int s)
{
if(mark[s]>)
{
num[s]=;
yb[s]=;
yd[s]=;
}
else
{
if(l[s]==r[s])
{
num[s]=;
yb[s]=;
yd[s]=;
}
else
{
if(yd[*s]==&&yb[*s+]==)
num[s]=num[*s]+num[*s+]-;
else num[s]=num[*s]+num[*s+];
yb[s]=yb[*s];
yd[s]=yd[*s+];
}
}
} void update(int tl,int tr,int sign,int s)
{
if(tl==l[s]&&tr==r[s])
{
mark[s]+=sign;
up(s);
return ;
}
int mid=(l[s]+r[s])/;
if(tr<=mid) update(tl,tr,sign,*s);
else if(tl>mid) update(tl,tr,sign,*s+);
else
{
update(tl,mid,sign,*s);
update(mid+,tr,sign,*s+);
}
up(s);
} void fuc()
{
int cnt=;
int mix=,mxx=-; for(int i=;i<n;i++)
{
line[cnt].k=g[i].y1;
line[cnt].b=g[i].x1;
line[cnt].d=g[i].x2-;//注意可能为负的
line[cnt].flag=;
cnt++; mix=min(mix,g[i].x1);
mxx=max(mxx,g[i].x2-); line[cnt].k=g[i].y2;
line[cnt].b=g[i].x1;
line[cnt].d=g[i].x2-;
line[cnt].flag=-;
cnt++;
}
sort(line,line+cnt,cmp);
build(mix,mxx,);
for(int i=;i<cnt-;i++)
{
update(line[i].b,line[i].d,line[i].flag,);
ans += (line[i+].k-line[i].k)**num[];
}
} int main()
{
while(scanf("%d",&n)!=EOF)
{
for(int i=;i<n;i++)
{
scanf("%d%d%d%d",&g[i].x1,&g[i].y1,&g[i].x2,&g[i].y2);
g[i].x1+=;
g[i].x2+=;
g[i].y1+=;
g[i].y2+=;
}
if(n==)
{
printf("0\n");
continue;
}
ans=;
fuc();
//
for(int i=;i<n;i++)
{
swap(g[i].x1,g[i].y1);
swap(g[i].x2,g[i].y2);
}
fuc();
printf("%d\n",ans);
}
return ;
}