【算法学习】扫描线

合集 - 【算法学习】(35)

1.【算法学习】排序2024-08-272.【算法学习】Manacher 马拉车2024-08-293.【算法学习】KMP 算法2024-09-014.LCA 最近公共祖先（树链和倍增）这次真有树链了！！！2024-08-305.线段覆盖问题2024-08-296.【算法学习】学换根dp有感2024-09-207.二分图最大匹配2024-09-178.【算法学习】01BFS2024-09-179.洛谷 P1892 [BOI2003] 团伙 种类并查集 扩展域并查集2024-09-1110.【算法学习】高斯消元法2024-09-1011.贝叶斯公式2024-09-0712.背包2024-09-0513.【算法学习】模拟退火2024-09-0514.【算法学习】基环树2024-10-2915.【算法学习】树链部分2024-10-2916.【算法学习】莫队2024-10-2917.【算法学习】分块九讲2024-10-2818.平衡树2024-10-2319.圆方树2024-10-2220.【算法学习】点分治2024-10-1121.公式2024-10-0922.【算法学习】笛卡尔树2024-10-0423.【算法学习】悬线法2024-10-0424.欧几里得算法与 EX2024-09-2325.【算法学习】逆元与求解2024-09-2326.【算法学习】费马定理2024-09-2327.三分2024-09-2328.裴蜀定理2024-09-2329.【算法学习】欧拉函数φ2024-09-2330.【算法学习】二维转一维问题2024-11-03

31.【算法学习】扫描线2024-11-01

32.【算法学习】矩阵乘法2024-10-3133.【算法学习】同余最短路2024-11-0434.【算法学习】组合数学2024-11-0535.【算法学习】反悔贪心2024-11-06

收起

这篇题解写的难以言喻，可能只有我能看的懂！

前言

虽然我觉得这个算法目前不太可能会考，但是我觉得挺有意思的，而且学个算法也挺好，我是为自己学的！！！

定义

扫描线可以求二维图形的面积，也可以求周长等多种用途……

P5490 【模板】扫描线 & 矩形面积并

这就是扫描线的模板题，我们有这么一个图形我们按照他的端点划分成许多个线段。

这个概括动图了扫描线的全过程。

我们将端点划分成线段，将线段抽象线段树的节点，我们用线段树维护线段的出现次数和线段长度。

然后我们将读入的矩阵下线段标为1，上线段标为-1，我们这一个节点（线段）的出现次数为0时是不能统计答案的。

我们按照高度排序是为了计算高度来计算面积，坐标都很大，我们需要离散化一下来缩小范围。

然后对于新来的线段我们要我们要更新目前整个线对应到矩阵上的长度，然后再对高度计算答案。

我们选择用左端点来对应这条线段，右端点也要参与计算进来。这里我们通过更新我们线扫到的长度。

然后就可以了，我全篇不知道在讲什么，就这样吧不管了。

极速版：

#include <bits/stdc++.h>
#define int long long
#define ls p<<1
#define rs p<<1|1 
#define fi first
#define se second
#define re register 
#define pir pair<int,int>
const int inf=1e9;
const int mod=998244353;
const int N=1e6+10;
using namespace std;
 
int n;
 
int dc[N];
 
struct ss{
	int l,r,h,mark;
}a[N];
 
bool cmp(ss g,ss h){
	return g.h<h.h;
}
 
int sum[N<<2];
int len[N<<2];
 
void pushup(int p,int l,int r){
	if(sum[p]){
		len[p]=dc[r+1]-dc[l];
	}
	else{
		len[p]=len[ls]+len[rs];
	}
}
 
void change(int p,int pl,int pr,int l,int r,int c){
	if(dc[pl]>=r||dc[pr+1]<=l){
		return;
	} 
	if(l<=dc[pl]&&dc[pr+1]<=r){
		sum[p]+=c;
		pushup(p,pl,pr);
		return;
	}
	int mid=(pl+pr)>>1;
	change(ls,pl,mid,l,r,c);
	change(rs,mid+1,pr,l,r,c);
	pushup(p,pl,pr); 
}
 
signed main(){
    ios::sync_with_stdio(false);
    cin.tie(nullptr); 
	
	cin>>n;
	
	for(int i=1;i<=n;i++){
		int x,y,xx,yy;
		cin>>x>>y>>xx>>yy;
		dc[i*2-1]=x;
		dc[i*2]=xx;
		a[i*2-1]={x,xx,y,1};
		a[i*2]={x,xx,yy,-1};
	}
	n<<=1;
	sort(a+1,a+n+1,cmp);
	sort(dc+1,dc+n+1);
	int tot=unique(dc+1,dc+n+1)-dc-1;
	int ans=0;
	for(int i=1;i<n;i++){
		change(1,1,tot-1,a[i].l,a[i].r,a[i].mark);
		ans+=len[1]*(a[i+1].h-a[i].h);
	}
	cout<<ans;
    return 0;
}