UVALive - 6887 Book Club (二分图求最大匹配||网络流求最大匹配)

传送门：https://icpcarchive.ecs.baylor.edu/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=4899

题目大意：

有n个人，m条有向边。对于一条边u->v表示u要从v得到一本书且仅一本（这一本数不一定是v原有的，也可以是v从别人那得到的），问在给出的边的条件下能不能使得每个人都得到一本书（非自己本人的）

题解：

目的就是让一个人可以能得到一本书，不管用什么的样的方式，只要符合条件就行,每个人都能有匹配就行了。所以只要求得最大匹配是否为n，判断一下。可以拆点建图，用二分图求或是用网络流求。

匈牙利算法（173ms）

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <queue>

using namespace std;
struct node
{
    int u,v;
    int nex;
}eage[40000];
int head[20000],cnt;
void ADD(int u,int v)
{
    eage[cnt].u=u;
    eage[cnt].v=v;
    eage[cnt].nex=head[u];
    head[u]=cnt++;
}
int used[20000],girl[20000];
bool dfs(int u)
{
    for(int i=head[u];i!=-1;i=eage[i].nex)
    {
        int v=eage[i].v;
        if(!used[v])
        {
            used[v]=1;
            if(girl[v]==-1||dfs(girl[v]))
            {
                girl[v]=u;
                return true;
            }
        }
    }
    return false;
}
int main()
{
    int n,m;
    while(~scanf("%d%d",&n,&m))
    {
        memset(head,-1,sizeof(head));
        cnt=0;
        int a,b;
        while(m--)
        {
            scanf("%d%d",&a,&b);
            ADD(a,b);
        }
        int all=0;
        memset(girl,-1,sizeof(girl));
        for(int i=0;i<n;i++)
        {
            memset(used,0,sizeof(used));
            if(dfs(i))all++;
        }
        if(all==n)printf("YES\n");
        else printf("NO\n");
    }
    return 0;
}

1. Hopcroft－Carp 算法(29ms)

#include <cstring>
#include <queue>
#include <cstdio>
#define INF 0x3f3f3f3f

using namespace std;
struct node
{
    int u,v,nex;
}eage[30000];
int head[20000],cnt;
void ADD(int u,int v)
{
    eage[cnt].u=u;
    eage[cnt].v=v;
    eage[cnt].nex=head[u];
    head[u]=cnt++;
}
int dis,Nx,Ny;
int Mx[20000],My[20000],dx[20000],dy[20000];
int vst[20000];
bool searchP()
{
    queue<int>q;
    dis=INF;
    memset(dx,-1,sizeof(dx));
    memset(dy,-1,sizeof(dy));
    for(int i=1;i<=Nx;i++)
    {
        if(Mx[i]==-1){q.push(i);
        dx[i]=0;}
    }
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
        if(dx[u]>dis)break;
        for(int i=head[u];i!=-1;i=eage[i].nex)
        {
            int v=eage[i].v;
            if(dy[v]==-1)
            {
                dy[v]=dx[u]+1;
                if(My[v]==-1)dis=dy[v];
                else
                {
                    dx[My[v]]=dy[v]+1;
                    q.push(My[v]);
                }
            }
        }
    }
    return dis!=INF;
}
bool DFS(int u)
{
    for(int i=head[u];i!=-1;i=eage[i].nex)
    {
        int v=eage[i].v;
        if(!vst[v]&&dy[v]==dx[u]+1)
        {
            vst[v]=1;
            if(My[v]!=-1&&dy[v]==dis)continue;
            if(My[v]==-1||DFS(My[v]))
            {
                My[v]=u;
                Mx[u]=v;
                return 1;
            }
        }
    }
    return 0;
}
int MaxMatch(int n)
{
    int res=0;
    memset(Mx,-1,sizeof(Mx));
    memset(My,-1,sizeof(Mx));
    Nx=Ny=n;
    while(searchP())
    {
        memset(vst,0,sizeof(vst));
        for(int i=1;i<=Nx;i++)
        if(Mx[i]==-1&&DFS(i))res++;
    }
    return res;
}
int main()
{
    int n,m;
    while(~scanf("%d%d",&n,&m))
    {
        int a,b;
        cnt=0;
        memset(head,-1,sizeof(head));
        while(m--)
        {
            scanf("%d%d",&a,&b);
            ADD(a+1,b+1);
        }
        int all=MaxMatch(n);
        if(all==n)printf("YES\n");
        else printf("NO\n");
    }
    return 0;
}

网络流（43ms）

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <queue>
#include <algorithm>
#include <cmath>
#define INF 0x3f3f3f3f
using namespace std;
typedef long long LL;
const int maxn=1e5+10;
int cas=1;
int d[maxn];///标记数组
int cur[maxn];
int start,tend;///源点和汇点
struct node///建边
{
    int to,cap,next;
}eage[maxn<<1];
int head[maxn];
bool vis[maxn];///标记
int cnt;
///建边函数
void ADD(int u,int v,int c)
{
    eage[cnt].to=v;
    eage[cnt].cap=c;
    eage[cnt].next=head[u];
    head[u]=cnt++;

    eage[cnt].to=u;
    eage[cnt].cap=0;
    eage[cnt].next=head[v];
    head[v]=cnt++;
}
bool BFS()
{
    memset(d,-1,sizeof(d));
    int Q[maxn*2];///队列
    int Thead=0,Ttail=0;
    Q[Ttail++]=start;
    d[start]=0;///标记
    while(Thead<Ttail)
    {
        //cout<<"sasadfasfdassgfafgdagfdd"<<endl;
        int x=Q[Thead];
        if(x==tend)return true;
        for(int i=head[x];i!=-1;i=eage[i].next)
        {
            int temp=eage[i].to;
            if(d[temp]==-1&&eage[i].cap>0)///没有标记，且可行流大于0
            {
                d[temp]=d[x]+1;
                Q[Ttail++]=temp;
            }
        }
        Thead++;
    }
    return false;///汇点是否成功标号，也就是说是否找到曾广路
}

int DFS(int x,int cap)
{
    //cout<<"saasdadasd"<<endl;
    if(x==tend)return cap;
    int flow = 0 , f;
    for(int i=head[x];i!=-1;i=eage[i].next)
    {
        int temp = eage[i].to;
        if(d[temp]==d[x]+1&&eage[i].cap)
        {
            f=DFS(temp,min(cap-flow,eage[i].cap));
            eage[i].cap -= f;
            eage[i^1].cap+=f;
            flow+=f;
            if(flow==cap)break;
        }
    }
    if(!flow)d[x]=-2;///防止重搜，剪枝。否则会超时
    return flow;
}
int maxflow()
{
    int flow=0,f;
    while(BFS())
    {
        //cout<<"sad"<<endl;
        while((f=DFS(start,INF))>0)
            flow+=f;
    }
    return flow;
}
int main()
{
    int n,m;
    while(~scanf("%d%d",&n,&m))
    {
        cnt=0;
        memset(head,-1,sizeof(head));
        int u,v;
        start=0;
        tend=2*n+1;
        for(int i=1;i<=n;i++)ADD(start,i,1);
        for(int i=n+1;i<=2*n;i++)ADD(i,tend,1);
        while(m--)
        {
            scanf("%d%d",&u,&v);///u<->v的流量是w(双向边)
            ADD(u+1,v+1+n,1);
        }
        int ans=maxflow();
        if(ans==n)printf("YES\n");
        else printf("NO\n");
    }
    return 0;
}