西南科技大学院赛 I题 上决╇ф的精确打击问题 【最大流好题】 建模问题

传送门
题意: 就是一个矩阵中有些点上面有一些鸟, 然后你可以对每一行或者每一列进行开枪, 会打死这一行或者这一列上面所有的鸟, 然后每一行和每一列都有一个花费, 如果在相应的行(列)开枪就必须付出相应的花费, 并且每只鸟自己也有一个花费, 问打死所有的鸟最小的花费是多少.

思路： 直接进行一个建图, 那就是建一个源点和汇点, 然后每一个行也是一个点, 每一个列也是一个点, 源点到每一行的点的容量就是对应行的花费， 列点到汇点的容量就是对应列的花费, 然后对于有鸟的点, 其所在的行和列标的点之间的建一条边, 边的容量是对应点上面的花费, 然后从源点到汇点跑一遍最大流就是答案…. (可以好好想想, 正确性是一定的)

AC Code

const int maxn = 1e4+5;
int head[maxn], cnt;
int n, m, s, t;
bool vis[maxn];
int d[maxn], cur[maxn];
struct node {
    int u, v;
    int cap, flow, next;
} e[maxn*10]; //因为是双向边  所以记得开二倍

void init() {
    cnt = 0;
    Fill(head, -1);
}

void add(int u, int v, int c, int f) {
    e[cnt] = node{u, v, c, f, head[u]};
    head[u] = cnt++;
}

bool BFS() {
    queue<int> q; q.push(s);
    Fill(vis,0); vis[s] = 1;
    d[s] = 0;
    while (!q.empty()) {
        int u = q.front(); q.pop();
        for (int i = head[u] ; ~i ; i = e[i].next) {
            int to = e[i].v;
            if (!vis[to] && e[i].cap > e[i].flow) { //只考虑残量网络的弧
                vis[to] = 1;
                d[to] = d[u] + 1;
                q.push(to);
            }
        }
    }
    return vis[t];
}

int dfs(int x, int a) {
    if (x == t || a == 0) return a;
    int flow = 0, f;
    int &i = cur[x]; // 这个不能少, 表示形参改变对应实参的改变.
    for (i = head[x] ; ~i ; i = e[i].next) { //从上次考虑的弧
        int to = e[i].v;
        if (d[to]== d[x] + 1 && (f = dfs(to, min(a, e[i].cap - e[i].flow))) > 0) {
            e[i].flow += f;
            e[i^1].flow -= f;
            flow += f;
            a -= f;
            if (a == 0) break;
        }
    }
    return flow;
}

int Dinic() {
    int maxflow = 0;
    while (BFS()) {
        Fill(cur, 0);
        maxflow += dfs(s, inf);
    }
    return maxflow;
}
void solve() {
    int k;
    scanf("%d%d%d", &n, &m, &k);
    s = 0, t = n*m + 1; init();
    for (int i = 1 ; i <= n ; i ++) {
        int x; scanf("%d", &x);
        add(s, i, x, 0);
        add(i, s, 0, 0);
    }
    for (int i = 1 ; i <= m ; i ++) {
        int x; scanf("%d", &x);
        add(i+n, t, x, 0);
        add(t, i+n, 0, 0);
    }
    for (int i = 1 ; i <= k ; i ++) {
        int x, y, w;
        scanf("%d%d%d", &x, &y, &w);
        add(x, y+n, w, 0);
        add(y+n, x, 0, 0);
    }
    printf("%d\n", Dinic());
}