浅谈树状数组与线段树:https://www.cnblogs.com/AKMer/p/9946944.html
题目传送门:http://acm.hdu.edu.cn/showproblem.php?pid=1828
扫描线求周长。如果你不知道什么是扫描线可以先去做这题:https://www.cnblogs.com/AKMer/p/9951300.html
每次扫描线上变化的长度就是周长的一部分。横着做一遍竖着做一遍就行了。
如果两条边\(pos\)相同应该先加边后删边,否则重叠部分会多算。
时间复杂度:\(O(nlogn)\)
空间复杂度:\(O(n)\)
代码如下:
#include <cmath>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long ll;
const int maxn=10005;
int n;
ll ans;
int tmp[2][maxn];
int read() {
int x=0,f=1;char ch=getchar();
for(;ch<'0'||ch>'9';ch=getchar())if(ch=='-')f=-1;
for(;ch>='0'&&ch<='9';ch=getchar())x=x*10+ch-'0';
return x*f;
}
struct line {
int pos,st,ed,mark;
line() {}
line(int _pos,int _st,int _ed,int _mark) {
pos=_pos,st=_st,ed=_ed,mark=_mark;
}
bool operator<(const line &a)const {
if(pos==a.pos)return mark>a.mark;//pos相同优先加边,这样就不会算重叠部分的长度
return pos<a.pos;
}
}p[2][maxn];
struct segment_tree {
int len[maxn<<2],cnt[maxn<<2];
void clear() {
memset(len,0,sizeof(len));
memset(cnt,0,sizeof(cnt));
}
void updata(int id,int p,int l,int r) {
if(cnt[p])len[p]=tmp[id][r+1]-tmp[id][l];
else if(l==r)len[p]=0;
else len[p]=len[p<<1]+len[p<<1|1];
}
void change(int id,int p,int l,int r,int L,int R,int v) {
if(L<=l&&r<=R) {
cnt[p]+=v;
updata(id,p,l,r);
return;
}
int mid=(l+r)>>1;
if(L<=mid)change(id,p<<1,l,mid,L,R,v);
if(R>mid)change(id,p<<1|1,mid+1,r,L,R,v);
updata(id,p,l,r);
}
}T;
int main() {
while(~scanf("%d",&n)) {
ans=0;
for(int i=1;i<=n;i++) {
int x1=read(),y1=read(),x2=read(),y2=read();
tmp[0][(i<<1)-1]=y1,tmp[0][i<<1]=y2;
tmp[1][(i<<1)-1]=x1,tmp[1][i<<1]=x2;
p[0][(i<<1)-1]=line(x1,y1,y2,1);
p[0][i<<1]=line(x2,y1,y2,-1);
p[1][(i<<1)-1]=line(y1,x1,x2,1);
p[1][i<<1]=line(y2,x1,x2,-1);
}for(int i=0;i<2;i++) {
T.clear();
sort(p[i]+1,p[i]+2*n+1);
sort(tmp[i]+1,tmp[i]+2*n+1);
int cnt=unique(tmp[i]+1,tmp[i]+2*n+1)-tmp[i]-1;
for(int j=1;j<=2*n;j++) {
p[i][j].st=lower_bound(tmp[i]+1,tmp[i]+cnt+1,p[i][j].st)-tmp[i];
p[i][j].ed=lower_bound(tmp[i]+1,tmp[i]+cnt+1,p[i][j].ed)-tmp[i];
int tmp=T.len[1];
T.change(i,1,1,cnt-1,p[i][j].st,p[i][j].ed-1,p[i][j].mark);
ans+=abs(T.len[1]-tmp);//差值就是周长的一部分
}
}
printf("%lld\n",ans);
}
return 0;
}