矩形周长并 poj 1177Picture/hdu 1828

同上篇一样，具体的去看论文吧

我没看过论文，不知道正解，但是我的代码能够A提，这就足够了

分两个方向分别扫一边

bug：

1.先处理入边，在处理出边

2.线段树query函数注意一下（代码种有标记）

#include <iostream>
#include <cstdio>
#include <algorithm>
using namespace std;
#define mid (L+R>>1)
#define lson L,mid,P<<1
#define rson mid,R,P<<1|1
const int M = 10010;
class segtree{
public:
int left[M<<2],right[M<<2];
int sum[M<<2],num[M<<2],s[M],n;

void init(int ss[],int k){
n=0;
for(int i=0;i<k;i++){
if(i==0||s
!=ss[i]){
s[++n]=ss[i];
}
}
__build(1,n,1);
}
void __build(int L,int R,int P){
left[P]=s[L],right[P]=s[R];
num[P]=0,sum[P]=0;
if(L+1==R)return ;
__build(lson);
__build(rson);
}
void pushup(int x,int L,int R){
if(num[x]){
sum[x]=right[x]-left[x];
}else{
if(L+1!=R){
sum[x]=sum[x<<1]+sum[x<<1|1];
}else{
sum[x]=0;
}
}
}
void __update(int L,int R,int P,int l,int r,int f){
if(l==left[P]&&r==right[P]){
num[P]+=f;
pushup(P,L,R);
return ;
}
if(l<right[P<<1]){
__update(lson,l,min(right[P<<1],r),f);
}
if(r>right[P<<1]){
__update(rson,max(l,right[P<<1]),r,f);
}
pushup(P,L,R);
}
int __query(int L,int R,int P,int l,int r){
if(sum[P]==right[P]-left[P]){//注意细节
return r-l;
}
if(l==left[P]&&r==right[P]){
return sum[P];
}
int ret=0;
if(l<right[P<<1]){
ret+=__query(lson,l,min(r,right[P<<1]));
}
if(r>right[P<<1]){
ret+=__query(rson,max(l,right[P<<1]),r);
}
return ret;
}
void update(int l,int r,int f){__update(1,n,1,l,r,f);}
int query(int l,int r){return __query(1,n,1,l,r);}
};
struct Line{
int x,y1,y2,f;
void set(int _x,int _y1,int _y2,int _f){
x=_x,y1=_y1,y2=_y2,f=_f;
}
}ll[M<<1],rr[M<<1];
int xx[M<<1],yy[M<<1],n;

bool cmp(Line a,Line b){
if(a.x!=b.x){
return a.x<b.x;
}else{
return a.f>b.f;
}
}
int main(){
while(cin>>n){
segtree *tree;
tree = new segtree;
for(int i=0,x1,y1,x2,y2;i<n;i++){
scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
ll[i<<1].set(x1,y1,y2,1);
ll[i<<1|1].set(x2,y1,y2,-1);
rr[i<<1].set(y1,x1,x2,1);
rr[i<<1|1].set(y2,x1,x2,-1);
yy[i<<1]=y1,yy[i<<1|1]=y2;
xx[i<<1]=x1,xx[i<<1|1]=x2;
}

sort(ll,ll+n*2,cmp);
sort(yy,yy+n*2);
tree->init(yy,n*2);
int sum=ll[0].y2-ll[0].y1;
tree->update(ll[0].y1,ll[0].y2,ll[0].f);
for(int i=1;i<n*2;i++){
if(ll[i].f==1){
sum+=(ll[i].y2-ll[i].y1-tree->query(ll[i].y1,ll[i].y2));
tree->update(ll[i].y1,ll[i].y2,ll[i].f);
}else{
tree->update(ll[i].y1,ll[i].y2,ll[i].f);
sum+=(ll[i].y2-ll[i].y1-tree->query(ll[i].y1,ll[i].y2));
}
}
delete tree;
tree = new segtree;

sort(rr,rr+n*2,cmp);
sort(xx,xx+n*2);
tree->init(xx,n*2);
sum+=rr[0].y2-rr[0].y1;
tree->update(rr[0].y1,rr[0].y2,rr[0].f);
for(int i=1;i<n*2;i++){
if(rr[i].f==1){
sum+=(rr[i].y2-rr[i].y1-tree->query(rr[i].y1,rr[i].y2));
tree->update(rr[i].y1,rr[i].y2,rr[i].f);
}else{
tree->update(rr[i].y1,rr[i].y2,rr[i].f);
sum+=(rr[i].y2-rr[i].y1-tree->query(rr[i].y1,rr[i].y2));
}
}
cout<<sum<<endl;
delete tree;
}
return 0;
}