2024江苏省赛 H. 完蛋，我被房产包围了 【费用流、分时图】

完蛋，我被房产包围了

$n\le 200,\sum n\le 10^4$

求出最大利润

思路

每个代理商每次买房狂潮只能卖出 $1$ 套房子，小红卖出一套房子贬值 $1$ 元，小绿卖出一套房子贬值 $\lceil \frac{a_i}{10}\rceil$，小蓝卖出一套房子不贬值。
因此我们应该尽可能让小蓝卖房子，其次再交给小红，最后没得选择才交给小绿，因为卖出房子总归是有利润的，只是多少的问题，所以最优的策略一定是尽可能卖更多数量的房子。

我们先将小红和小蓝占满，让他们每次狂潮都有房子卖，他们卖不完的房子再交给贬值最多的小绿，由于小绿优先选择价值最高的房子贩卖，所以我们把不需要卖的房子交给小绿，并不会影响我们最终的最大利润决策

通过上述的分析，我们不难发现：对于每一次 买房狂潮 ，我们应该尽可能跑满最大流，并且要最小化贬值，也就是最小费用最大流

我们可以这样建模：

1. 对于每个代理商 $j$，建立其在时刻 $i$ 的点：$id(i,j)$，这里总共有 $3n$ 个分时点
2. 对于每个代理商 $j$ 的每个时刻点，我们连边 $id(i-1,j)\to id(i,j)$，容量为 $\infty$，费用为 $0$，表示这个代理商手里持有的房子，随着时间推移而保留
3. 对于每次买房狂潮，我们对当前时刻 $i$ 的每个代理商 $j$ 连边：$id(i,j)\to T$，容量为 $1$，费用为 $0$，表示每个代理商在这个时刻最多贩卖一套房子
4. 对于每个房子 $x$，连边：$S\to x$，容量为 $1$，费用为 $0$，表示这个房子最多被卖 $1$ 次（这里新增 $O(n)$ 个点），连边：$x\to T$，容量为 $1$，费用为 $a_x$，表示这个房子不交给代理商，最后没有被卖出去的贬值，也即是一块钱利润都没有
5. 对于每个房子 $x$，我们在其分配的时刻 $i$，对每个代理商连边：
   $x\to id(i,0)$，容量为 $1$，费用为 $1$，表示通过小红卖这套房要贬值 $1$；
   $x\to id(i,1)$ ，容量为 $1$，费用为 $\lceil \frac{a_i}{10}\rceil$，表示通过小绿卖这套房要贬值 $\lceil \frac{a_i}{10}\rceil$；
   $x\to id(i,2)$ ，容量为 $1$，费用为 $0$，表示通过小蓝卖不贬值

最后我们用 $\sum a_i-(S\to T$ 的最小费用最大流) 就是答案了

#include<bits/stdc++.h>
#define fore(i,l,r)	for(int i=(int)(l);i<(int)(r);++i)
#define fi first
#define se second
#define endl '\n'
#define ull unsigned long long
#define ALL(v) v.begin(), v.end()
#define Debug(x, ed) std::cerr << #x << " = " << x << ed;

const int INF=0x3f3f3f3f;
const long long INFLL=1e18;

typedef long long ll;

struct MCF {
    struct Edge {
        int v, c, w; //边终点、容量、费用
        Edge(int v, int c, int w) : v(v), c(c), w(w) {}
    };
    const int n;
    std::vector<Edge> e;
    std::vector<std::vector<int>> g;
    std::vector<ll> h, dis;
    std::vector<int> pre;
    bool dijkstra(int s, int t) {
        dis.assign(n + 1, std::numeric_limits<ll>::max());
        pre.assign(n + 1, -1);
        std::priority_queue<std::pair<ll, int>, std::vector<std::pair<ll, int>>, std::greater<std::pair<ll, int>>> que;
        dis[s] = 0;
        que.emplace(0, s);
        while (!que.empty()) {
            ll d = que.top().first;
            int u = que.top().second;
            que.pop();
            if (dis[u] < d) continue;
            for (int i : g[u]) {
                int v = e[i].v;
                int c = e[i].c;
                int w = e[i].w;
                if (c > 0 && dis[v] > d + h[u] - h[v] + w) {
                    dis[v] = d + h[u] - h[v] + w;
                    pre[v] = i;
                    que.emplace(dis[v], v);
                }
            }
        }
        return dis[t] != std::numeric_limits<ll>::max();
    }
    MCF(int n) : n(n), g(n + 1) {}
    void addEdge(int u, int v, int c, int w) {
        g[u].push_back(e.size());
        e.emplace_back(v, c, w);
        g[v].push_back(e.size());
        e.emplace_back(u, 0, -w);
    }
    std::pair<int, ll> flow(int s, int t) {
        int flow = 0;
        ll cost = 0;
        h.assign(n + 1, 0);
        while (dijkstra(s, t)) {
            for (int i = 1; i <= n; ++i) h[i] += dis[i];
            int aug = std::numeric_limits<int>::max();
            for (int i = t; i != s; i = e[pre[i] ^ 1].v) aug = std::min(aug, e[pre[i]].c);
            for (int i = t; i != s; i = e[pre[i] ^ 1].v) {
                e[pre[i]].c -= aug;
                e[pre[i] ^ 1].c += aug;
            }
            flow += aug;
            cost += ll(aug) * h[t];
        }
        return std::make_pair(flow, cost);
    }
};

int main(){
    std::ios::sync_with_stdio(false);
    std::cin.tie(nullptr);
    std::cout.tie(nullptr);
    int t;
    std::cin >> t;
    while(t--){
        int n;
        std::cin >> n;
        MCF mcf(4 * n + 10);
        int tot = 3 * n + 1;
        int S = ++tot, T = ++tot;

        auto get_id = [&](int time, int j) -> int {
            return n * j + time;
        };

        int ans = 0;
        fore(i, 0, n){
            if(i > 0){
                fore(j, 0, 3){
                    int lst = get_id(i - 1, j), now = get_id(i, j);
                    mcf.addEdge(lst, now, INF, 0);
                }
            }

            int opt;
            std::cin >> opt;
            if(opt == 1){
                int w;
                std::cin >> w;
                ans += w;
                int id = ++tot;
                mcf.addEdge(S, id, 1, 0);
                mcf.addEdge(id, get_id(i, 0), 1, 1);
                mcf.addEdge(id, get_id(i, 1), 1, (w + 9) / 10);
                mcf.addEdge(id, get_id(i, 2), 1, 0);
                mcf.addEdge(id, T, 1, w);
            }
            else{
                fore(j, 0, 3){
                    int now = get_id(i, j);
                    mcf.addEdge(now, T, 1, 0);
                }
            }
        }

        ans -= mcf.flow(S, T).se;

        std::cout << ans << endl;
    }
    return 0;
}