P3128 [USACO15DEC] Max Flow P

//https://www.luogu.com.cn/problem/P3128 
#include <bits/stdc++.h>
using namespace std;
const int N=5e4+10,M=N*2;
int e[M],ne[M],h[N],idx;
int depth[N],fa[N][20];
int q[N];
int n,m;
int ans=0;
int d[N];
inline int read(){
	int x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-') f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
	return x*f;
}
void write(int x){
	if(x<0) putchar('-'),x=-x;
	if(x>9) write(x/10);
	putchar(x%10+'0');
}
void add(int a,int b){
    e[idx]=b,ne[idx]=h[a],h[a]=idx++;
}
void bfs(int root){
    memset(depth,0x3f,sizeof depth);
    depth[0]=0,depth[root]=1;
    int hh=0,tt=1;
    q[0]=root;
    while(hh<tt){
        int u=q[hh++];
        for(int i=h[u];~i;i=ne[i]){
            int v=e[i];
            if(depth[v]>depth[u]+1){
                depth[v]=depth[u]+1;
                q[tt++]=v;
                fa[v][0]=u;
                for(int k=1;k<=16;k++) fa[v][k]=fa[fa[v][k-1]][k-1];
            }
        }
    }
}
int dfs(int u,int father){
    int s=d[u];
    for(int i=h[u];~i;i=ne[i]){
        int j=e[i];
        if(j!=father){
            s+=dfs(j,u);
        }
    }
	if(s>ans) ans=s;
    return s;
}
int lca(int a,int b){
    if(depth[a]<depth[b]) swap(a,b);
    for(int k=16;k>=0;k--) if(depth[fa[a][k]]>=depth[b]) a=fa[a][k];
    if(a==b) return a;
    for(int k=16;k>=0;k--){
        if(fa[a][k]!=fa[b][k]){
            a=fa[a][k],b=fa[b][k];
        }
    }
    return fa[a][0];
}
int main(){
    n=read(),m=read();
    memset(h,-1,sizeof h);
    for(int i=0,a,b;i<n-1;i++){
        a=read(),b=read();
        add(a,b),add(b,a);
    }
    bfs(1);
    for(int i=0,a,b;i<m;i++){
        a=read(),b=read();
        d[a]++,d[b]++;
        int p=lca(a,b);
        d[p]--;
        if(fa[p][0]) d[fa[p][0]]--;
    }
    dfs(1,-1);
    write(ans);
    return 0;
}
posted @ 2025-05-10 10:17  九三青梧  阅读(8)  评论(0)    收藏  举报