【网络流#6】POJ 3041 Asteroids 二分图最大匹配 - 《挑战程序设计竞赛》例题

学习网络流中ing...作为初学者练习是不可少的~~~构图方法因为书上很详细了，所以就简单说一说

---

 

把光束作为图的顶点，小行星当做连接顶点的边，建图，由于 最小顶点覆盖 等于 二分图最大匹配 ，因此求二分图最大匹配即可。

邻接矩阵，DFS寻找增广路，匈牙利算法

邻接矩阵：复杂度O(n^3)

如果使用邻接表：复杂度O(n*m)

#include<cstdio>

#include<cstring>

#include<cmath>

#include<iostream>

#include<algorithm>

#include<set>

#include<map>

#include<stack>

#include<vector>

#include<queue>

#include<string>

#include<sstream>

#define eps 1e-9

#define ALL(x) x.begin(),x.end()

#define INS(x) inserter(x,x.begin())

#define FOR(i,j,k) for(int i=j;i<=k;i++)

using namespace std;

typedef long long LL;

int i,j,k,n,m,x,y,T,ans,big,cas,match[505],g[505][505],r,c;

bool flag,vis[505];



bool hungary(int u)//从u开始寻找增广路 

{

    for (int i=1;i<=n;i++)

    {

        if (g[u][i]&&!vis[i])//G数组存储邻接矩阵，match[i]表示当前结点i匹配的结点

        {

            vis[i]=1;

            if (match[i]==0||hungary(match[i]))

                return match[i]=u;

        }

    }

    return 0;

}



int main()

{

    scanf("%d%d",&n,&k);

    for (i=1;i<=k;i++)

    {

        scanf("%d%d",&r,&c);

        g[r][c]=1;

    }

    ans=0;



        /*匈牙利算法过程 - 寻找增广路*/

    for (i=1;i<=n;i++)

    {

        memset(vis,0,sizeof(vis));

        if (hungary(i)) ans++;

    }

    /*寻找结束*/



    printf("%d\n",ans);

    return 0;

}