poj_2594_匈牙利

题目描述：

有向图，无环路。求由多少个起点开始走路，能够走遍所有点。（所走路径中间的点可重复）。

 

解题思路:

联想到air raid那题的最小路径覆盖：如何找到最小的边，覆盖所能够覆盖到的点——这里直接利用匈牙利算法，找到的是不重复走点所能找到的最大匹配路径。
（因为匈牙利中一个点只能被一个前驱结点所link）。所以需要转化下图：把所有能到达的路都转化为边——转化方法类似迪杰斯特拉。这样就能够使用最小路径覆盖了。

 

代码：

#include <stdio.h>
#include <stdlib.h>
#define N 501

int c[N][N], D[N][N], visit[N], d[N], n, m, link[N];

int dfs(int v){
   int i;
   for(i=1;i<=n;i++){
      if(D[v][i] && !visit[i]){
         visit[i]=1;
         if(link[i]==-1 || dfs(link[i])){
            link[i] = v;
            return 1;
         }
      }
   }
   return 0;
}

main(){
   int i,j,k,x,y,cnt,f;
  
   scanf("%d %d",&n, &m);
   while(n!=0){
       memset(c,0,sizeof(c));
       for(i=1;i<=m;i++){
          scanf("%d %d",&x, &y);
          c[x][y] = 1;
       }
       // update edge
       memset(D,0,sizeof(D));
       for(i=1;i<=n;i++){ //更新每个点所能到达的其余点
          memset(visit,0,sizeof(visit));
          visit[i] = 1;
          for(j=1;j<=n;j++)
             d[j] = c[i][j];
         
          for(cnt=1;cnt<=n-1;cnt++){//每个点最多找n-1次
              f = 0;
              for(j=1;j<=n;j++){
                 if(!visit[j] && d[j]){ //找到，加入
                    D[i][j] = 1;
                    visit[j] = 1;
                    f = 1;
                    for(k=1;k<=n;k++){ //更新
                       if(!visit[k] && c[j][k]){
                          d[k] = 1;
                       }
                    }
                 }
              }
              if(f==0) //没有找到
                 break;
          }
         
          for(j=1;j<=n;j++) //i所能到达的点
             if(d[j])
                D[i][j] = 1;  
       }
       //debug
//       for(i=1;i<=n;i++){
//          for(j=1;j<=n;j++)
//             printf("%d ",D[i][j]);
//          printf("\n");
//       }
      
       for(i=1;i<=n;i++)
          link[i] = -1;
         
       cnt = 0;
       for(i=1;i<=n;i++){
          memset(visit,0,sizeof(visit));
          if(dfs(i)){
             cnt++;     
          }
       }
       printf("%d\n",n-cnt);
       scanf("%d %d",&n,&m);
   }
  
   system("pause");
   return 0;
}