我们发现,对于一个有向图,他若是形成了一个换。那么我们把他拆点,成二分图。
会发现这个二分图一定会形成一个完备匹配。
那么这题就明了了。我们建立二分图,其中边的权值取负数,然后跑一遍KM就可以了。
我邻接表用多了,没怎么用邻接矩阵。这里刚开始忘给矩阵付初始的值-inf了。。。。。。
代码如下
#include<iostream>
#include<algorithm>
#include<queue>
#include<cstring>
using namespace std;
const int max_n = 250;
const int inf = 1e9;
int nx,ny;
int linker[max_n],lx[max_n],ly[max_n],slack[max_n];
int visx[max_n],visy[max_n],w[max_n][max_n];
int DFS(int x){
visx[x]=1;
for(int y=1;y<=ny;y++){
if(visy[y])
continue;
int tmp=lx[x]+ly[y]-w[x][y];
if(tmp==0){
visy[y]=1;
if(linker[y]==-1 || DFS(linker[y])){
linker[y]=x;
return 1;
}
}else if(slack[y]>tmp){
slack[y]=tmp;
}
}
return 0;
}
int KM(){
int i,j;
memset(linker,-1,sizeof(linker));
memset(ly,0,sizeof(ly));
for(i=1;i<=nx;i++)
for(j=1,lx[i]=-inf;j<=ny;j++)
if(w[i][j]>lx[i])
lx[i]=w[i][j];
for(int x=1;x<=nx;x++){
for(i=1;i<=ny;i++)
slack[i]=inf;
while(1){
memset(visx,0,sizeof(visx));
memset(visy,0,sizeof(visy));
if(DFS(x))
break;
int d=inf;
for(i=1;i<=ny;i++)
if(!visy[i] && d>slack[i])
d=slack[i];
for(i=1;i<=nx;i++)
if(visx[i])
lx[i]-=d;
for(i=1;i<=ny;i++)
if(visy[i])
ly[i]+=d;
else
slack[i]-=d;
}
}
int res=0;
for(i=1;i<=ny;i++)
if(linker[i]!=-1)
res+=w[linker[i]][i];
return res;
}
int n,m;
int main(){
int T;scanf("%d",&T);
while (T--){
scanf("%d %d",&n,&m);
for (int i=1;i<=n;++i)for (int j=1;j<=n;++j)w[i][j]=-inf;
for (int i=1;i<=m;++i){
int u,v,c;
scanf("%d %d %d",&u,&v,&c);
w[u][v]=max(-c,w[u][v]);
}nx=ny=n;
printf("%d\n",-KM());
}
}
京公网安备 11010502036488号