扫描线模板

董晓的板子,算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);
}

  • 2
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值