参考博客:http://www.cnblogs.com/scau20110726/archive/2013/04/13/3018687.html
方法:
把矩形分成横线和竖线去处理,可知是完全相同的操作,我们来讲下怎么算出横线部分,竖线部分就是照搬即可。
将横线保存在一个表中,按横线所处的竖直位置排序(升序),另外每条横线带一个标记值,原矩形的下线为1,上线为-1(对应过去就是插入线段和删除线段)
从低到高扫描横线,没扫到一条横线就能计算出一部分横线值。计算方法是算出现在总区间的被覆盖的长度,然后求出与上一次的总区间的覆盖长度的差(即相减求绝对值),因为每次添加了一条线段,如果没有没有使总区间覆盖长度发生变化,说明这条线段其实在多边形的内部,被覆盖掉了,不能计算,只要能引起总区间长度发生变化的,说明该线段不被覆盖不被包含,要计算。
而竖线部分的做法是一样的,把竖线保存在一个表中,按水平位置排序(升序),每条横线带一个标记值,原矩形的左线为1,右线为-1,然后同样地操作
我用的是第一种方法,思路和上面的博客一样,只有一个地方需要注意。
注意:扫描线的排序方法 return a.h==b.h?a.status>b.status:a.h<b.h;
为什么这么排序?样例 :2 0 0 1 1 0 1 1 2 自己手动模拟一下这组样例就会明白了
代码:
///#include<bits/stdc++.h>
///#include<unordered_map>
///#include<unordered_set>
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<string>
#include<cmath>
#include<queue>
#include<bitset>
#include<set>
#include<stack>
#include<map>
#include<list>
#include<new>
#include<vector>
#define MT(a,b) memset(a,b,sizeof(a));
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const double pai=acos(-1.0);
const double E=2.718281828459;
const int MAXN=1e5;
const ll mod=1e9+7;
const int INF=0x3f3f3f3f;
struct node
{
ll l;
ll r;
ll h;
ll status;
bool friend operator<(node a,node b)
{
///重要的一步 样例 2 0 0 1 1 0 1 1 2
return a.h==b.h?a.status>b.status:a.h<b.h;
}
}line[10005];
ll x[10005];
struct xx
{
ll s;
ll sum;
}root[40005];
struct pp
{
ll x1;
ll y1;
ll x2;
ll y2;
}point[5005];
void pushup(int sign,int l,int r)
{
if(root[sign].s)
root[sign].sum=x[r+1]-x[l];
else if(l==r)
root[sign].sum=0;
else
root[sign].sum=root[sign<<1].sum+root[sign<<1|1].sum;
}
void updateroot(int sign,int l,int r,int a,int b,int f)
{
if(l==a&&r==b)
{
root[sign].s+=f;
pushup(sign,l,r);
return ;
}
int mid=(l+r)>>1;
if(b<=mid)
updateroot(sign<<1,l,mid,a,b,f);
else if(a>mid)
updateroot(sign<<1|1,mid+1,r,a,b,f);
else
{
updateroot(sign<<1,l,mid,a,mid,f);
updateroot(sign<<1|1,mid+1,r,mid+1,b,f);
}
pushup(sign,l,r);
}
int main()
{
ll n;
ll x1,y1,x2,y2;
while(scanf("%lld",&n)!=EOF)
{
memset(root,0,sizeof(root));
int sign=0;
for(int i=1;i<=n;i++)
{
scanf("%lld %lld %lld %lld",&x1,&y1,&x2,&y2);
point[i]=pp{x1,y1,x2,y2};
line[++sign]=node{x1,x2,y1,1};
x[sign]=x1;
line[++sign]=node{x1,x2,y2,-1};
x[sign]=x2;
}
sort(line+1,line+1+sign);
sort(x+1,x+1+sign);
int d=unique(x+1,x+1+sign)-(x+1);
ll len=abs(line[sign].r-line[sign].l);
ll pre=0;
for(int i=1;i<sign;i++)
{
int a=lower_bound(x+1,x+1+d,line[i].l)-x;
int b=lower_bound(x+1,x+1+d,line[i].r)-x-1;
updateroot(1,1,d-1,a,b,line[i].status);
len+=abs(root[1].sum-pre);
pre=root[1].sum;
}
memset(root,0,sizeof(root));
sign=0;
for(int i=1;i<=n;i++)
{
x1=point[i].x1;
y1=point[i].y1;
x2=point[i].x2;
y2=point[i].y2;
line[++sign]=node{y1,y2,x1,1};
x[sign]=y1;
line[++sign]=node{y1,y2,x2,-1};
x[sign]=y2;
}
sort(x+1,x+1+sign);
sort(line+1,line+1+sign);
d=unique(x+1,x+1+sign)-(x+1);
len+=abs(line[sign].r-line[sign].l);
pre=0;
for(int i=1;i<sign;i++)
{
int a=lower_bound(x+1,x+1+d,line[i].l)-x;
int b=lower_bound(x+1,x+1+d,line[i].r)-x-1;
updateroot(1,1,d-1,a,b,line[i].status);
len+=abs(root[1].sum-pre);
pre=root[1].sum;
}
printf("%lld\n",len);
}
return 0;
}