UOJ#80. 二分图最大权匹配 模板

#80. 二分图最大权匹配

• 描述
• 提交
• 自定义测试

从前一个和谐的班级，有 nlnl 个是男生，有 nrnr 个是女生。编号分别为 1,…,nl1,…,nl 和 1,…,nr1,…,nr。

有若干个这样的条件：第 vv 个男生和第 uu 个女生愿意结为配偶，且结为配偶后幸福程度为 ww。

请问这个班级里幸福程度之和最大是多少？

输入格式

第一行三个正整数，nl,nr,mnl,nr,m。

接下来 mm 行，每行三个整数 v,u,wv,u,w 表示第 vv 个男生和第 uu 个女生愿意结为配偶，且幸福程度为 ww。保证 1≤v≤nl1≤v≤nl，1≤u≤nr1≤u≤nr，保证同一对 v,uv,u 不会出现两次。

输出格式

第一行一个整数，表示幸福程度之和的最大值。

接下来一行 nlnl 个整数，描述一组最优方案。第 vv 个整数表示 vv 号男生的配偶的编号。如果 vv 号男生没配偶请输出 00。

样例一

input

2 2 3
1 1 100
1 2 1
2 1 1

output

100
1 0

限制与约定

1≤nl,nr≤4001≤nl,nr≤400，1≤m≤1600001≤m≤160000，1≤w≤1091≤w≤109。

时间限制：1s1s

空间限制：256MB256MB

下载

样例数据下载

http://blog.csdn.net/c20180630/article/details/71080521这篇讲的不错

 #include <iostream>
 #include <cstdio>
 #include <cstring>
 #include <cstdlib>
 #include <algorithm>
 #include <queue>
 #include <vector>
 #include <map>
 #include <string>
 #include <cmath>
 #define min(a, b) ((a) < (b) ? (a) : (b))
 #define max(a, b) ((a) > (b) ? (a) : (b))
 #define abs(a) ((a) < 0 ? (-1 * (a)) : (a))
 template<class T>
 inline void swap(T &a, T &b)
 {
     T tmp = a;a = b;b = tmp;
 }
 inline void read(long long &x)
 {
     x = ;char ch = getchar(), c = ch;
     while(ch < '' || ch > '') c = ch, ch = getchar();
     while(ch <= '' && ch >= '') x = x *  + ch - '', ch = getchar();
     if(c == '-') x = -x;
 }

 const long long INF = 0x3f3f3f3f3f3f3f3f;
 const long long MAXN =  + ;

 long long g[MAXN][MAXN], n1, n2, m, tmp1, tmp2, lab1[MAXN], lab2[MAXN], lk1[MAXN], lk2[MAXN], pre[MAXN], sla[MAXN], vis[MAXN];
 void fac(long long x)
 {
     memset(vis, , sizeof(vis)), memset(pre, , sizeof(pre)), memset(sla, 0x7f, sizeof(sla)), vis[] = ;
     long long y;
     do
     {
         y = ;
         for(long long i = ;i <= n2;++ i)
         {
             if(vis[i]) continue;
             if(lab1[x] + lab2[i] - g[x][i] < sla[i]) sla[i] = lab1[x] + lab2[i] - g[x][i], pre[i] = x;
             if(sla[i] < sla[y]) y = i;
         }
         long long d = sla[y];
         for(long long i = ;i <= n1;++ i) if(vis[lk1[i]]) lab1[i] -= d;
         for(long long i = ;i <= n2;++ i) if(vis[i]) lab2[i] += d; else sla[i] -= d;
         vis[y] = ;
     }while(x = lk2[y]);
     for(;y;swap(y, lk1[lk2[y] = pre[y]]));
 }
 long long KM()
 {
     for(long long i = ;i <= n1;++ i) fac(i);
     long long ans = ;
     for(long long i = ;i <= n1;++ i) ans += g[i][lk1[i]];
     return ans;
 }
 int main()
 {
     read(n1), read(n2), n2 = max(n2, n1), read(m);
     for(long long i = ;i <= m;++ i)read(tmp1), read(tmp2), read(g[tmp1][tmp2]), lab1[tmp1] = max(g[tmp1][tmp2], lab1[tmp1]);
     printf("%lld\n", KM());
     for(long long i = ;i <= n1;++ i) printf("%lld ", g[i][lk1[i]] ? lk1[i] : );
     return ;
 }

UPJ#80