匈牙利算法求最大匹配

#define MAXN 1000         // 实际问题时需要修改
int  mat[MAXN][MAXN];    // 邻接矩阵mat 的0行0列不用 
int  nx, ny;                         // 实际问题时矩阵的行列数 
int  fy[MAXN], matx[MAXN], maty[MAXN]; 
int  path( int  u ) 
{     int v; 
      for (v=1;v<=ny;v++) 
           if (mat[u][v] && fy[v]<0) 
           {         fy[v]=1; 
                     if (maty[v]<0 || path(maty[v])) 
                     {   matx[u]=v;    maty[v]=u;   return 1;   }
           } 
      return 0; 
} 


int MaximumMatch( ) 
{	int  i, ret=0;
	memset(matx,-1,sizeof(matx));
	memset(maty,-1,sizeof(maty));
	for (i=1;i<=nx;i++) if (matx[i]<0)
	 {	memset(fy,-1,sizeof(fy));
   		ret+=path(i);
	}
	return ret;
}