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

最新推荐文章于 2021-11-07 19:23:19 发布

CHN_JZ

最新推荐文章于 2021-11-07 19:23:19 发布

阅读量762

点赞数

本文链接：https://blog.csdn.net/CHN_JZ/article/details/78313026

版权

CHNJZ的OI学习总结同时被 3 个专栏收录

158 篇文章 1 订阅

订阅专栏

-------------各大网站-------------

137 篇文章 0 订阅

订阅专栏

-------------图论-------------

29 篇文章 0 订阅

订阅专栏

在艾泽拉斯大陆遭遇死亡之翼袭击的背景下，本文介绍了一种基于二分图匹配的算法来解决如何最高效地将各族人民引导至安全避难所的问题。通过使用Dinic算法优化匹配过程，确保在有限时间内所有种族都能通过狭窄的地精洞口进入地下避难所。

摘要由CSDN通过智能技术生成

题目描述

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

解题思路

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

我只想说出题人肯定~~吃屎了~~想让我们知道做二分图匹配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;
}