[APIO2016]烟火表演

2021-06-20 01:04

阅读:586

标签:image   color   getch   www   wap   stream   +=   int   描述   

题目描述

https://www.lydsy.com/JudgeOnline/problem.php?id=4585

题解

这题太神了。

我们可以先列出一个dp方程,dp[x][d]表示x节点到所有叶子的距离的d时的代价。

结论1:对于每个点来说,这个dp数组为二维平面上是一个下凸函数。

证明:对于叶子来说一定成立,在w[x]处为0,然后小于w[x]的部分斜率为-1,大与w[x]的部斜率为1。

对于非叶子节点 ,它的函数时有儿子们加起来的,也成立。

然后我们考虑一个点如何向父亲转移。

这是要分四种情况,假设当前节点的斜率为0的部分为L~R。

x

因为边权不能为负,所以我们只能从小的地方往大的地方转移。

技术分享图片

比如这个蓝点,它只能从它以及前面的地方转移,但从自己转移时最优的。

x>=L&&x

在L处答案为f(L)+w,没往右动一格代价会-1。

x>=L+w[u]&&x

这个相当于直接转移了,没有代价。

x>=R+w[u] f‘(x)=f(x)+(x-R)-w

相当于是走过了,会产生代价。

我们发现转移大概长这样(继续盗图)。

技术分享图片

我们把左边的点向上动一段后插入斜率为-1的线,再把斜率为0的部分向右平移,最后面是斜率为1的线。

然后这种操作就可以维护了。

思路大概就是只维护拐点,用一个可并堆,每次把右边的部分弹掉。

结论2:每次合并到一个非叶子节点时,斜率为0的线右边有n个点,n为该点的儿子数。

证明:因为最后那一块斜率一定为n(n个斜率为1的直线相加),然后往前经过一个拐点,斜率会-1。

然后把LR取出来做平移,再插回去,最后一直合并到1。

结论3:每经过一个拐点,斜率-1(其实它和上面的结论一个意思)。

然后就可以利用最后的拐点直接算答案了。

怎么算呢?我们把斜率>=0的直线弹掉,把前面的所有直线斜率+1,那么每条直线都会产生deltax的代价,总共有xn的代价,那我们就在答案里+xn,此时最后一条直线斜率为0,把它删掉。

然后一直做,最后得到的是f(0)-f(min),f(0)就是所有边权之和,那么f(min)就可以求出来了。

代码

#include
#include#define N 600002
using namespace std;
typedef long long ll;
int deep[N],d[N],n,m,fa[N],T[N],tot;
ll sum,w[N];
inline ll rd(){
    ll x=0;char c=getchar();bool f=0;
    while(!isdigit(c)){if(c==-)f=1;c=getchar();}
    while(isdigit(c)){x=(x1)+(x3)+(c^48);c=getchar();}
    return f?-x:x;
}
struct tr{
    ll v;int l,r;
}tr[N];
int merge(int x,int y){
    if(!x||!y)return x^y;
    if(tr[x].vtr[y].v)swap(x,y);
    tr[x].r=merge(tr[x].r,y);
    if(deep[tr[x].l]deep[tr[x].r])swap(tr[x].l,tr[x].r);
    deep[x]=deep[tr[x].r]+1;
    return x;
}
inline int pop(int x){return merge(tr[x].l,tr[x].r);}
int main(){
    n=rd();m=rd();
    for(int i=2;ii){
        fa[i]=rd();w[i]=rd();sum+=w[i];d[fa[i]]++;
    }
    for(int i=n+m;i>=2;--i){
        ll l=0,r=0;
        if(in){
            while(--d[i])T[i]=pop(T[i]);
            l=tr[T[i]].v;T[i]=pop(T[i]);
            r=tr[T[i]].v;T[i]=pop(T[i]);
        }
        tr[++tot].v=l+w[i];tr[++tot].v=r+w[i];
        T[i]=merge(T[i],merge(tot-1,tot));
        T[fa[i]]=merge(T[fa[i]],T[i]);
    }
    while(d[1]--)T[1]=pop(T[1]);
    while(T[1]){sum-=tr[T[1]].v;T[1]=pop(T[1]);}
    coutsum;
    return 0;
}

 

[APIO2016]烟火表演

标签:image   color   getch   www   wap   stream   +=   int   描述   

原文地址:https://www.cnblogs.com/ZH-comld/p/10269171.html


评论


亲,登录后才可以留言!