作爲一道模板題讓我卡了一個月……
對於線段樹+離散化新手而言這實在是太難了……
有關離散化:
可以查看這一篇文章:https://www.jianshu.com/p/9347659dcf18
一種縮小數據範圍的小技巧,在本題中用於縮短線段樹的區間大小
基本思路
想象一條直線,從左往右掃,我們要計算的就是直線掃過的面積。
很容易可以想出把整個圖形按照每一條豎線切割,分割成幾個矩形。
具體做法可以看代碼。
#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;
long long tag[1600005],tree[1600005],cnt,ans,n,x,x2,y,y2,line[1600005];
struct data
{
long long x,l,r,flag;
}edge[1600005];
bool comp(data a,data b)
{
if (a.x!=b.x)return a.x<b.x;
else return a.l<b.l;
}
/*void push_down(int root,int l,int r)
{
int lson=root<<1;
int rson=root<<1|1;
int mid=(l+r)>>1;
tag[lson]+=tag[root];
if (tag[lson]) tree[lson]=line[mid]-line[l];
tag[rson]+=tag[root];
if (tag[rson]) tree[rson]=line[r]-line[mid];
tag[root]=0;
}*/
void add(int root,int l,int r,long long al,long long ar,int num)
{
if (al<=line[l]&&line[r]<=ar)
{
tag[root]+=num;
if (!tag[root]) tree[root]=tree[root<<1]+tree[root<<1|1];
if (tag[root]) tree[root]=line[r]-line[l];
return ;
//tag存儲節點代表的區間被覆蓋的次數
//tree存儲節點代表的區間被覆蓋的長度
//值得注意的是,即使當前區間沒有被完全覆蓋,tree[root]也不一定爲0。
}
if (l+1==r) return ;
if (ar<=line[l]||al>=line[r]) return ;//加入等號是因爲若只有端點相同,實質上也是沒有交集的。
//push_down(root,l,r);
//無需pushdown的原因爲:對於任意一個區間,完全覆蓋它的豎線一定是成對的(因爲有左右兩邊),即有進有出,不可能會出現負值。
//並且對於任意區間的子區間,子區間是否被覆蓋或者被覆蓋的長度對父區間毫無影響
//如果加入pushdown,會導致超時。
int mid=(l+r)>>1;
if (al<line[mid]) add(root<<1,l,mid,al,ar,num);
if (ar>line[mid]) add(root<<1|1,mid,r,al,ar,num);
if (tag[root]==0) tree[root]=tree[root<<1]+tree[root<<1|1];
else tree[root]=line[r]-line[l];
}
int main()
{
scanf("%lld",&n);
for (int i=1;i<=n;i++)
{
scanf("%lld%lld%lld%lld",&x,&y,&x2,&y2);
edge[++cnt].x=x;
edge[cnt].l=y;
edge[cnt].r=y2;
edge[cnt].flag=1;
line[cnt]=y;
edge[++cnt].x=x2;
edge[cnt].l=y;
edge[cnt].r=y2;
edge[cnt].flag=-1;
line[cnt]=y2;
//edge記錄豎線的位置
//line記錄每一個點的位置,方便離散化處理
}
sort(edge+1,edge+1+cnt,comp);
sort(line+1,line+1+cnt);
int newend=unique(line+1,line+1+cnt)-line-1;//去重函數,此時newend恰好是line最後一個元素的下標
for (int i=1;i<cnt;i++)
{
add(1,1,newend,edge[i].l,edge[i].r,edge[i].flag);
ans+=tree[1]*(edge[i+1].x-edge[i].x);//整個區間被覆蓋的總長度*到下一條豎邊的距離,得到這一部分的面積
//最後一條豎邊不需要處理
}
printf("%lld",ans);
return 0;
}