• 匈牙利算法

  

  創新互聯公司專注于企業全網整合營銷推廣、網站重做改版、納雍網站定制設計、自適應品牌網站建設、H5技術、商城建設、集團公司官網建設、外貿網站制作、高端網站制作、響應式網頁設計等建站業務，價格優惠性價比高，為納雍等各大城市提供網站開發制作服務。

  struct edge{int u,v;edge *next;}*head[N],e[N];
  void add(int u,int v){
  edge *p=&e[cnt++];
  p->u=u;p->v=v;p->next=head[u];head[u]=p;
  }
  bool dfs(int u){
  for(edge *p=head[u];p;p=p->next){
      if(!vis[p->v]){
          vis[p->v]=1;
          if(!y[p->v]||dfs(y[p->v])){// 如果u->v 可以連標記 或者回溯到y[v]連的連其他的路
              x[u]=p->v;y[p->v]=u;//找增廣路
              return true;
          }
      }
  }return false;
  }

  大匹配=最小點集
  大獨立集=最小邊集(最小路徑覆蓋)=|V|-大匹配

• KM算法

#include<bits/stdc++.h>
#define me(a,x) memset(a,x,sizeof(a))
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const int mod=1e9+7;
const int N=2e3+5;
const int MAX=0x7fffffff;
const int MIN=0x80000000;
int nx,ny;
int w[N][N],lx[N],ly[N];
bool sx[N],sy[N];
int match[N][N];
bool dfs(int u){
    sx[u]=true;
    for(int v=0;v<ny;v++){
        int d=lx[u]+ly[v]-w[u][v];
        if(!sy[v]&&!d){
            sy[v]=true;
            if(match[v]==-1||dfs(match[v])){
                match[v]=u;return true;
            }
        }
    }return false;
}
int KM(){
    for(int i=0;i<nx;i++){
        ly[i]=lx[i]=0;
        for(int j=0;j<ny;j++)lx[i]=max(lx[i],w[i][j]);
    }
    me(match,-1);
    for(int u=0;u<nx;u++){
       while(true){
            me(sx,0),me(sy,0);
            if(dfs(u))break;
            int dis=MAX;
            for(int i=0;i<nx;i++){
                if(sx[i]){
                    for(int j=0;j<ny;j++){
                        if(!sy[j]){
                            dis=min(dis,lx[i]+ly[j]-w[i][j]);
                        }
                    }
                }
            }
            if(dis==0)continue;
            if(dis==MAX)return -1;
            for(int i=0;i<nx;i++){
                if(sx[i])lx[i]-=dis;
            }
            for(int i=0;i<ny;i++){
                if(sy[i])ly[i]+=dis;
            }
        }
    }
    int sum=0;
    for(int i=0;i<ny;i++){
        if(match[i]>=0)sum+=w[match[i]][i];
    }return sum;
}

#include<bits/stdc++.h>
#define me(a,x) memset(a,x,sizeof(a))
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const int mod=1e9+7;
const int N=2e3+5;
const int MAX=0x7fffffff;
const int MIN=0x80000000;
int nx,ny;
int w[N][N],lx[N],ly[N];
bool sx[N],sy[N];
int match[N][N];
int slack[N];
bool dfs(int u){
    sx[u]=true;
    for(int v=0;v<ny;v++){
        int d=lx[u]+ly[v]-w[u][v];
        if(!sy[v]&&!d){
            sy[v]=true;
            if(match[v]==-1||dfs(match[v])){
                match[v]=u;return true;
            }
        }else if(!sy[v]&&d)slack[i]=min(slack[i],d);
    }return false;
}
int KM(){
    for(int i=0;i<nx;i++){
        ly[i]=lx[i]=0;
        for(int j=0;j<ny;j++)lx[i]=max(lx[i],w[i][j]);
    }
    me(match,-1);
    for(int u=0;u<nx;u++){
        for(int i=0;i<N;i++)slack[i]=MAX;
       while(true){
            me(sx,0),me(sy,0);
            if(dfs(u))break;
            int dis=MAX;
            for(int i=0;i<ny;i++){
                if(!sy[i])dis=min(dis,slack[i]);
            }
            if(dis==0)continue;
            if(dis==MAX)return -1;
            for(int i=0;i<nx;i++){
                if(sx[i])lx[i]-=dis;
            }
            for(int i=0;i<ny;i++){
                if(sy[i])ly[i]+=dis;
                else slack[i]-=dis;
            }
        }
    }
    int sum=0;
    for(int i=0;i<ny;i++){
        if(match[i]>=0)sum+=w[match[i]][i];
    }return sum;
}

另外有需要云服務器可以了解下創新互聯scvps.cn，海內外云服務器15元起步，三天無理由+7*72小時售后在線，公司持有idc許可證，提供“云服務器、裸金屬服務器、高防服務器、香港服務器、美國服務器、虛擬主機、免備案服務器”等云主機租用服務以及企業上云的綜合解決方案，具有“安全穩定、簡單易用、服務可用性高、性價比高”等特點與優勢，專為企業上云打造定制，能夠滿足用戶豐富、多元化的應用場景需求。

網頁標題：二分圖匹配問題-創新互聯
標題鏈接：http://www.yijiale78.com/article28/dodijp.html

成都網站建設公司_創新互聯，為您提供網站營銷、網站設計、網站建設、網站制作、品牌網站制作、網站排名

廣告

聲明：本網站發布的內容（圖片、視頻和文字）以用戶投稿、用戶轉載內容為主，如果涉及侵權請盡快告知，我們將會在第一時間刪除。文章觀點不代表本網站立場，如需處理請聯系客服。電話：028-86922220；郵箱：631063699@qq.com。內容未經允許不得轉載，或轉載時需注明來源： 創新互聯