AC日记——【模板】二分图匹配 洛谷 P3386

题目背景

二分图

题目描述

给定一个二分图，结点个数分别为n,m，边数为e，求二分图最大匹配数

输入输出格式

输入格式：

第一行，n,m,e

第二至e+1行，每行两个正整数u,v，表示u,v有一条连边

输出格式：

共一行，二分图最大匹配

输入输出样例

输入样例#1：

1 1 1
1 1

输出样例#1：

1

说明

n,m<=1000，1<=u<=n，1<=v<=m

因为数据有坑，可能会遇到v>m的情况。请把v>m的数据自觉过滤掉。

算法：二分图匹配

思路：

　　二分图模板；

来，上代码：

#include <queue>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>

#define maxn 1005
#define INF 0x7fffffff

using namespace std;

struct EdgeType {
    int v,f,e;
};
struct EdgeType edge[maxn*maxn*];

int cnt,deep[maxn<<],ans,e;
int n,m,head[maxn<<],s=,t=(maxn<<)-;

char Cget;

inline void in(int &now)
{
    now=,Cget=getchar();
    while(Cget>''||Cget<'') Cget=getchar();
    while(Cget>=''&&Cget<='')
    {
        now=now*+Cget-'';
        Cget=getchar();
    }
}

bool bfs()
{
    for(int i=s;i<=t;i++) deep[i]=-;
    queue<int>que;deep[s]=,que.push(s);
    while(!que.empty())
    {
        int now=que.front();que.pop();
        for(int i=head[now];i;i=edge[i].e)
        {
            if(edge[i].f>&&deep[edge[i].v]<)
            {
                deep[edge[i].v]=deep[now]+;
                if(edge[i].v==t) return true;
                que.push(edge[i].v);
            }
        }
    }
    return false;
}

int flowing(int now,int flow)
{
    if(now==t||flow<=) return flow;
    int oldflow=;
    for(int i=head[now];i;i=edge[i].e)
    {
        if(edge[i].f<=||deep[edge[i].v]!=deep[now]+) continue;
        int pos=flowing(edge[i].v,min(edge[i].f,flow));
        if(pos>)
        {
            flow-=pos;
            oldflow+=pos;
            edge[i].f-=pos;
            edge[i^].f+=pos;
            if(flow==) return oldflow;
        }
    }
    if(oldflow==) deep[now]=-;
    return oldflow;
}

int main()
{
    in(n),in(m),in(e);
    for(int i=;i<=n;i++)
    {
        edge[++cnt].v=i,edge[cnt].f=,edge[cnt].e=head[s],head[s]=cnt;
        edge[++cnt].v=s,edge[cnt].f=,edge[cnt].e=head[i],head[i]=cnt;
    }
    for(int i=+n;i<=m+n;i++)
    {
        edge[++cnt].v=t,edge[cnt].f=,edge[cnt].e=head[i],head[i]=cnt;
        edge[++cnt].v=i,edge[cnt].f=,edge[cnt].e=head[t],head[t]=cnt;
    }
    int u,v;
    while(e--)
    {
        in(u),in(v);v+=n;
        edge[++cnt].v=v,edge[cnt].f=,edge[cnt].e=head[u],head[u]=cnt;
        edge[++cnt].v=u,edge[cnt].f=,edge[cnt].e=head[v],head[v]=cnt;
    }
    while(bfs()) ans+=flowing(s,INF);
    cout<<ans;
    return ;
}