题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=3488
思路：KM算法求最小权。。。
ps:任何完全图的二分匹配其实就是n个环的并
View Code
1 #include<iostream> 2 const int MAXN=222; 3 const int inf=1<<30; 4 using namespace std; 5 int n,m; 6 int map[MAXN][MAXN]; 7 int lx[MAXN],ly[MAXN]; 8 int match[MAXN]; 9 bool visitx[MAXN],visity[MAXN]; 10 11 int Hungary(int u){ 12 visitx[u]=true; 13 for(int i=1;i<=n;i++){ 14 if(!visity[i]&&lx[u]+ly[i]==map[u][i]){ 15 visity[i]=true; 16 if(match[i]==-1||Hungary(match[i])){ 17 match[i]=u; 18 return true; 19 } 20 } 21 } 22 return false; 23 } 24 25 26 void KM_prefect_match(){ 27 int tmp; 28 for(int i=1;i<=n;i++){ 29 lx[i]=-inf; //初始化应该赋为无穷小 30 } 31 memset(ly,0,sizeof(ly)); 32 for(int i=1;i<=n;i++){ 33 for(int j=1;j<=n;j++){ 34 lx[i]=max(lx[i],map[i][j]); 35 } 36 } 37 for(int i=1;i<=n;i++) 38 { 39 while(1){ 40 memset(visitx,false,sizeof(visitx)); 41 memset(visity,false,sizeof(visity)); 42 if(Hungary(i)) 43 break; 44 else { 45 tmp=inf; 46 for(int j=1;j<=n;j++)if(visitx[j]){ //x在交错树中 47 for(int k=1;k<=n;k++){ 48 if(!visity[k]&&tmp>lx[j]+ly[k]-map[j][k]){ 49 tmp=lx[j]+ly[k]-map[j][k]; 50 } 51 } 52 } 53 for(int j=1;j<=n;j++){ 54 if(visitx[j]) 55 lx[j]-=tmp; 56 if(visity[j]) 57 ly[j]+=tmp; 58 } 59 } 60 } 61 } 62 } 63 64 65 int main(){ 66 int _case; 67 scanf("%d",&_case); 68 while(_case--){ 69 scanf("%d%d",&n,&m); 70 memset(match,-1,sizeof(match)); 71 for(int i=1;i<=n;i++){ 72 for(int j=1;j<=n;j++){ 73 map[i][j]=-inf; 74 } 75 } 76 for(int i=1;i<=m;i++){ 77 int x,y,w; 78 scanf("%d%d%d",&x,&y,&w); 79 if(map[x][y]<-w){ 80 map[x][y]=-w; 81 } 82 } 83 KM_prefect_match(); 84 int ans=0; 85 for(int i=1;i<=n;i++){ 86 ans+=map[match[i]][i]; 87 } 88 printf("%d\n",-ans); 89 } 90 return 0; 91 }
题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=3435
View Code
1 #include<iostream> 2 const int MAXN=1010; 3 const int inf=1<<30; 4 using namespace std; 5 int n,m; 6 int map[MAXN][MAXN]; 7 int match[MAXN]; 8 int lx[MAXN],ly[MAXN]; 9 bool visitx[MAXN],visity[MAXN]; 10 11 int Hungary(int u){ 12 visitx[u]=true; 13 for(int i=1;i<=n;i++){ 14 if(!visity[i]&&lx[u]+ly[i]==map[u][i]){ 15 visity[i]=true; 16 if(match[i]==-1||Hungary(match[i])){ 17 match[i]=u; 18 return true; 19 } 20 } 21 } 22 return false; 23 } 24 25 26 void KM_prefect_match(){ 27 int tmp; 28 for(int i=1;i<=n;i++){ 29 lx[i]=-inf; 30 } 31 memset(ly,0,sizeof(ly)); 32 for(int i=1;i<=n;i++){ 33 for(int j=1;j<=n;j++){ 34 lx[i]=max(lx[i],map[i][j]); 35 } 36 } 37 for(int i=1;i<=n;i++) 38 { 39 while(1){ 40 memset(visitx,false,sizeof(visitx)); 41 memset(visity,false,sizeof(visity)); 42 if(Hungary(i)) 43 break; 44 else { 45 tmp=inf; 46 for(int j=1;j<=n;j++)if(visitx[j]){ 47 for(int k=1;k<=n;k++){ 48 if(!visity[k]&&tmp>lx[j]+ly[k]-map[j][k]){ 49 tmp=lx[j]+ly[k]-map[j][k]; 50 } 51 } 52 } 53 for(int j=1;j<=n;j++){ 54 if(visitx[j]) 55 lx[j]-=tmp; 56 if(visity[j]) 57 ly[j]+=tmp; 58 } 59 } 60 } 61 } 62 } 63 64 int main(){ 65 int T,_case=1; 66 scanf("%d",&T); 67 while(T--){ 68 scanf("%d%d",&n,&m); 69 memset(match,-1,sizeof(match)); 70 for(int i=1;i<=n;i++){ 71 for(int j=1;j<=n;j++){ 72 map[i][j]=-inf; 73 } 74 } 75 for(int i=1;i<=m;i++){ 76 int x,y,w; 77 scanf("%d%d%d",&x,&y,&w); 78 if(map[x][y]<-w){ 79 map[x][y]=map[y][x]=-w; 80 } 81 } 82 KM_prefect_match(); 83 printf("Case %d: ",_case++); 84 int ans=0; 85 int flag=false; 86 for(int i=1;i<=n;i++){ 87 if(match[i]==-1||map[match[i]][i]==-inf){ 88 flag=true; 89 break; 90 } 91 ans+=map[match[i]][i]; 92 } 93 if(flag){ 94 printf("NO\n"); 95 }else 96 printf("%d\n",-ans); 97 } 98 return 0; 99 }