二分匹配之最大权值匹配算法---KM模板

百科：http://baike.baidu.com/link?url=vbM3H4XmfrsWfP-epdlR2sVKSNzOq4hXnWDqm5uo8fd7VWsF2SmhDV35XyVUDvVjvrtf42RUITJuNCHn-7_x6K

大神总结：http://www.cnblogs.com/skyming/archive/2012/02/18/2356919.html

代码：

 #include<stdio.h>
 #include<string.h>
 const int N=,INF=0x3f3f3f3f;
 int lx[N],ly[N],vx[N],vy[N],slack[N],match[N];
 int a[N][N];
 int nx,ny;
 int dfs(int u)
 {
     vx[u]=;
     for(int i=;i<=ny;i++)
     {
         if(vy[i])    continue;
         int t=lx[u]+ly[i]-a[u][i];
         if(t==)
         {
             vy[i]=;
             if(match[i]==-||dfs(match[i]))
             {
                 match[i]=u;
                 return ;
             }
         }
         else if(slack[i]>t)
             slack[i]=t;
     }
     return ;
 }
 int KM()
 {
     int i,j,d;
     memset(ly,,sizeof(ly));
     memset(match,-,sizeof(match));
     for(i=,lx[i]=-INF;i<=nx;i++)
         for(j=;j<=ny;j++)
             if(a[i][j]>lx[i])
                 lx[i]=a[i][j];
     for(i=;i<=nx;i++)
     {
         memset(slack,0x3f,sizeof(slack));
         while()
         {
             memset(vx,,sizeof(vx));
             memset(vy,,sizeof(vy));
             if(dfs(i))    break;
             d=INF;
             for(j=;j<=ny;j++)
                 if(!vy[j]&&slack[j]<d)
                     d=slack[j];
             //if(d==INF)    break;//该点找不到任何匹配
             for(j=;j<=nx;j++)
                 if(vx[j])
                     lx[j]-=d;
             for(j=;j<=ny;j++)
             {
                 if(vy[j])
                     ly[j]+=d;
                 else
                     slack[j]-=d;
             }
         }
     }
     int ans=,sum=;
     for(i=;i<=ny;i++)
         if(match[i]>-&&a[match[i]][i])
             ans+=a[match[i]][i],sum++;
     printf("%d++\n",sum);//匹配数
     return ans;
 }
 int main()
 {
     int m,u,v,w,n,i;
     scanf("%d",&nx);
     scanf("%d",&ny);
     scanf("%d",&m);
     memset(a,,sizeof(a));
     while(m--)
     {
         scanf("%d%d%d",&u,&v,&w);
         a[u][v]=a[v][u]=w;
     }
     printf("%d\n",KM());
     return ;
 }