经验首页 前端设计 程序设计 Java相关 移动开发 数据库/运维 软件/图像 大数据/云计算 其他经验
当前位置:技术经验 » 程序设计 » C++ » 查看文章
BZOJ1927: [Sdoi2010]星际竞速(最小费用最大流 最小路径覆盖)
来源:cnblogs  作者:自为风月马前卒  时间:2018/12/7 9:33:20  对本文有异议

题意

题目链接

Sol

看完题不难想到最小路径覆盖,但是带权的咋做啊?qwqqq

首先冷静思考一下:最小路径覆盖 = \(n - \text{二分图最大匹配数}\)

为什么呢?首先最坏情况下是用\(n\)条路径去覆盖(就是\(n\)个点),根据二分图的性质,每个点只能有一个和他配对,这样就保证了,每多出一个匹配,路径数就会\(-1\)

扩展到有边权的图也是同理的,\(i\)表示二分图左侧的点,\(i'\)表示二分图右侧的点,对于两点\(u, v\),从\(u\)\(v'\)\((1, w_i)\)的边(前面是流量,后面是费用)

接下来从\(S\)\(i\)\((1, 0)\)的边,从\(i'\)\(T\)\((1, 0)\)的边,从\(S\)\(i'\)\((1, A_i)\)的边

跑最小费用最大流即可

  1. #include<bits/stdc++.h>
  2. #define chmin(x, y) (x = x < y ? x : y)
  3. #define chmax(x, y) (x = x > y ? x : y)
  4. using namespace std;
  5. const int MAXN = 1e6 + 10, INF = 1e9 + 10;
  6. inline int read() {
  7. char c = getchar(); int x = 0, f = 1;
  8. while(c < '0' || c > '9') {if(c == '-') f = -1; c = getchar();}
  9. while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
  10. return x * f;
  11. }
  12. int N, M, S, T, dis[MAXN], vis[MAXN], Pre[MAXN], ansflow, anscost, A[MAXN];
  13. struct Edge {
  14. int u, v, f, w, nxt;
  15. }E[MAXN];
  16. int head[MAXN], num = 0;
  17. inline void add_edge(int x, int y, int f, int w) {
  18. E[num] = (Edge) {x, y, f, w, head[x]}; head[x] = num++;
  19. }
  20. inline void AddEdge(int x, int y, int f, int w) {
  21. add_edge(x, y, f, w); add_edge(y, x, 0, -w);
  22. }
  23. bool SPFA() {
  24. queue<int> q; q.push(S);
  25. memset(dis, 0x3f, sizeof(dis));
  26. memset(vis, 0, sizeof(vis));
  27. dis[S] = 0;
  28. while(!q.empty()) {
  29. int p = q.front(); q.pop(); vis[p] = 0;
  30. for(int i = head[p]; ~i; i = E[i].nxt) {
  31. int to = E[i].v;
  32. if(E[i].f && dis[to] > dis[p] + E[i].w) {
  33. dis[to] = dis[p] + E[i].w; Pre[to] = i;
  34. if(!vis[to]) vis[to] = 1, q.push(to);
  35. }
  36. }
  37. }
  38. return dis[T] <= INF;
  39. }
  40. void F() {
  41. int canflow = INF;
  42. for(int i = T; i != S; i = E[Pre[i]].u) chmin(canflow, E[Pre[i]].f);
  43. for(int i = T; i != S; i = E[Pre[i]].u) E[Pre[i]].f -= canflow, E[Pre[i] ^ 1].f += canflow;
  44. anscost += canflow * dis[T];
  45. }
  46. void MCMF() {
  47. while(SPFA()) F();
  48. }
  49. int main() {
  50. memset(head, -1, sizeof(head));
  51. N = read(); M = read(); S = N * 2 + 2, T = N * 2 + 3;
  52. for(int i = 1; i <= N; i++) A[i] =read(), AddEdge(S, i, 1, 0), AddEdge(i + N, T, 1, 0), AddEdge(S, i + N, 1, A[i]);
  53. for(int i = 1; i <= M; i++) {
  54. int x = read(), y = read(), v = read();
  55. if(x > y) swap(x, y);
  56. AddEdge(x, y + N, 1, v);
  57. }
  58. MCMF();
  59. printf("%d\n", anscost);
  60. return 0;
  61. }
 友情链接:直通硅谷  点职佳  北美留学生论坛

本站QQ群:前端 618073944 | Java 606181507 | Python 626812652 | C/C++ 612253063 | 微信 634508462 | 苹果 692586424 | C#/.net 182808419 | PHP 305140648 | 运维 608723728

W3xue 的所有内容仅供测试,对任何法律问题及风险不承担任何责任。通过使用本站内容随之而来的风险与本站无关。
关于我们  |  意见建议  |  捐助我们  |  报错有奖  |  广告合作、友情链接(目前9元/月)请联系QQ:27243702 沸活量
皖ICP备17017327号-2 皖公网安备34020702000426号