hdu5452 Minimum Cut（lca+树上差分）

题意：

给定n个点m条边的无权无向图G，
给定图中的一棵生成树T，
要求找到一个G的割，满足恰好割掉了T的一条边，
问满足条件的最小割是多少。

数据范围：n<=2e4，m<=2e5

解法：

考虑枚举T被割掉的边，因为T是一棵树，割掉一条边之后就变为两个连通块，
那么最小割显然就是非树边中连接两个连通块的边的数量+1，加1是因为还要算上这条树边。
那么问题变为对于每条树边，计算有多少条非树边连接两端的连通块。

反过来想，考虑每条非树边对树边的贡献，
如果存在这样一条非树边（图中红色部分）：

那么有贡献的树边就是（图中绿色的边）：

因为当选择删掉这些绿色的边时，红色的这条边连接两个连通块。

那么解决方法很显然了，对于每条非树边(x,v)，
将树中x到v的每条边都加上1，
最后答案就是树的最小边权+1

因为修改的是边权，差分的是点权，还用了边权转点权的技巧。

code：

#include <bits/stdc++.h>
using namespace std;
const int maxm=2e5+10;
const int maxd=20;
int head[maxm],nt[maxm<<1],to[maxm<<1],tot;
int f[maxm][25];
int num[maxm];
int d[maxm];
int n,m;
int ans;
void add(int x,int y){
    tot++;nt[tot]=head[x];head[x]=tot;to[tot]=y;
}
void dfs(int x,int fa){
    for(int j=1;j<=maxd;j++){
        f[x][j]=f[f[x][j-1]][j-1];
    }
    for(int i=head[x];i;i=nt[i]){
        int v=to[i];
        if(v==fa)continue;
        d[v]=d[x]+1;
        f[v][0]=x;
        dfs(v,x);
    }
}
void dfs2(int x,int fa){
    for(int i=head[x];i;i=nt[i]){
        int v=to[i];
        if(v==fa)continue;
        dfs2(v,x);
        num[x]+=num[v];
        ans=min(ans,num[v]);
    }
}
int lca(int a,int b){
    if(d[a]<d[b])swap(a,b);//把a定为更深的
    for(int i=maxd;i>=0;i--){//爬到同一高度
        if(d[f[a][i]]>=d[b]){
            a=f[a][i];
        }
    }
    if(a==b)return a;
    for(int i=maxd;i>=0;i--){
        if(f[a][i]!=f[b][i]){
            a=f[a][i];
            b=f[b][i];
        }
    }
    return f[a][0];
}
void init(){
    for(int i=1;i<=n;i++){
        head[i]=0;
        num[i]=0;
        d[i]=0;
    }
    for(int j=0;j<=maxd;j++){
        for(int i=1;i<=n;i++){
            f[i][j]=0;
        }
    }
    tot=0;
}
signed main(){
    int T;scanf("%d",&T);
    int cas=1;
    while(T--){
        init();
        scanf("%d%d",&n,&m);
        for(int i=1;i<n;i++){
            int a,b;scanf("%d%d",&a,&b);
            add(a,b);
            add(b,a);
        }
        d[1]=1;
        dfs(1,1);
        for(int i=n;i<=m;i++){
            int a,b;scanf("%d%d",&a,&b);
            int lc=lca(a,b);
            num[a]++;num[b]++;
            num[lc]-=2;
        }
        ans=1e9;
        dfs2(1,1);
        printf("Case #%d: ",cas++);
        printf("%d\n",ans+1);
    }
    return 0;
}

---