[二分+最大流]51 Nod——1757 大灾变

题目描述

死亡之翼降临了！艾泽拉斯大陆的子民们必须逃出他的魔爪！
艾泽拉斯的结构是一棵树，这棵树上的一些节点是地精建造的通往地下避难所的洞口。
除了这些洞口之外，树上的每个节点上都有一个种族，每个种族通过树上的一条边都需要一个单位时间。
因为地精比较矮小，所以洞口很窄，每个单位时间只能让一个种族通过，但是一个单位时间内的一个节点上可以存在多个种族。
地精们需要你求出最少需要多少单位时间才能让所有种族躲进地下避难所。
【注意题目有修改，洞口不一定是叶子节点】

解题思路

这题目解法比较明显，我不想多说。

我只想说出题人肯定吃屎了想让我们知道做二分图匹配Dinic比匈牙利更优。

#include<cstdio>
#include<cstring>
using namespace std;
const int maxn=2005,maxm=400005,INF=2147483647,maxt=100005;
int n,L,R,m,mid,cst[maxn][maxn],a[maxn],h[maxn][maxn];
int tot,son[maxm],nxt[maxm],time,lnk[maxt];
int que[maxt],cap[maxm],hed,til,flow[maxm],dis[maxt],B[maxt],s,t;
bool b[maxn],vis[maxt];
void DFS(int g,int x,int fa,int dis){
    cst[g][x]=dis;
    for (int i=1;i<=h[x][0];i++) if (h[x][i]!=fa) DFS(g,h[x][i],x,dis+1);
}
void add(int x,int y,int c){
    nxt[++tot]=lnk[x];lnk[x]=tot;son[tot]=y;cap[tot]=c;flow[tot]=0;
    nxt[++tot]=lnk[y];lnk[y]=tot;son[tot]=x;cap[tot]=0;flow[tot]=0;
}
int min(int x,int y){if (x<y) return x;return y;}
bool BFS(){
    memset(vis,0,sizeof(vis));
    hed=0,til=1;que[1]=s;vis[s]=1;
    while(hed!=til){
        int x=que[++hed];
        for (int j=lnk[x];j!=-1;j=nxt[j])
        if (!vis[son[j]]&&cap[j]>flow[j]){
            que[++til]=son[j];vis[son[j]]=1;
            dis[son[j]]=dis[x]+1;
        }
    }
    return vis[t];
}
int DFS(int x,int mi){
    if (x==t||mi==0) return mi;
    int num=0;
    for (int &j=B[x];j!=-1;j=nxt[j])
    if (dis[x]+1==dis[son[j]]){
        int now=DFS(son[j],min(mi,cap[j]-flow[j]));
        if (now){
            flow[j]+=now;flow[j^1]-=now;
            num+=now;mi-=now;
            if (!mi) break;
        }
    }
    return num;
}
int Dinic(){
    int num=0;
    while(BFS()){
        memcpy(B,lnk,sizeof(B));
        num+=DFS(0,INF);
    }
    return num;
}
bool check(int x){
    s=0;t=n+m*x+1;
    memset(lnk,-1,sizeof(lnk));tot=-1;
    for (int i=1;i<=n;i++) if (!b[i]) add(s,i,1);
    for (int i=1;i<=m;i++){
        for (int j=1;j<=x;j++) add(n+(i-1)*x+j,t,1);
        for (int j=1;j<x;j++) add(n+(i-1)*x+j,n+(i-1)*x+j+1,INF);
    }
    for (int i=1;i<=n;i++) if (!b[i])
    for (int j=1;j<=m;j++) if (cst[i][a[j]]<=x) add(i,n+(j-1)*x+cst[i][a[j]],1);
    //printf("%d %d\n",Dinic(),x);
    return Dinic()==n-m;
}
int main(){
    freopen("exam.in","r",stdin);
    freopen("exam.out","w",stdout);
    scanf("%d%d",&n,&m);
    for (int i=1;i<n;i++){
        int x,y;scanf("%d%d",&x,&y);
        h[x][++h[x][0]]=y;
        h[y][++h[y][0]]=x;
    }
    for (int i=1;i<=m;i++) scanf("%d",&a[i]),b[a[i]]=1;
    for (int i=1;i<=n;i++) DFS(i,i,0,0);
    L=1,R=n;
    while(L<=R){
        mid=L+(R-L>>1);
        if (check(mid)) R=mid-1;else L=mid+1;
    }
    printf("%d\n",L);
    return 0;
}