zoukankan      html  css  js  c++  java
  • [BZOJ1877][SDOI2009]晨跑[最大流+费用流]

    天数最多 长度最小 天数是流量,长度是费用

    每个点拆成两个点限流1,就能保证只走一次 然后跑费用流

    #include <bits/stdc++.h>
    using namespace std;
    #define MP make_pair
    #define pb push_back
    #define read2(a, b) (read(a), read(b))
    #define read3(a, b, c) (read(a), read(b), read(c))
    #define lop(i,a,b) for(register int i = (a); i <= (b); ++i)
    #define dlop(i,a,b) for(register int i = (a); i >= (b); --i)
    #define eps (1e-7)
    #define fir first
    #define sec second
        
    const int inf = 0x3f3f3f3f-1;
    const int MAXN = 4e2+7;
    const int MAXM = 8e4+7;
    typedef long long LL;
    typedef long double LF;
    typedef pair<int,int> pii;
    typedef unsigned long long ULL;
    typedef unsigned int uint;
    template<class T> void read(T & x) {
      register int c = getchar(), f = 1;x = 0;
      while(!isdigit(c)) {if (c == '-') f = -f;c = getchar();}
      while(isdigit(c)) x = x * 10 + c - '0', c = getchar();
      x *= f;
    }
        
    int n, m, k, maxflow, dep[MAXN], s, t, cur[MAXN], head[MAXN], ans, dis[MAXN];
    bool inq[MAXN];
        
    struct Edge {
      int u, v, w, c, next;
    } G[MAXM]; int tot = 1;
    inline void add(int u, int v, int w, int c) {
      // cout << u << ' ' << v << ' ' << w << endl;
      G[++tot] = (Edge){u, v, w, c, head[u]}; head[u] = tot;
      G[++tot] = (Edge){v, u, 0, -c, head[v]}; head[v] = tot;
    } 
    bool spfa(int s, int t) {
      memset(dis, 0x3f, sizeof dis);
      memset(inq, 0, sizeof inq);
      queue<int>q; q.push(s); dis[s] = 0, inq[s] = 1;
      while(!q.empty()) {
        int u = q.front();
        inq[u] = 0;
        q.pop();
        for(int i = head[u]; i; i = G[i].next) {
          int v = G[i].v, w = G[i].w, c = G[i].c;
          if (dis[v] > dis[u] + c && w) {
            dis[v] = dis[u] + c;
            cur[v] = i;
            if (!inq[v]) q.push(v), inq[v] = 1;
          }
        }
      }
      return dis[t] < inf; 
    }
      
    void update(int s, int t) {
      int i, x = inf;
      for(i = cur[t]; i; x = min(x, G[i].w), i = cur[G[i].u]);
      for(i = cur[t]; i; G[i].w -= x, G[i^1].w += x, ans += x * G[i].c, i = cur[G[i].u]);
      maxflow += x;
    }
     
    void EK(int s, int t) {
      while(spfa(s, t)) update(s, t);
    }
     
    int main(void) {
      // freopen("data.in", "r", stdin);
      read2(n, m);
      for(int i = 2; i < n; ++i) add(i, i+n, 1, 0);
      for(int u, v, w, i = 1; i <= m; ++i) {
        read3(u, v, w); add(u+n, v, 1, w);
      }
      EK(1+n, n);
      cout << maxflow << ' ' << ans;
      return 0;
    }
    
  • 相关阅读:
    一个优秀的博主 有助于前端进阶
    JS异步的理解以及promise
    node-sass出错
    vuecli脚手架搭建
    GitHub上传项目方法
    ES6转ES5
    webpack---安装与卸载
    CSS之BFC详解
    Vue.js学习 --环境搭建
    python is和==的区别
  • 原文地址:https://www.cnblogs.com/storz/p/10191475.html
Copyright © 2011-2022 走看看