【Luogu P5490】掃描線

Luogu P5490

作爲一道模板題讓我卡了一個月……

對於線段樹+離散化新手而言這實在是太難了……

有關離散化:

可以查看這一篇文章: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;
}
發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章