bzoj 1877 最小费用流

思路：挺裸的费用流，拆拆点就好啦。

#include<bits/stdc++.h>
#define LL long long
#define fi first
#define se second
#define mk make_pair
#define PII pair<int, int>
#define y1 skldjfskldjg
#define y2 skldfjsklejg
using namespace std;

const int N =  + ;
const int M = 5e5 + ;
const int inf = 0x3f3f3f3f;
const LL INF = 0x3f3f3f3f3f3f3f3f;
const int mod = ;

int n, m, edgenum, S, T, head[N], dist[N], pre[N];
bool vis[N];

struct Edge {
    int from, to, cap, flow, cost, next;
} edge[M];

void init() {
    edgenum = ;
    memset(head, -, sizeof(head));
}

void addEdge(int u, int v, int w, int c) {
    Edge E1 = {u, v, w, , c, head[u]};
    edge[edgenum] = E1;
    head[u] = edgenum++;
    Edge E2 = {v, u, , , -c, head[v]};
    edge[edgenum] = E2;
    head[v] = edgenum++;
}

bool SPFA(int s, int t) {
    queue<int> Q;
    memset(dist, INF, sizeof(dist));
    memset(vis, false, sizeof(vis));
    memset(pre, -, sizeof(pre));
    dist[s] = ; vis[s] = true; Q.push(s);
    while(!Q.empty()) {
        int u = Q.front(); Q.pop(); vis[u] = false;
        for(int i = head[u]; i != -; i = edge[i].next) {
            Edge E = edge[i];
            if(dist[E.to] > dist[u] + E.cost && E.cap > E.flow) {
                dist[E.to] = dist[u] + E.cost;
                pre[E.to] = i;
                if(!vis[E.to]) {
                    vis[E.to] = true;
                    Q.push(E.to);
                }
            }
        }
    }
    return pre[t] != -;
}

void MCMF(int s, int t, int &cost, int &flow) {
    flow = ; cost = ;
    while(SPFA(s, t)) {
        int Min = INF;
        for(int i = pre[t]; i != -; i = pre[edge[i^].to]) {
            Edge E = edge[i];
            Min = min(Min, E.cap - E.flow);
        }
        for(int i = pre[t]; i != -; i = pre[edge[i^].to]) {
            edge[i].flow += Min;
            edge[i^].flow -= Min;
            cost += edge[i].cost * Min;
        }
        flow += Min;
    }
}
int main() {
    init();
    scanf("%d%d", &n, &m);
    S = , T =  * n;
    for(int i = ; i < n; i++) addEdge(i, i + n, , );
    addEdge(,  + n, inf, ); addEdge(n, n + n, inf, );
    for(int i = ; i <= m; i++) {
        int u, v, c;
        scanf("%d%d%d", &u, &v, &c);
        addEdge(u + n, v, , c);
    }

    int cost, flow;
    MCMF(S, T, cost, flow);
    printf("%d %d\n", flow, cost);
    return ;
}

/*
*/