bzoj3924 幻想乡战略游戏【树链剖分+线段树】

解题思路：

太蒟蒻了，不会动态点分治，只能链剖暴力移中心了，不过还好bzoj数据水。
很明显是求待修改下的带权重心，于是链剖建线段树维护每个点的权值，注意到一个子树在该线段树上编号是连续的（以前竟然没发现，惭愧……），所以以u为根的子树的权值和就是sum（pos[u]~pos[u]+size[u]-1）。
每次修改后暴力往更优的方向走就行了。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<cmath>
#include<algorithm>
#include<ctime>
#include<vector>
#include<queue>
#define ll long long
using namespace std;

int getint()
{
    int i=0,f=1;char c;
    for(c=getchar();(c!='-')&&(c<'0'||c>'9');c=getchar());
    if(c=='-')c=getchar(),f=-1;
    for(;c>='0'&&c<='9';c=getchar())i=(i<<3)+(i<<1)+c-'0';
    return i*f;
}

const int N=100005;
int n,m,totsize,rt=1;
int ecnt,first[N],nxt[N<<1],to[N<<1],w[N<<1];
int tot,fa[N],son[N],dep[N],top[N],pos[N],size[N],tr[N<<2];
ll ans;

void add(int x,int y,int z)
{
    nxt[++ecnt]=first[x],first[x]=ecnt,to[ecnt]=y,w[ecnt]=z;
}

void dfs1(int u)
{
    size[u]=1;
    for(int e=first[u];e;e=nxt[e])
    {
        int v=to[e];
        if(v==fa[u])continue;
        fa[v]=u,dep[v]=dep[u]+w[e];
        dfs1(v);
        size[u]+=size[v];
        if(size[v]>size[son[u]])son[u]=v;
    }
}

void dfs2(int u,int f)
{
    pos[u]=++tot,top[u]=f;
    if(son[u])dfs2(son[u],f);
    for(int e=first[u];e;e=nxt[e])
    {
        int v=to[e];
        if(v!=fa[u]&&v!=son[u])
            dfs2(v,v);
    }
}

int lca(int x,int y)
{
    while(top[x]!=top[y])
        dep[top[x]]>dep[top[y]]?x=fa[top[x]]:y=fa[top[y]];
    return dep[x]<dep[y]?x:y;
}

int dis(int x,int y)
{
    return dep[x]+dep[y]-2*dep[lca(x,y)];
}

void modify(int k,int l,int r,int p,int v)
{
    if(l==r)
    {
        tr[k]+=v;
        return;
    }
    int mid=l+r>>1;
    if(p<=mid)modify(k<<1,l,mid,p,v);
    else modify(k<<1|1,mid+1,r,p,v);
    tr[k]=tr[k<<1]+tr[k<<1|1];
}

int query(int k,int l,int r,int x,int y)
{
    if(x<=l&&r<=y)return tr[k];
    int mid=l+r>>1;
    if(y<=mid)return query(k<<1,l,mid,x,y);
    else if(x>mid)return query(k<<1|1,mid+1,r,x,y);
    else return query(k<<1,l,mid,x,mid)+query(k<<1|1,mid+1,r,mid+1,y);
}

ll calc(int u,int v)
{
    int t1,t2;
    if(u==fa[rt])
    {
        t1=query(1,1,n,pos[rt],pos[rt]+size[rt]-1);
        t2=totsize-t1;
    }
    else
    {
        t2=query(1,1,n,pos[u],pos[u]+size[u]-1);
        t1=totsize-t2;
    }
    return ans+1ll*(t1-t2)*v;
}

void move(int u)
{
    for(int e=first[u];e;e=nxt[e])
    {
        int v=to[e];
        ll t=calc(v,w[e]);
        if(t<ans)
        {
            ans=t,rt=v;
            move(rt);break;
        }
    }
}

int main()
{
    //freopen("lx.in","r",stdin);
    int x,y,z;
    n=getint(),m=getint();
    for(int i=1;i<n;i++)
    {
        x=getint(),y=getint(),z=getint();
        add(x,y,z),add(y,x,z);
    }
    dfs1(1);
    dfs2(1,1);
    while(m--)
    {
        x=getint(),y=getint();
        ans+=1ll*dis(x,rt)*y;
        totsize+=y;
        modify(1,1,n,pos[x],y);
        move(rt);
        cout<<ans<<'\n';
    }
}