線段樹+掃描線【p1884】[Usaco12FEB]過度種植(銀)Overplanting …
阿新 • • 發佈:2018-10-31
Description
在一個笛卡爾平面座標系裡(則X軸向右是正方向,Y軸向上是正方向),有\(N(1<=N<=1000)\)個矩形,第i個矩形的左上角座標是\((x1, y1)\),右下角座標是\((x2,y2)\)。問這\(N\)個矩形所覆蓋的面積是多少?注意:被重複覆蓋的區域的面積只算一次。
Input
第一行,一個整數N。 \((1<=N<=1000)\)。
接下來有\(N\)行,每行描述一個矩形的資訊,分別是矩形的\(x1、y1、x2、y2\)。
其中 \(−10^8<=x1,y1,x2,y2<=10^8\)。
Ouput
一個整數,被N個矩形覆蓋的區域的面積。
難得遇到一個裸的掃描線的題,竟然沒切掉 emmm.
看到\(x,y\)的座標範圍,離散化就好了!
沒有一遍切,竟然是沒開\(long \ \ long\)!!!
太難受了,關於這個的話就不多BB,網上講解很多.
大家可以去搜一下。(貌似NOIP不會考,暫且學了)
將來有時間寫講解好了 qwq.
程式碼
#include<cstdio> #include<algorithm> #include<iostream> #define int long long #define R register using namespace std; const int gz=10086; inline void in(int &x) { int f=1;x=0;char s=getchar(); while(!isdigit(s)){if(s=='-')f=-1;s=getchar();} while(isdigit(s)){x=x*10+s-'0';s=getchar();} x*=f; } struct cod { int l,r,h; int f; bool operator <(const cod&a)const { return h<a.h; } }edge[gz]; struct tre { int l,r,s; int len; }tr[gz]; #define ls o<<1 #define rs o<<1|1 int x[gz],n,tot; void build(R int o,R int l,R int r) { tr[o].l=l;tr[o].r=r; if(l==r)return; R int mid=(l+r)>>1; build(ls,l,mid); build(rs,mid+1,r); } inline void up(R int o) { if(tr[o].s) tr[o].len=x[tr[o].r+1]-x[tr[o].l]; else if(tr[o].l==tr[o].r) tr[o].len=0; else tr[o].len=tr[ls].len+tr[rs].len; } void change(R int o,R int l,R int r,R int del) { if(tr[o].l==l and tr[o].r==r) { tr[o].s+=del; up(o); return; } R int mid=(tr[o].l+tr[o].r)>>1; if(r<=mid) change(ls,l,r,del); else if(l>mid) change(rs,l,r,del); else change(ls,l,mid,del),change(rs,mid+1,r,del); up(o); } signed main() { in(n); for(R int i=1;i<=n;i++) { R int x1,x2,y1,y2; in(x1),in(y1),in(x2),in(y2); edge[++tot].l=x1;edge[tot].r=x2;edge[tot].f=-1; edge[tot].h=y1;x[tot]=x1; edge[++tot].l=x1;edge[tot].r=x2;edge[tot].f=1; edge[tot].h=y2;x[tot]=x2; } sort(edge+1,edge+tot+1); sort(x+1,x+tot+1); int new_n=1; for(R int i=2;i<=tot;i++) if(x[new_n]!=x[i])x[++new_n]=x[i]; build(1,1,new_n); int ans=0; for(R int i=1;i<=tot;i++) { R int l=lower_bound(x+1,x+new_n+1,edge[i].l)-x; R int r=lower_bound(x+1,x+new_n+1,edge[i].r)-x-1; change(1,l,r,edge[i].f); ans+=(edge[i+1].h-edge[i].h)*tr[1].len; } printf("%lld",ans); }