Scaffolding

最新推荐文章于 2022-11-12 15:52:23 发布

StaroForgin

最新推荐文章于 2022-11-12 15:52:23 发布

阅读量253

点赞数 1

分类专栏： ------图论------ # 树 # 二叉搜索树文章标签：算法动态规划

本文链接：https://blog.csdn.net/Tan_tan_tann/article/details/120047362

版权

------图论------ 同时被 3 个专栏收录

101 篇文章 0 订阅

订阅专栏

树

57 篇文章 0 订阅

订阅专栏

二叉搜索树

7 篇文章 0 订阅

订阅专栏

这篇博客探讨了一种使用Treap数据结构解决在有限材料情况下搭建多层建筑的最小运输次数优化问题。文章详细介绍了如何通过广义笛卡尔树的优化方法，将问题转化为Treap的构造和动态规划转移，最终实现线性时间复杂度的解决方案。

摘要由CSDN通过智能技术生成

Scaffolding

题解

考虑我们怎样放才是最优的放法。
很明显，当我们上楼后是不能带着材料下来的，所以我们带上楼的材料最多也只是刚好够放完上面一层，剩下的材料我们都在下一层中分完后再拿到上面去。
但我们该怎么分层呢？
很明显，对于一个完整的矩形我们是可以将它完整地分为一层的，我们对于这种完整的矩形我们在哪个位置上楼并不会影响我们的答案。
所以，我们要做的是将这样的柱状图分为多个完整的矩形处理。
观察到 $n\leqslant 10^5$ ，广义笛卡尔树的 $O\left(n^2\right)$ 暴力建图是明显不行的，考虑用 $t r e a p$ 的方法建图，这部分是 $O\left(n\right)$ 的。

在建图的同时，我们可以在树上进行 $d p$ 转移。
由于我们建完上面的楼后可能会剩一部分材料，而这部分材料只能拿到下面一层，所以我们肯定是从从叶子往根 $d p$ ，某个节点有多的材料，人给它父亲。
我们记节点 $i$ 代表的矩阵的大小为 $sum_{i}$ ， $tot_{i}$ 表示 $i$ 的子树的总矩阵大小， $f_{i}$ 表示建完点 $i$ 的子树需要的最少运送次数， $g_{i}$ 表示建完 $i$ 的子树可以给它的祖先节点留下多少材料。
容易得到 $d p$ 转移式：
$f_{u}=\sum_{v\in son_{u}}f_{v}+\max(0,\left\lceil\frac{sum_{u}-\sum_{v\in son_{u}}g_{v}}{m}\right\rceil)$
$g_{u}=mf_{u}-tot_{u}$
这些转移都可以在遍历整棵树时进行，树在我们输入的时候就可以建好了。

时间复杂度 $O\left(n\right)$ 。

源码

#include<cstdio>
#include<cmath>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
#define MAXN 100005
#define lowbit(x) (x&-x)
#define reg register
#define pb push_back
#define mkpr make_pair
#define fir first
#define sec second
typedef long long LL;
#define int LL
typedef unsigned long long uLL;
const int INF=0x3f3f3f3f;
const int mo=1e9+7;
const int inv2=499122177;
const int jzm=2333;
const int zero=10000;
const int lim=30000;
const int orG=3,invG=332748118;
const double Pi=acos(-1.0);
const double eps=1e-5;
typedef pair<int,int> pii;
template<typename _T>
_T Fabs(_T x){return x<0?-x:x;}
template<typename _T>
void read(_T &x){
	_T f=1;x=0;char s=getchar();
	while(s>'9'||s<'0'){if(s=='-')f=-1;s=getchar();}
	while('0'<=s&&s<='9'){x=(x<<3)+(x<<1)+(s^48);s=getchar();}
	x*=f;
}
template<typename _T>
void print(_T x){if(x<0){x=(~x)+1;putchar('-');}if(x>9)print(x/10);putchar(x%10+'0');}
LL gcd(LL a,LL b){return !b?a:gcd(b,a%b);}
int add(int x,int y,int p){return x+y<p?x+y:x+y-p;}
void Add(int &x,int y,int p){x=add(x,y,p);}
int qkpow(int a,int s,int p){int t=1;while(s){if(s&1LL)t=1ll*a*t%p;a=1ll*a*a%p;s>>=1LL;}return t;}
int n,m,h[MAXN],lson[MAXN],rson[MAXN],sta[MAXN],siz[MAXN],stak;
LL f[MAXN],g[MAXN],sum[MAXN];
void sakura(int rt,int minn){
	int sumg=0;siz[rt]++;
	if(lson[rt])sakura(lson[rt],h[rt]),siz[rt]+=siz[lson[rt]],
		f[rt]+=f[lson[rt]],sumg+=g[lson[rt]],sum[rt]+=sum[lson[rt]];
	if(rson[rt])sakura(rson[rt],h[rt]),siz[rt]+=siz[rson[rt]],
		f[rt]+=f[rson[rt]],sumg+=g[rson[rt]],sum[rt]+=sum[rson[rt]];
	int tmp=siz[rt]*(h[rt]-minn)-sumg;sum[rt]+=(h[rt]-minn)*siz[rt]; 
	if(tmp>0)f[rt]+=(tmp+m-1LL)/m;
	g[rt]=1ll*m*f[rt]-sum[rt];
}
signed main(){
	while(scanf("%lld %lld",&n,&m)!=EOF){
		for(int i=1;i<=n;i++){
			read(h[i]);bool flag=0;
			while(stak&&h[sta[stak]]>h[i])stak--,flag=1;
			if(sta[stak])lson[sta[stak]]=i;
			if(flag)rson[i]=sta[stak+1];sta[++stak]=i;
		}
		sakura(sta[1],0);
		printf("%lld\n",f[sta[1]]);stak=0;
		for(int i=1;i<=n;i++)lson[i]=siz[i]=rson[i]=f[i]=g[i]=sum[i]=0;
	}
	return 0;
}