【题解】Catering World Finals 2015 上下界费用流

Prelude

传送到Codeforces：0.0

---

Solution

板子题，在这里贴个板子。
这题面是smg？题面中有说每个点只能经过一次吗？是我瞎了吗？
因为这WA on test 27一个小时，烦死了，浪费时间。

---

Code

#include 
#include 
#include 
#include 
#include 

using namespace std;
const int MAXN = 1000;
const int MAXM = 100000;
const int INF = 0x3f3f3f3f;
int _w;

int read() {
    int x;
    _w = scanf( "%d", &x );
    return x;
}

namespace MCMF {
    struct Edge {
        int u, v, c, f, w;
        Edge() {}
        Edge( int u, int v, int c, int f, int w ):
            u(u), v(v), c(c), f(f), w(w) {}
    };
    
    int n, m, s, t;
    int head[MAXN], nxt[MAXM];
    Edge edge[MAXM];
    
    void init( int _n ) {
        n = _n, m = 0;
        for( int i = 0; i < n; ++i )
            head[i] = -1;
    }
    void adde( int u, int v, int c, int w ) {
        edge[m] = Edge(u, v, c, 0, w);
        nxt[m] = head[u], head[u] = m++;
        edge[m] = Edge(v, u, 0, 0, -w);
        nxt[m] = head[v], head[v] = m++;
    }
    
    queue q;
    int dis[MAXN], inq[MAXN], res[MAXN], from[MAXN];
    
    bool spfa() {
        for( int i = 0; i < n; ++i )
            dis[i] = INF, inq[i] = 0;
        dis[s] = 0, inq[s] = 1, q.push(s), res[s] = INF;
        while( !q.empty() ) {
            int u = q.front(); q.pop();
            inq[u] = 0;
            for( int i = head[u]; ~i; i = nxt[i] ) {
                const Edge &e = edge[i];
                if( e.c > e.f && dis[u] + e.w < dis[e.v] ) {
                    dis[e.v] = dis[u] + e.w;
                    from[e.v] = i;
                    res[e.v] = min( res[u], e.c-e.f );
                    if( !inq[e.v] )
                        inq[e.v] = 1, q.push(e.v);
                }
            }
        }
        return dis[t] != INF;
    }
    void augment() {
        int f = res[t], u = t;
        while( u != s ) {
            int i = from[u];
            edge[i].f += f;
            edge[i^1].f -= f;
            u = edge[i].u;
        }
    }
    int solve( int _s, int _t ) {
        s = _s, t = _t;
        int cost = 0;
        while( spfa() ) {
            cost += res[t] * dis[t];
            augment();
        }
        return cost;
    }
}

int n, k, g[MAXN][MAXN];
int s, t, ss, tt, nid;
int in[MAXN], out[MAXN];

void adde( int u, int v, int l, int r, int w ) {
    MCMF::adde(ss, v, l, 0);
    MCMF::adde(u, tt, l, 0);
    MCMF::adde(u, v, r-l, w);
}

void solve() {
    s = nid++, t = nid++, ss = nid++, tt = nid++;
    for( int i = 1; i <= n; ++i )
        in[i] = nid++, out[i] = nid++;
    MCMF::init(nid);
    adde(t, s, 0, INF, 0);
    adde(s, out[1], 0, k, 0);
    for( int i = 2; i <= n; ++i ) {
        adde(out[i], t, 0, INF, 0);
        adde(in[i], out[i], 1, 1, 0); // 每个点只能经过一次
    }
    for( int i = 1; i <= n; ++i )
        for( int j = i+1; j <= n; ++j )
            adde(out[i], in[j], 0, INF, g[i][j]);
    printf( "%d\n", MCMF::solve(ss, tt) );
}

int main() {
    n = read()+1, k = read();
    for( int i = 1; i <= n; ++i )
        for( int j = i+1; j <= n; ++j )
            g[i][j] = read();
    solve();
    return 0;
}