【算法学习】悬线法

合集 - 【算法学习】(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-04

23.【算法学习】悬线法2024-10-04

24.欧几里得算法与 EX2024-09-2325.【算法学习】逆元与求解2024-09-2326.【算法学习】费马定理2024-09-2327.三分2024-09-2328.裴蜀定理2024-09-2329.【算法学习】欧拉函数φ2024-09-2330.【算法学习】二维转一维问题2024-11-0331.【算法学习】扫描线2024-11-0132.【算法学习】矩阵乘法2024-10-3133.【算法学习】同余最短路2024-11-0434.【算法学习】组合数学2024-11-0535.【算法学习】反悔贪心2024-11-06

收起

简单介绍

学习笔记

悬线法，相当于有一个限高绳，向左向右找到不低于这个高度的左右边界。

例题

SP1805 例题

分类讨论：

• 当 \(i=1\)，到达边界停止。

• 当 \(a_i>a_{i-1}\)，低于高度，停止拓展。

• 当 \(a_i<=a_{i-1}\)，可以扩展，直接继承 \(l_i=l_{l_{i-1}}\)。

相同的求右端点，最后求最大矩阵面积就是左右边界距离乘高度。

点击查看代码

#include <bits/stdc++.h>
#define ll long long
using namespace std;
const int N=1e5+10;
int n;
 
int a[N];
int l[N],r[N];
 
int main(){
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    while(1){
    	cin>>n;
    	if(n==0){
    		break;
		}
		for(int i=1;i<=n;i++){
			cin>>a[i];
			l[i]=i;
			r[i]=i;
		}
		for(int i=1;i<=n;i++){
			while(l[i]-1>=1&&a[i]<=a[l[i]-1]){
				l[i]=l[l[i]-1];
			}
		}
		for(int i=n;i>=1;i--){
			while(r[i]+1<=n&&a[i]<=a[r[i]+1]){
				r[i]=r[r[i]+1];
			}
		}
		ll ans=0;
		for(int i=1;i<=n;i++){
			ans=max(ans,(ll)(r[i]-l[i]+1)*a[i]);
		}
		cout<<ans<<"\n";
	}
    return 0;
}

再给一个例题，P4147 玉蟾宫。

看到这题是不是还是没啥思路，我们这题也用悬线法，对每层只要一直是 F 就增加高度否则归零，然后每层左右更新边界，计算面积即可。

配图就比较形象了。

同样代码再来一次

#include <bits/stdc++.h>
#define ll long long
using namespace std;
const int N=1e5+10;
int n,m;
int a[N];
int l[N],r[N];
int ans;
int main(){
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    cin>>n>>m;
    for(int i=1;i<=n;i++){
    	for(int j=1;j<=m;j++){
    		l[j]=r[j]=j;
		}
		char c;
		for(int j=1;j<=m;j++){
			cin>>c;
			if(c=='R'){
				a[j]=0;
			}
			else if(c=='F'){
				a[j]++;
			}
		}
		for(int j=1;j<=m;j++){
			while(l[j]!=1&&a[j]<=a[l[j]-1]){
				l[j]=l[l[j]-1];
			}
		}
		for(int j=m;j>=1;j--){
			while(r[j]!=m&&a[j]<=a[r[j]+1]){
				r[j]=r[r[j]+1];
			}
		}
		for(int j=1;j<=m;j++){
			ans=max(ans,a[j]*(r[j]-l[j]+1));
		}
	}
    cout<<ans*3;
    return 0;
}