【POJ2594】Treasure Exploration 宝藏勘探

Description
　　你曾经看过的任何关于宝藏勘探的书吗？你有没有看到过任何关于宝藏勘探的电影？你曾经探讨过宝藏吗？如果你从来没有这种经历，你将永远不知道会宝藏探索带给您是多么的有趣。最近，一家名为EUC（探索未知的公司）计划，以探索一个不为人知的地方在火星上，这里被认为有丰富的宝藏。对于快速发展的技术和火星上恶劣的环境，EUC送了一些机器人去探索宝藏。简单的说，我们现有一个由N点(编号从1到n)来表示勘探的地方，用一些单向道路连接着的火星地图。这意味着通过这些道路，机器人只能从一端走向另一端，但不能返回。在这个地图上没有圈。该机器人可以通过火箭传送到火星上任何一个地方。着陆后，机器人可以通过连接着的道路探索一些点，您应该注意到的是两个不同的机器人探索的道路，可能包含一些相同点。由于财政的原因，EUC想用最少的机器人去探索火星上的各个点。作为一个ICPCer，具有良好的编程技巧，请你帮助EUC？
Input
　　输入文件包含多组测试数据，对于每组测试数据第一行包含两个整数n(1<=n<=500)和m(0<=m<=5000)，分别表示火星地图上的点数和边数，接下来的m行，每行包含两个用空格分开的整数A和B，表示一条从A到B的路径(0<=A,B<=n)，输入文件以一行两个0 0作为结束标记。
Output
　　输出文件对于每一组测试数据输出一行表示需要的最少机器人的数目。
Sample Input
1 0
2 1
1 2
2 0
0 0
Sample Output
1
1
2

---

求最小路径覆盖
网络流不会
二分图吧匈牙利算法
注意 一个点可以走多次
用传递闭包把间接相连的点改为直接相连

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#define N 505
using namespace std;

int g[N][N],vst[N],my[N],n,m;
int dfs(int i)
{
    for(int j=1;j<=n;j++)
    {
        if(g[i][j]&&!vst[j])
        {
        vst[j]=1;
        if(!my[j]||dfs(my[j]))
        {
            my[j]=i;
            return 1;
        }
        }
    }
    return 0;
}
void Hungary()
{
    int ans=0;
    memset(my,0,sizeof(my));
    for(int i=1;i<=n;i++)
    {
        memset(vst,0,sizeof(vst));
        if(dfs(i)) ans++;
    }
    printf("%d\n",n-ans);
}
int main()
{
    int i,j,k,x,y;
   cin>>n>>m;
   while(n!=0||m!=0)
   {
    memset(g,0,sizeof(g));
      for(i=1;i<=m;i++) cin>>x>>y,g[x][y]=1;
      for(i=1;i<=n;i++)
      for(j=1;j<=n;j++)
      if(g[j][i])
      for(k=1;k<=n;k++)
      if(g[i][k]) g[j][k]=1;
      Hungary();
      cin>>n>>m;
   }
   return 0;
}