ZJOI2015 幻想乡战略游戏 动态点分治_树链剖分_未调完

Description

傲娇少女幽香正在玩一个非常有趣的战略类游戏，本来这个游戏的地图其实还不算太大，幽香还能管得过来，但是不知道为什么现在的网游厂商把游戏的地图越做越大，以至于幽香一眼根本看不过来，更别说和别人打仗了。 在打仗之前，幽香现在面临一个非常基本的管理问题需要解决。 整个地图是一个树结构，一共有n块空地，这些空地被n-1条带权边连接起来，使得每两个点之间有一条唯一的路径将它们连接起来。在游戏中，幽香可能在空地上增加或者减少一些军队。同时，幽香可以在一个空地上放置一个补给站。 如果补给站在点u上，并且空地v上有dv个单位的军队，那么幽香每天就要花费dv×dist(u,v)的金钱来补给这些军队。由于幽香需要补给所有的军队，因此幽香总共就要花费为Sigma(Dv*dist(u，v),其中1<=V<=N)的代价。其中dist(u,v)表示u个v在树上的距离（唯一路径的权和）。 因为游戏的规定，幽香只能选择一个空地作为补给站。在游戏的过程中，幽香可能会在某些空地上制造一些军队，也可能会减少某些空地上的军队，进行了这样的操作以后，出于经济上的考虑，幽香往往可以移动他的补给站从而省一些钱。但是由于这个游戏的地图是在太大了，幽香无法轻易的进行最优的安排，你能帮帮她吗？ 你可以假定一开始所有空地上都没有军队。

题解：

最近写动态点分治快自闭了QAQ.......

有时间再来写一下题解.

其实我不会告诉你我还没调出来呢QAQ

Code:

// luogu-judger-enable-o2
// luogu-judger-enable-o2
#include <bits/stdc++.h>
#define setIO(s) freopen(s".in","r",stdin) // ,freopen(s".out","w",stdout)
#define maxn 300000
#define N 300010
#define ll long long
#define inf 0x7f7f7f
using namespace std;
int n,m,cnt,sn,root;
int hd[N],to[N],nx[N],val[N];
void add(int u,int v,int c)
{
    nx[++cnt] = hd[u],hd[u] = cnt,to[cnt] = v,val[cnt] = c;
}
namespace heavyedge
{
    int fa[N],hson[N],siz[N],top[N],dep[N];
    ll dis[N];
    void dfs1(int u,int ff)
    {
        fa[u] = ff;
        for(int i = hd[u]; i ; i = nx[i])
        {
            if(to[i] == ff) continue;
            dep[to[i]] = dep[u] + 1,dis[to[i]] = dis[u] + val[i];
            dfs1(to[i],u);
            siz[u] += siz[to[i]];
            if(siz[to[i]] > siz[hson[u]]) hson[u] = to[i];
        }
    }
    void dfs2(int u,int tp)
    {
        top[u] = tp;
        if(hson[u]) dfs2(hson[u],tp);
        for(int i = hd[u]; i ; i = nx[i])
        {
            if(to[i] == hson[u] || to[i] == fa[u]) continue;
            dfs2(to[i],to[i]);
        }
    }
    int LCA(int a,int b)
    {
        while(top[a] != top[b]) dep[top[a]] < dep[top[b]] ? b = fa[top[b]]: a = fa[top[a]];
        return dep[a] < dep[b] ? a : b;
    }
    int main()
    {
        dfs1(1,0),dfs2(1,1);
        return 0;
    }
};
ll Dis(int a,int b)
{
    return heavyedge::dis[a] + heavyedge::dis[b] - (heavyedge::dis[heavyedge::LCA(a,b)] << 1);
}
int siz[N],vis[N],Fa[N],hto[N],f[N],tmp;
ll sumv[N],sumd[N];
void GetRoot(int u,int ff)
{
    siz[u] = 1,f[u] = 0;
    for(int i = hd[u]; i ; i = nx[i])
    {
        if(vis[to[i]] || to[i] == ff) continue;
        GetRoot(to[i],u);
        siz[u] += siz[to[i]];
        f[u] = max(f[u],siz[to[i]]);
    }
    f[u] = max(f[u],sn - siz[u]);
    if(f[root] > f[u]) root = u;
}
void dfs(int u)
{
    vis[u] = 1;
    for(int i = hd[u]; i ; i = nx[i])
    {
        if(vis[to[i]]) continue;
        root = 0, sn = siz[to[i]], GetRoot(to[i],u), hto[i] = root,Fa[root] = u, dfs(root);
    }
}
#define fax(i) (i + n)
void Update(int x,int w)
{
    sumv[x] += w;
    for(int i = x; Fa[i]; i = Fa[i])
    {
        ll dis = Dis(x, Fa[i]);
        sumv[Fa[i]] += w;
        sumd[Fa[i]] += 1ll*dis * w;
        sumd[fax(i)] += 1ll*dis * w;
    }
}
ll calc(int x)
{
    ll res = sumd[x];
    for(int i = x; Fa[i]; i = Fa[i])
    {
        ll dis = Dis(x, Fa[i]);
        res += 1ll*(sumv[Fa[i]] - sumv[i]) * dis;
        res += sumd[Fa[i]] - sumd[fax(i)];
    }
    return res;
}
ll Query(int u)
{
    ll tmp = calc(u);
    for(int i = hd[u]; i ; i = nx[i])
    {
        if(calc(to[i]) < tmp) return Query(to[i]);
    }
    return tmp;
}
int main()
{

    scanf("%d%d",&n,&m);
    for(int i = 1,a,b,c;i < n; ++i) scanf("%d%d%d",&a,&b,&c),add(a,b,c),add(b,a,c);
    heavyedge :: main();
    f[0] = inf, sn = n, root = 0, GetRoot(1,0),tmp = root,dfs(root), root = tmp;
    int u,e;
    while(m--)
    {
        scanf("%d%d",&u,&e),Update(u,e),printf("%lld\n",Query(root));
    }
    return 0;
}