董晓的板子,算len的时候r要+1,进去修改的时候r要-1,数组其实要开16*n,因为n会扩大两倍。
面积...
#include <iostream>
#include <cstdio>
#include <algorithm>
using namespace std;
#define ls u<<1
#define rs u<<1|1
const int N=200005;
struct line{ //扫描线
int x1,x2,y;
int tag; //入边:+1,出边:-1
bool operator<(line &t){return y<t.y;}
}L[N];
struct tree{ //线段树
int l,r;
int cnt,len; //区间覆盖次数和覆盖长度
}tr[N*8];
int X[N]; //X坐标
void build(int u,int l,int r){ //建树
tr[u]={l,r,0,0};
if(l==r) return;
int mid=(l+r)>>1;
build(ls,l,mid);
build(rs,mid+1,r);
}
void pushup(int u){
int l=tr[u].l, r=tr[u].r; //r → X[r+1]
if(tr[u].cnt) tr[u].len=X[r+1]-X[l];
else tr[u].len=tr[ls].len+tr[rs].len;
}
void change(int u,int l,int r,int tag){ //区修
if(l>tr[u].r || r<tr[u].l)return;
if(l<=tr[u].l && tr[u].r<=r){
tr[u].cnt+=tag;
pushup(u); //这里会爆4倍空间
return;
}
change(ls,l,r,tag);
change(rs,l,r,tag);
pushup(u);
}
int main(){
int n,x1,x2,y1,y2; scanf("%d",&n);
for(int i=1; i<=n; i++){
scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
L[i]={x1,x2,y1,1};
L[n+i]={x1,x2,y2,-1};
X[i]=x1; X[n+i]=x2;
}
n*=2;
sort(L+1,L+n+1); //扫描线排序
sort(X+1,X+n+1); //X坐标排序
int s=unique(X+1,X+n+1)-X-1; //去重
build(1,1,s-1);
long long ans=0;
for(int i=1; i<n; i++){
int l=lower_bound(X+1,X+s+1,L[i].x1)-X;
int r=lower_bound(X+1,X+s+1,L[i].x2)-X;
change(1,l,r-1,L[i].tag); //x2 → r-1
ans+=1ll*tr[1].len*(L[i+1].y-L[i].y);
}
printf("%lld\n",ans);
}
周长....
#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
#define ls u<<1
#define rs u<<1|1
#define N 10005
struct line{ //扫描线
int l,r,y;
int tag; //入边:+1, 出边:-1
bool operator<(line &t){
return y==t.y ? tag>t.tag : y<t.y;
} //入边优先,两矩形相切等价于相交
}L[N];
struct tree{ //线段树
int l,r;
int cnt,len; //区间覆盖次数和覆盖长度
int sum; //区间覆盖竖边条数
bool lcover,rcover;//左右端点是否被覆盖
}tr[N*8];
int X[N]; //X坐标
void build(int u,int l,int r){ //建树
tr[u].l=l, tr[u].r=r;
if(l==r) return;
int mid=(l+r)>>1;
build(ls,l,mid);
build(rs,mid+1,r);
}
void pushup(int u){ //上传
int l=tr[u].l,r=tr[u].r;
if(tr[u].cnt){
tr[u].len=X[r+1]-X[l];
tr[u].sum=2;
tr[u].lcover=tr[u].rcover=true;
}
else{
tr[u].len=tr[ls].len+tr[rs].len;
tr[u].sum=tr[ls].sum+tr[rs].sum;
tr[u].lcover=tr[ls].lcover;
tr[u].rcover=tr[rs].rcover;
if(tr[ls].rcover && tr[rs].lcover)
tr[u].sum-=2; //左右区间连续
}
}
void change(int u,int l,int r,int tag){ //区修
if(l>tr[u].r || r<tr[u].l) return;
if(l<=tr[u].l && tr[u].r<=r){
tr[u].cnt+=tag;
pushup(u);
return;
}
change(ls,l,r,tag);
change(rs,l,r,tag);
pushup(u);
}
int main(){
int n,x1,y1,x2,y2,last=0;
scanf("%d",&n);
for(int i=1; i<=n; i++){
scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
L[i]={x1,x2,y1,1};
L[n+i]={x1,x2,y2,-1};
X[i]=x1; X[n+i]=x2;
}
n*=2;
sort(L+1,L+n+1);
sort(X+1,X+n+1);
int s=unique(X+1,X+n+1)-X-1;
build(1,1,s-1);
int res=0;
for(int i=1; i<n; i++){
int l=lower_bound(X+1,X+s+1,L[i].l)-X;
int r=lower_bound(X+1,X+s+1,L[i].r)-X;
change(1,l,r-1,L[i].tag);
res+=abs(tr[1].len-last); //横边
last=tr[1].len;
res+=tr[1].sum*(L[i+1].y-L[i].y); //竖边
}
res+=L[n].r-L[n].l; //最后一条扫描线
printf("%d\n",res);
}