匈牙利算法

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(tǒng)[v]連的連其他的路
x[u]=p->v;y[p->v]=u;//找增廣路
return true;
}
}
}return false;
}
大匹配=最小點集
大獨立集=最小邊集(最小路徑覆蓋)=|V|-大匹配
#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;
}另外有需要云服務(wù)器可以了解下創(chuàng)新互聯(lián)scvps.cn,海內(nèi)外云服務(wù)器15元起步,三天無理由+7*72小時售后在線,公司持有idc許可證,提供“云服務(wù)器、裸金屬服務(wù)器、高防服務(wù)器、香港服務(wù)器、美國服務(wù)器、虛擬主機、免備案服務(wù)器”等云主機租用服務(wù)以及企業(yè)上云的綜合解決方案,具有“安全穩(wěn)定、簡單易用、服務(wù)可用性高、性價比高”等特點與優(yōu)勢,專為企業(yè)上云打造定制,能夠滿足用戶豐富、多元化的應(yīng)用場景需求。
網(wǎng)頁標題:二分圖匹配問題-創(chuàng)新互聯(lián)
標題鏈接:http://www.chinadenli.net/article28/dodijp.html
成都網(wǎng)站建設(shè)公司_創(chuàng)新互聯(lián),為您提供網(wǎng)站營銷、網(wǎng)站設(shè)計、網(wǎng)站建設(shè)、網(wǎng)站制作、品牌網(wǎng)站制作、網(wǎng)站排名
聲明:本網(wǎng)站發(fā)布的內(nèi)容(圖片、視頻和文字)以用戶投稿、用戶轉(zhuǎn)載內(nèi)容為主,如果涉及侵權(quán)請盡快告知,我們將會在第一時間刪除。文章觀點不代表本網(wǎng)站立場,如需處理請聯(lián)系客服。電話:028-86922220;郵箱:631063699@qq.com。內(nèi)容未經(jīng)允許不得轉(zhuǎn)載,或轉(zhuǎn)載時需注明來源: 創(chuàng)新互聯(lián)
猜你還喜歡下面的內(nèi)容