匈牙利演算法求最大匹配

#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;}