[gotoac]二分图最大匹配hungary & 二分图最佳完美匹配KM(邻接表&邻接矩阵) beta1

struct hungary{
int n,m;            //X集和Y集的点数
int match[M];       //匹配
bool vis[M];        //是否在增广路上
vector<int>g[M];    //邻接表
void init(int _n,int _m){
n=_n,m=_m,clr(match,-1,m);
for(int i=0;i<n;i++) g[i].clear();
}
void insert(int u,int v){g[u].push_back(v);}
bool augment(int u){//增广
for(int i=0,v;i<g[u].size();i++){
if(!vis[v=g[u][i]]){
vis[v]=1;
if(match[v]==-1 || augment(match[v])){
match[v]=u;
return true;
}
}
}
return false;
}
int maxMatch(){
int res=0;
for(int i=0;i<n;i++){
clr(vis,0,m);
if(augment(i)) res++;
}
return res;
}
};

struct KM{//邻接表
typedef int type;
struct edge{
int from,to;
type w;
edge(int u,int v,type _w){from=u,to=v,w=_w;}
};
int n,m;            // X/Y集的点数
int match[M];       // 匹配
type lx[M],ly[M],d; // X/Y项标,松弛量
bool vx[M],vy[M];   // X/Y集的点是否在已经访问
vector<edge>edg;    // 边表
vector<int>g[M];    // 邻接表
void init(int _n,int _m){
n=_n,m=_m;
clr(match,-1,m),edg.clear();
for(int i=0;i<n;i++) g[i].clear();
}
void insert(int u,int v,type w){
edg.push_back(edge(u,v,w));
g[u].push_back(edg.size()-1);
}
bool augment(int u){// 增广
vx[u]=1;
for(int i=0;i<g[u].size();i++){
edge& e=edg[g[u][i]];
if(vy[e.to]) continue;
if(lx[u]+ly[e.to]==e.w){
vy[e.to]=1;
if(match[e.to]==-1 || augment(match[e.to])){
match[e.to]=u;
return true;
}
}else getmin(d,lx[u]+ly[e.to]-e.w);
}
return false;
}
type maxWeight(){
int i,j;
for(i=0;i<m;i++) ly[i]=0;
for(i=0;i<n;i++){
lx[i]=-INF;
for(j=0;j<g[i].size();j++)
getmax(lx[i],edg[g[i][j]].w);
}
for(i=0;i<n;i++){
for(;;){
clr(vx,0,n),clr(vy,0,m);
d=-1;
if(augment(i)) break;
for(j=0;j<n;j++) if(vx[j]) lx[j]-=d;
for(j=0;j<m;j++) if(vy[j]) ly[j]+=d;
}
}
type res=0;
for(i=0;i<n;i++) res+=lx[i];
for(i=0;i<m;i++) res+=ly[i];
return res;
}
};

struct KM{//邻接矩阵
typedef int type;
int n,m,match[M];          // X/Y集的点数,匹配
type lx[M],ly[M],w[M][M],d;// X/Y项标,松弛量
bool vx[M],vy[M];          // X/Y集的点是否在已经访问
void init(int _n,int _m){
n=_n,m=_m,clr(match,-1,m);
for(int i=0;i<n;i++) for(int j=0;j<m;j++) w[i][j]=-INF;
}
bool augment(int u){//增广
vx[u]=1;
for(int v=0;v<m;v++){
if(vy[v]) continue;
if(lx[u]+ly[v]==w[u][v]){
vy[v]=1;
if(match[v]==-1 || augment(match[v])){
match[v]=u;
return true;
}
}else getmin(d,lx[u]+ly[v]-w[u][v]);
}
return false;
}
type maxWeight(){
int i,j;
for(i=0;i<m;i++) ly[i]=0;
for(i=0;i<n;i++){
lx[i]=-INF;
for(j=0;j<m;j++)
getmax(lx[i],w[i][j]);
}
for(i=0;i<n;i++){
for(;;){
clr(vx,0,n),clr(vy,0,m);
d=-1;
if(augment(i)) break;
for(j=0;j<n;j++) if(vx[j]) lx[j]-=d;
for(j=0;j<m;j++) if(vy[j]) ly[j]+=d;
}
}
type res=0;
for(i=0;i<n;i++) res+=lx[i];
for(i=0;i<m;i++) res+=ly[i];
return res;
}
};