题目:http://acm.hdu.edu.cn/showproblem.php?pid=1828
题目大意:就是给你n个矩形,每个矩形由左上和右上两个点确定,问你他们的周长并是多少?
思路:线段树+扫描线。同样需要左闭右开处理,稍微和面积并有点不同,底线的长度求法和面积并一样,关键是还要-上次的底线长度,求这个值abs,然后这就是
update当前这条line之后新增的横的长度。然后就是求竖的那种线有几根的问题,如果某个区间被某几条线完全覆盖,那么我们就设它的
竖线有两根,然后在定义两个数组 lbd,rbd,表示这个区间的左右端点存不存在竖线(因为是左闭右开,右端点是-1的,其实是l、r+1 存不存在竖线),
push_up的时候,如果左儿子的rbd存在,右儿子的lbd存在,就要 两个numseg值加起来-2,因为这两条线重合了。具体的可以看代码里的 push_up,我感觉
那里是关键。
代码是直接从面积那里并改过来的。
题目大意:就是给你n个矩形,每个矩形由左上和右上两个点确定,问你他们的周长并是多少?
思路:线段树+扫描线。同样需要左闭右开处理,稍微和面积并有点不同,底线的长度求法和面积并一样,关键是还要-上次的底线长度,求这个值abs,然后这就是
update当前这条line之后新增的横的长度。然后就是求竖的那种线有几根的问题,如果某个区间被某几条线完全覆盖,那么我们就设它的
竖线有两根,然后在定义两个数组 lbd,rbd,表示这个区间的左右端点存不存在竖线(因为是左闭右开,右端点是-1的,其实是l、r+1 存不存在竖线),
push_up的时候,如果左儿子的rbd存在,右儿子的lbd存在,就要 两个numseg值加起来-2,因为这两条线重合了。具体的可以看代码里的 push_up,我感觉
那里是关键。
代码是直接从面积那里并改过来的。
代码如下:
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
#define lson rt<<1,l,mid
#define rson rt<<1|1,mid+1,r
const int MAXN = 5005<<1;
struct Line
{
int l,r,h;
int flag;
Line(){}
Line(int a,int b,int c,int d)
{
l = a;
r = b;
h = c;
flag = d;
}
bool operator < (const Line &tmp) const
{
return h < tmp.h;
}
} line[MAXN];
int x[MAXN];
int Bin(int l,int r,int key)
{
while(l <= r)
{
int mid = (l+r)>>1;
if(x[mid] == key)
{
return mid;
}
else if(x[mid] < key)
l = mid+1;
else r = mid-1;
}
return -1;
}
int sum[MAXN<<2],cnt[MAXN<<2],numseg[MAXN<<2],lbd[MAXN<<2],rbd[MAXN<<2];
void push_up(int rt,int l,int r)
{
if(cnt[rt])
{
sum[rt] = x[r+1]-x[l];
numseg[rt] = 2;
lbd[rt] = rbd[rt] = 1;
}
else if(l == r) sum[rt] = numseg[rt] = lbd[rt] = rbd[rt] = 0;
else
{
sum[rt] = sum[rt<<1]+sum[rt<<1|1];
numseg[rt] = numseg[rt<<1]+numseg[rt<<1|1];
if(lbd[rt<<1|1] && rbd[rt<<1]) numseg[rt] -= 2;
lbd[rt] = lbd[rt<<1];
rbd[rt] = rbd[rt<<1|1];
}
}
void update(int rt,int l,int r,int a,int b,int c)
{
if(a <= l && b >= r)
{
cnt[rt] += c;
push_up(rt,l,r);
return ;
}
int mid = l+r>>1;
if(a <= mid) update(lson,a,b,c);
if(b > mid) update(rson,a,b,c);
push_up(rt,l,r);
}
void build(int rt,int l,int r)
{
cnt[rt] = sum[rt] = numseg[rt] = lbd[rt] = rbd[rt] = 0;
if(l == r)
return ;
int mid = l+r>>1;
build(lson);
build(rson);
}
int main()
{
int n;
while(~scanf("%d",&n))
{
int a,b,c,d;
int tot = 0;
for(int i = 0;i < n;i++)
{
scanf("%d%d%d%d",&a,&b,&c,&d);
x[tot] = a;
line[tot++] = Line(a,c,b,1);
x[tot] = c;
line[tot++] = Line(a,c,d,-1);
}
sort(x,x+tot);
sort(line,line+tot);
int m = unique(x,x+tot)-x;
build(1,0,m-1);
int ans = 0;
int last = 0;
for(int i = 0;i < tot;i++)
{
int l = Bin(0,m-1,line[i].l);
int r = Bin(0,m-1,line[i].r);
if(l >= r) continue;
//printf("l = %d,r = %d,m = %d\n",l,r,m);
update(1,0,m-1,l,r-1,line[i].flag);
ans += numseg[1]*(i == tot-1 ? 0 : (line[i+1].h-line[i].h));
ans += abs(sum[1]-last);
last = sum[1];
//printf("%d,%d,ans = %d\n",numseg[1],sum[1],ans);
}
printf("%d\n",ans);
}
return 0;
}
/*
2
10 10 20 20
15 15 25 25
*/