luogu P3994 高速公路

背景：

逃了物理竞赛，机房就$2$人。

题目传送门：

https://www.luogu.org/problemnew/show/P3994

题意：

有一棵树，其中第$i$个点到第$j$个点的贡献为$dis(i,j)\ast P_i+Q_i$.若$i$是$j$的祖先，则保证$P_i\le P_j$。

思路：

显然是一个在树上维护的斜率优化的$DP$。
对于一条根到叶子节点的路径的路径，它的斜率一定是单调的。那就好办了，因为它们就可以放在一起维护了。而对于不在同一棵子树内的儿子节点，回溯时注意单调队列内的点被删除的放回即可。

设$f_i$表示$i$点到$1$号根节点的路径的最小值。
则$f_i={\sum }_{j\in i.ancestor}{f}_j+P_i\ast dis(i,j)+Q_i$。
由于$i,j$在同一条链上，因此$dis(i,j)=dis(i,1)-dis(j-1)$。
设$di{s}_i=di{s}_{1,i}$，则：$f_i=f_j+{\sum }_{j.ancestor}{P}_i(di{s}_i-di{s}_j)+Q_i$。
得到斜率方程：$\frac{f_j-f_k}{di{s}_j-di{s}_k}<p_i$。
套进去即可。

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#define LL long long
using namespace std;
	int n,m,len=0,head,tail;
	struct node{int x,y,z,next;} a[2000010];
	int last[1000010];
	LL p[1000010],q[1000010],dis[1000010],f[1000010];
	int que[1000010];
LL calc1(int x,int y)
{
	return dis[y]-dis[x];
}
LL calc2(int x,int y)
{
	return f[y]-f[x];
}
double calc(int x,int y)
{
	return (double)(calc2(x,y))/(double)(calc1(x,y));
}
void ins(int x,int y,int z)
{
	a[++len]=(node){x,y,z,last[x]}; last[x]=len;
}
int last_head[1000010],last_tail[1000010],tmp[1000010],pos[1000010];
void dfs(int x,int fa)
{
	last_head[x]=head,last_tail[x]=tail;
	while(head<tail&&calc(que[head],que[head+1])<=(double)p[x]) head++;
	f[x]=f[que[head]]+p[x]*calc1(que[head],x)+q[x];
	while(head<=tail&&calc(que[tail],x)<=calc(que[tail-1],que[tail])) tail--;
	tail++;
	tmp[x]=que[tail];
	pos[x]=tail;
	que[tail]=x;
	for(int i=last[x];i;i=a[i].next)
	{
		int y=a[i].y;
		if(y==fa) continue;
		dis[y]=dis[x]+a[i].z;
		dfs(y,x);
	}
	que[pos[x]]=tmp[x];
	head=last_head[x],tail=last_tail[x];
}
int main()
{
	int x,z;
	scanf("%d",&n);
	for(int i=2;i<=n;i++)
	{
		scanf("%d %d %lld %lld",&x,&z,&p[i],&q[i]);
		ins(i,x,z),ins(x,i,z);
	}
	head=1;
	tail=0;
	que[0]=0;
	dfs(1,0);
	for(int i=2;i<=n;i++)
		printf("%lld\n",f[i]);
}