P3128 Max Flow P(LCA,树上差分)

[USACO15DEC]Max Flow P - 洛谷

给定一棵树，有k次操作,每次u 到 v 路径上的点权值加一。求权值最大的点。

思路：

将区间操作变为点操作

u的点权加一，表示u到root的所有点权值加一，那么u到v的路径点权加一可转化为：

w[u]++;

w[v]++;

w[lca(u,v)]--;

w[fa[lca(u,v)]--; //fa[lca(u,v)为lca(u,v)的父节点

该点的最终权值即为子树权值和其本身之和，dfs求解即可。

Code:

//LCA 树上差分
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cstdio>
#include<cmath>
#include<vector>
#include<queue>
#include<stack>
#include<map>
#include<unordered_map>
using namespace std;
#define x first
#define y second
#define PII pair<int,int>
#define V vector<int>
#define endl "\n"
#define init(array) memset(array,0,sizeof array)
typedef long long ll;
typedef unsigned long long llu;
const int INF=0x3f3f3f3f;
const int N=5e4+10,M=2e5+10;
int n,m;
int h[N],e[M],ne[M],idx;
int depth[N];
int fa[N][16];
int w[N];
void add(int a,int b)
{
    e[idx]=b,ne[idx]=h[a],h[a]=idx++;
}
void bfs()
{
    memset(depth,0x3f,sizeof depth);
    queue<int> q;
    depth[1]=1,depth[0]=0;
    q.push(1);
    while(q.size())
    {
        int t=q.front();
        q.pop();
        for(int i=h[t];~i;i=ne[i])
        {
            int j=e[i];
            if(depth[j]>depth[t]+1)
            {
                fa[j][0]=t;
                for(int k=1;k<=15;k++)
                    fa[j][k]=fa[fa[j][k-1]][k-1];
                depth[j]=depth[t]+1;
                q.push(j);
            }
        }
    }
}
int lca(int a,int b)
{
    if(depth[a]<depth[b]) swap(a,b);
    for(int k=15;k>=0;k--)
    {
        if(depth[fa[a][k]]>=depth[b])
            a=fa[a][k];  
    }
    if(a==b) return a;
    for(int k=15;k>=0;k--)
    {
        if(fa[a][k]!=fa[b][k])
        {
            a=fa[a][k];
            b=fa[b][k];
        }
    }
    return fa[a][0];
}
void dfs(int u,int fa)
{
    for(int i=h[u];~i;i=ne[i])
    {
        int j=e[i];
        if(j==fa) continue;
        dfs(j,u);
        w[u]+=w[j];
    }
}
int main()
{
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    cout.tie(nullptr);
    memset(h,-1,sizeof h);
    cin>>n>>m;
    for(int i=1;i<n;i++)
    {
        int a,b;
        cin>>a>>b;
        add(a,b);
        add(b,a);
    }
    bfs();
    while(m--)
    {
        int a,b;
        cin>>a>>b;
        int tmp=lca(a,b);
        w[a]++;
        w[b]++;
        w[tmp]--;
        w[fa[tmp][0]]--;
    }
    dfs(1,0);
    int res=0;
    for(int i=1;i<=n;i++) res=max(res,w[i]);
    cout<<res;
    return 0;
}