[Noi2014]购票

来自FallDream 的博客，未经允许，请勿转载，谢谢。

---

今年夏天，NOI在SZ市迎来了她30周岁的生日。来自全国 n 个城市的OIer们都会从各地出发，到SZ市参加这次盛会。

       全国的城市构成了一棵以SZ市为根的有根树，每个城市与它的父亲用道路连接。为了方便起见，我们将全国的 n 个城市用 1 到 n 的整数编号。其中SZ市的编号为 1。对于除SZ市之外的任意一个城市 v，我们给出了它在这棵树上的父亲城市 fv  以及到父亲城市道路的长度 sv。

从城市 v 前往SZ市的方法为：选择城市 v 的一个祖先 a，支付购票的费用，乘坐交通工具到达 a。再选择城市 a 的一个祖先 b，支付费用并到达 b。以此类推，直至到达SZ市。

对于任意一个城市 v，我们会给出一个交通工具的距离限制 lv。对于城市 v 的祖先 a，只有当它们之间所有道路的总长度不超过 lv  时，从城市 v 才可以通过一次购票到达城市 a，否则不能通过一次购票到达。对于每个城市 v，我们还会给出两个非负整数 pv,qv  作为票价参数。若城市 v 到城市 a 所有道路的总长度为 d，那么从城市 v 到城市 a 购买的票价为 dpv+qv。

每个城市的OIer都希望自己到达SZ市时，用于购票的总资金最少。你的任务就是，告诉每个城市的OIer他们所花的最少资金是多少。

n<=2*10^5

dp式子比较容易列出来 f[i]=min f[j]+q[i]+p[i]*(d[i]-d[j)

考虑斜率优化，发现j比k优，当且仅当p<(f[j]-f[k])/(d[j]-d[k]) 所以考虑维护凸包

可以直接线段树+无旋treap 大力维护凸包

也可以cdq分治上树，每次像点分治那样做，然后用包含根的那个子树更新其它子树。

更新时候按照深度(用于更新的点直接放深度，否则放深度减去limit)排序从下往上做

复杂度nlogn^2

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
#define MN 200000
#define INF 2000000000
using namespace std;
inline ll read()
{
    ll x =  , f = ; char ch = getchar();
    while(ch < '' || ch > ''){ if(ch == '-') f = -;  ch = getchar();}
    while(ch >= '' && ch <= ''){x = x *  + ch - '';ch = getchar();}
    return x * f;
}
 
ll q[MN+],p[MN+],l[MN+],dep[MN+],f[MN+];
int n,head[MN+],cnt=,top,size[MN+],b[MN+],fa[MN+],mx[MN+],Mx[MN+];
struct edge{int to,next;ll w;}e[MN*+];
inline double Getslop(int x,int y){return !y?-INF:(double)(f[y]-f[x])/(dep[y]-dep[x]);}
struct data{int x,kind;ll dep;}a[MN+];
struct CH
{
    int s[MN+],tp;
    void clear() {tp=;}
    void ins(int x)
    {
        while(tp>&&Getslop(s[tp],s[tp-])<Getslop(x,s[tp])) --tp;
        s[++tp]=x;
    }
    int query(ll p)
    {
        int l=,r=tp,mid,ans;s[tp+]=;
        while(l<=r)
        {
            mid=l+r>>;
            if(Getslop(s[mid],s[mid+])<(double)p) ans=mid,r=mid-;
            else l=mid+;
        }
        return s[ans];
    }
}H;
inline void ins(int f,int t,ll w){e[++cnt]=(edge){t,head[f],w};head[f]=cnt;}
 
void Pre(int x)
{
    size[x]=;mx[x]=;
    for(int i=head[x];i;i=e[i].next)
        dep[e[i].to]=dep[x]+e[i].w,Pre(e[i].to),size[x]+=size[e[i].to],mx[x]=max(mx[x],size[e[i].to]);
}
 
void GetRt(int x,int tot,int &rt)
{
    size[x]=;mx[x]=;
    for(int i=head[x];i;i=e[i].next)
        if(!b[e[i].to]) GetRt(e[i].to,tot,rt),size[x]+=size[e[i].to],mx[x]=max(mx[x],size[e[i].to]);
    Mx[x]=max(mx[x],tot-size[x]);
    if(size[x]>&&Mx[x]<Mx[rt]) rt=x;
}
 
void Dfs(int x)
{
    a[++top]=(data){x,,dep[x]-l[x]};
    for(int i=head[x];i;i=e[i].next)
        if(!b[e[i].to]) Dfs(e[i].to);
}
bool cmp(data x,data y){return x.dep==y.dep?x.kind>y.kind:x.dep>y.dep;}
ll GetAns(int x,int y){return f[y]+p[x]*(dep[x]-dep[y])+q[x];}
void Solve(int x,int Sz)
{
    if(Sz<=) return;int rt=;
    GetRt(x,Sz,rt);
    for(int i=head[rt];i;i=e[i].next) b[e[i].to]=;
    Solve(x,Sz-size[rt]+);
    a[top=]=(data){rt,,dep[rt]};
    for(int t=fa[rt];t!=fa[x];t=fa[t]) a[++top]=(data){t,,dep[t]};
    for(int i=head[rt];i;i=e[i].next) Dfs(e[i].to);
    sort(a+,a+top+,cmp);H.clear();
    for(int i=;i<=top;++i)
        if(a[i].kind) H.ins(a[i].x);
        else if(H.tp) f[a[i].x]=min(f[a[i].x],GetAns(a[i].x,H.query(p[a[i].x])));
    for(int i=head[rt];i;i=e[i].next) Solve(e[i].to,size[e[i].to]);
} 
 
int main()
{
    n=read();read();
    memset(f,,sizeof(f));
    for(int i=;i<=n;++i)
    {
        fa[i]=read();ll w=read();
        p[i]=read(),q[i]=read();l[i]=read();
        ins(fa[i],i,w);
    }
    Pre();f[]=;Mx[]=INF;Solve(,size[]);
    for(int i=;i<=n;++i) printf("%lld\n",f[i]);
    return ;
}