Minimum Cut HDU - 5452

http://acm.hdu.edu.cn/showproblem.php?pid=5452

考虑枚举每一条树边 拿掉这条边后会把树分成两个子图 还需要删除的就是两个端点在不同子图上的非树边

所以就计算每条非树边对树边的影响 这些树边也就是非树边两端点在树上确定的那条路径上的边 考虑树剖链更新 但是会T

可以发现更新和查询是分开的 所以利用差分数组的思想 对两点++ 对两点的lca-=2 然后后序遍历一遍 每个点把所有子节点的值都拿上来 最后扫一遍数组就可以了

#include <bits/stdc++.h>
using namespace std;
#define N 0x3f3f3f3f

struct node1
{
    int u;
    int v;
};

struct node2
{
    int v;
    int next;
};

node1 pre[200010];
node2 edge[40010];
int dp[20010][20];
int first[20010],fa[20010],deep[20010],val[20010];
int n,m,num,sum;

void addedge(int u,int v)
{
    edge[num].v=v;
    edge[num].next=first[u];
    first[u]=num++;
}

void dfsI(int cur)
{
    int i,v,w;
    for(i=first[cur];i!=-1;i=edge[i].next)
    {
        v=edge[i].v;
        if(v!=fa[cur])
        {
            dp[v][0]=cur;
            fa[v]=cur,deep[v]=deep[cur]+1;
            dfsI(v);
        }
    }
    return;
}

void solve()
{
    int i,j;
    dp[1][0]=0;
    fa[1]=0,deep[1]=1;
    dfsI(1);
    for(j=1;(1<<j)<=n;j++)
    {
        for(i=1;i<=n;i++)
        {
            dp[i][j]=dp[dp[i][j-1]][j-1];
        }
    }
    return;
}

int getlca(int u,int v)
{
    int i;
    if(deep[u]<deep[v]) swap(u,v);
    for(i=log2(n);i>=0;i--)
    {
        if(deep[dp[u][i]]>=deep[v])
        {
            u=dp[u][i];
        }
    }
    if(u==v) return u;
    for(i=log2(n);i>=0;i--)
    {
        if(dp[u][i]!=dp[v][i])
        {
            u=dp[u][i];
            v=dp[v][i];
        }
    }
    return dp[u][0];
}

void dfsII(int cur)
{
    int i,v;
    for(i=first[cur];i!=-1;i=edge[i].next)
    {
        v=edge[i].v;
        if(v!=fa[cur])
        {
            dfsII(v);
            val[cur]+=val[v];
        }
    }
}

int main()
{
    int t,cas,i,ans,lca;
    scanf("%d",&t);
    for(cas=1;cas<=t;cas++)
    {
        scanf("%d%d",&n,&m);
        memset(first,-1,sizeof(first));
        num=0;
        for(i=1;i<=n-1;i++)
        {
            scanf("%d%d",&pre[i].u,&pre[i].v);
            addedge(pre[i].u,pre[i].v);
            addedge(pre[i].v,pre[i].u);
        }
        for(;i<=m;i++) scanf("%d%d",&pre[i].u,&pre[i].v);
        solve();

        for(i=1;i<=m;i++)
        {
            if(deep[pre[i].u]<deep[pre[i].v]) swap(pre[i].u,pre[i].v);
        }

        memset(val,0,sizeof(val));
        for(i=n;i<=m;i++)
        {
            lca=getlca(pre[i].u,pre[i].v);
            val[lca]-=2,val[pre[i].u]++,val[pre[i].v]++;
        }

        sum=0;
        dfsII(1);
        ans=N;
        for(i=2;i<=n;i++) ans=min(ans,val[i]);

        printf("Case #%d: %d\n",cas,ans+1);
    }
    return 0;
}

/*
1
8 10
1 2
2 3
3 4
4 6
6 7
7 8
8 5

1 4
6 5
5 1
*/