P2050 [NOI2012]美食节 动态加边加点

修车数据加强版 需要动态加边加点

#include<bits/stdc++.h>
using namespace std;
const int INF = 0x7f7f7f7f;
const int MAXN = , MAXM = ;
int Head[MAXN], cur[MAXN], lev[MAXN], to[MAXM << ], nxt[MAXM << ], f[MAXM << ], mono[MAXM << ], ed = , S, T;
int pre[MAXN];
int t[][];
int n, m, sum, now;
bool exist[MAXN];
void addedge(int u, int v, int cap, int val) {
        to[++ed] = v;
        nxt[ed] = Head[u];
        Head[u] = ed;
        f[ed] = cap;
        mono[ed] = val;
        to[++ed] = u;
        nxt[ed] = Head[v];
        Head[v] = ed;
        f[ed] = ;
        mono[ed] = - * val;
        return;
}
bool BFS() {
        int u;
        queue<int>q;
        memset(exist, false, sizeof(exist));
        memset(lev, , sizeof(lev));
        lev[S] = pre[S] = ;
        q.push(S);
        while (q.size()) {
                u = q.front();
                q.pop();
                exist[u] = false;
                for (int i = Head[u]; i; i = nxt[i])
                        if (f[i] && lev[u] + mono[i] < lev[to[i]]) {
                                lev[to[i]] = lev[u] + mono[i];
                                pre[to[i]] = i;
                                if (!exist[to[i]]) {
                                        exist[to[i]] = true;
                                        q.push(to[i]);
                                }
                        }
        }
        memcpy(cur, Head, sizeof(Head));
        return lev[T] != INF;
}
int DFS(int u, int maxf) {
        if (u == T || !maxf) {
                return maxf;
        }
        exist[u] = true;
        int cnt = ;
        for (int &i = cur[u], tem; i; i = nxt[i])
                if (f[i] && lev[u] + mono[i] == lev[to[i]]) {
                        if (exist[to[i]]) {
                                continue;
                        }
                        tem = DFS(to[i], min(f[i], maxf));
                        maxf -= tem;
                        f[i] -= tem;
                        f[i ^ ] += tem;
                        cnt += tem;
                        if (!maxf) {
                                break;
                        }
                }
        if (!cnt) {
                lev[u] = - * INF;
        }
        exist[u] = false;
        return cnt;
}
int Augment() {
        int delta = INF;
        for (int i = pre[T]; i; i = pre[to[i ^ ]])
                if (f[i] < delta) {
                        delta = f[i];
                }
        for (int i = pre[T]; i; i = pre[to[i ^ ]]) {
                f[i] -= delta;
                f[i ^ ] += delta;
        }
        return delta * lev[T];
}
int MCMF() {
        int ans = ;
        memset(exist, false, sizeof(exist));
        while (BFS())
                //ans+=DFS(S,INF)*lev[T];
        {
                int j = to[pre[T] ^ ];
                addedge(j + , T, , );
                j++;
                for (int i = ; i <= n; i++) {
                        int x = (j - n - ) / sum + ;
                        int y = j - n - (x - ) * sum;
                        addedge(i, j, , y * t[x][i]);
                }
                ans += Augment();
        }
        return ans;
}
int main() {
        sum = ;
        scanf("%d %d", &n, &m);
        S = ;
        for (int i = ; i <= n; i++) {
                scanf("%d", &now);
                sum += now;
                addedge(S, i, now, );
        }
        T = n + sum * m + ;
        for (int i = n + ; i <= T - ; i += sum) {
                addedge(i, T, , );
        }
        for (int i = ; i <= n; i++) {
                for (int j = ; j <= m; j++) {
                        scanf("%d", &t[j][i]);
                }
        }
        for (int i = ; i <= n; i++) {
                for (int j = n + ; j <= T - ; j += sum) {
                        int x = (j - n - ) / sum + ;
                        int y = j - n - (x - ) * sum;
                        addedge(i, j, , y * t[x][i]);
                }
        }
        int ans = MCMF();
        cout << ans << endl;
        return ;
}